Tempo and Mode of Transposable Element Activity in Drosophila

The evolutionary dynamics of transposable element (TE) insertions have been of continued interest since TE activity has important implications for genome evolution and adaptation. Here, we infer the transposition dynamics of TEs by comparing their abundance in natural D. melanogaster and D. simulans populations. Sequencing pools of more than 550 South African flies to at least 320-fold coverage, we determined the genome wide TE insertion frequencies in both species. We suggest that the predominance of low frequency insertions in the two species (>80% of the insertions have a frequency <0.2) is probably due to a high activity of more than 58 families in both species. We provide evidence for 50% of the TE families having temporally heterogenous transposition rates with different TE families being affected in the two species. While in D. melanogaster retrotransposons were more active, DNA transposons showed higher activity levels in D. simulans. Moreover, we suggest that LTR insertions are mostly of recent origin in both species, while DNA and non-LTR insertions are older and more frequently vertically transmitted since the split of D. melanogaster and D. simulans. We propose that the high TE activity is of recent origin in both species and a consequence of the demographic history, with habitat expansion triggering a period of rapid evolution.     

The Drosophila Gene for Antizyme Requires Ribosomal Frameshifting for Expression and Contains an Intronic Gene for snRNP Sm D3 on the Opposite Strand

Previously, a Drosophila melanogaster sequence with high homology to the sequence for mammalian antizyme (ornithine decarboxylase antizyme) was reported. The present study shows that homology of this coding sequence to its mammalian antizyme counterpart also extends to a 5′ open reading frame (ORF) which encodes the amino-terminal part of antizyme and overlaps the +1 frame (ORF2) that encodes the carboxy-terminal three-quarters of the protein. Ribosomes shift frame from the 5′ ORF to ORF2 with an efficiency regulated by polyamines. At least in mammals, this is part of an autoregulatory circuit. The shift site and 23 of 25 of the flanking nucleotides which are likely important for efficient frameshifting are identical to their mammalian homologs. In the reverse orientation, within one of the introns of the Drosophila antizyme gene, the gene for snRNP Sm D3 is located. Previously, it was shown that two closely linked P-element transposon insertions caused the gutfeeling phenotype of embryonic lethality and aberrant neuronal and muscle cell differentiation. The present work shows that defects in either snRNP Sm D3 or antizyme, or both, are likely causes of the phenotype.     

Drosophila Embryogenesis Scales Uniformly across Temperature in Developmentally Diverse Species

Temperature affects both the timing and outcome of animal development, but the detailed effects of temperature on the progress of early development have been poorly characterized. To determine the impact of temperature on the order and timing of events during Drosophila melanogaster embryogenesis, we used time-lapse imaging to track the progress of embryos from shortly after egg laying through hatching at seven precisely maintained temperatures between 17.5°C and 32.5°C. We employed a combination of automated and manual annotation to determine when 36 milestones occurred in each embryo. D. melanogaster embryogenesis takes 33 hours at 17.5°C, and accelerates with increasing temperature to a low of 16 hours at 27.5°C, above which embryogenesis slows slightly. Remarkably, while the total time of embryogenesis varies over two fold, the relative timing of events from cellularization through hatching is constant across temperatures. To further explore the relationship between temperature and embryogenesis, we expanded our analysis to cover ten additional Drosophila species of varying climatic origins. Six of these species, like D. melanogaster, are of tropical origin, and embryogenesis time at different temperatures was similar for them all. D. mojavensis, a sub-tropical fly, develops slower than the tropical species at lower temperatures, while D. virilis, a temperate fly, exhibits slower development at all temperatures. The alpine sister species D. persimilis and D. pseudoobscura develop as rapidly as tropical flies at cooler temperatures, but exhibit diminished acceleration above 22.5°C and have drastically slowed development by 30°C. Despite ranging from 13 hours for D. erecta at 30°C to 46 hours for D. virilis at 17.5°C, the relative timing of events from cellularization through hatching is constant across all species and temperatures examined here, suggesting the existence of a previously unrecognized timer controlling the progress of embryogenesis that has been tuned by natural selection as each species diverges.     

Population-Genomic Analysis Identifies a Low Rate of Global Adaptive Fixation in the Proteins of the Cyclical Parthenogen Daphnia magna

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared with Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative ‘arms race’ genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction.     

Drosophila Interspecific Hybrids Phenocopy piRNA-Pathway Mutants

The Piwi-interacting RNA (piRNA) pathway defends the germline of animals from the deleterious activity of selfish transposable elements (TEs) through small-RNA mediated silencing. Adaptation to novel invasive TEs is proposed to occur by incorporating their sequences into the piRNA pool that females produce and deposit into their eggs, which then propagates immunity against specific TEs to future generations. In support of this model, the F1 offspring of crosses between strains of the same Drosophila species sometimes suffer from germline derepression of paternally inherited TE families, caused by a failure of the maternal strain to produce the piRNAs necessary for their regulation. However, many protein components of the Drosophila piRNA pathway exhibit signatures of positive selection, suggesting that they also contribute to the evolution of host genome defense. Here we investigate piRNA pathway function and TE regulation in the F1 hybrids of interspecific crosses between D. melanogaster and D. simulans and compare them with intraspecific control crosses of D. melanogaster. We confirm previous reports showing that intraspecific crosses are characterized by derepression of paternally inherited TE families that are rare or absent from the maternal genome and piRNA pool, consistent with the role of maternally deposited piRNAs in shaping TE silencing. In contrast to the intraspecific cross, we discover that interspecific hybrids are characterized by widespread derepression of both maternally and paternally inherited TE families. Furthermore, the pattern of derepression of TE families in interspecific hybrids cannot be attributed to their paucity or absence from the piRNA pool of the maternal species. Rather, we demonstrate that interspecific hybrids closely resemble piRNA effector-protein mutants in both TE misregulation and aberrant piRNA production. We suggest that TE derepression in interspecific hybrids largely reflects adaptive divergence of piRNA pathway genes rather than species-specific differences in TE-derived piRNAs.     

An Age-of-Allele Test of Neutrality for Transposable Element Insertions

How natural selection acts to limit the proliferation of transposable elements (TEs) in genomes has been of interest to evolutionary biologists for many years. To describe TE dynamics in populations, previous studies have used models of transposition–selection equilibrium that assume a constant rate of transposition. However, since TE invasions are known to happen in bursts through time, this assumption may not be reasonable. Here we propose a test of neutrality for TE insertions that does not rely on the assumption of a constant transposition rate. We consider the case of TE insertions that have been ascertained from a single haploid reference genome sequence. By conditioning on the age of an individual TE insertion allele (inferred by the number of unique substitutions that have occurred within the particular TE sequence since insertion), we determine the probability distribution of the insertion allele frequency in a population sample under neutrality. Taking models of varying population size into account, we then evaluate predictions of our model against allele frequency data from 190 retrotransposon insertions sampled from North American and African populations of Drosophila melanogaster. Using this nonequilibrium neutral model, we are able to explain ∼80% of the variance in TE insertion allele frequencies based on age alone. Controlling for both nonequilibrium dynamics of transposition and host demography, we provide evidence for negative selection acting against most TEs as well as for positive selection acting on a small subset of TEs. Our work establishes a new framework for the analysis of the evolutionary forces governing large insertion mutations like TEs, gene duplications, or other copy number variants.     

Behavioral and Genomic Sensory Adaptations Underlying the Pest Activity of Drosophila suzukii

Studying how novel phenotypes originate and evolve is fundamental to the field of evolutionary biology as it allows us to understand how organismal diversity is generated and maintained. However, determining the basis of novel phenotypes is challenging as it involves orchestrated changes at multiple biological levels. Here, we aim to overcome this challenge by using a comparative species framework combining behavioral, gene expression, and genomic analyses to understand the evolutionary novel egg-laying substrate-choice behavior of the invasive pest species Drosophila suzukii. First, we used egg-laying behavioral assays to understand the evolution of ripe fruit oviposition preference in D. suzukii compared with closely related species D. subpulchrella and D. biarmipes as well as D. melanogaster. We show that D. subpulchrella and D. biarmipes lay eggs on both ripe and rotten fruits, suggesting that the transition to ripe fruit preference was gradual. Second, using two-choice oviposition assays, we studied how D. suzukii, D. subpulchrella, D. biarmipes, and D. melanogaster differentially process key sensory cues distinguishing ripe from rotten fruit during egg-laying. We found that D. suzukii’s preference for ripe fruit is in part mediated through a species-specific preference for stiff substrates. Last, we sequenced and annotated a high-quality genome for D. subpulchrella. Using comparative genomic approaches, we identified candidate genes involved in D. suzukii’s ability to seek out and target ripe fruits. Our results provide detail to the stepwise evolution of pest activity in D. suzukii, indicating important cues used by this species when finding a host, and the molecular mechanisms potentially underlying their adaptation to a new ecological niche.     

High genetic and epigenetic variation of transposable elements: Potential drivers to rapid adaptive evolution for the noxious invasive weed Mikania micrantha

Why invasive species can rapidly adapt to novel environments is a puzzling question known as the genetic paradox of invasive species. This paradox is explainable in terms of transposable elements (TEs) activity, which are theorized to be powerful mutational forces to create genetic variation. Mikania micrantha, a noxious invasive weed, in this sense provides an excellent opportunity to test the explanation. The genetic and epigenetic variation of 21 invasive populations of M. micrantha in southern China have been examined by using transposon display (TD) and transposon methylation display (TMD) techniques to survey 12 TE superfamilies. Our results showed that M. micrantha populations maintained an almost equally high level of TE‐based genetic and epigenetic variation and they have been differentiated into subpopulations genetically and epigenetically. A similar positive spatial genetic and epigenetic structure pattern was observed within 300 m. Six and seven TE superfamilies presented significant genetic and epigenetic isolation by distance (IBD) pattern. In total, 59 genetic and 86 epigenetic adaptive TE loci were identified. Of them, 51 genetic and 44 epigenetic loci were found to correlate with 25 environmental variables (including precipitation, temperature, vegetation coverage, and soil metals). Twenty‐five transposon‐inserted genes were sequenced and homology‐based annotated, which are found to be involved in a variety of molecular and cellular functions. Our research consolidates the importance of TE‐associated genetic and epigenetic variation in the rapid adaptation and invasion of M. micrantha.     

Endogenous siRNAs and piRNAs derived from transposable elements and genes in the malaria vector mosquito Anopheles gambiae

Background

The siRNA and piRNA pathways have been shown in insects to be essential for regulation of gene expression and defence against exogenous and endogenous genetic elements (viruses and transposable elements). The vast majority of endogenous small RNAs produced by the siRNA and piRNA pathways originate from repetitive or transposable elements (TE). In D. melanogaster, TE-derived endogenous siRNAs and piRNAs are involved in genome surveillance and maintenance of genome integrity. In the medically relevant malaria mosquito Anopheles gambiae TEs constitute 12-16% of the genome size. Genetic variations induced by TE activities are known to shape the genome landscape and to alter the fitness in An. gambiae.

Results

Here, using bioinformatics approaches we analyzed the small RNA data sets from 6 libraries formally reported in a previous study and examined the expression of the mixed germline/somatic siRNAs and piRNAs produced in adult An. gambiae females. We characterized a large population of TE-derived endogenous siRNAs and piRNAs, which constitutes 56-60% of the total siRNA and piRNA reads in the analysed libraries. Moreover, we identified a number of protein coding genes producing gene-specific siRNAs and piRNAs that were generally expressed at much lower levels than the TE-associated small RNAs. Detailed sequence analysis revealed that An. gambiae piRNAs were produced by both “ping-pong” dependent (TE-associated piRNAs) and independent mechanisms (genic piRNAs). Similarly to D. melanogaster, more than 90% of the detected piRNAs were produced from TE-associated clusters in An. gambiae. We also found that biotic stress as blood feeding and infection with Plasmodium parasite, the etiological agent of malaria, modulated the expression levels of the endogenous siRNAs and piRNAs in An. gambiae.

Conclusions

We identified a large and diverse set of the endogenously derived siRNAs and piRNAs that share common and distinct aspects of small RNA expression across insect species, and inferred their impact on TE and gene activity in An. gambiae. The TE-specific small RNAs produced by both the siRNA and piRNA pathways represent an important aspect of genome stability and genetic variation, which might have a strong impact on the evolution of the genome and vector competence in the malaria mosquitoes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1436-1) contains supplementary material, which is available to authorized users.     

Testing the palindromic target site model for DNA transposon insertion using the Drosophila melanogaster P-element

Understanding the molecular mechanisms that influence transposable element target site preferences is a fundamental challenge in functional and evolutionary genomics. Large-scale transposon insertion projects provide excellent material to study target site preferences in the absence of confounding effects of post-insertion evolutionary change. Growing evidence from a wide variety of prokaryotes and eukaryotes indicates that DNA transposons recognize staggered-cut palindromic target site motifs (TSMs). Here, we use over 10 000 accurately mapped P-element insertions in the Drosophila melanogaster genome to test predictions of the staggered-cut palindromic target site model for DNA transposon insertion. We provide evidence that the P-element targets a 14-bp palindromic motif that can be identified at the primary sequence level, which predicts the local spacing, hotspots and strand orientation of P-element insertions. Intriguingly, we find that the although P-element destroys the complete 14-bp target site upon insertion, the terminal three nucleotides of the P-element inverted repeats complement and restore the original TSM, suggesting a mechanistic link between transposon target sites and their terminal inverted repeats. Finally, we discuss how the staggered-cut palindromic target site model can be used to assess the accuracy of genome mappings for annotated P-element insertions.     

The P-M Hybrid Dysgenesis Cline in Eastern Australian Drosophila melanogaster: Discrete P, Q and M Regions Are Nearly Contiguous

The dramatic latitudinal cline in P-M hybrid dysgenesis characteristics along the east coast of Australia is not smooth. Tests of recent collections of Drosophila melanogaster from the southeastern coast define the previously described cline as comprising three discrete, apparently contiguous regions of P, Q and M phenotypes, respectively. Northern populations from Cairns (16.9°SLat) to Ourimbah (33.4°SLat) are phenotypically P; populations from Wollongong (34.4°SLat) to Eden (37.1°SLat) are Q; and populations from Genoa (37.5°SLat) to Cygnet (43.2°SLat) are M. The decline in P activity from northern Queensland (55-60% gonadal dysgenesis (GD) in cross A) to mid-New South Wales (20-30% GD in cross A) is gradual; proceeding south, there then is a sharp drop to Q populations (<10% GD in crosses A and A*). This drop in P activity occurs in only 150 km, across the urban and suburban area of Sydney. Q populations are then found south to Eden, but Genoa, only about 50 km further southeast, is clearly M (48% GD in cross A*), as are two populations further south. The two discontinuities in the P-M cline do not correspond to obvious climatic differences along the coast, nor to obvious barriers to dispersal of D. melanogaster. The cline has apparently not moved between 1983 and 1985-1986.     

A scenario for the hobo transposable element invasion, deduced from the structure of natural populations of Drosophila melanogaster using tandem TPE repeats

Temporal surveys of hobo transposable elements in natural populations reveal a historical pattern suggesting a recent world-wide invasion of D. melanogaster by these transposons, perhaps following a recent horizontal transfer. To clarify the dynamics of hobo elements in natural populations, and thus to provide further data for our understanding of the hobo invasion, TPE tandem repeats, observed in the polymorphic S region of the element, were used as molecular markers. The number of TPE repeats was studied in 101 current populations from around the world, and in 63 strains collected in the past. This revealed a geographical distribution which seems to have been stable since the beginning of the 1960s. This distribution is compatible with a number of hypotheses for the dynamics of hobo elements. We propose a scenario based on an invasion in two stages: first, a complete invasion by elements with three TPE repeats, followed by the beginning of a new invasion involving hobo elements with five or seven repeats.     

Gene-by-Temperature Interactions and Candidate Plasticity Genes for Morphological Traits in Drosophila melanogaster

Understanding the genetic architecture of any quantitative trait requires identifying the genes involved in its expression in different environmental conditions. This goal can be achieved by mutagenesis screens in genetically tractable model organisms such as Drosophila melanogaster. Temperature during ontogenesis is an important environmental factor affecting development and phenotypic variation in holometabolous insects. In spite of the importance of phenotypic plasticity and genotype by environment interaction (GEI) for fitness related traits, its genetic basis has remained elusive. In this context, we analyzed five different adult morphological traits (face width, head width, thorax length, wing size and wing shape) in 42 co-isogenic single P-element insertional lines of Drosophila melanogaster raised at 17°C and 25°C. Our analyses showed that all lines differed from the control for at least one trait in males or females at either temperature. However, no line showed those differences for all traits in both sexes and temperatures simultaneously. In this sense, the most pleiotropic candidate genes were CG34460, Lsd-2 and Spn. Our analyses also revealed extensive genetic variation for all the characters mostly indicated by strong GEIs. Further, our results indicate that GEIs were predominantly explained by changes in ranking order in all cases suggesting that a moderate number of genes are involved in the expression of each character at both temperatures. Most lines displayed a plastic response for at least one trait in either sex. In this regard, P-element insertions affecting plasticity of a large number of traits were associated to the candidate genes Btk29A, CG43340, Drak and jim. Further studies will help to elucidate the relevance of these genes on the morphogenesis of different body structures in natural populations of D. melanogaster.     

Sex and tissue‐specific evolution of developmental plasticity in Drosophila melanogaster

Developmental plasticity influences the size of adult tissues in insects. Tissues can have unique responses to environmental perturbation during development; however, the prevalence of within species evolution of tissue‐specific developmental plasticity remains unclear. To address this, we studied the effects of temperature and nutrition on wing and femur size in D. melanogaster populations from a temperate and tropical region. Wings were more sensitive to temperature, while wings and femurs were equally responsive to nutrition in both populations and sexes. The temperate population was larger under all conditions, except for femurs of starved females. In line with this, we observed greater femur size plasticity in response to starvation in temperate females, leading to differences in sexual dimorphism between populations such that the slope of the reaction norm of sexual dimorphism in the tropical population was double that of the temperate population. Lastly, we observed a significant trend for steeper slopes of reaction norms in temperate than in tropical females, but not in males. These findings highlight that plasticity divergence between populations can evolve heterogeneously across sexes and tissues and that nutritional plasticity can alter sexual dimorphism in D. melanogaster.     

Therapeutic effects of Typha elephantina leave’s extract against paracetamol induced renal injury in rabbits

Present study focuses on ameliorative potential of Typha elephantina leave’s aqueous (TE.AQ) extract against Paracetamol (PCM) induced toxicity in rabbits. We fed the male rabbits with 300 mg PCM in alone and in combination with TE.AQ at different doses i.e. (100, 200 and 300 mg/kg body weight) or silymarin (100 mg/kg) daily for 21 days. PCM in alone significantly (P < 0.5) increased serum urea, uric acid, creatinine, total protein, albumin, globulin and blood urea nitrogen. Serum sodium, potassium and magnesium level were high. The glutathione, radical scavenging activity and Thiobarbituric acid reactive substances were significantly reduced. Treatment with TE.AQ at dose rate 300 mg/kg body weight and Silymarin significantly ameliorated all the parameters when compared with PCM administered group. The 100 and 200 mg of TE.AQ showed no significant effects. The histopathological examination confirmed the therapeutic potential of TE.AQ. These results established the presence of natural antioxidants in Typha elephantina leaves.     

EXPERIMENTS ON TOLERATION OF TEMPERATURE BY DROSOPHILA

1. Most wild stocks of Drosophila melanogaster can be bred indefinitely on banana agar at a temperature of 31°C. There is no relation between the geographical origin of these stocks and their ability to tolerate this temperature. 2. A single wild stock has been found which will breed for only one generation at temperatures above 29°C. The offspring hatched at 31°C. will breed normally at 24°C. This difference from other wild stocks is apparently genetic, but its genetic basis has not yet been worked out. 3. The mutant stocks of D. melanogaster tested by us will breed for only one generation at 31°C. and their offspring at this temperature are also fertile at 24°C. This condition is apparently a physiological effect of the presence of any of the mutant genes in a homozygous condition. 4. Similar tests indicate that wild stocks of D. virilis and Chymomyza procnemis will breed at 31°C., while D. simulans, D. immigrans, and D. funebris will not. The last two species are northern forms not commonly found in the tropics. 5. Both male and female flies from mutant stocks hatched at 31°C. produce offspring at this temperature if mated to flies hatched at 24°C. Their germ cells are therefore capable of development, and the cause of their failure to develop at 31°C. when inbred must lie either in the failure of the germ cells to reach each other or in the fertilization process itself.