Justin P. Blumenstiel

Identification of EMS-induced mutations in Drosophila melanogaster by whole-genome sequencing.

Next-generation methods for rapid whole-genome sequencing enable the identification of single-base-pair mutations in Drosophila by comparing a chromosome bearing a new mutation to the unmutagenized sequence. To validate this approach, we sought to identify the molecular lesion responsible for a recessive EMS-induced mutation affecting egg shell morphology by using Illumina next-generation sequencing. After obtaining sufficient sequence from larvae that were homozygous for either wild-type or mutant chromosomes, we obtained high-quality reads for base pairs composing ∼70% of the third chromosome of both DNA samples. We verified 103 single-base-pair changes between the two chromosomes. Nine changes were nonsynonymous mutations and two were nonsense mutations. One nonsense mutation was in a gene, encore, whose mutations produce an egg shell phenotype also observed in progeny of homozygous mutant mothers. Complementation analysis revealed that the chromosome carried a new functional allele of encore, demonstrating that one round of next-generation sequencing can identify the causative lesion for a phenotype of interest. This new method of whole-genome sequencing represents great promise for mutant mapping in flies, potentially replacing conventional methods.

Evolutionary dynamics of transposable elements in a small RNA world

Transposable elements (TEs) are selfish elements that cause harmful mutations, contribute to the structure of regulatory networks and shape the architecture of genomes. Natural selection against their harmful effects has long been considered the dominant force limiting their spread. It is now clear that a genome defense system of RNA-mediated silencing also plays a crucial role in limiting TE proliferation. A full understanding of TE evolutionary dynamics must consider how these forces jointly determine their proliferation within genomes. Here I consider these forces from two perspectives - dynamics within populations and evolutionary games within the germline. The analysis of TE dynamics from these two perspectives promises to provide new insight into their role in evolution.

Genome complexity reduction for SNP genotyping analysis

Efficient single nucleotide polymorphism (SNP) genotyping methods are necessary to accomplish many current gene discovery goals. A crucial element in large-scale SNP genotyping is the number of individual biochemical reactions that must be performed. An efficient method that can be used to simultaneously amplify a set of genetic loci across a genome with high reliability can provide a valuable tool for large-scale SNP genotyping studies. In this paper we describe and characterize a method that addresses this goal. We have developed a strategy for reducing genome complexity by using degenerate oligonucleotide primer (DOP)-PCR and applied this strategy to SNP genotyping in three complex eukaryotic genomes; human, mouse, and Arabidopsis thaliana. Using a single DOP-PCR primer, SNP loci spread throughout a genome can be amplified and accurately genotyped directly from a DOP-PCR product mixture. DOP-PCRs are extremely reproducible. The DOP-PCR method is transferable to many species of interest. Finally, we describe an in silico approach that can effectively predict the SNP loci amplified in a given DOP-PCR, permitting the design of an efficient set of reactions for large-scale, genome-wide SNP studies.

A Whole-Chromosome Analysis of Meiotic Recombination in Drosophila melanogaster

Although traditional genetic assays have characterized the pattern of crossing over across the genome in Drosophila melanogaster, these assays could not precisely define the location of crossovers. Even less is known about the frequency and distribution of noncrossover gene conversion events. To assess the specific number and positions of both meiotic gene conversion and crossover events, we sequenced the genomes of male progeny from females heterozygous for 93,538 X chromosomal single-nucleotide and InDel polymorphisms. From the analysis of the 30 F1 hemizygous X chromosomes, we detected 15 crossover and 5 noncrossover gene conversion events. Taking into account the nonuniform distribution of polymorphism along the chromosome arm, we estimate that most oocytes experience 1 crossover event and 1.6 gene conversion events per X chromosome pair per meiosis. An extrapolation to the entire genome would predict approximately 5 crossover events and 8.6 conversion events per meiosis. Mean gene conversion tract lengths were estimated to be 476 base pairs, yielding a per nucleotide conversion rate of 0.86 × 10−5 per meiosis. Both of these values are consistent with estimates of conversion frequency and tract length obtained from studies of rosy, the only gene for which gene conversion has been studied extensively in Drosophila. Motif-enrichment analysis revealed a GTGGAAA motif that was enriched near crossovers but not near gene conversions. The low-complexity and frequent occurrence of this motif may in part explain why, in contrast to mammalian systems, no meiotic crossover hotspots have been found in Drosophila.

Molecular evolution under increasing transposable element burden in Drosophila: A speed limit on the evolutionary arms race

BackgroundGenome architecture is profoundly influenced by transposable elements (TEs), and natural selection against their harmful effects is a critical factor limiting their spread. Genome defense by the piRNA silencing pathway also plays a crucial role in limiting TE proliferation. How these two forces jointly determine TE abundance is not well understood. To shed light on the nature of factors that predict TE success, we test three distinct hypotheses in the Drosophila genus. First, we determine whether TE abundance and relaxed genome-wide purifying selection on protein sequences are positively correlated. This serves to test the hypothesis that variation in TE abundance in the Drosophila genus can be explained by the strength of natural selection, relative to drift, acting in parallel against mildly deleterious non-synonymous mutations. Second, we test whether increasing TE abundance is correlated with an increased rate of amino-acid evolution in genes encoding the piRNA machinery, as might be predicted by an evolutionary arms race model. Third, we test whether increasing TE abundance is correlated with greater codon bias in genes of the piRNA machinery. This is predicted if increasing TE abundance selects for increased efficiency in the machinery of genome defense.ResultsSurprisingly, we find neither of the first two hypotheses to be true. Specifically, we found that genome-wide levels of purifying selection, measured by the ratio of non-synonymous to synonymous substitution rates (ω), were greater in species with greater TE abundance. In addition, species with greater TE abundance have greater levels of purifying selection in the piRNA machinery. In contrast, it appears that increasing TE abundance has primarily driven adaptation in the piRNA machinery by increasing codon bias.ConclusionsThese results indicate that within the Drosophila genus, a historically reduced strength of selection relative to drift is unlikely to explain patterns of increased TE success across species. Other factors, such as ecological exposure, are likely to contribute to variation in TE abundances within species. Furthermore, constraints on the piRNA machinery may temper the evolutionary arms race that would drive increasing rates of evolution at the amino acid level. In the face of these constraints, selection may act primarily by improving the translational efficiency of the machinery of genome defense through efficient codon usage.

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.

The Formation of the Central Element of the Synaptonemal Complex May Occur by Multiple Mechanisms: The Roles of the N- and C-Terminal Domains of the Drosophila C(3)G Protein in Mediating Synapsis and Recombination

In Drosophila melanogaster oocytes, the C(3)G protein comprises the transverse filaments (TFs) of the synaptonemal complex (SC). Like other TF proteins, such as Zip1p in yeast and SCP1 in mammals, C(3)G is composed of a central coiled-coil-rich domain flanked by N- and C-terminal globular domains. Here, we analyze in-frame deletions within the N- and C-terminal regions of C(3)G in Drosophila oocytes. As is the case for Zip1p, a C-terminal deletion of C(3)G fails to attach to the lateral elements of the SC. Instead, this C-terminal deletion protein forms a large cylindrical polycomplex structure. EM analysis of this structure reveals a polycomplex of concentric rings alternating dark and light bands. However, unlike both yeast and mammals, all three proteins deleted for N-terminal regions completely abolished both SC and polycomplex formation. Both the N- and C-terminal deletions significantly reduce or abolish meiotic recombination similarly to c(3)G null homozygotes. To explain these data, we propose that in Drosophila the N terminus, but not the C-terminal globular domain, of C(3)G is critical for the formation of antiparallel pairs of C(3)G homodimers that span the central region and thus for assembly of complete TFs, while the C terminus is required to affix these homodimers to the lateral elements.

Whole-Genome Analysis of Individual Meiotic Events in Drosophila melanogaster Reveals That Noncrossover Gene Conversions Are Insensitive to Interference and the Centromere Effect

A century of genetic analysis has revealed that multiple mechanisms control the distribution of meiotic crossover events. In Drosophila melanogaster, two significant positional controls are interference and the strongly polar centromere effect. Here, we assess the factors controlling the distribution of crossovers (COs) and noncrossover gene conversions (NCOs) along all five major chromosome arms in 196 single meiotic divisions to generate a more detailed understanding of these controls on a genome-wide scale. Analyzing the outcomes of single meiotic events allows us to distinguish among different classes of meiotic recombination. In so doing, we identified 291 NCOs spread uniformly among the five major chromosome arms and 541 COs (including 52 double crossovers and one triple crossover). We find that unlike COs, NCOs are insensitive to the centromere effect and do not demonstrate interference. Although the positions of COs appear to be determined predominately by the long-range influences of interference and the centromere effect, each chromosome may display a different pattern of sensitivity to interference, suggesting that interference may not be a uniform global property. In addition, unbiased sequencing of a large number of individuals allows us to describe the formation of de novo copy number variants, the majority of which appear to be mediated by unequal crossing over between transposable elements. This work has multiple implications for our understanding of how meiotic recombination is regulated to ensure proper chromosome segregation and maintain genome stability.

Corolla Is a Novel Protein That Contributes to the Architecture of the Synaptonemal Complex of Drosophila

In most organisms the synaptonemal complex (SC) connects paired homologs along their entire length during much of meiotic prophase. To better understand the structure of the SC, we aim to identify its components and to determine how each of these components contributes to SC function. Here, we report the identification of a novel SC component in Drosophila melanogaster female oocytes, which we have named Corolla. Using structured illumination microscopy, we demonstrate that Corolla is a component of the central region of the SC. Consistent with its localization, we show by yeast two-hybrid analysis that Corolla strongly interacts with Cona, a central element protein, demonstrating the first direct interaction between two inner-synaptonemal complex proteins in Drosophila. These observations help provide a more complete model of SC structure and function in Drosophila females.

piRNAs Are Associated with Diverse Transgenerational Effects on Gene and Transposon Expression in a Hybrid Dysgenic Syndrome of D. virilis.

Sexual reproduction allows transposable elements (TEs) to proliferate, leading to rapid divergence between populations and species. A significant outcome of divergence in the TE landscape is evident in hybrid dysgenic syndromes, a strong form of genomic incompatibility that can arise when (TE) family abundance differs between two parents. When TEs inherited from the father are absent in the mothers genome, TEs can become activated in the progeny, causing germline damage and sterility. Studies in Drosophila indicate that dysgenesis can occur when TEs inherited paternally are not matched with a pool of corresponding TE silencing PIWI-interacting RNAs (piRNAs) provisioned by the female germline. Using the D. virilis syndrome of hybrid dysgenesis as a model, we characterize the effects that divergence in TE profile between parents has on offspring. Overall, we show that divergence in the TE landscape is associated with persisting differences in germline TE expression when comparing genetically identical females of reciprocal crosses and these differences are transmitted to the next generation. Moreover, chronic and persisting TE expression coincides with increased levels of genic piRNAs associated with reduced gene expression. Combined with these effects, we further demonstrate that gene expression is idiosyncratically influenced by differences in the genic piRNA profile of the parents that arise though polymorphic TE insertions. Overall, these results support a model in which early germline events in dysgenesis establish a chronic, stable state of both TE and gene expression in the germline that is maintained through adulthood and transmitted to the next generation. This work demonstrates that divergence in the TE profile is associated with diverse piRNA-mediated transgenerational effects on gene expression within populations.