BACKGROUND
in sexually reproducing eukaryotes, successful meiosis ensures faithful transmission of a haploid set of chromosomes to the next generation.
problems arising during meiosis can lead to meiotic arrest, chromosomal nondisjunction, and infertility.
a key step in meiosis is the close alignment of homologous chromosomes, a process known as synapsis.
synapsis is typically essential for establishing meiotic crossovers and a specialized, tripartite protein structure known as the synaptonemal complex  forms the foundation for synapsis .
the sc has been cytologically observed across eukaryotes and the molecular components have been characterized in a range of model organisms including arabidopsis thaliana, caenorhabditis elegans, drosophila melanogaster, saccharomyces cerevisiae, mus musculus, and several species of hydra  <cit> .
across this diverse group of eukaryotes the sc maintains, with some exceptions, evolutionary conservation in both structure as a tripartite complex and function in meiotic recombination and synapsis  <cit> .
the sc consists of three main parts in most eukaryotes: the lateral elements , the transverse filaments  and the central element  .
two les run along the length of each pair of sister chromatids and directly interact with the meiotic cohesin complex.
the tfs extend out from the les, resembling rungs of a ladder connecting the juxtaposed chromosomes.
the ce is a solid visible element in the center of the tfs and secures them in place.
some eukaryotes lack an observable sc including schizosaccharomyces pombe and aspergillus nidulans .
in the case of s. pombe, the sc may have been replaced by thin thread-like structures known as the linear elements  <cit> .
d. melanogaster males also lack the sc.
this coincides with the fact that d. melanogaster males also have no meiotic recombination.
these observations indicate that other mechanisms can ensure proper chromosome segregation in the absence of the sc.
despite the strong structural conservation across eukaryotes, the proteins that comprise the sc are strikingly varied  <cit> .
based on the fact that several eukaryote lineages lack the sc , several authors have theorized that the sc has evolved independently multiple times  <cit> .
however, a recent analysis  <cit>  found that m. musculus sc proteins formed monophyletic groups with orthologs in metazoans ranging from cnidarians to humans.
this supports a hypothesis of a single sc origin in at least all metazoans.
the sc of the ecdysozoa  appears substantially different from the sc in other metazoans.
sc components from species such as d. melanogaster and c. elegans show low conservation outside arthropods and nematodes, respectively.
the reason for such lack of conservation of sc components is unknown  <cit> .
several sc proteins have been identified and characterized in d. melanogaster.
five such proteins are included in this study.
em studies in d. melanogaster females indicate the sc is similar in structure to other eukaryotes  <cit>  and all five proteins are contained within the tripartite structure  .
ord and cm have been identified as two of the le proteins in drosophila  <cit> .
ord localizes to the chromosome arms during early prophase i and is necessary for chromosome segregation, loading of the cohesin complex on the chromosomal axis, normal levels of meiotic recombination, and sc stability  <cit> .
its role in crossing over is not entirely understood as recombination is not completely eliminated in ord mutants and there is an increased amount of dsb repair via the sister chromatid.
this suggests that ord suppresses sister chromatid exchange  <cit> .
cm is also a component of the le and is responsible for chromosome core formation  <cit> , sc-dependent meiotic dsb repair, and assembling a continuous ce  <cit> .
the n-terminus of cm lies within the inner region of the le and the c-terminus is assumed to face the central region  <cit> .
so far, cg is the only known drosophila tf protein  <cit> .
like other tf proteins, it has globular n- and c-terminal domains and an internal coiled-coil central domain  <cit> .
cg forms into parallel dimers with the n-terminal globular domains extending into the ce and the c-terminal domains are anchored to the le  <cit> .
cg is necessary for synapsis, conversion of dsbs into crossovers  <cit>  and perhaps gene conversion  <cit> .
finally, the ce is comprised of two other proteins along with the cg n-termini, corona and corolla.
corona, commonly referred to as cona, is a pillar-like protein that aligns outside of the dense ce  <cit> .
cona promotes dsb maturation into crossovers and synapsis does not occur in cona mutants  <cit> .
additionally, cona both co-localizes with cg and stabilizes cg polycomplexes  <cit> .
corolla is also localized within the ce and interacts with cona  <cit> .
thought to be comprised of coiled-coil domains much like cg, it is also essential for sc function and recombination.
all of these proteins have roles exclusive to female meiosis except for ord, which also functions in sister-chromatid cohesion in meiosis i and ii and is necessary for gametogenesis in both drosophila sexes  <cit> .fig.
<dig> a model of the synaptonemal complex in drosophila.
this model is adapted from lake & hawley,  <dig>  <cit>  and collins et al.
<dig>  <cit> .
so far five genes have been found coding sc components and their proteins localized within the structure: ord, cm, cg, corolla, and cona
two hypotheses have been proposed to explain the lack of conservation of the sc: genetic drift and positive selection.
a high rate of evolutionary drift in protein evolution in caenorhabditis and drosophila has been proposed to explain the evolution of the lamin proteins  <cit>  and ribosomal proteins in ecdysozoa  <cit>  as well as olfactory genes in drosophila  <cit> .
low levels of purifying selection on drosophila sc components would allow it to diverge at a high rate resulting in little conservation.
low levels of purifying selection might be expected if the major function of the sc was simply to hold homologs together at a proper distance.
under this scenario, there may be few selective constraints on the particular amino acids that function primarily as structural spacers within the sc.
alternatively, positive selection may contribute to the rapid evolution of sc components.
many studies have demonstrated that reproductive proteins evolve rapidly .
in fact, population genetic analyses in d. melanogaster and close relatives have previously revealed that ord shows a significant deviation from neutrality in d. simulans, with more non-synonymous fixations than expected  <cit> .
recurrent meiotic drive and selection to ameliorate this conflict has been proposed to drive positive selection in meiosis genes .
we aimed to perform a molecular evolutionary analysis of the sc proteins in drosophila to determine what forces may be driving the high rate of evolution of these proteins.
using the genomic sequence data available for different drosophila species and d. melanogaster population data, we aimed to test the null hypothesis that divergence in sc proteins is effectively neutral.
in addition, we sought to test the hypothesis that patterns of molecular evolution in sc components are uniform across the genus.
finally, we examined available d. melanogaster population data to determine if any deviations from neutrality have occurred in recent time, which would be consistent with ongoing positive selection.
methods
ortholog search
the amino acid sequences of cm , cg , cona , corolla  and ord  in d. melanogaster were acquired from flybase  <dig>   <cit> .
an additional sc component, solo, was not examined due to the fact that it is an alternative splice variant of vasa, which is known to play a role in pirna biogenesis  <cit> .
these were used in a tblastn  <cit>  homolog search in  <dig> available genomes of drosophila species with a liberal cutoff of e =  <dig> .
this liberal cutoff was chosen to ensure detection of highly divergent orthologs that were subjected to further validation.
d. melanogaster, d. sechellia, d. yakuba, d. erecta, d. ananassae, d. pseudoobscura, d. willistoni, d. virilis, d. mojavensis, and d. grimshawi genomes were obtained from flybase  <cit> .
the genomes for d. ficusphila , d. eugracilis , d. biarmipes , d. takahashii , d. elegans , d. bipectinata , and d. miranda  were obtained from ncbi.
the genome of d. simulans was obtained from the andolfatto lab server  <cit>  and the d. mauritiana genome was obtained from the schlötterer lab server  <cit> .
to identify highly divergent orthologs, an additional tblastn search was performed using the most diverged protein sequence captured in the original tblastn search.
these results were combined with results from blastp searches of annotated proteins using the d. melanogaster protein sequence.
finally, we included additional ortholog searches with hmmer  <dig> b <dig>  <cit>  and phylomedb v <dig>  <cit>  as well as orthologs listed in orthodb v <dig>  <cit> .
this combined approach allowed us to obtain a broad list of candidate orthologs for each of the five sc components.
orthology was then evaluated for candidates by using a reciprocal best blast hits approach with tblastn.
in all cases where orthology was determined the second reciprocal blast hit e-value was substantially worse than the ortholog e-value.
in addition, synteny for orthologs was evaluated , though it should be noted that there is substantial gene shuffling within muller elements across the genus  <cit> .
sequence retrieval
upon identification of orthologs, sequences from annotated and un-annotated genomes were extracted using identical approaches to limit biases that might arise from using gene annotations only from annotated genomes.
dna sequences 3000 bp upstream and downstream of identified orthologous sequences were first extracted.
these were analyzed with fgenesh+, a hidden markov model protein-based gene predictor used to identify the open reading frames in un-annotated dna sequence using a known protein sequence as a guide  <cit> .
we included  <dig> nucleotides of upstream and downstream flanking sequence to ensure that parts of the open reading frame not originally identified in tblastn were included.
the d. melanogaster amino acid sequence was used as the guide.
sequence alignments and drosophila phylogeny
sequence alignments were generated using coding sequences  obtained with fgenesh+ from each species using both translational mafft  <cit>  and translational muscle  <cit>  in geneious v <dig>   <cit>  with default parameters.
sequence alignments were also generated using codon-based prank  <cit>  based on a pre-determined phylogenetic tree  with the “-f” option allowing insertions.
these three alignment programs were used to evaluate sensitivity of results to alignment procedure.
concatenated alignments of sc sequences  were used to generate phylogenetic trees required for prank alignment and other analyses.
phylogenetic analysis was performed using the cipres science gateway v <dig>  with raxml-hpc blackbox using default parameters and a gtr model with  <dig> bootstrap iterations  <cit> .
the tree topologies produced by concatenated mafft and muscle alignments were identical to each other.
the sc gene tree topology also matched the known phylogeny for the drosophila species used in this analysis  <cit> .
molecular evolutionary analysis
the global omega  value, often referred to as the global dn/ds estimate, is a measure of the average selective pressure acting on a gene across an entire phylogeny  <cit> .
global ω for each alignment was calculated using hyphy with a gtr model  <cit>  and also with the one-ratio model f3x <dig> codon model  in the codeml program of paml v <dig>   <cit> .
both analyses made use of the tree topology obtained from phylogenetic analysis described above.
global ω estimates were obtained using all available orthologs, a smaller subset of  <dig> species within the melanogaster group , and an even smaller subset of six species within the melanogaster subgroup .
estimates were obtained at different levels of divergence to account for potential problems that might occur in the alignment of highly diverged protein sequence.
to quantify heterogeneity in selection pressure, alignments were analyzed with ga branch using a gtr model of nucleotide substitution  <cit>  and the previously described phylogenetic tree.
analysis was performed using datamonkey, the hyphy web server  <cit> .
ga branch uses a genetic algorithm and the akaike information criterion to identify the best fitting model for the number of branch ω classes.
this allows one to evaluate evidence for heterogeneity in ω across the tree.
a model-averaged probability of positive selection  on any of these branches is used to test whether positive selection has occurred.
an analysis of ω was also performed in paml  <cit>  by comparing two different codon based models of evolution.
a likelihood ratio-test was performed to compare a model allowing a beta-distributed value of global ω ranging from zero to one  to a model that also included an additional class of codons with ω greater than one .
both of these models were run with the f3xf <dig> codon model using the nucleotide frequencies at each codon separately and the phylogenetic tree constructed above.
tests of neutrality using polymorphism and divergence
while codon models of molecular evolution provide insight into long-term patterns of selection acting on protein coding sequence, population genetic analyses allow for tests of neutrality in more recent time.
mcdonald-kreitman  tests of neutrality were performed using polymorphism data from two d. melanogaster populations and d. simulans and d. yakuba reference genomes served as outgroups.
the drosophila genetic reference panel v <dig>   <cit>  provided dna sequences from  <dig> d. melanogaster isofemale lines collected from a population in raleigh, north carolina.
in addition,  <dig> genomes from the drosophila population genomics project v <dig>   <cit>  from  <dig> separate populations in sub-sahara africa were used.
sc gene sequences were collected using blast with d. melanogaster reference genes as the query.
blast was performed locally in geneious.
gaps in the alignment were removed and mk tests were performed online with the standardized and generalized mk test website  <cit> .
polarized mk tests were also performed using d. yakuba sequences to polarize lineage-specific substitutions.
in addition, gammamap  <cit>  was used to identify particular codons within the sc genes of d. melanogaster that have likely been fixed by positive selection.
a challenge of the mk test is that polymorphic sites are treated equally and allele frequencies are not taken into account.
in contrast, gammamap utilizes population and divergence data fully.
under a codon model of evolution, polymorphism and divergence data are used to estimate the distribution of fitness effects  for new mutations and substitutions.
gammamap estimates the γ parameter for each codon along the length of the gene.
γ is the population-scaled selection coefficient, γ = 2pnes, where p is the ploidy level, ne is the effective population size, and s is the fitness advantage of a derived allele relative to the ancestral allele if the derived amino acid differs from the ancestral allele.
evidence for positive selection driving an amino acid substitution in d. melanogaster was deemed significant if the probability of γ greater than  <dig> was greater than  <dig>  in d. melanogaster following wilson et al.
.
in addition, dnasp  <dig> .
<dig>  <cit>  was used to estimate average pairwise differences within each gene  and we compared these to the average pairwise site differences for other meiosis genes previously measured  <cit> .
tajima’s d was also calculated in dnasp  <cit> .
haplotype structure was illustrated with phylogenetic trees built using upgma, a hierarchal clustering method  <cit> , in geneious  <dig> .
<dig>  <cit> .
RESULTS
distant orthologs of drosophila sc components are elusive using diverse search methods
we assembled a list of candidate orthologs of sc components in d. melanogaster using blast, the hmmer search tool  <cit> , and by consulting databases of listed orthologous genes including phylomedb and orthodb .
orthologs were validated using the reciprocal best blast hit approach and hits were consistent with prior ortholog annotations.
only cm and ord orthologs could be identified in all drosophila species and further outside the genus .
the le gene sequences were identified in every drosophila species by tblastn and in several closely related diptera species using blastp against annotated proteins .
these include bactrocera cucurbitae , b. dorsilas , ceratitis capitata , musca domestica  and glossina morsitans morsitans .
the remaining three sc components, cg, corolla, and cona, could be identified in all species of drosophila with annotated genomes using blastp .
the one exception is that cona could not be identified within d. willistoni .
none of the tf and ce gene sequences could be identified outside of the drosophila genus.
these results suggest that the tf and ce proteins are less conserved than those comprising the le.
sc genes are evolving quickly and according to position within the sc
hyphy and paml were used to calculate global ω with sequences obtained from the tblastn search.
orthologs that were only identified with blastp could not be reasonably aligned.
thus, the orthologs of cg in d. willistoni and the drosophila subgenus and orthologs of cona in d. ananassae, d. bipectinata, d. willistoni, and the drosophila subgenus were not included in the molecular evolutionary analyses .
to account for possible issues with alignment quality for divergent sequences, we generated alignments with mafft, muscle, and prank.
the global ω estimates were robust to the three alignment methods .
to account for long divergence times between many of the drosophila species, global ω was also estimated across three different scales of divergence.
we selected a subset of  <dig> species within the melanogaster group  and an even smaller set of six species within the melanogaster subgroup .
global ω estimates were similar across different scales of divergence and different alignment methods .
the global ω of each sc component was higher than the median ω for each gene ontology  category in drosophila  <cit> .
the majority of genes within drosophila have ω estimates less than  <dig>   <cit>  and only two go categories have a median ω greater than  <dig>    <cit> .
ord has the lowest ω amongst all the sc genes at ~  <dig>  which is twice as high as the median ω for odorant binding genes and greater than the reported value for seminal fluid proteins  in the d. melanogaster species group  <cit> .fig.
<dig> the ω ratio of each sc gene is high and robust to alignment and divergence level.
ω was estimated in paml using a gtr nucleotide substitution model with 95 % confidence intervals.
the ratio remains relatively consistent for each alignment program used  and the sampling from the phylogeny.
ω estimates are reflected in spatial orientation of the protein within the complex respective to the chromosome; ω increases as distance from the chromosome increases.
the dotted line indicates the median ω of odorant binding genes, the go term with the largest median ω, from the  <dig> drosophila genomes consortium  <cit>
there is an apparent relationship between position within the sc and ω. although the le component ord is evolving at more than twice the average genome-wide rate ratio, it has the lowest value of ω in the sc .
cona is evolving with the highest rate ratio  and the global ω estimate is even higher within the species of the melanogaster subgroup .
the estimates of ω increase as a function of position within the sc: lateral element components evolve the slowest, central element components evolve the fastest, and cg, which functions as a transverse filament, evolves at an intermediate rate.
because we have only characterized five proteins, there is little power in a test for significance in this relationship.
however, it is worth noting that this result is robust to different time scales of analysis.
evolutionary rate ratio variation and signatures of positive selection
we further tested for heterogeneity in ω estimates across the genus.
ga branch  <cit>  uses a genetic algorithm to estimate and evaluate evidence for multiple classes of ω within a phylogenetic context using the akaike information criteria.
it further tests a model for averaged probability for ω >  <dig> for each branch.
results from ga branch indicate that the evolutionary rate of sc components has varied considerably .
cg and corolla have the fewest evolutionary rate ratio classes , ord had the most , and cm and cona both have four rate ratio classes .
there was support for positive selection  on at least one branch in every sc-coding gene except cm.
corolla had the highest ω estimate in any of the ga branch analyses.
corolla also demonstrated a strong signature of positive selection on the branch containing d. biarmipes and d. takahashii and also the branch prior to the split between d. eugracilis and the melanogaster subgroup .
cona shows the most branches with signatures of positive selection .
the le protein ord has the lowest global ω but shows multiple branches with high probabilities of positive selection within the obscura group and prior to the d. eugracilis and melanogaster subgroup divergence.
along with the fact that ord had the most ω rate classes, this suggests that the evolution of ord is highly variable even amongst sc components.
it should be noted that since alignment of divergent sequences can be challenging, ω estimates on deep internal branches might not be precise.
however, rate ratio variation and significant evidence for positive selection are clearly evident on terminal branches.
in particular, for each gene, support for the highest ω class on the phylogeny is evident on at least one terminal or near terminal branch.fig.
<dig> ga branch analysis of the drosophila phylogenetic tree reveals heterogeneity of evolutionary rates for each sc gene.
a supported rate ratio classes correspond to branch colors.
numbers on the branches present the posterior probability that a gene has evolved under positive selection along that particular branch.
the phylogenetic trees correspond to the sequences of ord, cm, cg, corolla, and cona used in the molecular evolution analyses.
b evolutionary rate covariation analysis.
evolutionary rate covariation values  are above the diagonal.
values closer to  <dig> indicate higher levels of covariation.
p-values are below the diagonal
given this rate ratio heterogeneity, we sought to evaluate whether changes in ω estimates tended to co-occur among sc components.
this would be the case if structural changes in one sc component drove structural changes in other sc components.
a simple test for a correlation between branch ω estimates of different sc components must control for shared demographic changes that influence all proteins in the genome.
therefore, we employed the method of evolutionary rate covariation  <cit> .
clear, alignable orthologs of cona are found in the fewest number of species and cona was not included in this analysis, limiting this analysis to four sc components.
we find significant evidence that ω estimates are correlated between ord and corolla and also ord and cm .
cg shows no significant evidence of evolutionary rate co-variation with any other component, even though it interacts with both the lateral element and the central element.
evidence for positive selection across the genus was evaluated using the m <dig> vs. m <dig> test in paml.
two models of evolution were compared using a likelihood ratio test; a model with beta-distributed ω values less than one  and the same model with an additional class of codons with ω values greater than one   <cit> .
a significant likelihood test indicates a signature of positive selection.
positive selection is evident in corolla and this result is robust to both alignment procedure and sampling across different levels of divergence .
ga branch also identified at least one branch with evidence of positive selection within each of the three levels of divergence.
cg also demonstrated evidence for positive selection within the drosophila genus and melanogaster group but none was detected within the six species in the melanogaster subgroup.
this is consistent with results from ga branch that only identified branches with ω estimates near one outside of this clade.
in contrast, ord showed significant evidence for positive selection in the melanogaster subgroup and nowhere else.
the likelihood ratio tests and ga branch both suggest that while ord is the most conserved of the sc components, positive selection intermittently contributes to its divergence.
no signatures of positive selection were detected in cm and cona.
for cm, this is consistent with results from ga branch.
however, the failure to reject a model of neutral evolution in cona stands in contrast to the positive selection detected on multiple branches by ga branch.
this may be explained by the fact that the cona coding sequence is much shorter and multiple branches were identified to be very conserved in ga branch.
under these circumstances, global paml analysis of cona may have reduced power to detect a class of codons with ω greater than one.table  <dig>
p-values of a likelihood ratio test between a model of variable selection pressures with no positive selection  and the same model with positive selection  for each sc component*
<dig> e-03
<dig> e-03
<dig> e-02
<dig>
<dig> e-07
<dig> e-03
<dig> e-16
<dig> e-03
<dig> e-03
<dig>
<dig> e-07
<dig> e-08
<dig> e-03
<dig> e-09
<dig> e-07
<dig>
<dig> e-04
<dig> e-03
*significant p-values in bold
the results of ga branch and paml complement each other and detect positive selection in most of the sc components.
both agree that cm shows no sign of positive selection anywhere in the phylogeny or across different divergence times.
the tf protein cg does show signatures of positive selection outside of the melanogaster subgroup in both tests.
likewise, corolla shows evidence of positive selection throughout the drosophila phylogeny across different time scales of divergence.
despite having the lowest calculated ω, ord shows strong a signature of positive selection within the melanogaster subgroup.
polymorphism and divergence in the d. melanogaster subgroup
to characterize the forces that have shaped the evolution of sc components in more recent time, we turn to readily available population data for d. melanogaster.
we used the second drosophila population genomics project african survey of  <dig> genomes from  <dig> african d. melanogaster populations  <cit>  as well as  <dig> genomes made available by the drosophila genetic reference panel, a sampling of inbred lines from raleigh, north carolina  <cit> .
we performed a series of mcdonald-kreitman  tests  <cit>  using d. simulans sequences as an outgroup to test neutrality in divergence of sc components.
to account for deleterious recessive polymorphisms that are retained at low frequencies, we removed singletons, doubletons, and tripletons.
additionally, the mk test can be used to calculate an alpha parameter – the proportion of substitutions that are positively selected  <cit> .
a negative alpha value indicates the fixation or segregation of deleterious mutations within the gene.
polarized mk tests were also performed with the d. yakuba sequence as an outgroup.
the mk test revealed evidence for deviation from neutrality in some, but not all, sc components.
using population genetic data from d. melanogaster and d. simulans as an outgroup, an unpolarized mk test does not localize signatures of deviation from neutrality to a certain branch.
polarizing fixations on the d. melanogaster branch with d. yakuba as an additional outgroup allows one to determine whether the deviation from neutrality occurred on the d. melanogaster lineage.
across all tests, we find no evidence for recent selection in cm and cona , consistent with molecular evolutionary analyses.
in contrast to its overall slowest ω estimate, but consistent with paml results in the d. melanogaster subgroup , ord is the only gene found to deviate from neutrality in both the polarized and unpolarized mk tests , supporting previous results  <cit> .
positive alpha values from the polarized mk test indicate recent positive selection in d. melanogaster.
evidence for positive selection was found for cg and cona in the unpolarized test using african populations only.
however, polarized tests that examine fixations on the d. melanogaster lineage fail to reject neutrality for cg and cona.
thus, the signature of positive selection in cg and cona can be attributed to changes on the d. simulans lineage.table  <dig> mcdonald-kreitman tests  detecting deviation from neutrality within two population samples of d. melanogaster for all sc components*
ord
<dig>
<dig>
<dig>
p
<dig>
<dig>
<dig>
cm
p
cg
<dig>
p
<dig>
corolla
p
cona
<dig>
p
<dig>
*significant p-values in bold
adetects deviation within d. melanogaster and d. simulans
bdetects deviation within d. melanogaster exclusively
further investigation revealed d. simulans was more highly divergent when compared to both d. melanogaster in four sc components , with ord being the exception.
cg and cona both show an excess of non-synonymous divergence within d. simulans .
thus, the results of the mk tests for cg and cona can be explained by an excess level of non-synonymous divergence on the d. simulans lineage.
this observation is also made in the ga branch analysis .
though not significant, both cm and cona show a similar pattern of increased non-synonymous divergence in d. simulans.
pooling polarized fixations in every sc gene revealed significantly more non-synonymous fixations in d. simulans than d. melanogaster .table  <dig> fisher’s exact tests reveal an increase of non-synonymous substitutions on the d. simulans lineage*
d. mel
d. sim
d. mel
d. sim
ord
cm
cg
<dig>
<dig>
corolla
cona
<dig>
<dig>
<dig>
<dig>
*significant p-values in bold
the mk test is inadequate for identifying the codons that have been fixed positive selection.
we therefore complemented the mk approach using gammamap  <cit>  to estimate the γ selection coefficient for each codon.
similar to the mk test, gammamap utilizes both polymorphism and divergence data.
however, it also makes use of frequency data to estimate the strength of selection that has acted individual codons.
the selection coefficient is expressed in terms of γ, which is equal to 2pnes, twice the product of the effect population size multiplied by the ploidy level and the selection coefficient.
in accordance to wilson et al.
2011  <cit> , we used the probability of γ >  <dig> being 50 % or greater as a cutoff for a significant signature of positive selection  <cit> .
since we were using polymorphism data from d. melanogaster, we did not perform estimation of γ in d. simulans.
overall, signatures of positive selection on the d. melanogaster lineage are demonstrated for all sc proteins across the entire length, with the exception of cona.
the distribution of putative selection effects were similar using data from two subpopulations of d. melanogaster , though more codon variants were deemed significant for evidence of positively selection using data from the north american subpopulations compared to african populations.
for example, results for corolla using african data provide no significant evidence for recent positive selection at the 50 % threshold, in contrast to results using north american data.
this is likely an effect of recent demographic history in north america .
additionally, corolla sequences contain many low-frequency segregating alleles that are potentially deleterious.
using dgrp data, no codons in cm were identified to be under significant positive selection while there were six noted in using population data from africa .
overall, many of the same codons estimated to be putatively positively selected using data from one population were also were also found using data from the other population.
ord and corolla show evidence of weak positive selection in specific regions, specifically between codons  <dig> and  <dig> in ord and between codons  <dig> and  <dig> in corolla .
evidence for selection was also concentrated in cm between codons  <dig> and  <dig>  but using data from africa, these sites were not above our threshold of 50 % probability of γ >  <dig>  while there were many codons identified to be under significant positive selection in cg , codons under positive selection appeared dispersed along the length of the coding sequence.
cona showed no particular codons under selection in both d. melanogaster samples despite having the highest calculated global ω. this coincides with the failure to detect deviation from neutrality in the polarized mk test  and a drastic reduction of ω in d. melanogaster according to ga branch .fig.
<dig> gammamap reveals the posterior probability for positive selection coefficient at each codon using population data from the dpgp.
in concordance with wilson et al.
<dig>  a codon is under significant signature of selection when the posterior probabilities of selection  are greater than  <dig>  .
vertical bars illustrate minor allele frequencies in d. melanogaster  and the substitutions are the circular dots.
the colors correspond to d. melanogaster non-synonymous  and synonymous  variants as well as d. simulans non-synonymous  and synonymous  variants.
estimated number of selected codons is indicated in the upper right of each plot
finally, pairwise nucleotide polymorphism  was calculated for each sc gene.
overall, there is a similar level of nucleotide diversity in every sc component when compared to π genome-wide and mean π for meiosis genes reported in anderson, et al.
<dig>  <cit> .
the one exception was for corolla .
corolla estimates of synonymous π are considerably lower in both north america and africa.
considering tajima’s d, only corolla demonstrated a strong negative value , possibly an indication of ongoing positive selection within corolla.
a sliding window analysis of π and tajima’s d reveal that the central region of corolla,  <dig> to  <dig> nucleotides downstream of the start codon, is almost entirely lacking polymorphism save one doubleton in the african populations  and two singletons within north carolina .
in north carolina populations, 250 bp sliding windows within this region reveal gene regions where π =  <dig> .
flanking this central region, polymorphism increases and tajima’s d is negative as many of the site-wise differences can be attributed to singletons, doubletons, and tripletons.
haplotype structure within corolla is illustrated with dendrograms constructed using upgma  <cit> .
a region of possible recurrent selection shows a higher proportion of individuals carrying a single haplotype with no diversity .
crucially, within this span, there are  <dig> base pairs that are completely monomorphic in both africa and north carolina.
flanking this region, there is an increase of diversity and fewer individuals carry the haplotype with no diversity .
this pattern was also observed in the north carolina population .
strikingly, within the 178 bp monomorphic span, there are eight non-synonymous substitutions and one synonymous substitution between d. melanogaster and d. simulans with ω estimated to be  <dig> .
this also corresponds to the region identified with gammamap with the highest density of codons characterized by the highest probability that γ >  <dig> .
this suggests that ongoing positive selection has driven rapid and recurrent change in the protein coding sequence of corolla.
the low levels of nucleotide diversity within corolla in d. melanogaster can not be attributed to strong purifying selection since ka/ks values, another indicator of selective pressure, between d. melanogaster and d. simulans are high .
in the african populations, the genomic region including corolla has reduced polymorphism compared to flanking regions .
however, the signature is less clear within the north carolina population  possibly due to overall less nucleotide diversity in the dgrp sequences in comparison to the dgpg sequences.
this pattern of reduced polymorphism in a 3 kb region is weaker than other signatures of recent positive selection in d. melanogaster .
this may indicate that this pattern of reduced polymorphism in corolla may be a remnant of positive selection that is not as recent or as strong as other examples of recent positive selection.fig.
<dig> sliding window estimates of pairwise diversity and tajima’s d reveal potential signature of positive selection resulting in loss of haplotype diversity.
pairwise diversity  and tajima’s d measured in 250 bp windows along the length of corolla within the dpgp sequences.
introns are indicated in gray bars.
dendrograms constructed using a hky model of upgma between nucleotides 1- <dig> , 701- <dig> , and 1301- <dig>  downstream of the translation start site
discussion
the sc has been identified across diverse eukaryotes with only a few rare exceptions  <cit> .
homologous protein components of the sc can be found in metazoans ranging from mammals to hydra, indicating that the sc is very likely present at the origin of animals.
however, these metazoan sc components are very difficult to detect in ecdysozoa, including d. melanogaster and c. elegans, despite the fact that em studies identify the sc to be structurally similar.
two hypotheses exist for the presence of the sc in the ecdysozoans: either there has been non-homologous replacement of the sc or an extreme amount of divergence in sc homologs from other lineages.
in support of the hypothesis that a high rate of divergence explains lack of apparent sc protein homology between ecdysozoa and other metazoans, we presented evidence that the sc is evolving very rapidly within the drosophila genus.
importantly, there is a relationship between the estimated global ω estimates for each protein and the ability to identify orthologs in divergent taxa.
only two genes, ord and cm, were identified outside of the drosophila genus.
these both comprise the lateral element, interact with chromatin, and their ω estimates are the lowest.
in contrast, cg, corolla, and cona have higher ω estimates and ortholog identification was more difficult in divergent taxa.
therefore, it is reasonable to conclude that the failure to identify orthologs for sc components outside of the drosophila genus is due to their fast rate of evolution, not necessarily by de novo origination within drosophila  <cit> .
such rapid sequence divergence between orthologs may also suggest that sequence identity is not essential for structural integrity of the sc, despite many drosophila-specific sc components sharing remarkable functional homology with sc components in other eukaryotes.
further resolution of this question may require additional approaches to orthology detection that incorporate structural information and ancestral state reconstruction.
alternatively, proteomic analysis of the sc in species outside the drosophila genus may also identify orthologs that this analysis did not.
we further demonstrate that rapid divergence of sequence identity is not effectively neutral and can in part be explained by prevalent and recurrent positive selection within the drosophila species examined.
using ga branch, we find that sc evolution is not uniform as originally hypothesized.
we provide evidence for a range of ω estimates that have significantly fluctuated across time.
ga branch analysis indicated that cona, a component of the ce, had the greatest number of branches with evidence of positive selection.
across the full phylogeny and also the melanogaster group, a comparison of m <dig> and m <dig> models in paml identified the strongest signatures of positive selection in corolla, also a component of the ce ; this same gene also posed a challenge for ortholog detection outside of the genus.
in contrast, ord, a component of the le, was estimated to have the lowest global ω across the genus and a strong signature of positive selection was observed only when examining the six species within the melanogaster group.
we found an increased ω for sc components that do not directly interact with chromatin: components of the ce have the highest ω estimates, components of the le have the lowest and cg, which comprises the transverse filament, has an intermediate estimate.
a higher rate of evolution for ce proteins in drosophila is concordant with the observation that ce components are more dynamic across metazoans compared to other components  <cit> .
from a structural perspective, the chromatin interaction required of the le may constrain the rate of evolution.
however, ce proteins likely interact with a variety of other meiotic proteins.
therefore, a higher rate of evolution in ce proteins may be partly driven by changes in these interactions.
as the sc is so conserved across eukaryotes, what can explain recurrent positive selection of the sc in drosophila?
as previously mentioned, sc components are highly divergent in both drosophila and caenorhabditis.
since both of these genera are in the ecdysozoa, there may be a shared cause of rapid sc divergence within these two lineages.
one shared cause may be the fact that both d. melanogaster and c. elegans have dsb-independent synapsis.
this may lead to reduced constraint on sc components, though it is hard to see how this would lead to recurrent positive selection.
alternatively, there may be different underlying causes for rapid divergence in these two lineages.
there are several features of meiosis that make these lineages unique.
caenorhabditis species have holocentric chromosomes with complete crossover interference.
drosophila males lack both the sc and meiotic recombination.
thus, multiple forces may independently contribute to the high rate of sc protein evolution in these two lineages.
one possibility is that the rapid evolution and positive selection in sc proteins of drosophila is driven by an interaction between the sex-specific nature of the sc and the rapid turnover of centromeric sequences caused by recurrent bouts of meiotic drive.
previous studies have suggested that sex-specific function can relax selective constraint on a gene and allow it to diverge more freely.
this has been proposed to explain the higher divergence of maternally expressed genes such as bicoid .
all of the sc proteins studied have no phenotypic effect in males when mutant, with the exception of ord which also plays a role in sister chromatid cohesion in the first and second division of meiosis in both sexes  <cit> .
this additional burden of constraint required by being functional in both sexes may explain why ord has the lowest ω value among the sc genes examined.
because the sc is expressed only in females, it may be particularly influenced by rapid evolution of centromeric sequences driven by meiotic drive.
in contrast to male meiosis where all four meiotic products become functional gametes, only one of four meiotic products becomes the egg pronucleus, with the remaining three forming polar bodies.
therefore, strong selection in female meiosis can favor a centromere that is biased to enter the pronucleus over an opposing centromere.
a centromeric variant that strongly distorts meiosis in its favor will sweep through the population even though it may convey deleterious effects such as interfering in male spermatogenesis .
this form of competition has been proposed drive rapid evolution of centromeric sequences .
rapid evolution of centromeric sequences arising from centromere drive has also been proposed to explain signatures of positive selection on centromere-associated proteins such as the centromeric variant of histone h <dig>  <cit> .
sc components also have specialized functions at centromeres.
across diverse organisms, early centromeric associations are mediated by components of the sc  <cit> .
for example, in budding yeast, the tf protein zip <dig> is required for early centromere coupling  <cit> , though not through formation of the sc per se  <cit> .
in drosophila, sc components have the unique property of mediating centromere pairing in mitotically dividing germ cells  <cit> .
additionally, the drosophila sc is essential for centromere synapsis within the chromocenter  <cit>  where the sc is first assembled prior to assembling along the length of the chromosome arms.
finally, across diverse organisms, the sc persists in centromeric regions long after sc disassembly from the euchromatin  <cit> .
this persistence likely facilitates proper chromosomal segregation  <cit> .
due to these multiple functions at the centromere, and as has been proposed for centromeric histones  <cit> , positive selection in sc components may be driven by the need to accommodate rapid turnover of centromere sequences driven by bouts of centromere drive in female meiosis.
this signal may be enhanced by the sex-specific nature of the sc in drosophila.
additional support for this hypothesis lies in the conservation of cm when compared to other sc components.
our analyses showed few signs of positive selection in cm beyond its high global ω, which was higher than ord.
in the studies of sc centromere association, cm mutants either showed partially reduced centromere clustering  <cit>  or no effect  <cit> .
cm may show a weaker signature of positive selection compared to other sc components because it has a limited role in centromeric clustering.
it is also worth noting that the sc plays a crucial role in establishing the landscape of recombination in meiosis.
recent studies have shown that selection may act to modify recombination landscapes as a means to reduce the cost of female meiotic drive, particularly by modulating recombination rates near centromeres  <cit> .
previous studies have also shown that the centromere can vary significantly in its effects on local recombination in closely related species of the d. melanogaster group  <cit> .
overall, we propose that positive selection may jointly arise from the role that sc components have at rapidly evolving centromeres and modulation of recombination rates in these regions.
a combination of these forces, along with sex-specificity, may play an important role in driving rapid evolution of this highly conserved structure in drosophila.
CONCLUSIONS
the sc shows little sequence conservation across eukaryotes despite its conserved function in meiotic segregation and recombination.
the genes comprising the drosophila sc show almost no apparent homology when compared to sc components in other model organisms.
we have determined that the sc components in drosophila are evolving rapidly and their ω estimates are higher than observed for most genes.
we conclude that this can be partly explained by positive selection detected in nearly every sc gene.
this contrasts to our understanding of the sc as a conserved structure necessary for fertility.
we propose that the combination of the female-exclusive function of the sc within drosophila, its role in meiotic recombination, and its interaction with centromeres is driving the rapid evolution of the sc within drosophila.
availability of data and material
the datasets supporting the conclusions of this article are included within this article and its additional files  <dig>   <dig> and  <dig>
additional files
additional file 1: supplementary tables and figures.
tables summarizing orthology search and detection, likelihood values for paml tests, and population genetic parameters estimated from the data.
figures describing ω calculated from hyphy, ga branch using muscle- and prank-aligned sequences, the gammamap output for the dgrp, and per-site differences and tajima’s d estimates along the length of corolla for the dpgp.
additional file 2: aligned sequences of drosophila sc genes.
text file containing the aligned sequences of each sc gene from the respective drosophila species.
alignments were performed by mafft, muscle, and prank.
additional file 3: aligned d. melanogaster population sequences from the drosophila population genomics project .
sequences were aligned to d. simulans and d. yakuba.
sequences that did not align or were missing data were removed prior to the mk tests.
additional file 4: aligned d. melanogaster population sequences from the drosophila genetics reference panel .
sequences were aligned to d. simulans and d. yakuba.
sequences that did not align or were missing data were removed prior to the mk tests.
abbreviations
cecentral element
dfedistribution of fitness effects
dgrpdrosophila genomics reference panel
dpgpdrosophila population genomics project
lelateral element
mk testsmcdonald-kreitman tests
scsynaptonemal complex
tftransverse filament
competing interests
i declare that i have no significant competing financial, professional or personal interests that might have influenced the performance or presentation of the work described in this manuscript.
author’s contributions
jpb and lwh conceived the study and design.
lwh performed the analysis.
jpb and lwh wrote the manuscript.
both authors read and approved the final manuscript.

 in sexually reproducing eukaryotes, successful meiosis ensures faithful transmission of a haploid set of chromosomes to the next generation. problems arising during meiosis can lead to meiotic arrest, chromosomal nondisjunction, and infertility. a key step in meiosis is the close alignment of homologous chromosomes, a process known as synapsis. synapsis is typically essential for establishing meiotic crossovers and a specialized, tripartite protein structure known as the synaptonemal complex forms the foundation for synapsis. the sc has been cytologically observed across eukaryotes and the molecular components have been characterized in a range of model organisms including arabidopsis thaliana, caenorhabditis elegans, drosophila melanogaster, saccharomyces cerevisiae, mus musculus, and several species of hydra. across this diverse group of eukaryotes the sc maintains, with some exceptions, evolutionary conservation in both structure as a tripartite complex and function in meiotic recombination and synapsis. the sc consists of three main parts in most eukaryotes: the lateral elements, the transverse filaments and the central element. two les run along the length of each pair of sister chromatids and directly interact with the meiotic cohesin complex. the tfs extend out from the les, resembling rungs of a ladder connecting the juxtaposed chromosomes. the ce is a solid visible element in the center of the tfs and secures them in place. some eukaryotes lack an observable sc including schizosaccharomyces pombe and aspergillus nidulans. in the case of s. pombe, the sc may have been replaced by thin thread-like structures known as the linear elements. d. melanogaster males also lack the sc. this coincides with the fact that d. melanogaster males also have no meiotic recombination. these observations indicate that other mechanisms can ensure proper chromosome segregation in the absence of the sc. despite the strong structural conservation across eukaryotes, the proteins that comprise the sc are strikingly varied. based on the fact that several eukaryote lineages lack the sc, several authors have theorized that the sc has evolved independently multiple times. however, a recent analysis found that m. musculus sc proteins formed monophyletic groups with orthologs in metazoans ranging from cnidarians to humans. this supports a hypothesis of a single sc origin in at least all metazoans. the sc of the ecdysozoa appears substantially different from the sc in other metazoans. sc components from species such as d. melanogaster and c. elegans show low conservation outside arthropods and nematodes, respectively. the reason for such lack of conservation of sc components is unknown. several sc proteins have been identified and characterized in d. melanogaster. five such proteins are included in this study. em studies in d. melanogaster females indicate the sc is similar in structure to other eukaryotes and all five proteins are contained within the tripartite structure. ord and cm have been identified as two of the le proteins in drosophila. ord localizes to the chromosome arms during early prophase i and is necessary for chromosome segregation, loading of the cohesin complex on the chromosomal axis, normal levels of meiotic recombination, and sc stability. its role in crossing over is not entirely understood as recombination is not completely eliminated in ord mutants and there is an increased amount of dsb repair via the sister chromatid. this suggests that ord suppresses sister chromatid exchange. cm is also a component of the le and is responsible for chromosome core formation, sc-dependent meiotic dsb repair, and assembling a continuous ce. the n-terminus of cm lies within the inner region of the le and the c-terminus is assumed to face the central region. so far, cg is the only known drosophila tf protein. like other tf proteins, it has globular n- and c-terminal domains and an internal coiled-coil central domain. cg forms into parallel dimers with the n-terminal globular domains extending into the ce and the c-terminal domains are anchored to the le. cg is necessary for synapsis, conversion of dsbs into crossovers and perhaps gene conversion. finally, the ce is comprised of two other proteins along with the cg n-termini, corona and corolla. corona, commonly referred to as cona, is a pillar-like protein that aligns outside of the dense ce. cona promotes dsb maturation into crossovers and synapsis does not occur in cona mutants. additionally, cona both co-localizes with cg and stabilizes cg polycomplexes. corolla is also localized within the ce and interacts with cona. thought to be comprised of coiled-coil domains much like cg, it is also essential for sc function and recombination. all of these proteins have roles exclusive to female meiosis except for ord, which also functions in sister-chromatid cohesion in meiosis i and ii and is necessary for gametogenesis in both drosophila sexes.fig. a model of the synaptonemal complex in drosophila. this model is adapted from lake & hawley, and collins et al.. so far five genes have been found coding sc components and their proteins localized within the structure: ord, cm, cg, corolla, and cona two hypotheses have been proposed to explain the lack of conservation of the sc: genetic drift and positive selection. a high rate of evolutionary drift in protein evolution in caenorhabditis and drosophila has been proposed to explain the evolution of the lamin proteins and ribosomal proteins in ecdysozoa as well as olfactory genes in drosophila. low levels of purifying selection on drosophila sc components would allow it to diverge at a high rate resulting in little conservation. low levels of purifying selection might be expected if the major function of the sc was simply to hold homologs together at a proper distance. under this scenario, there may be few selective constraints on the particular amino acids that function primarily as structural spacers within the sc. alternatively, positive selection may contribute to the rapid evolution of sc components. many studies have demonstrated that reproductive proteins evolve rapidly. in fact, population genetic analyses in d. melanogaster and close relatives have previously revealed that ord shows a significant deviation from neutrality in d. simulans, with more non-synonymous fixations than expected. recurrent meiotic drive and selection to ameliorate this conflict has been proposed to drive positive selection in meiosis genes. we aimed to perform a molecular evolutionary analysis of the sc proteins in drosophila to determine what forces may be driving the high rate of evolution of these proteins. using the genomic sequence data available for different drosophila species and d. melanogaster population data, we aimed to test the null hypothesis that divergence in sc proteins is effectively neutral. in addition, we sought to test the hypothesis that patterns of molecular evolution in sc components are uniform across the genus. finally, we examined available d. melanogaster population data to determine if any deviations from neutrality have occurred in recent time, which would be consistent with ongoing positive selection. methods ortholog search the amino acid sequences of cm, cg, cona, corolla and ord in d. melanogaster were acquired from flybase. an additional sc component, solo, was not examined due to the fact that it is an alternative splice variant of vasa, which is known to play a role in pirna biogenesis. these were used in a tblastn homolog search in available genomes of drosophila species with a liberal cutoff of e = . this liberal cutoff was chosen to ensure detection of highly divergent orthologs that were subjected to further validation. d. melanogaster, d. sechellia, d. yakuba, d. erecta, d. ananassae, d. pseudoobscura, d. willistoni, d. virilis, d. mojavensis, and d. grimshawi genomes were obtained from flybase. the genomes for d. ficusphila, d. eugracilis, d. biarmipes, d. takahashii, d. elegans, d. bipectinata, and d. miranda were obtained from ncbi. the genome of d. simulans was obtained from the andolfatto lab server and the d. mauritiana genome was obtained from the schlötterer lab server. to identify highly divergent orthologs, an additional tblastn search was performed using the most diverged protein sequence captured in the original tblastn search. these results were combined with results from blastp searches of annotated proteins using the d. melanogaster protein sequence. finally, we included additional ortholog searches with hmmer b and phylomedb v as well as orthologs listed in orthodb v. this combined approach allowed us to obtain a broad list of candidate orthologs for each of the five sc components. orthology was then evaluated for candidates by using a reciprocal best blast hits approach with tblastn. in all cases where orthology was determined the second reciprocal blast hit e-value was substantially worse than the ortholog e-value. in addition, synteny for orthologs was evaluated, though it should be noted that there is substantial gene shuffling within muller elements across the genus. sequence retrieval upon identification of orthologs, sequences from annotated and un-annotated genomes were extracted using identical approaches to limit biases that might arise from using gene annotations only from annotated genomes. dna sequences 3000 bp upstream and downstream of identified orthologous sequences were first extracted. these were analyzed with fgenesh+, a hidden markov model protein-based gene predictor used to identify the open reading frames in un-annotated dna sequence using a known protein sequence as a guide. we included nucleotides of upstream and downstream flanking sequence to ensure that parts of the open reading frame not originally identified in tblastn were included. the d. melanogaster amino acid sequence was used as the guide. sequence alignments and drosophila phylogeny sequence alignments were generated using coding sequences obtained with fgenesh+ from each species using both translational mafft and translational muscle in geneious v with default parameters. sequence alignments were also generated using codon-based prank based on a pre-determined phylogenetic tree with the “-f” option allowing insertions. these three alignment programs were used to evaluate sensitivity of results to alignment procedure. concatenated alignments of sc sequences were used to generate phylogenetic trees required for prank alignment and other analyses. phylogenetic analysis was performed using the cipres science gateway v with raxml-hpc blackbox using default parameters and a gtr model with bootstrap iterations. the tree topologies produced by concatenated mafft and muscle alignments were identical to each other. the sc gene tree topology also matched the known phylogeny for the drosophila species used in this analysis. molecular evolutionary analysis the global omega value, often referred to as the global dn/ds estimate, is a measure of the average selective pressure acting on a gene across an entire phylogeny. global ω for each alignment was calculated using hyphy with a gtr model and also with the one-ratio model f3x codon model in the codeml program of paml v. both analyses made use of the tree topology obtained from phylogenetic analysis described above. global ω estimates were obtained using all available orthologs, a smaller subset of species within the melanogaster group, and an even smaller subset of six species within the melanogaster subgroup. estimates were obtained at different levels of divergence to account for potential problems that might occur in the alignment of highly diverged protein sequence. to quantify heterogeneity in selection pressure, alignments were analyzed with ga branch using a gtr model of nucleotide substitution and the previously described phylogenetic tree. analysis was performed using datamonkey, the hyphy web server. ga branch uses a genetic algorithm and the akaike information criterion to identify the best fitting model for the number of branch ω classes. this allows one to evaluate evidence for heterogeneity in ω across the tree. a model-averaged probability of positive selection on any of these branches is used to test whether positive selection has occurred. an analysis of ω was also performed in paml by comparing two different codon based models of evolution. a likelihood ratio-test was performed to compare a model allowing a beta-distributed value of global ω ranging from zero to one to a model that also included an additional class of codons with ω greater than one. both of these models were run with the f3xf codon model using the nucleotide frequencies at each codon separately and the phylogenetic tree constructed above. tests of neutrality using polymorphism and divergence while codon models of molecular evolution provide insight into long-term patterns of selection acting on protein coding sequence, population genetic analyses allow for tests of neutrality in more recent time. mcdonald-kreitman tests of neutrality were performed using polymorphism data from two d. melanogaster populations and d. simulans and d. yakuba reference genomes served as outgroups. the drosophila genetic reference panel v provided dna sequences from d. melanogaster isofemale lines collected from a population in raleigh, north carolina. in addition, genomes from the drosophila population genomics project v from separate populations in sub-sahara africa were used. sc gene sequences were collected using blast with d. melanogaster reference genes as the query. blast was performed locally in geneious. gaps in the alignment were removed and mk tests were performed online with the standardized and generalized mk test website. polarized mk tests were also performed using d. yakuba sequences to polarize lineage-specific substitutions. in addition, gammamap was used to identify particular codons within the sc genes of d. melanogaster that have likely been fixed by positive selection. a challenge of the mk test is that polymorphic sites are treated equally and allele frequencies are not taken into account. in contrast, gammamap utilizes population and divergence data fully. under a codon model of evolution, polymorphism and divergence data are used to estimate the distribution of fitness effects for new mutations and substitutions. gammamap estimates the γ parameter for each codon along the length of the gene. γ is the population-scaled selection coefficient, γ = 2pnes, where p is the ploidy level, ne is the effective population size, and s is the fitness advantage of a derived allele relative to the ancestral allele if the derived amino acid differs from the ancestral allele. evidence for positive selection driving an amino acid substitution in d. melanogaster was deemed significant if the probability of γ greater than was greater than in d. melanogaster following wilson et al.. in addition, dnasp. was used to estimate average pairwise differences within each gene and we compared these to the average pairwise site differences for other meiosis genes previously measured. tajima’s d was also calculated in dnasp. haplotype structure was illustrated with phylogenetic trees built using upgma, a hierarchal clustering method, in geneious..  distant orthologs of drosophila sc components are elusive using diverse search methods we assembled a list of candidate orthologs of sc components in d. melanogaster using blast, the hmmer search tool, and by consulting databases of listed orthologous genes including phylomedb and orthodb. orthologs were validated using the reciprocal best blast hit approach and hits were consistent with prior ortholog annotations. only cm and ord orthologs could be identified in all drosophila species and further outside the genus. the le gene sequences were identified in every drosophila species by tblastn and in several closely related diptera species using blastp against annotated proteins. these include bactrocera cucurbitae, b. dorsilas, ceratitis capitata, musca domestica and glossina morsitans morsitans. the remaining three sc components, cg, corolla, and cona, could be identified in all species of drosophila with annotated genomes using blastp. the one exception is that cona could not be identified within d. willistoni. none of the tf and ce gene sequences could be identified outside of the drosophila genus. these results suggest that the tf and ce proteins are less conserved than those comprising the le. sc genes are evolving quickly and according to position within the sc hyphy and paml were used to calculate global ω with sequences obtained from the tblastn search. orthologs that were only identified with blastp could not be reasonably aligned. thus, the orthologs of cg in d. willistoni and the drosophila subgenus and orthologs of cona in d. ananassae, d. bipectinata, d. willistoni, and the drosophila subgenus were not included in the molecular evolutionary analyses. to account for possible issues with alignment quality for divergent sequences, we generated alignments with mafft, muscle, and prank. the global ω estimates were robust to the three alignment methods. to account for long divergence times between many of the drosophila species, global ω was also estimated across three different scales of divergence. we selected a subset of species within the melanogaster group and an even smaller set of six species within the melanogaster subgroup. global ω estimates were similar across different scales of divergence and different alignment methods. the global ω of each sc component was higher than the median ω for each gene ontology category in drosophila. the majority of genes within drosophila have ω estimates less than and only two go categories have a median ω greater than. ord has the lowest ω amongst all the sc genes at ~  which is twice as high as the median ω for odorant binding genes and greater than the reported value for seminal fluid proteins in the d. melanogaster species group.fig. the ω ratio of each sc gene is high and robust to alignment and divergence level. ω was estimated in paml using a gtr nucleotide substitution model with 95 % confidence intervals. the ratio remains relatively consistent for each alignment program used and the sampling from the phylogeny. ω estimates are reflected in spatial orientation of the protein within the complex respective to the chromosome; ω increases as distance from the chromosome increases. the dotted line indicates the median ω of odorant binding genes, the go term with the largest median ω, from the drosophila genomes consortium there is an apparent relationship between position within the sc and ω. although the le component ord is evolving at more than twice the average genome-wide rate ratio, it has the lowest value of ω in the sc. cona is evolving with the highest rate ratio and the global ω estimate is even higher within the species of the melanogaster subgroup. the estimates of ω increase as a function of position within the sc: lateral element components evolve the slowest, central element components evolve the fastest, and cg, which functions as a transverse filament, evolves at an intermediate rate. because we have only characterized five proteins, there is little power in a test for significance in this relationship. however, it is worth noting that this result is robust to different time scales of analysis. evolutionary rate ratio variation and signatures of positive selection we further tested for heterogeneity in ω estimates across the genus. ga branch uses a genetic algorithm to estimate and evaluate evidence for multiple classes of ω within a phylogenetic context using the akaike information criteria. it further tests a model for averaged probability for ω >  for each branch. results from ga branch indicate that the evolutionary rate of sc components has varied considerably. cg and corolla have the fewest evolutionary rate ratio classes, ord had the most, and cm and cona both have four rate ratio classes. there was support for positive selection on at least one branch in every sc-coding gene except cm. corolla had the highest ω estimate in any of the ga branch analyses. corolla also demonstrated a strong signature of positive selection on the branch containing d. biarmipes and d. takahashii and also the branch prior to the split between d. eugracilis and the melanogaster subgroup. cona shows the most branches with signatures of positive selection. the le protein ord has the lowest global ω but shows multiple branches with high probabilities of positive selection within the obscura group and prior to the d. eugracilis and melanogaster subgroup divergence. along with the fact that ord had the most ω rate classes, this suggests that the evolution of ord is highly variable even amongst sc components. it should be noted that since alignment of divergent sequences can be challenging, ω estimates on deep internal branches might not be precise. however, rate ratio variation and significant evidence for positive selection are clearly evident on terminal branches. in particular, for each gene, support for the highest ω class on the phylogeny is evident on at least one terminal or near terminal branch.fig. ga branch analysis of the drosophila phylogenetic tree reveals heterogeneity of evolutionary rates for each sc gene. a supported rate ratio classes correspond to branch colors. numbers on the branches present the posterior probability that a gene has evolved under positive selection along that particular branch. the phylogenetic trees correspond to the sequences of ord, cm, cg, corolla, and cona used in the molecular evolution analyses. b evolutionary rate covariation analysis. evolutionary rate covariation values are above the diagonal. values closer to indicate higher levels of covariation. p-values are below the diagonal given this rate ratio heterogeneity, we sought to evaluate whether changes in ω estimates tended to co-occur among sc components. this would be the case if structural changes in one sc component drove structural changes in other sc components. a simple test for a correlation between branch ω estimates of different sc components must control for shared demographic changes that influence all proteins in the genome. therefore, we employed the method of evolutionary rate covariation. clear, alignable orthologs of cona are found in the fewest number of species and cona was not included in this analysis, limiting this analysis to four sc components. we find significant evidence that ω estimates are correlated between ord and corolla and also ord and cm. cg shows no significant evidence of evolutionary rate co-variation with any other component, even though it interacts with both the lateral element and the central element. evidence for positive selection across the genus was evaluated using the m vs. m test in paml. two models of evolution were compared using a likelihood ratio test; a model with beta-distributed ω values less than one and the same model with an additional class of codons with ω values greater than one. a significant likelihood test indicates a signature of positive selection. positive selection is evident in corolla and this result is robust to both alignment procedure and sampling across different levels of divergence. ga branch also identified at least one branch with evidence of positive selection within each of the three levels of divergence. cg also demonstrated evidence for positive selection within the drosophila genus and melanogaster group but none was detected within the six species in the melanogaster subgroup. this is consistent with results from ga branch that only identified branches with ω estimates near one outside of this clade. in contrast, ord showed significant evidence for positive selection in the melanogaster subgroup and nowhere else. the likelihood ratio tests and ga branch both suggest that while ord is the most conserved of the sc components, positive selection intermittently contributes to its divergence. no signatures of positive selection were detected in cm and cona. for cm, this is consistent with results from ga branch. however, the failure to reject a model of neutral evolution in cona stands in contrast to the positive selection detected on multiple branches by ga branch. this may be explained by the fact that the cona coding sequence is much shorter and multiple branches were identified to be very conserved in ga branch. under these circumstances, global paml analysis of cona may have reduced power to detect a class of codons with ω greater than one.table p-values of a likelihood ratio test between a model of variable selection pressures with no positive selection and the same model with positive selection for each sc component* e-03 e-03 e-02 e-07 e-03 e-16 e-03 e-03 e-07 e-08 e-03 e-09 e-07 e-04 e-03 *significant p-values in bold the results of ga branch and paml complement each other and detect positive selection in most of the sc components. both agree that cm shows no sign of positive selection anywhere in the phylogeny or across different divergence times. the tf protein cg does show signatures of positive selection outside of the melanogaster subgroup in both tests. likewise, corolla shows evidence of positive selection throughout the drosophila phylogeny across different time scales of divergence. despite having the lowest calculated ω, ord shows strong a signature of positive selection within the melanogaster subgroup. polymorphism and divergence in the d. melanogaster subgroup to characterize the forces that have shaped the evolution of sc components in more recent time, we turn to readily available population data for d. melanogaster. we used the second drosophila population genomics project african survey of genomes from african d. melanogaster populations as well as genomes made available by the drosophila genetic reference panel, a sampling of inbred lines from raleigh, north carolina. we performed a series of mcdonald-kreitman tests using d. simulans sequences as an outgroup to test neutrality in divergence of sc components. to account for deleterious recessive polymorphisms that are retained at low frequencies, we removed singletons, doubletons, and tripletons. additionally, the mk test can be used to calculate an alpha parameter – the proportion of substitutions that are positively selected. a negative alpha value indicates the fixation or segregation of deleterious mutations within the gene. polarized mk tests were also performed with the d. yakuba sequence as an outgroup. the mk test revealed evidence for deviation from neutrality in some, but not all, sc components. using population genetic data from d. melanogaster and d. simulans as an outgroup, an unpolarized mk test does not localize signatures of deviation from neutrality to a certain branch. polarizing fixations on the d. melanogaster branch with d. yakuba as an additional outgroup allows one to determine whether the deviation from neutrality occurred on the d. melanogaster lineage. across all tests, we find no evidence for recent selection in cm and cona, consistent with molecular evolutionary analyses. in contrast to its overall slowest ω estimate, but consistent with paml results in the d. melanogaster subgroup, ord is the only gene found to deviate from neutrality in both the polarized and unpolarized mk tests, supporting previous results. positive alpha values from the polarized mk test indicate recent positive selection in d. melanogaster. evidence for positive selection was found for cg and cona in the unpolarized test using african populations only. however, polarized tests that examine fixations on the d. melanogaster lineage fail to reject neutrality for cg and cona. thus, the signature of positive selection in cg and cona can be attributed to changes on the d. simulans lineage.table mcdonald-kreitman tests detecting deviation from neutrality within two population samples of d. melanogaster for all sc components* ord p cm p cg p corolla p cona p *significant p-values in bold adetects deviation within d. melanogaster and d. simulans bdetects deviation within d. melanogaster exclusively further investigation revealed d. simulans was more highly divergent when compared to both d. melanogaster in four sc components, with ord being the exception. cg and cona both show an excess of non-synonymous divergence within d. simulans. thus, the results of the mk tests for cg and cona can be explained by an excess level of non-synonymous divergence on the d. simulans lineage. this observation is also made in the ga branch analysis. though not significant, both cm and cona show a similar pattern of increased non-synonymous divergence in d. simulans. pooling polarized fixations in every sc gene revealed significantly more non-synonymous fixations in d. simulans than d. melanogaster.table fisher’s exact tests reveal an increase of non-synonymous substitutions on the d. simulans lineage* d. mel d. sim d. mel d. sim ord cm cg corolla cona *significant p-values in bold the mk test is inadequate for identifying the codons that have been fixed positive selection. we therefore complemented the mk approach using gammamap to estimate the γ selection coefficient for each codon. similar to the mk test, gammamap utilizes both polymorphism and divergence data. however, it also makes use of frequency data to estimate the strength of selection that has acted individual codons. the selection coefficient is expressed in terms of γ, which is equal to 2pnes, twice the product of the effect population size multiplied by the ploidy level and the selection coefficient. in accordance to wilson et al. 2011 , we used the probability of γ >  being 50 % or greater as a cutoff for a significant signature of positive selection. since we were using polymorphism data from d. melanogaster, we did not perform estimation of γ in d. simulans. overall, signatures of positive selection on the d. melanogaster lineage are demonstrated for all sc proteins across the entire length, with the exception of cona. the distribution of putative selection effects were similar using data from two subpopulations of d. melanogaster, though more codon variants were deemed significant for evidence of positively selection using data from the north american subpopulations compared to african populations. for example, results for corolla using african data provide no significant evidence for recent positive selection at the 50 % threshold, in contrast to results using north american data. this is likely an effect of recent demographic history in north america. additionally, corolla sequences contain many low-frequency segregating alleles that are potentially deleterious. using dgrp data, no codons in cm were identified to be under significant positive selection while there were six noted in using population data from africa. overall, many of the same codons estimated to be putatively positively selected using data from one population were also were also found using data from the other population. ord and corolla show evidence of weak positive selection in specific regions, specifically between codons and in ord and between codons and in corolla. evidence for selection was also concentrated in cm between codons and but using data from africa, these sites were not above our threshold of 50 % probability of γ >  while there were many codons identified to be under significant positive selection in cg, codons under positive selection appeared dispersed along the length of the coding sequence. cona showed no particular codons under selection in both d. melanogaster samples despite having the highest calculated global ω. this coincides with the failure to detect deviation from neutrality in the polarized mk test and a drastic reduction of ω in d. melanogaster according to ga branch.fig. gammamap reveals the posterior probability for positive selection coefficient at each codon using population data from the dpgp. in concordance with wilson et al. a codon is under significant signature of selection when the posterior probabilities of selection are greater than. vertical bars illustrate minor allele frequencies in d. melanogaster and the substitutions are the circular dots. the colors correspond to d. melanogaster non-synonymous and synonymous variants as well as d. simulans non-synonymous and synonymous variants. estimated number of selected codons is indicated in the upper right of each plot finally, pairwise nucleotide polymorphism was calculated for each sc gene. overall, there is a similar level of nucleotide diversity in every sc component when compared to π genome-wide and mean π for meiosis genes reported in anderson, et al.. the one exception was for corolla. corolla estimates of synonymous π are considerably lower in both north america and africa. considering tajima’s d, only corolla demonstrated a strong negative value, possibly an indication of ongoing positive selection within corolla. a sliding window analysis of π and tajima’s d reveal that the central region of corolla, to nucleotides downstream of the start codon, is almost entirely lacking polymorphism save one doubleton in the african populations and two singletons within north carolina. in north carolina populations, 250 bp sliding windows within this region reveal gene regions where π = . flanking this central region, polymorphism increases and tajima’s d is negative as many of the site-wise differences can be attributed to singletons, doubletons, and tripletons. haplotype structure within corolla is illustrated with dendrograms constructed using upgma. a region of possible recurrent selection shows a higher proportion of individuals carrying a single haplotype with no diversity. crucially, within this span, there are base pairs that are completely monomorphic in both africa and north carolina. flanking this region, there is an increase of diversity and fewer individuals carry the haplotype with no diversity. this pattern was also observed in the north carolina population. strikingly, within the 178 bp monomorphic span, there are eight non-synonymous substitutions and one synonymous substitution between d. melanogaster and d. simulans with ω estimated to be. this also corresponds to the region identified with gammamap with the highest density of codons characterized by the highest probability that γ > . this suggests that ongoing positive selection has driven rapid and recurrent change in the protein coding sequence of corolla. the low levels of nucleotide diversity within corolla in d. melanogaster can not be attributed to strong purifying selection since ka/ks values, another indicator of selective pressure, between d. melanogaster and d. simulans are high. in the african populations, the genomic region including corolla has reduced polymorphism compared to flanking regions. however, the signature is less clear within the north carolina population possibly due to overall less nucleotide diversity in the dgrp sequences in comparison to the dgpg sequences. this pattern of reduced polymorphism in a 3 kb region is weaker than other signatures of recent positive selection in d. melanogaster. this may indicate that this pattern of reduced polymorphism in corolla may be a remnant of positive selection that is not as recent or as strong as other examples of recent positive selection.fig. sliding window estimates of pairwise diversity and tajima’s d reveal potential signature of positive selection resulting in loss of haplotype diversity. pairwise diversity and tajima’s d measured in 250 bp windows along the length of corolla within the dpgp sequences. introns are indicated in gray bars. dendrograms constructed using a hky model of upgma between nucleotides 1-, 701-, and 1301- downstream of the translation start site discussion the sc has been identified across diverse eukaryotes with only a few rare exceptions. homologous protein components of the sc can be found in metazoans ranging from mammals to hydra, indicating that the sc is very likely present at the origin of animals. however, these metazoan sc components are very difficult to detect in ecdysozoa, including d. melanogaster and c. elegans, despite the fact that em studies identify the sc to be structurally similar. two hypotheses exist for the presence of the sc in the ecdysozoans: either there has been non-homologous replacement of the sc or an extreme amount of divergence in sc homologs from other lineages. in support of the hypothesis that a high rate of divergence explains lack of apparent sc protein homology between ecdysozoa and other metazoans, we presented evidence that the sc is evolving very rapidly within the drosophila genus. importantly, there is a relationship between the estimated global ω estimates for each protein and the ability to identify orthologs in divergent taxa. only two genes, ord and cm, were identified outside of the drosophila genus. these both comprise the lateral element, interact with chromatin, and their ω estimates are the lowest. in contrast, cg, corolla, and cona have higher ω estimates and ortholog identification was more difficult in divergent taxa. therefore, it is reasonable to conclude that the failure to identify orthologs for sc components outside of the drosophila genus is due to their fast rate of evolution, not necessarily by de novo origination within drosophila. such rapid sequence divergence between orthologs may also suggest that sequence identity is not essential for structural integrity of the sc, despite many drosophila-specific sc components sharing remarkable functional homology with sc components in other eukaryotes. further resolution of this question may require additional approaches to orthology detection that incorporate structural information and ancestral state reconstruction. alternatively, proteomic analysis of the sc in species outside the drosophila genus may also identify orthologs that this analysis did not. we further demonstrate that rapid divergence of sequence identity is not effectively neutral and can in part be explained by prevalent and recurrent positive selection within the drosophila species examined. using ga branch, we find that sc evolution is not uniform as originally hypothesized. we provide evidence for a range of ω estimates that have significantly fluctuated across time. ga branch analysis indicated that cona, a component of the ce, had the greatest number of branches with evidence of positive selection. across the full phylogeny and also the melanogaster group, a comparison of m and m models in paml identified the strongest signatures of positive selection in corolla, also a component of the ce; this same gene also posed a challenge for ortholog detection outside of the genus. in contrast, ord, a component of the le, was estimated to have the lowest global ω across the genus and a strong signature of positive selection was observed only when examining the six species within the melanogaster group. we found an increased ω for sc components that do not directly interact with chromatin: components of the ce have the highest ω estimates, components of the le have the lowest and cg, which comprises the transverse filament, has an intermediate estimate. a higher rate of evolution for ce proteins in drosophila is concordant with the observation that ce components are more dynamic across metazoans compared to other components. from a structural perspective, the chromatin interaction required of the le may constrain the rate of evolution. however, ce proteins likely interact with a variety of other meiotic proteins. therefore, a higher rate of evolution in ce proteins may be partly driven by changes in these interactions. as the sc is so conserved across eukaryotes, what can explain recurrent positive selection of the sc in drosophila? as previously mentioned, sc components are highly divergent in both drosophila and caenorhabditis. since both of these genera are in the ecdysozoa, there may be a shared cause of rapid sc divergence within these two lineages. one shared cause may be the fact that both d. melanogaster and c. elegans have dsb-independent synapsis. this may lead to reduced constraint on sc components, though it is hard to see how this would lead to recurrent positive selection. alternatively, there may be different underlying causes for rapid divergence in these two lineages. there are several features of meiosis that make these lineages unique. caenorhabditis species have holocentric chromosomes with complete crossover interference. drosophila males lack both the sc and meiotic recombination. thus, multiple forces may independently contribute to the high rate of sc protein evolution in these two lineages. one possibility is that the rapid evolution and positive selection in sc proteins of drosophila is driven by an interaction between the sex-specific nature of the sc and the rapid turnover of centromeric sequences caused by recurrent bouts of meiotic drive. previous studies have suggested that sex-specific function can relax selective constraint on a gene and allow it to diverge more freely. this has been proposed to explain the higher divergence of maternally expressed genes such as bicoid. all of the sc proteins studied have no phenotypic effect in males when mutant, with the exception of ord which also plays a role in sister chromatid cohesion in the first and second division of meiosis in both sexes. this additional burden of constraint required by being functional in both sexes may explain why ord has the lowest ω value among the sc genes examined. because the sc is expressed only in females, it may be particularly influenced by rapid evolution of centromeric sequences driven by meiotic drive. in contrast to male meiosis where all four meiotic products become functional gametes, only one of four meiotic products becomes the egg pronucleus, with the remaining three forming polar bodies. therefore, strong selection in female meiosis can favor a centromere that is biased to enter the pronucleus over an opposing centromere. a centromeric variant that strongly distorts meiosis in its favor will sweep through the population even though it may convey deleterious effects such as interfering in male spermatogenesis. this form of competition has been proposed drive rapid evolution of centromeric sequences. rapid evolution of centromeric sequences arising from centromere drive has also been proposed to explain signatures of positive selection on centromere-associated proteins such as the centromeric variant of histone h. sc components also have specialized functions at centromeres. across diverse organisms, early centromeric associations are mediated by components of the sc. for example, in budding yeast, the tf protein zip is required for early centromere coupling , though not through formation of the sc per se. in drosophila, sc components have the unique property of mediating centromere pairing in mitotically dividing germ cells. additionally, the drosophila sc is essential for centromere synapsis within the chromocenter where the sc is first assembled prior to assembling along the length of the chromosome arms. finally, across diverse organisms, the sc persists in centromeric regions long after sc disassembly from the euchromatin. this persistence likely facilitates proper chromosomal segregation. due to these multiple functions at the centromere, and as has been proposed for centromeric histones, positive selection in sc components may be driven by the need to accommodate rapid turnover of centromere sequences driven by bouts of centromere drive in female meiosis. this signal may be enhanced by the sex-specific nature of the sc in drosophila. additional support for this hypothesis lies in the conservation of cm when compared to other sc components. our analyses showed few signs of positive selection in cm beyond its high global ω, which was higher than ord. in the studies of sc centromere association, cm mutants either showed partially reduced centromere clustering or no effect. cm may show a weaker signature of positive selection compared to other sc components because it has a limited role in centromeric clustering. it is also worth noting that the sc plays a crucial role in establishing the landscape of recombination in meiosis. recent studies have shown that selection may act to modify recombination landscapes as a means to reduce the cost of female meiotic drive, particularly by modulating recombination rates near centromeres. previous studies have also shown that the centromere can vary significantly in its effects on local recombination in closely related species of the d. melanogaster group. overall, we propose that positive selection may jointly arise from the role that sc components have at rapidly evolving centromeres and modulation of recombination rates in these regions. a combination of these forces, along with sex-specificity, may play an important role in driving rapid evolution of this highly conserved structure in drosophila.  the sc shows little sequence conservation across eukaryotes despite its conserved function in meiotic segregation and recombination. the genes comprising the drosophila sc show almost no apparent homology when compared to sc components in other model organisms. we have determined that the sc components in drosophila are evolving rapidly and their ω estimates are higher than observed for most genes. we conclude that this can be partly explained by positive selection detected in nearly every sc gene. this contrasts to our understanding of the sc as a conserved structure necessary for fertility. we propose that the combination of the female-exclusive function of the sc within drosophila, its role in meiotic recombination, and its interaction with centromeres is driving the rapid evolution of the sc within drosophila. availability of data and material the datasets supporting the conclusions of this article are included within this article and its additional files and additional files additional file 1: supplementary tables and figures. tables summarizing orthology search and detection, likelihood values for paml tests, and population genetic parameters estimated from the data. figures describing ω calculated from hyphy, ga branch using muscle- and prank-aligned sequences, the gammamap output for the dgrp, and per-site differences and tajima’s d estimates along the length of corolla for the dpgp. additional file 2: aligned sequences of drosophila sc genes. text file containing the aligned sequences of each sc gene from the respective drosophila species. alignments were performed by mafft, muscle, and prank. additional file 3: aligned d. melanogaster population sequences from the drosophila population genomics project. sequences were aligned to d. simulans and d. yakuba. sequences that did not align or were missing data were removed prior to the mk tests. additional file 4: aligned d. melanogaster population sequences from the drosophila genetics reference panel. sequences were aligned to d. simulans and d. yakuba. sequences that did not align or were missing data were removed prior to the mk tests. abbreviations cecentral element dfedistribution of fitness effects dgrpdrosophila genomics reference panel dpgpdrosophila population genomics project lelateral element mk testsmcdonald-kreitman tests scsynaptonemal complex tftransverse filament competing interests i declare that i have no significant competing financial, professional or personal interests that might have influenced the performance or presentation of the work described in this manuscript. author’s contributions jpb and lwh conceived the study and design. lwh performed the analysis. jpb and lwh wrote the manuscript. both authors read and approved the final manuscript. 
the synaptonemal complex is a highly conserved meiotic structure that functions to pair homologs and facilitate meiotic recombination in most eukaryotes. the sc is required for meiotic recombination in drosophila and absence of these proteins leads to reduced crossing over and chromosomal nondisjunction. despite the conserved nature of the sc and the key role that these five proteins have in meiosis in d. melanogaster, they display little apparent sequence conservation outside the genus. to identify factors that explain this lack of apparent conservation, we performed a molecular evolutionary analysis of these genes across the drosophila genus. sc components within drosophila show little apparent sequence homology to those identified in other model organisms due to their rapid evolution. sc gene sequences have a higher dn/ds rate ratio than the genome wide average and this can in part be explained by the action of positive selection in almost every sc component. finally, using population genetic approaches, we demonstrate that positive selection on sc components is ongoing. we also propose that positive selection in the sc is driven by its sex-specificity combined with its role in facilitating both recombination and centromere clustering in the face of recurrent bouts of drive in female meiosis. drosophilasynaptonemal complexpositive selectionhttp://dx.doi.org/ /100000001national science foundationmcb-1022165nsf mcb-1413532blumenstiel justin p. issue-copyright-statement© the author 2016 first, we propose that a high rate of evolution can be partly explained by low purifying selection on protein components whose function is to simply hold chromosomes together. it further appears that ω estimates for the five sc components are in accordance with their physical position within the sc. for the five sc components, gene sequence similarity declines rapidly with increasing phylogenetic distance and only ord and cm are identifiable outside of the drosophila genus. five drosophila sc proteins have been identified and localized within the complex: cg, cm, cona, ord, and the newly identified corolla. we propose that the drosophila sc is evolving rapidly due to two combined effects.
 the synaptonemal complex is a highly conserved meiotic structure that functions to pair homologs and facilitate meiotic recombination in most eukaryotes. five drosophila sc proteins have been identified and localized within the complex: cg, cm, cona, ord, and the newly identified corolla. the sc is required for meiotic recombination in drosophila and absence of these proteins leads to reduced crossing over and chromosomal nondisjunction. despite the conserved nature of the sc and the key role that these five proteins have in meiosis in d. melanogaster, they display little apparent sequence conservation outside the genus. to identify factors that explain this lack of apparent conservation, we performed a molecular evolutionary analysis of these genes across the drosophila genus.  for the five sc components, gene sequence similarity declines rapidly with increasing phylogenetic distance and only ord and cm are identifiable outside of the drosophila genus. sc gene sequences have a higher dn/ds rate ratio than the genome wide average and this can in part be explained by the action of positive selection in almost every sc component. across the genus, there is significant variation in ω for each protein. it further appears that ω estimates for the five sc components are in accordance with their physical position within the sc. components interacting with chromatin evolve slowest and components comprising the central elements evolve the most rapidly. finally, using population genetic approaches, we demonstrate that positive selection on sc components is ongoing.  sc components within drosophila show little apparent sequence homology to those identified in other model organisms due to their rapid evolution. we propose that the drosophila sc is evolving rapidly due to two combined effects. first, we propose that a high rate of evolution can be partly explained by low purifying selection on protein components whose function is to simply hold chromosomes together. we also propose that positive selection in the sc is driven by its sex-specificity combined with its role in facilitating both recombination and centromere clustering in the face of recurrent bouts of drive in female meiosis. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords drosophilasynaptonemal complexpositive selectionhttp://dx.doi.org/ /100000001national science foundationmcb-1022165nsf mcb-1413532blumenstiel justin p. issue-copyright-statement© the author 2016 