Search

The rates of transcription and transposition of retrotransposons vary between lines of Drosophila melanogaster. We have studied the genetics of differences in copia retrotransposon activity by quantitative trait loci (QTL) mapping. Ninety-eight recombinant inbred lines were constructed from two parental lines exhibiting a 10-fold difference in copia transcript level and a 100-fold difference in transposition rate. The lines were scored for 126 molecular markers, copia transcript level, and rate of copia transposition. Transcript level correlated with copia copy number, and the difference in copia copy number Between parental lines accounted for 45.1% of copia transcript-level difference. Most of the remaining difference was accounted for by two transcript-level QTL mapping to cytological positions 27B-30D and 50F-57C on the second chromosome, which accounted for 11.5 and 30.4%, respectively. copia transposition rate was controlled By interacting QTL mapping to the region 27B-48D on the second and 61A-65A and 97D-100A on the third chromosome. The genes controlling copia transcript level are thus not necessarily those involved in controlling copia transposition rate. Segregation of modifying genes, rather than mutations, might explain the variability in copia retrotransposon activity between lines.

The magnitude of segregating variation for bristle number in Drosophila melanogaster exceeds that predicted from models of mutation-selection balance. To evaluate the hypothesis that genotype-environment interaction (GEI) maintains variation for bristle number in nature, we quantified the extent of GEI for abdominal and sternopleural bristles among 98 recombinant inbred lines, derived from two homozygous laboratory strains, in three temperature environments. There was considerable GEl for both bristle traits, which was mainly attributable to changes in rank order of line means. We conducted a genome-wide screen for quantitative trait loci (QTLs) affecting bristle number in each sex and temperature environment, using a dense (3.2-cM) marker map of polymorphic insertion sites of roo transposable elements. Nine sternopleural and 11 abdominal bristle number QTLs were detected. Significant GEI was exhibited by 14 QTLs, but there was heterogeneity among QTLs in their sensitivity to thermal and sexual environments. To further evaluate the hypothesis that GEl maintains variation for bristle number, we require estimates of allelic effects across environments at genetic loci affecting the traits. This level of resolution may be achievable for Drosophila bristle number because candidate loci affecting bristle development often map to the same location as bristle number QTLs.

Single P-element mutagenesis was used to construct 1094 lines with P[lArB] inserts on all three major chromosomes in an isogenic background previously free of P elements. The effects of insertions on bristle number and on viability were assessed by comparison to 392 control lines. The variance and effects of P-element inserts on bristle number and viability were larger than those inferred from spontaneous mutations. The distributions of effects on bristle number were symmetrical and highly leptokurtic, such that a few inserts with large effects caused most of the increase in variance. The distribution of effects on viability were negatively skewed and platykurtic. On average, the effects of P-element insertions on bristle number were partly recessive and on viability were completely recessive. P-element inserts with large effects on bristle number tended to have reduced viability, but the correlation between the absolute value of the effects on bristle number and on viability was not strong. Fifty P-element inserts tagging quantitative trait loci (QTLs) with large effects on bristle number were mapped cytogenetically. Two P-element-induced scabrous alleles and five extramacrochaetae alleles were generated. Single P-element mutagenesis is a powerful method for identifying QTLs at the level of genetic locus.

The technique of divergent artificial selection for about 150 generations and arbitrary mating for 180 generations is used for determining the relationship between abdominal bristle number and sternopleural bristle number in Drosophila melanogaster. The bristle number for replicate selection lines is found to be larger than that of the lines maintained without artificial selection. The unselected lines exhibit an average fitness level which is higher than that of the selection lines.

An average divergence, following 125 generations of selection of 12.0 abdominal and 8.2 sternopleural bristles from the accumulation of new mutations influencing bristle number, resulted from replicated divergent artificial selection for abdominal and sternopleural responses from an inbred strain of Drosophila melanogaster. Extremely asymmetrical responses to selection were evident. Low abdominal and high sternopleural bristle numbers exhibited greater numbers.

Abstract Background Understanding animal development and physiology at a molecular-biological level has been advanced by the ability to determine at high resolution the repertoire of mRNA molecules by whole transcriptome resequencing. This includes the ability to detect and quantify rare abundance transcripts and isoform-specific mRNA variants produced from a gene. The sex hierarchy consists of a pre-mRNA splicing cascade that directs the production of sex-specific transcription factors that specify nearly all sexual dimorphism. We have used deep RNA sequencing to gain insight into how the Drosophila sex hierarchy generates somatic sex differences, by examining gene and transcript isoform expression differences between the sexes in adult head tissues. Results Here we find 1,381 genes that differ in overall expression levels and 1,370 isoform-specific transcripts that differ between males and females. Additionally, we find 512 genes not regulated downstream of transformer that are significantly more highly expressed in males than females. These 512 genes are enriched on the × chromosome and reside adjacent to dosage compensation complex entry sites, which taken together suggests that their residence on the × chromosome might be sufficient to confer male-biased expression. There are no transcription unit structural features, from a set of features, that are robustly significantly different in the genes with significant sex differences in the ratio of isoform-specific transcripts, as compared to random isoform-specific transcripts, suggesting that there is no single molecular mechanism that generates isoform-specific transcript differences between the sexes, even though the sex hierarchy is known to include three pre-mRNA splicing factors. Conclusions We identify thousands of genes that show sex-specific differences in overall gene expression levels, and identify hundreds of additional genes that have differences in the abundance of isoform-specific transcripts. No transcription unit structural feature was robustly enriched in the sex-differentially expressed transcript isoforms. Additionally, we found that many genes with male-biased expression were enriched on the × chromosome and reside adjacent to dosage compensation entry sites, suggesting that differences in sex chromosome composition contributes to dimorphism in gene expression. Taken together, this study provides new insight into the molecular underpinnings of sexual differentiation.

Abstract Background RNA-seq is revolutionizing the way we study transcriptomes. mRNA can be surveyed without prior knowledge of gene transcripts. Alternative splicing of transcript isoforms and the identification of previously unknown exons are being reported. Initial reports of differences in exon usage, and splicing between samples as well as quantitative differences among samples are beginning to surface. Biological variation has been reported to be larger than technical variation. In addition, technical variation has been reported to be in line with expectations due to random sampling. However, strategies for dealing with technical variation will differ depending on the magnitude. The size of technical variance, and the role of sampling are examined in this manuscript. Results In this study three independent Solexa/Illumina experiments containing technical replicates are analyzed. When coverage is low, large disagreements between technical replicates are apparent. Exon detection between technical replicates is highly variable when the coverage is less than 5 reads per nucleotide and estimates of gene expression are more likely to disagree when coverage is low. Although large disagreements in the estimates of expression are observed at all levels of coverage. Conclusions Technical variability is too high to ignore. Technical variability results in inconsistent detection of exons at low levels of coverage. Further, the estimate of the relative abundance of a transcript can substantially disagree, even when coverage levels are high. This may be due to the low sampling fraction and if so, it will persist as an issue needing to be addressed in experimental design even as the next wave of technology produces larger numbers of reads. We provide practical recommendations for dealing with the technical variability, without dramatic cost increases.

Abstract Background Allele-specific expression (ASE) assays can be used to identify cis, trans, and cis-by-trans regulatory variation. Understanding the source of expression variation has important implications for disease susceptibility, phenotypic diversity, and adaptation. While ASE is commonly measured via relative fluorescence at a SNP, next generation sequencing provides an opportunity to measure ASE in an accurate and high-throughput manner using read counts. Results We introduce a Solexa-based method to perform large numbers of ASE assays using only a single lane of a Solexa flowcell. In brief, transcripts of interest, which contain a known SNP, are PCR enriched and barcoded to enable multiplexing. Then high-throughput sequencing is used to estimate allele-specific expression using sequencing counts. To validate this method, we measured the allelic bias in a dilution series and found high correlations between measured and expected values (r>0.9, p < 0.001). We applied this method to a set of 5 genes in a Drosophila simulans parental mix, F1 and introgression and found that for these genes the majority of expression divergence can be explained by cis-regulatory variation. Conclusion We present a new method with the capacity to measure ASE for large numbers of assays using as little as one lane of a Solexa flowcell. This will be a valuable technique for molecular and population genetic studies, as well as for verification of genome-wide data sets.

Abstract Background The pea aphid, Acyrthosiphon pisum, is an emerging genomic model system for studies of polyphenisms, bacterial symbioses, host-plant specialization, and the vectoring of plant viruses. Here we provide estimates of nucleotide diversity and linkage disequilibrium (LD) in native (European) and introduced (United States) populations of the pea aphid. Because introductions can cause population bottlenecks, we hypothesized that U.S. populations harbor lower levels of nucleotide diversity and higher levels of LD than native populations. Results We sampled four non-coding loci from 24 unique aphid clones from the U. S. (12 from New York and 12 from California) and 24 clones from Europe (12 alfalfa and 12 clover specialists). For each locus, we sequenced approximately 1 kb from two amplicons spaced ~10 kb apart to estimate both short range and longer range LD. We sequenced over 250 kb in total. Nucleotide diversity averaged 0.6% across all loci and all populations. LD decayed slowly within ~1 kb but reached much lower levels over ~10 kb. Contrary to our expectations, neither LD nor nucleotide diversity were significantly different between native and introduced populations. Conclusion Both introduced and native populations of pea aphids exhibit low levels of nucleotide diversity and moderate levels of LD. The introduction of pea aphids to North America has not led to a detectable reduction of nucleotide diversity or increase in LD relative to native populations.

Abstract Background A molecular process based genotype-to-phenotype map will ultimately enable us to predict how genetic variation among individuals results in phenotypic alterations. Building such a map is, however, far from straightforward. It requires understanding how molecular variation re-shapes developmental and metabolic networks, and how the functional state of these networks modifies phenotypes in genotype specific way. We focus on the latter problem by describing genetic variation in transcript levels of genes in the InR/TOR pathway among 72 Drosophila melanogaster genotypes. Results We observe tight co-variance in transcript levels of genes not known to influence each other through direct transcriptional control. We summarize transcriptome variation with factor analyses, and observe strong co-variance of gene expression within the dFOXO-branch and within the TOR-branch of the pathway. Finally, we investigate whether major axes of transcriptome variation shape phenotypes expected to be influenced through the InR/TOR pathway. We find limited evidence that transcript levels of individual upstream genes in the InR/TOR pathway predict fly phenotypes in expected ways. However, there is no evidence that these effects are mediated through the major axes of downstream transcriptome variation. Conclusion In summary, our results question the assertion of the 'sparse' nature of genetic networks, while validating and extending candidate gene approaches in the analyses of complex traits.

Abstract Background Co-expression of genes that physically cluster together is a common characteristic of eukaryotic transcriptomes. This organization of transcriptomes suggests that coordinated evolution of gene expression for clustered genes may also be common. Clusters where expression evolution of each gene is not independent of their neighbors are important units for understanding transcriptome evolution. Results We used a common microarray platform to measure gene expression in seven closely related species in the Drosophila melanogaster subgroup, accounting for confounding effects of sequence divergence. To summarize the correlation structure among genes in a chromosomal region, we analyzed the fraction of variation along the first principal component of the correlation matrix. We analyzed the correlation for blocks of consecutive genes to assess patterns of correlation that may be manifest at different scales of coordinated expression. We find that expression of physically clustered genes does evolve in a coordinated manner in many locations throughout the genome. Our analysis shows that relatively few of these clusters are near heterochromatin regions and that these clusters tend to be over-dispersed relative to the rest of the genome. This suggests that these clusters are not the byproduct of local gene clustering. We also analyzed the pattern of co-expression among neighboring genes within a single Drosophila species: D. simulans. For the co-expression clusters identified within this species, we find an under-representation of genes displaying a signature of recurrent adaptive amino acid evolution consistent with previous findings. However, clusters displaying co-evolution of expression among species are enriched for adaptively evolving genes. This finding points to a tie between adaptive sequence evolution and evolution of the transcriptome. Conclusion Our results demonstrate that co-evolution of expression in gene clusters is relatively common among species in the D. melanogaster subgroup. We consider the possibility that local regulation of expression in gene clusters may drive the connection between adaptive sequence and coordinated gene expression evolution.

A study investigates survival in a high-oxygen environment at two adult densities using a set of recombinant inbred lines (RILs) isolated from a natural population of Drosophila melanogaster as well as the genetic correlations and quantitative trait loci (QTLs) for oxidative stress survival and early-age egg productions. The QTLs for oxidative stress survival are relatively abundant and found across the second and third chromosomes, suggesting that the genetic architecture of the trait is quite complex, but there is little evidence that average early-age production and oxidative stress survival are correlated.