Thermal evolution of gene expression profiles in Drosophila subobscura

Background: Despite its pervasiveness, the genetic basis of adaptation resulting in variation directly or indirectly related to temperature (climatic) gradients is poorly understood. By using 3-fold replicated laboratory thermal stocks covering much of the physiologically tolerable temperature range for the temperate (i.e., cold tolerant) species Drosophila subobscura we have assessed whole-genome transcriptional responses after three years of thermal adaptation, when the populations had already diverged for inversion frequencies, pre-adult life history components, and morphological traits. Total mRNA from each population was compared to a reference pool mRNA in a standard, highly replicated two-colour competitive hybridization experiment using cDNA microarrays.Results: A total of 306 (6.6%) cDNA clones were identified as 'differentially expressed' (following a false discovery rate correction) after contrasting the two furthest apart thermal selection regimes (i.e., 13°C vs . 22°C), also including four previously reported candidate genes for thermotolerance in Drosophila (Hsp26, Hsp68, Fst, and Treh). On the other hand, correlated patterns of gene expression were similar in cold- and warm-adapted populations. Analysis of functional categories defined by the Gene Ontology project point to an overrepresentation of genes involved in carbohydrate metabolism, nucleic acids metabolism and regulation of transcription among other categories. Although the location of differently expressed genes was approximately at random with respect to chromosomes, a physical mapping of 88 probes to the polytene chromosomes of D. subobscura has shown that a larger than expected number mapped inside inverted chromosomal segments.Conclusion: Our data suggest that a sizeable number of genes appear to be involved in thermal adaptation in Drosophila, with a substantial fraction implicated in metabolism. This apparently illustrates the formidable challenge to understanding the adaptive evolution of complex trait variation. Furthermore, some clustering of genes within inverted chromosomal sections was detected. Disentangling the effects of inversions will be obviously required in any future approach if we want to identify the relevant candidate genes.

Precision and scaling in morphogen gradient read-out.

Morphogen gradients infer cell fate as a function of cellular position. Experiments in Drosophila embryos have shown that the Bicoid (Bcd) gradient is precise and exhibits some degree of scaling. We present experimental results on the precision of Bcd target genes for embryos with a single, double or quadruple dose of bicoid demonstrating that precision is highest at mid-embryo and position dependent, rather than gene dependent. This confirms that the major contribution to precision is achieved already at the Bcd gradient formation. Modeling this dynamic process, we investigate precision for inter-embryo fluctuations in different parameters affecting gradient formation. Within our modeling framework, the observed precision can only be achieved by a transient Bcd profile. Studying different extensions of our modeling framework reveals that scaling is generally position dependent and decreases toward the posterior pole. Our measurements confirm this trend, indicating almost perfect scaling except for anterior most expression domains, which overcompensate fluctuations in embryo length.

Short- and long-term changes in chromosomal inversion polymorphism and global warming: Drosophila subobscura from the Balkans

The chromosomal inversion polymorphism of Drosophila subobscura is adaptive to environmental changes. The population of Petnica, Serbia, was chosen to analyze short- and long-term changes in this polymorphism. Short-term changes were studied in the samples collected in May, June, and August of 1995. The inversion polymorphism varied over these months, although various interpretations are possible. To analyze long-term changes, samples obtained in May 1995 and May 2010 were compared. The frequency of the 'cold' adapted inversions (Ast, Jst, Ust, Est, and Ost) decreased and that of the 'warm' adapted inversions (A2, J1, U1+2, and O3+4) increased, from 1995 to 2010. These changes are consistent with the general increase in temperature recorded in Petnica for the same period. Finally, the possible response of chromosomal polymorphism to global warming was analyzed at the regional level (Balkan peninsula). This polymorphism depends on the ecological conditions of the populations, and the changes observed appear to be consistent with global warming expectations. Natural selection seems to be the main mechanism responsible for the evolution of this chromosomal polymorphism.

Short- and long-term changes in chromosomal inversion polymorphism and global warming: Drosophila subobscura from the Balkans

The chromosomal inversion polymorphism of Drosophila subobscura is adaptive to environmental changes. The population of Petnica, Serbia, was chosen to analyze short- and long-term changes in this polymorphism. Short-term changes were studied in the samples collected in May, June, and August of 1995. The inversion polymorphism varied over these months, although various interpretations are possible. To analyze long-term changes, samples obtained in May 1995 and May 2010 were compared. The frequency of the 'cold' adapted inversions (Ast, Jst, Ust, Est, and Ost) decreased and that of the 'warm' adapted inversions (A2, J1, U1+2, and O3+4) increased, from 1995 to 2010. These changes are consistent with the general increase in temperature recorded in Petnica for the same period. Finally, the possible response of chromosomal polymorphism to global warming was analyzed at the regional level (Balkan peninsula). This polymorphism depends on the ecological conditions of the populations, and the changes observed appear to be consistent with global warming expectations. Natural selection seems to be the main mechanism responsible for the evolution of this chromosomal polymorphism.

Bacterial diversity of soil under eucalyptus assessed by 16S rDNA sequencing analysis

Studies on the impact of Eucalyptus spp. on Brazilian soils have focused on soil chemical properties and isolating interesting microbial organisms. Few studies have focused on microbial diversity and ecology in Brazil due to limited coverage of traditional cultivation and isolation methods. Molecular microbial ecology methods based on PCR amplified 16S rDNA have enriched the knowledge of soils microbial biodiversity. The objective of this work was to compare and estimate the bacterial diversity of sympatric communities within soils from two areas, a native forest (NFA) and an eucalyptus arboretum (EAA). PCR primers, whose target soil metagenomic 16S rDNA were used to amplify soil DNA, were cloned using pGEM-T and sequenced to determine bacterial diversity. From the NFA soil 134 clones were analyzed, while 116 clones were analyzed from the EAA soil samples. The sequences were compared with those online at the GenBank. Phylogenetic analyses revealed differences between the soil types and high diversity in both communities. Soil from the Eucalyptus spp. arboretum was found to have a greater bacterial diversity than the soil investigated from the native forest area.

Robustness of Drosophila early embryo and wing imaginal disc development

Within a developing organism, cells require information on where they are in order to differentiate into the correct cell-type. Pattern formation is the process by which cells acquire and process positional cues and thus determine their fate. This can be achieved by the production and release of a diffusible signaling molecule, called a morphogen, which forms a concentration gradient: exposure to different morphogen levels leads to the activation of specific signaling pathways. Thus, in response to the morphogen gradient, cells start to express different sets of genes, forming domains characterized by a unique combination of differentially expressed genes. As a result, a pattern of cell fates and specification emerges.Though morphogens have been known for decades, it is not yet clear how these gradients form and are interpreted in order to yield highly robust patterns of gene expression. During my PhD thesis, I investigated the properties of Bicoid (Bcd) and Decapentaplegic (Dpp), two morphogens involved in the patterning of the anterior-posterior axis of Drosophila embryo and wing primordium, respectively. In particular, I have been interested in understanding how the pattern proportions are maintained across embryos of different sizes or within a growing tissue. This property is commonly referred to as scaling and is essential for yielding functional organs or organisms. In order to tackle these questions, I analysed fluorescence images showing the pattern of gene expression domains in the early embryo and wing imaginal disc. After characterizing the extent of these domains in a quantitative and systematic manner, I introduced and applied a new scaling measure in order to assess how well proportions are maintained. I found that scaling emerged as a universal property both in early embryos (at least far away from the Bcd source) and in wing imaginal discs (across different developmental stages). Since we were also interested in understanding the mechanisms underlying scaling and how it is transmitted from the morphogen to the target genes down in the signaling cascade, I also quantified scaling in mutant flies where this property could be disrupted. While scaling is largely conserved in embryos with altered bcd dosage, my modeling suggests that Bcd trapping by the nuclei as well as pre-steady state decoding of the morphogen gradient are essential to ensure precise and scaled patterning of the Bcd signaling cascade. In the wing imaginal disc, it appears that as the disc grows, the Dpp response expands and scales with the tissue size. Interestingly, scaling is not perfect at all positions in the field. The scaling of the target gene domains is best where they have a function; Spalt, for example, scales best at the position in the anterior compartment where it helps to form one of the anterior veins of the wing. Analysis of mutants for pentagone, a transcriptional target of Dpp that encodes a secreted feedback regulator of the pathway, indicates that Pentagone plays a key role in scaling the Dpp gradient activity.

Idiosyncratic evolution of maternal effects in response to juvenile malnutrition in Drosophila.

Maternal effects often affect fitness traits, but there is little experimental evidence pertaining to their contribution to response to selection imposed by novel environments. We studied the evolution of maternal effects in Drosophila populations selected for tolerance to chronic larval malnutrition. To this end, we performed pairwise reciprocal F1 crosses between six selected (malnutrition tolerant) populations and six unselected control populations and assessed the effect of cross direction on larval growth and developmental rate, adult weight and egg-to-adult viability expressed under the malnutrition regime. Each pair of reciprocal crosses revealed large maternal effects (possibly including cytoplasmic genetic effects) on at least one trait, but the magnitude, sign and which traits were affected varied among populations. Thus, maternal effects contributed significantly to the response to selection imposed by the malnutrition regime, but these changes were idiosyncratic, suggesting a rugged adaptive landscape. Furthermore, although the selected populations evolved both faster growth and higher viability, the maternal effects on growth rate and viability were negatively correlated across populations. Thus, genes mediating maternal effects can evolve to partially counteract the response to selection mediated by the effects of alleles on their own carriers' phenotype, and maternal effects may contribute to evolutionary trade-offs between components of offspring fitness.

Vanishing chromosomal inversion clines in Drosophila subobscura from Chile: Is behavioral thermoregulation to blame?

Chromosomal inversion clines paralleling the long-standing ones in native Palearctic populations of Drosophila subobscura evolved swiftly after this species invaded the Americas in the late 1970s and early 1980s. However, the new clines did not consistently continue to converge on the Old World baseline. Our recent survey of Chilean populations of D. subobscura shows that inversion clines have faded or even changed sign with latitude. Here, we investigate the hypothesis that this fading of inversion clines might be due to the Bogert effect, namely, that flies' thermoregulatory behavior has eventually compensated for environmental variation in temperature, thus buffering selection on thermal-related traits. We show that latitudinal divergence in thermal preference (T-p) has evolved in Chile for females, with higher-latitude flies having a lower mean T-p. Plastic responses in T-p also lessen latitudinal thermal variation because flies developed at colder temperatures prefer warmer microclimates. Our results are consistent with the idea that active behavioral thermoregulation might buffer environmental variation and reduce the potential effect of thermal selection on other traits as chromosomal arrangements.

Gene flow and gene flux shape evolutionary patterns of variation in Drosophila subobscura

Gene flow (defined as allele exchange between populations) and gene flux (defined as allele exchange during meiosis in heterokaryotypic females) are important factors decreasing genetic differentiation between populations and inversions. Many chromosomal inversions are under strong selection and their role in recombination reduction enhances the maintenance of their genetic distinctness. Here we analyze levels and patterns of nucleotide diversity, selection and demographic history, using 37 individuals of Drosophila subobscura from Mount Parnes (Greece) and Barcelona (Spain). Our sampling focused on two frequent O-chromosome arrangements that differ by two overlapping inversions (OST and O3+4), which are differentially adapted to the environment as observed by their opposing latitudinal clines in inversion frequencies. The six analyzed genes (Pif1A, Abi, Sqd, Yrt, Atpa and Fmr1) were selected for their location across the O-chromosome and their implication in thermal adaptation. Despite the extensive gene flux detected outside the inverted region, significant genetic differentiation between both arrangements was found inside it. However, high levels of gene flow were detected for all six genes when comparing the same arrangement among populations. These results suggest that the adaptive value of inversions is maintained, regardless of the lack of genetic differentiation within arrangements from different populations, and thus favors the Local Adaptation hypothesis over the Coadapted Genome hypothesis as the basis of the selection acting on inversions in these populations.

Inversion polymorphism in two Serbian natural populations of Drosophila subobscura: Analysis of long-term changes

To study whether inversions (or arrangements) by themselves or karyotypes are the main global warming adaptive target of natural selection, two Drosophila subobscura Serbian populations (Apatin and Petnica) were re analyzed using different statistical approaches. Both populations were sampled in an approximately 15 years period: Apatin in 1994 and 2008 + 2009 and Petnica in 1995 and 2010. For all chromosomes, the four collections studied were in Hardy-Weinberg equilibrium. Thus, it seems that inversions (or arrangements) combined at random to constitute populations" karyotypes. However, there were differences in karyotypic fre quencies along the years, although they were significant only for Apatin population. It is possible to conclude that inversions (or arrangements) are likely the target of natural selection, because they presented long term changes, but combine at random to generate the corresponding karyotypic combinations. As a consequence, the frequencies of karyotypes also change along time.

Family Size Evolution in Drosophila Chemosensory Gene Families: A Comparative Analysis with a Critical Appraisal of Methods

Gene turnover rates and the evolution of gene family sizes are important aspects of genome evolution. Here, we use curated sequence data of the major chemosensory gene families from Drosophila-the gustatory receptor, odorant receptor, ionotropic receptor, and odorant-binding protein families-to conduct a comparative analysis among families, exploring different methods to estimate gene birth and death rates, including an ad hoc simulation study. Remarkably, we found that the state-of-the-art methods may produce very different rate estimates, which may lead to disparate conclusions regarding the evolution of chemosensory gene family sizes in Drosophila. Among biological factors, we found that a peculiarity of D. sechellia's gene turnover rates was a major source of bias in global estimates, whereas gene conversion had negligible effects for the families analyzed herein. Turnover rates vary considerably among families, subfamilies, and ortholog groups although all analyzed families were quite dynamic in terms of gene turnover. Computer simulations showed that the methods that use ortholog group information appear to be the most accurate for the Drosophila chemosensory families. Most importantly, these results reveal the potential of rate heterogeneity among lineages to severely bias some turnover rate estimation methods and the need of further evaluating the performance of these methods in a more diverse sampling of gene families and phylogenetic contexts. Using branch-specific codon substitution models, we find further evidence of positive selection in recently duplicated genes, which attests to a nonneutral aspect of the gene birth-and-death process.

Gene expression following induction of regeneration in Drosophila wing imaginal discs

Background: Regeneration is the ability of an organism to rebuild a body part that has been damaged or amputated, and can be studied at the molecular level using model organisms. Drosophila imaginal discs, which are the larval primordia of adult cuticular structures, are capable of undergoing regenerative growth after transplantation and in vivo culture into the adult abdomen. Results: Using expression profile analyses, we studied the regenerative behaviour of wing discs at 0, 24 and 72 hours after fragmentation and implantation into adult females. Based on expression level, we generated a catalogue of genes with putative role in wing disc regeneration, identifying four classes: 1) genes with differential expression within the first 24 hours; 2) genes with differential expression between 24 and 72 hours; 3) genes that changed significantly in expression levels between the two time periods; 4) genes with a sustained increase or decrease in their expression levels throughout regeneration. Among these genes, we identified members of the JNK and Notch signalling pathways and chromatin regulators. Through computational analysis, we recognized putative binding sites for transcription factors downstream of these pathways that are conserved in multiple Drosophilids, indicating a potential relationship between members of the different gene classes. Experimental data from genetic mutants provide evidence of a requirement of selected genes in wing disc regeneration. Conclusions: We have been able to distinguish various classes of genes involved in early and late steps of the regeneration process. Our data suggests the integration of signalling pathways in the promoters of regulated genes.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background: The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians. Results: In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern. Conclusions: The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.