Micael Reis

Resolving the phylogenetic relationships and evolutionary history of the Drosophila virilis group using multilocus data

The Drosophila virilis group is one of the major lineages of Drosophila previously recognised and it has been used as a model for different types of studies. It comprises 13 species whose phylogenetic relationships are not well resolved. In the present study, six nuclear genes (Adh, fused, Gpdh, NonA, CG9631 and CG7219) and the mitochondrial ribosomal RNA genes (12S-16S) have been used to estimate the evolutionary tree of the group using different methods of phylogenetic reconstruction. Different competing evolutionary hypotheses have also been compared using the Approximately Unbiased test to further evaluate the robustness of the inferred trees. Results are, in general, consistent with previous studies in recovering the four major lineages of the group (D. virilis phylad, Drosophila montana subphylad, Drosophila kanekoi subphylad and Drosophila littoralis subphylad), although D. kanekoi, D. littoralis and Drosophila ezoana are here inferred to be more closely related to the D. virilis phylad than to the D. montana subphylad. The age of the crown group, estimated with a Bayesian method that assumes a relaxed molecular clock, is placed in the late Miocene (∼ 10 Mya). The oldest lineages also appeared during this period (∼ 7.5 to ∼ 8.9 Mya), while the ages of the basal nodes of the montana subphylad and the virilis phylad are located in the early Pliocene (∼ 4.9 and ∼ 4.1 Mya). Major cladogenesis events correlate to geological and palaeoclimatic occurrences that most likely affected the freshwater and deciduous forests where these species are found. The inferred biogeographical history of the group, based on the statistical dispersal-vicariance analysis, indicates that the last common ancestor of the group had a Holarctic distribution from which the North American and the Eurasian lineages evolved as a result of a vicariant event.

A comparative study of the short term cold resistance response in distantly related Drosophila species: the role of regucalcin and frost.

The molecular basis of short term cold resistance (indexed as chill-coma recovery time) has been mostly addressed in D. melanogaster, where candidate genes (Dca (also known as smp-30) and Frost (Fst)) have been identified. Nevertheless, in Drosophila, the ability to tolerate short term exposure to low temperatures evolved several times independently. Therefore, it is unclear whether variation in the same candidate genes is also responsible for short term cold resistance in distantly related Drosophila species. It should be noted that Dca is a candidate gene for cold resistance in the Sophophora subgenus only, since there is no orthologous gene copy in the Drosophila subgenus. Here we show that, in D. americana (Drosophila subgenus), there is a north-south gradient for a variant at the 5′ non-coding region of regucalcin (a Dca-like gene; in D. melanogaster the proteins encoded by the two genes share 71.9% amino acid identities) but in our D. americana F2 association experiment there is no association between this polymorphism and chill-coma recovery times. Moreover, we found no convincing evidence that this gene is up-regulated after cold shock in both D. americana and D. melanogaster. Size variation in the Fst PEST domain (putatively involved in rapid protein degradation) is observed when comparing distantly related Drosophila species, and is associated with short term cold resistance differences in D. americana. Nevertheless, this effect is likely through body size variation. Moreover, we show that, even at two hours after cold shock, when up-regulation of this gene is maximal in D. melanogaster (about 48 fold expression change), in D. americana this gene is only moderately up-regulated (about 3 fold expression change). Our work thus shows that there are important differences regarding the molecular basis of cold resistance in distantly related Drosophila species.

Drosophila americana as a Model Species for Comparative Studies on the Molecular Basis of Phenotypic Variation

Understanding the molecular basis of within and between species phenotypic variation is one of the main goals of Biology. In Drosophila, most of the work regarding this issue has been performed in D. melanogaster, but other distantly related species must also be studied to verify the generality of the findings obtained for this species. Here, we make the case for D. americana, a species of the virilis group of Drosophila that has been diverging from the model species, D. melanogaster, for approximately 40 Myr. To determine the suitability of this species for such studies, polymorphism and recombination estimates are presented for D. americana based on the largest nucleotide sequence polymorphism data set so far analyzed (more than 100 data sets) for this species. The polymorphism estimates are also compared with those obtained from the comparison of the genome assembly of two D. americana strains (H5 and W11) here reported. As an example of the general utility of these resources, we perform a preliminary study on the molecular basis of lifespan differences in D. americana. First, we show that there are lifespan differences between D. americana populations from different regions of the distribution range. Then, we perform five F2 association experiments using markers for 21 candidate genes previously identified in D. melanogaster. Significant associations are found between polymorphism at two genes (hep and Lim3) and lifespan. For the F2 association study involving the two sequenced strains (H5 and W11), we identify amino acid differences at Lim3 and Hep that could be responsible for the observed changes in lifespan. For both genes, no large gene expression differences were observed between the two strains.

An old bilbo-like non-LTR retroelement insertion provides insight into the relationship of species of the virilis group

In Drosophila, at the population and species level, fixation of a TE insertion is an unlikely fate. Of the few reported fixations at the species level most involve non-LTR retroelements. In this work we report the fixation of a non-LTR retroelement in five species (Drosophila littoralis, Drosophila virilis, Drosophila lummei, Drosophila americana and Drosophila novamexicana) of the virilis group of Drosophila. In most species, this TE insertion is being lost through the accumulation of small deletions, but there is also evidence for the accumulation of large deletions. In the americana lineage an insertion of about 900 bp of the non-LTR retroelement is a marker for the Xc inversion. This insertion is, at most, 80 kb away from the basal Xc inversion breakpoint. The presence of a bilbo-like element in D. littoralis but not in D. kanekoi, suggests that D. littoralis is more closely related to species of the virilis phylad than to species of the montana phylad, which is in contrast with the traditional view. Nevertheless, the phylogenetic analyses here performed using a 7 gene dataset suggest that D. littoralis is indeed more closely related to species of the virilis phylad than to species of the montana phylad. The re-evaluation of the phylogenetic relationship of the species of the virilis group, under the assumption of a relaxed molecular clock, results in an estimated age of the bilbo-like element insertion of at least 7.5 Mya.

Drosophila Genes That Affect Meiosis Duration Are among the Meiosis Related Genes That Are More Often Found Duplicated

Using a phylogenetic approach, the examination of 33 meiosis/meiosis-related genes in 12 Drosophila species, revealed nine independent gene duplications, involving the genes cav, mre11, meiS332, polo and mtrm. Evidence is provided that at least eight out of the nine gene duplicates are functional. Therefore, the rate at which Drosophila meiosis/meiosis-related genes are duplicated and retained is estimated to be 0.0012 per gene per million years, a value that is similar to the average for all Drosophila genes. It should be noted that by using a phylogenetic approach the confounding effect of concerted evolution, that is known to lead to overestimation of the duplication and retention rate, is avoided. This is an important issue, since even in our moderate size sample, evidence for long-term concerted evolution (lasting for more than 30 million years) was found for the meiS332 gene pair in species of the Drosophila subgenus. Most striking, in contrast to theoretical expectations, is the finding that genes that encode proteins that must follow a close stoichiometric balance, such as polo, mtrm and meiS332 have been found duplicated. The duplicated genes may be examples of gene neofunctionalization. It is speculated that meiosis duration may be a trait that is under selection in Drosophila and that it has different optimal values in different species.

The Drosophila melanogaster methuselah gene: a novel gene with ancient functions.

The Drosophila melanogaster G protein-coupled receptor gene, methuselah (mth), has been described as a novel gene that is less than 10 million years old. Nevertheless, it shows a highly specific expression pattern in embryos, larvae, and adults, and has been implicated in larval development, stress resistance, and in the setting of adult lifespan, among others. Although mth belongs to a gene subfamily with 16 members in D. melanogaster, there is no evidence for functional redundancy in this subfamily. Therefore, it is surprising that a novel gene influences so many traits. Here, we explore the alternative hypothesis that mth is an old gene. Under this hypothesis, in species distantly related to D. melanogaster, there should be a gene with features similar to those of mth. By performing detailed phylogenetic, synteny, protein structure, and gene expression analyses we show that the D. virilis GJ12490 gene is the orthologous of mth in species distantly related to D. melanogaster. We also show that, in D. americana (a species of the virilis group of Drosophila), a common amino acid polymorphism at the GJ12490 orthologous gene is significantly associated with developmental time, size, and lifespan differences. Our results imply that GJ12490 orthologous genes are candidates for developmental time and lifespan differences in Drosophila in general.

Drosophila americana Diapausing Females Show Features Typical of Young Flies

Diapause is a period of arrested development which is controlled physiologically, preprogrammed environmentally and characterized by metabolic depression that can occur during any stage of insect development. Nevertheless, in the genus Drosophila, diapause is almost always associated with the cessation of ovarian development and reproductive activity in adult females. In this work, we show that, in D. americana (a temperate species of the virilis group), diapause is a genetically determined delay in ovarian development that is triggered by temperature and/or photoperiod. Moreover, we show that in this species diapause incidence increases with latitude, ranging from 13% in the southernmost to 91% in the northernmost range of the distribution. When exposed to diapause inducing conditions, both diapausing and non-diapausing females show a 10% increase in lifespan, that is further increased by 18.6% in diapausing females, although senescence is far from being negligible. ActinD1 expression levels suggest that diapausing females are biologically much younger than their chronological age, and that the fly as a whole, rather than the ovarian development alone, which is phenotypically more evident, is delayed by diapause. Therefore, diapause candidate genes that show expression levels that are compatible with flies younger than their chronological age may not necessarily play a role in reproductive diapause and in adaptation to seasonally varying environmental conditions.

Niche evolution and thermal adaptation in the temperate species Drosophila americana

The study of ecological niche evolution is fundamental for understanding how the environment influences species geographical distributions and their adaptation to divergent environments. Here, we present a study of the ecological niche, demographic history and thermal performance (locomotor activity, developmental time and fertility/viability) of the temperate species Drosophila americana and its two chromosomal forms. Temperature is the environmental factor that contributes most to the species and chromosomal forms ecological niches, although precipitation is also important in the model of the southern populations. The past distribution model of the species predicts a drastic reduction in the suitable area for the distribution of the species during the last glacial maximum (LGM), suggesting a strong bottleneck. However, DNA analyses did not detect a bottleneck signature during the LGM. These contrasting results could indicate that D. americana niche preference evolves with environmental change, and thus, there is no evidence to support niche conservatism in this species. Thermal performance experiments show no difference in the locomotor activity across a temperature range of 15 to 38 °C between flies from the north and the south of its distribution. However, we found significant differences in developmental time and fertility/viability between the two chromosomal forms at the models optimal temperatures for the two forms. However, results do not indicate that they perform better for the traits studied here in their respective optimal niche temperatures. This suggests that behaviour plays an important role in thermoregulation, supporting the capacity of this species to adapt to different climatic conditions across its latitudinal distribution.

The Drosophila melanogaster Muc68E Mucin Gene Influences Adult Size, Starvation Tolerance, and Cold Recovery

Mucins have been implicated in many different biological processes, such as protection from mechanical damage, microorganisms, and toxic molecules, as well as providing a luminal scaffold during development. Nevertheless, it is conceivable that mucins have the potential to modulate food absorption as well, and thus contribute to the definition of several important phenotypic traits. Here we show that the Drosophila melanogaster Muc68E gene is 40- to 60-million-yr old, and is present in Drosophila species of the subgenus Sophophora only. The central repeat region of this gene is fast evolving, and shows evidence for repeated expansions/contractions. This and/or frequent gene conversion events lead to the homogenization of its repeats. The amino acid pattern P[ED][ED][ST][ST][ST] is found in the repeat region of Muc68E proteins from all Drosophila species studied, and can occur multiple times within a single conserved repeat block, and thus may have functional significance. Muc68E is a nonessential gene under laboratory conditions, but Muc68E mutant flies are smaller and lighter than controls at birth. However, at 4 d of age, Muc68E mutants are heavier, recover faster from chill-coma, and are more resistant to starvation than control flies, although they have the same percentage of lipids as controls. Mutant flies have enlarged abdominal size 1 d after chill-coma recovery, which is associated with higher lipid content. These results suggest that Muc68E has a role in metabolism modulation, food absorption, and/or feeding patterns in larvae and adults, and under normal and stress conditions. Such biological function is novel for mucin genes.

Genes Belonging to the Insulin and Ecdysone Signaling Pathways Can Contribute to Developmental Time, Lifespan and Abdominal Size Variation in Drosophila americana

Even within a single genus, such as Drosophila, cases of lineage-specific adaptive evolution have been found. Therefore, the molecular basis of phenotypic variation must be addressed in more than one species group, in order to infer general patterns. In this work, we used D. americana, a species distantly-related to D. melanogaster, to perform an F2 association study for developmental time (DT), chill-coma recovery time (CRT), abdominal size (AS) and lifespan (LS) involving the two strains (H5 and W11) whose genomes have been previously sequenced. Significant associations were found between the 43 large indel markers developed here and DT, AS and LS but not with CRT. Significant correlations are also found between DT and LS, and between AS and LS, that might be explained by variation at genes belonging to the insulin and ecdysone signaling pathways. Since, in this F2 association study a single marker, located close to the Ecdysone receptor (EcR) gene, explained as much as 32.6% of the total variation in DT, we performed a second F2 association study, to determine whether large differences in DT are always due to variation in this genome region. No overlapping signal was observed between the two F2 association studies. Overall, these results illustrate that, in D. americana, pleiotropic genes involved in the highly-conserved insulin and ecdysone signaling pathways are likely responsible for variation observed in ecologically relevant phenotypic traits, although other genes are also involved.