I am an evolutionary biologist interested in sex chromosome evolution, sexual selection and speciation.
Caenorhabditis
Please check the work of students
Salt conditioning in hermaphrodite and gonochoristic species
Hybrid incompatibility in between-species crosses
Rana temporaria
Within population variation in Y chromosome differentiation level
No evidence that Y-chromosome differentiation affects male fitness in a Swiss population of common frogs
Male frogs in the same population with Y chromosomes of varying differentiation level from the X have no significant difference in morphometrics, amplexus success and fathering success.
The standard theory of sex chromosome evolution predicts that they will accumulate mutations that are beneficial for one sex and deleterious to the other. Such genes are termed sexually antagonistic.
The reason sex chromosomes attract sexually antagonistic genes is that they spend the majority of their time in one sex. For example, the Y is only found in males, and the X spends twice the time in females (XX) compared to males (XY) in mammals. By localising on sex chromosomes, sexually antagonistic genes can therefore maximise the positive effects to the sex in which they occur most often and minimise the negative effects the same alleles have in the opposite sex.
The common frog presents a rare system to test this idea, because populations can be found in which Y chromosomes of varying differentiation to the X exist. We worked on a population with 3 major Y-haplotypes: one which is fully differentiated across its length, one which is only differentiated close to a candidate sex-determining gene, and one which is indistinguishable from the X (XX males). These haplotypes represent different portions of the genome that spend their time in males; the absence of differentiation indicates regular recombination. According to the standard sex chromosome evolution theory, one might expect the fully differentiated haplotype to have accumulated more male beneficial alleles than the other haplotypes found in males.
To do this, we measured frogs in their natural environment during breeding season. It is a short 2-3 week window in which they all travel and meet in a pond at high latitude. They start breeding as soon as some of the pond is not frozen so they can enter it. My first encounter with the field site was in negative temperatures and wind. I was lucky I was not the one to enter the pond in these conditions to catch the frogs.
We looked for biological differences between the Y haplotypes by measuring male fitness. We did this directly, by comparing the proportion of offspring fathered by males with the fully differentiated Y compared to the proportion of free-swimming males and, similarly, the proportion of males with a fully differentiated Y that were found in amplexus with females, compared to free-swimming males. We also measured male size to check whether different Y haplotypes are associated with body size differences.
One can never prove a negative result, and it is possible that larger sampling or a different measurement would identify a difference between the haplotypes. For example, maybe the system is unstable, and variation is only possible from gene flow which introduces lost variation from nearby populations. Or perhaps there are performance differences between the haplotypes in different environments, such as different temperatures, or in their development time, which our study was not particularly good at capturing as it would require sampling over many years.
However, there is a good reason for the standard sex chromosome evolution model not to apply to frogs, which is that sex-chromosome recombination is different than what the model assumes. Frog recombination depends on phenotypic sex and is extremely limited in males to all chromosome ends, including the sex chromosomes. Consequently, a male-determining mutation will immediately restrict recombination to almost the whole chromosome carrying the male-determining allele. This would allow it to capture a male-beneficial mutation and is likely the reason for the spread of a new male-determining mutation. However, it would make it difficult for additional male-beneficial alleles to accumulate, because they can mostly arrive by mutation. Rare recombination events are possible through rare sex-reversal cases when XY individuals develop as females and thus have a female-like chromosome-wide recombination pattern which allows recombination on the Y. The rare sex-reversed individuals, which have been observed in nature, represent a fountain of youth for the Y chromosome, and are likely the source of the semi-differentiated Y haplotypes. Nevertheless, the commonly experiences lack of recombination throughout most of the Y length means that the deleterious effects of no recombination, which are responsible for sex chromosome degeneration according to the standard model, are immediately present in young sex chromosomes in frogs.
It will be interesting to model the evolutionary dynamics of the different forces affecting sex chromosome evolution in systems with achiasmatic meiosis, or meiosis that is extremely restricted to chromosome ends, as in frogs. This type of meiotic recombination is also found in many fishes and Drosophila. It is possible that the sex chromosomes of these species are primarily affected by deleterious mutations over a larger parameter space than they are affected by sexually antagonistic mutations. In this light, it is perhaps not surprising that different frog Y haplotypes do not have obvious differences in fitness.
The published papers is available here.
Mimulus guttatus
Quantitative effects of gene expression (eQTL) and the effects of chromosomal inversions
This project was supported by a National Science Foundation grant to John Kelly.
Chromosomal inversions are known to have strong phenotypic effects, for example in determining life histories, in many organisms. However the quantitative importance of chromosomal rearrangements, which include inversions and translocations, is not known, because only rearrangements of known effects have been historically studied.
We aim to advance our understanding of the role of chromosomal rearrangements in evolution by studying genome-wide inversions at the level of the individual, the population and the species.
To do so, we are sequencing and assembling multiple genomes at varying degrees of phylogenetic distance from each other in order to characterise the accumulation of chromosomal rearrangements over evolutionary time. We are simultaneously obtaining full transcriptome data from the same individuals to quantify the chromosomal rearrangement effects in gene expression.
The organism of choice is the Mimulus species complex, because it harbours known inversions with phenotypic effects. In addition, Mimulus species have small genomes and high rates of divergence between populations and species, while retaining genetic compatibility, making them ideal to study chromosomal rearrangements.
Silene
The projects were initiated by Lynda Delph.
QTL for traits under sex specific selection in between population crosses
The level of sexual dimorphism depends on the environment where Silene latifolia grows. We are looking for QTL based on a cross between plants from dry and wet environments, with a special interest to QTL that localise on the sex chromosomes.
I am analysing ≈30 traits from a large F2 cross of plants from Spain and Croatia. Some of the traits are involved in both adaptation to a dry/wet environment and reproductive fitness, in an antagonistic way. For example, fast growth and flowering is generally beneficial for males, but deleterious in a dry environment. Consequently, it is sexually antagonistic (has different fitness optima for the two sexes) in the less stressful Croatian environment.
Such sexually antagonistic traits are predicted to occur on the recombining part of the sex chromosome, because they are expected to benefit from a non-random association with sex, while remaining capable to responding to environmental change.
We are constructing a genetic map and performing QTL analysis to test this prediction.
RNAseq of experimentally evolved plants to male and female traits
Lines with large and small flowers (good for females and males, respecively) have been generated by artificial selection. Gene expression differences between those lines should reveal the first genes that respond to such sexually antagonistic selection.
Plants from the same population have been selected for 4 generations on flower size. Both large and small flower size lines have been generated.
Flower size has different optima yet differs between the sexes, with large flowers being beneficial to females. We are investigating the changes in gene expression associated with selection in one sex at a time.
We are particularly interested in the genomic location of genes that respond to sex-specific selection.
YY Silene latifolia are inviable
No YY plants are produced when the Y is wild type, confirming the Y is significantly degenerated.
The Y chromosome of many organisms occurs as a shadow of its former self, in a degenerated state. Silene latifolia is one of the main plant species that has provided extensive evidence that the Y chromosome degenerates. Under the microscope the Y looks larger than its previously identical partner and is enriched in transposable elements and other repeats, which are also visible cytogenetically. Sequencing studies have found that about 50% of the genes on the Y have become non-functional by mutations and there have probably been two separate events that stopped the recombination between the X and the Y, that have resulted in characteristic patterns of strata of different levels of X and Y divergence.
Such events are expected to have biological consequences. One way to assess them is to generate and observe individuals with a YY genotype. This has been previously done using mutated Y chromosomes, and the effect depended on the Y variant used.
At the end of the last field season, Lynda Delph came to my office holding a dry fruit from an otherwise normal male (1/1151 flowers was female). We quickly became very excited because the offspring of this fruit come from a XY x XY cross. 25% of them should have the YY genotype, allowing us to assess the ability of a wild type Y chromosome to replace the function of the X in plant development.
We used PCR to score DNA sequences that were previously found to only occur on the X or the Y chromosome to assign the offspring of our cross into the XX, XY and YY genotypes. We did not see any YY individuals, and this was supported statistically i.e., it is very unlikely we would not see any given our sample size. Our data confirm that YY offspring do not develop. This suggests that there are genes on the X that are required for generating a normal plant, i.e. that the Y chromosome is too degenerated to support normal development.
Furthermore, it is likely that this developmental failure affects the production of ovules, because the rare fruit that we used had very few ovules, and a normal proportion of aborted seeds.
Read the full paper here.
Mercurialis annua
Almost every plant species with separate male and female individuals represents a novel evolution of sex chromosomes. The study of different plant species is therefore a great opportunity to test the theoretical predictions about how sex chromosomes evolve and differentiate.
One such prediction involves sexually antagonistic selection: genes that are beneficial for one sex but deleterious to the other should accumulate on the sex chromosomes. The reasoning is that, unlike all other chromosomes, genes on the sex chromosomes spend more time on one sex than the other: genes on the Y are always found in males, while genes on the X are found in females 66% of the time. Therefore a gene that is deleterious to females and beneficial to males will confer the benefit 100% of the time if it is on the Y, and a similar logic applies to genes with the opposite fitness effects and the X chromosome (but there it is 66% of the time, which is still better than 50%, which is the value for genes on autosomes). Such sexually antagonistic genes are expected to accumulate on the sex chromosomes and select for a reduction of recombination between them and the sex determining region. In the short term those sex chromosomes will outcompete their recombining ancestors. However in the long-term this has negative consequences as they cannot remove deleterious mutations that evolve on them.
We see evidence of the deleterious mutation accumulation in the form of genetic degeneration (gene loss, invasion of repeats, size reduction of the non recombining Y chromosome) in many species with old sex chromosomes such as mammals, birds and Drosophila. In order to study the early events in sex chromosome evolution we need a system with young sex chromosomes.
The Mercurialis annua genome
Mercurialis annua is an annual plant distributed in Europe and the middle east. It is part of a species complex with varying ploidy levels. The diploid plant has separate sexes (dioecious), and the male flowers are easily distinguished from the female flowers as they grow on long stalks (peduncles) which help with pollen spreading. Other species are monoecious (male and female flowers in the same plant), while others are androdioecious (male and monoecious plants in the same population). There are also perenial relatives which are all dioecious.
Previous crosses between the species have established that maleness follows Mendelian inheritance, suggesting there is a Y chromosome.
We used a sex-specific genetic map based on a de novo assembled transcriptome to identify the non-recombining region in the young Y chromosome of the annual mercury Mercurialis annua. The sex chromosomes are homomorphic (indistinguishable from each other) yet 1/3 of the Y is not recombining.
The non-recombining region seems special in terms of its sex biased genes: It has more female-biased genes (feminisation in terms of gene number), while simultaneously it has high male biased expression for few genes (masculinisation in terms of gene expression extent). This could be the sign of ongoing sexual conflict.
The recombination suppression seems to be recent as we only found one clear case of a non-functional Y gene copy. Comparisons with the closely related species M. huetii showed that the recombination suppression happened after the species split.
Read the full paper here.
zaba
Double Y dosage reduces fertility, but not inviability in Mercurialis annua
Plants with 2 Y chromosomes suffer infertility, but not inviability, perhaps because XX individuals that produced male flowers had an advantage in establishing populations in the past. Or YY plants have abnormal gene dosage.
The continuous growth of plants is based on embryonic tissue they contain. Unlike animals, they make developmental pathway choices as they grow, through the differentiation of pluripotent cells into the different parts of the plant. This regeneration potential allows plants to grow from cuttings and reproduce clonally. When combined with the limited differentiation between the sexes, the regeneration ability of plants makes it possible for a normal male plant to produce female flowers, occasionally or in exceptional circumstances.
Mercurialis annua is an annual herb with male and female individuals and is an interesting plant model for studying sex chromosome evolution. While it has a XY sex chromosome system, the sex chromosomes are homomorphic and of limited differentiation, making them interesting for informing early sex-chromosome evolution theory. Males can be induced to produce female flowers by spraying them with a feminising hormone, or by repeatedly cutting the top of the plant. We used both approaches to generate female flowers and studied the phenotypic effects of having 2 Y chromosomes in the resulting offspring.
Since the female flowers originated in a XY individual, their progeny represent a XY x XY cross. We developed a molecular genetic assay to distinguish the X and Y chromosomes and genotyped the offspring. Both the genotyping and the sex ratio in the first generation suggested that YY individuals were perfectly viable. We also measured various phenotypes such as height and flower mass and did not find any major differences between XY and YY males.
We might expect the XY and YY males to differ because the Y chromosome is subject to opposing evolutionary forces. On one hand it is expected to accumulate male beneficial mutations, because the genes on the Y are always present in males. We might therefore expect a super-male phenotype from YY individuals. On the other hand, the Y has degenerated somewhat, which might result in lower fitness of YY individuals, if some genes on the Y are not working well, and a X copy of a gene is required to compensate for the degeneration of Y copies. In addition, there is an unusual gene dosage of YY plants whose ratios of Y: X or Y : Autosome genes have not been tested by selection, since these plants are very rare in nature. Like most untested changes to something that works, these dosage effects would likely work less well than the progenitor state, and result in lower fitness.
Our study of the sex ratio in the offspring of the XY and YY males and direct observations of pollen from the YY plants suggest that they are infertile, at least partially. It is surprising that fertility was affected before viability, because the Y is usually considered to be specialising in male function. For example genes on the Y improving pollen fertility are likely to be under strong selection to maintain good function, which would have selected against any mutations of genes on the Y that negatively affect pollen. Nevertheless, pollen infertility could have been caused by the dosage effects mentioned above.
We also came up with a selective explanation for our data, relating to the ability of plants to produce flowers of the opposite to their genotypic sex. This is termed sex inconstancy, and it has been found in field populations of M. annua. It is possible that the X has been under strong selection to maintain genes involved in pollen function in the recent past. For example it would be an advantage for XX individuals to occasionally produce functional male flowers, because it would allow long-distance dispersing XX individuals to establish a population on their own, without relying on long dispersing pollen for fertilisation.
In fact there is a hexaploid lineage of M. annua (in the Iberian Peninsula) in which some populations have hermaphrodite flowers, and some populations have males and hermaphrodites. This situation can be explained through a metapopulation model where hermaphrodite individuals can establish populations on their own by selfing. Once they build to a higher density they can be invaded by males, who outcompete hermaphrodites in male function. Diploid inconstant male individuals might have worked in a similar fashion in the past as they are functionally hermaphrodite even though they have separate male and female flowers. For example sex inconstancy might have conferred a fitness advantage at the end of the last ice age when new locations became available for colonisation. This would have generated strong selection to maintain pollen fertility function on the X.
More experiments will be required to understand the reasons behind the poor pollen production of YY individuals. Their initially unexpected infertility and viability has opened up new understanding in how the modular function of the plant body plan may have influenced the evolution of their sex chromosomes.
Drosophila pseudoobscura
The projects were a collaboration between the labs of Mike Ritchie, Rhonda Snook and Andy Cossins.
Experimental evolution under high and low sexual selection
Sexual selection leads to the evolution of traits of one sex based on their preference by the other. It is known to lead to rapid divergence and often results in extreme phenotypes. It can be a major force in the diversification of populations and generation of genetic diversity. The effects of sexual selection on phenotypic diversity have been the easiest to study while the effects on the underlying genetics have been elusive. By studying the genetic responses to sexual selection manipulation we hope to understand their relevance to the generation of genetic diversity and, ultimately, biodiversity.
Lines of the fruit fly Drosophila pseudoobscura have been evolving in the lab under different strengths of sexual selection, controlled through sex ratio manipulation, for about 150 generations. Populations experiencing a male-biased sex ratio experience greater male competition and stronger female choice on male characteristics. Previous studies have identified various phenotypic changes that affect mating success, including both courtship behaviour, such as courtship intensity and courtship song, and reproductive physiology characters, such as accessory gland size and effects of the males on the remating interval of females.
The use of next generation sequencing allows to study the underlying genetics behind such phenotypes. We hope to obtain information on the types of changes associated with responses to sexual selection. Such information would answer exciting evolutionary questions such as:
Whether the same phenotypic response is based on the same underlying genetics when it evolves multiple times independently.
Whether the X chromosome has disproportionately more changes than other chromosomes, as suggested by theory.
Whether changes in gene expression are associated with changes in gene sequence.
Reproductive tissue mating response after experimental sexual selection
Experimentally elevated sexual selection increases male manipulation of the female post-mating response, but also makes females poised in readiness for mating and more resistant to male manipulation.
We used tissue-specific transcriptomes of reproductive tissues from both sexes. Working with specific tissues is beneficial because it minimises the effects of allometry from the observed response, it can detect gene regulation evolution even if, for example, testes have evolved to be larger under high sexual selection. One disadvantage of using specific tissues from a small insect is that there are fewer cells that provide the RNA, so the dissections required to obtain sufficient RNA took a whole week. We had 4 2-hour blocks of dissections per day.
Some dissections involved mated females that had to be prepared first. This made for a very busy week, so that the following week of (lots of) RNA extractions seemed relatively mild in comparison. With so focused physical work events blur and it is difficult to recall the experience. Nevertheless, I remember we were late for the xmas party because we wanted to finish within the year, and that we looked at an inverted microscope after tissue disruption with metal beads, to make sure we had successfully pulverised the very small reproductive tissues, including the hard spermathecae. I also have a vague memory of delivering the samples in ice in person at the Centre for Genomics Research in Liverpool, which prepared and sequenced the libraries.
All this work was before we started analysing the data, which took the most time. Here I will focus on two of our results, both of which involve interactions. But first, a brief note on the experimental design.
High sexual selection females do not require mating to activate some gene expression
Our experimental design allowed us to look for an interaction between the effect of sexual selection (the overall difference between high and low sexual selection females i.e. E vs M contrast) and the effect of mating (i.e. virgin vs mated contrast). Identifying interactions is quite a standard statistical practice, the main issue is interpreting the result and understanding the biology behind an interaction. One of the most interesting results of the paper is the strong interaction we found between mating status and sexual selection history.
To visualise the interaction, we plotted the difference in expression between virgin and mated females, separately (on the x and y axes) for the E and M females.
It is quite striking that most of the genes significant for the interaction are in one direction (green points): they are genes expressed in virgin high sexual selection females (E), but only expressed after mating in monogamous females (M).
The result supports previous phenotypic studies suggesting that the high sexual selection females are primed for mating. This might be expected, since they have evolved in an environment where they have the opportunity to mate with different males. The fact that they do not need mating to activate some genes that monogamous females only activate after mating suggests the possibility that there may be a cost to having them expressed as virgins. Their functions are related to the immune system and egg production which makes sense for the female reproductive tract and further support our interpretation that E females are poised for mating.
A similar analysis in ovaries found many genes with a significant main effect of the selection treatment, few genes changing expression based on mating, but no genes significant for an interaction, suggesting the ovariole genes of the two types of female respond the same way to mating. One explanation is that the 6 h after mating in which we sampled for RNA did not change much in terms of gene expression in the ovaries. So perhaps the main conclusion is that it was a very good call to separately analyse the two female reproductive tissues as they show very different results.
High sexual selection males manipulate the female mating response more - but high sexual selection females seem more resistant to this manipulation
We were able to study the effect of different types of male in the female mating response by analysing females of both selection treatments mated with males they had coevolved with and with males they had not. This is a complicated interaction so we plot separately the two types of female.
Comparing the number of significant points between the two graphs illustrates that E females do not change the expression of as many genes as M females do, regardless of mating partner. This also supports the suggestion that E females are primed for mating, since they activate genes prior to mating, compared to the M females.
Let us now focus on the two types of male (E or M). Notice that the same type of male has a different colour in the two plots, for example high sexual selection males affect the expression of red genes for a) but blue genes for c). Once you realise this, it is quite clear that E males cause more genes to change expression between virgin and mated M females (a - red genes, 155+103) than E females (c - blue genes 66+73).
The results lead to two intriguing interpretations:
- High sexual selection males manipulate the female postmating response more than low sexual selection males
- Simultaneously, high sexual selection females are more resistant to this manipulation than low sexual selection females.
For the full story, including results on virgin tissues (where male accessory glands and female reproductive tracts had more gene expression changes between E and M treatments than testes and ovaries) and comparisons with proteomic work, see the paper.
Drosophila montana
Selection analysis and QTL for traits associated with reproductive isolation within species in Drosophila montana
The project was a collaboration between the labs of Mike Ritchie, Anneli Hoikkala, Roger Butlin and Jon Slate.
Sexual selection leads to the evolution of traits of one sex based on their preference by the other. It is known to lead to rapid divergence and often results in extreme phenotypes. Sexual selection is a major force in the evolution of populations and generates genetic diversity. Its effects on phenotypic diversity have been the easiest to study while the effects on the underlying genetics have been elusive. In this project we hope to shed light to genetic responses to sexual selection and understand their relevance to the generation of genetic diversity and, ultimately, biodiversity.
Drosophila montana belongs to the virilis group of Drosophila. The male courtship song has been shown to be necessary for successful courtship as well as the target of sexual selection. We are performing a QTL study of parameters of the male song known to influence female choice, as well as female choice itself. We are using wild-caught populations to minimise the effects of keeping flies in the lab in the observed variation.
One long-term aim of the project is to determine whether intrapopulation variation in courtship characters is due to candidate genes, usually identified from interpopulation crosses and mutations studied in lab populations. Alternatively, some of this variation could be due to novel QTLs.
Another long-term aim of the project is to assess whether the same regions of the genome are involved in behavioural reproductive isolation a) in different populations of the same species and b) in closely related species. We hope that such comparative studies will increase our understanding on how repeatable evolution is and whether some outcomes can be predicted.
Podisma pedestris
My PhD was at Queen Mary University of London with Richard Nichols.
I worked on the Podisma pedestris hybrid zone, trying to understand the cause of hybrid inviability. Selection against hybrids is one of the two forces that maintain the hybrid zone (the other being dispersal of the pure races into the zone), therefore attributing a particular genomic sequence to it would be of great significance in understanding the zone dynamics.
Ribosomal DNA cytogenetics across a hybrid zone
A difference in rDNA sequence had previously been reported between the two sides of the hybrid zone. I took a cytogenetics approach (fluorescent in-situ hybridisation) to directly visualise the rDNA loci in various grasshopper populations. These data were combined with data at the sequence level from Irene Keller, a postdoc in the lab.
We found that rDNA does not follow concerted evolution (Keller et al., 2006) in P. pedestris. Unlike most other organisms, many copies of the rDNA sequence are differ than other copies, and many seem not be be functional. We also found very high within-species variability in the number and location of rDNA loci, though the broad distribution of loci has a biogeographic signal (Veltsos et al., 2009). We could identify 5 distinct subpopulations, based on the rDNA loci segregating in them. Two of these populations, on either side of the hybrid zone, have no loci in common, so that some of their F2 hybrid offspring would not have any rDNA loci, due to simple mendelian segregation. Nevertheless there seems to be rDNA sequence introgression throughout the hybrid zone (Keller et al., 2008), so rDNA does not seem to be the main driver of hybrid inviability.
Simulation of a sex chromosome hybrid zone
The two hybridising races have a sex chromosome difference: the ancestral race has XO males while the derived race has a Neo-XY system, caused by a chromosomal fusion of an autosome and the X chromosome. The emerging Neo-Y chromosome is restricted to males but not involved in sex determination (as is the case in Drosophila).
I studied the effect this fusion may have in the hybrid zone using computer simulations. The idea is that the sexually antagonistic effects (good for males, bad for females) that quickly arise on a Neo-Y chromosome may generate some of the hybrid inviability observed in the hybrid zone. For example the Neo-Y can enter hybrid females, which is likely to be deleterious.
The simulations show that selection on the Y chromosome affects the ancestral, unfused X chromosome. Selection against the Y in hybrid females, also selects against the unfused X (the fused X chromosome already has a homologous part to the Y, leaving no space for it), while selection for the Y in males prevents the Y from loss, as would happen in a single panmictic populations. The overall effect is that if the Y is sexually antagonistic, selection on the Neo-Y will lead to the establishment of the whole XY system throughout the species range (Veltsos et al., 2008), leaving the genetic variation on the autosomes to neutrally admix. The model provides a way for new sex chromosome systems to establish throughout a species range without the species going through an extreme bottleneck.