Anoja Perera

Aneuploidy Underlies Rapid Adaptive Evolution of Yeast Cells Deprived of a Conserved Cytokinesis Motor

The ability to evolve is a fundamental feature of biological systems, but the mechanisms underlying this capacity and the evolutionary dynamics of conserved core processes remain elusive. We show that yeast cells deleted of MYO1, encoding the only myosin II normally required for cytokinesis, rapidly evolved divergent pathways to restore growth and cytokinesis. The evolved cytokinesis phenotypes correlated with specific changes in the transcriptome. Polyploidy and aneuploidy were common genetic alterations in the best evolved strains, and aneuploidy could account for gene expression changes due directly to altered chromosome stoichiometry as well as to downstream effects. The phenotypic effect of aneuploidy could be recapitulated with increased copy numbers of specific regulatory genes in myo1Delta cells. These results demonstrate the evolvability of even a well-conserved process and suggest that changes in chromosome stoichiometry provide a source of heritable variation driving the emergence of adaptive phenotypes when the cell division machinery is strongly perturbed.

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.

Comparison of Commercially Available Target Enrichment Methods for Next-Generation Sequencing

Isolating high-priority segments of genomes greatly enhances the efficiency of next-generation sequencing (NGS) by allowing researchers to focus on their regions of interest. For the 2010-11 DNA Sequencing Research Group (DSRG) study, we compared outcomes from two leading companies, Agilent Technologies (Santa Clara, CA, USA) and Roche NimbleGen (Madison, WI, USA), which offer custom-targeted genomic enrichment methods. Both companies were provided with the same genomic sample and challenged to capture identical genomic locations for DNA NGS. The target region totaled 3.5 Mb and included 31 individual genes and a 2-Mb contiguous interval. Each company was asked to design its best assay, perform the capture in replicates, and return the captured material to the DSRG-participating laboratories. Sequencing was performed in two different laboratories on Genome Analyzer IIx systems (Illumina, San Diego, CA, USA). Sequencing data were analyzed for sensitivity, specificity, and coverage of the desired regions. The success of the enrichment was highly dependent on the design of the capture probes. Overall, coverage variability was higher for the Agilent samples. As variant discovery is the ultimate goal for a typical targeted sequencing project, we compared samples for their ability to sequence single-nucleotide polymorphisms (SNPs) as a test of the ability to capture both chromosomes from the sample. In the targeted regions, we detected 2546 SNPs with the NimbleGen samples and 2071 with Agilents. When limited to the regions that both companies included as baits, the number of SNPs was ∼1000 for each, with Agilent and NimbleGen finding a small number of unique SNPs not found by the other.

Incomplete penetrance and phenotypic variability characterize Gdf6-attributable oculo-skeletal phenotypes

Proteins of the bone morphogenetic protein (BMP) family are known to have a role in ocular and skeletal development; however, because of their widespread expression and functional redundancy, less progress has been made identifying the roles of individual BMPs in human disease. We identified seven heterozygous mutations in growth differentiation factor 6 (GDF6), a member of the BMP family, in patients with both ocular and vertebral anomalies, characterized their effects with a SOX9-reporter assay and western analysis, and demonstrated comparable phenotypes in model organisms with reduced Gdf6 function. We observed a spectrum of ocular and skeletal anomalies in morphant zebrafish, the latter encompassing defective tail formation and altered expression of somite markers noggin1 and noggin2. Gdf6(+/-) mice exhibited variable ocular phenotypes compatible with phenotypes observed in patients and zebrafish. Key differences evident between patients and animal models included pleiotropic effects, variable expressivity and incomplete penetrance. These data establish the important role of this determinant in ocular and vertebral development, demonstrate the complex genetic inheritance of these phenotypes, and further understanding of BMP function and its contributions to human disease.

Piwi Is Required in Multiple Cell Types to Control Germline Stem Cell Lineage Development in the Drosophila Ovary

The piRNA pathway plays an important role in maintaining genome stability in the germ line by silencing transposable elements (TEs) from fly to mammals. As a highly conserved piRNA pathway component, Piwi is widely expressed in both germ cells and somatic cells in the Drosophila ovary and is required for piRNA production in both cell types. In addition to its known role in somatic cap cells to maintain germline stem cells (GSCs), this study has demonstrated that Piwi has novel functions in somatic cells and germ cells of the Drosophila ovary to promote germ cell differentiation. Piwi knockdown in escort cells causes a reduction in escort cell (EC) number and accumulation of undifferentiated germ cells, some of which show active BMP signaling, indicating that Piwi is required to maintain ECs and promote germ cell differentiation. Simultaneous knockdown of dpp, encoding a BMP, in ECs can partially rescue the germ cell differentiation defect, indicating that Piwi is required in ECs to repress dpp. Consistent with its key role in piRNA production, TE transcripts increase significantly and DNA damage is also elevated in the piwi knockdown somatic cells. Germ cell-specific knockdown of piwi surprisingly causes depletion of germ cells before adulthood, suggesting that Piwi might control primordial germ cell maintenance or GSC establishment. Finally, Piwi inactivation in the germ line of the adult ovary leads to gradual GSC loss and germ cell differentiation defects, indicating the intrinsic role of Piwi in adult GSC maintenance and differentiation. This study has revealed new germline requirement of Piwi in controlling GSC maintenance and lineage differentiation as well as its new somatic function in promoting germ cell differentiation. Therefore, Piwi is required in multiple cell types to control GSC lineage development in the Drosophila ovary.

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.

Sequencing of the Dutch elm disease fungus genome using the Roche/454 GS-FLX Titanium System in a comparison of multiple genomics core facilities.

As part of the DNA Sequencing Research Group of the Association of Biomolecular Resource Facilities, we have tested the reproducibility of the Roche/454 GS-FLX Titanium System at five core facilities. Experience with the Roche/454 system ranged from <10 to >340 sequencing runs performed. All participating sites were supplied with an aliquot of a common DNA preparation and were requested to conduct sequencing at a common loading condition. The evaluation of sequencing yield and accuracy metrics was assessed at a single site. The study was conducted using a laboratory strain of the Dutch elm disease fungus Ophiostoma novo-ulmi strain H327, an ascomycete, vegetatively haploid fungus with an estimated genome size of 30-50 Mb. We show that the Titanium System is reproducible, with some variation detected in loading conditions, sequencing yield, and homopolymer length accuracy. We demonstrate that reads shorter than the theoretical minimum length are of lower overall quality and not simply truncated reads. The O. novo-ulmi H327 genome assembly is 31.8 Mb and is comprised of eight chromosome-length linear scaffolds, a circular mitochondrial conti of 66.4 kb, and a putative 4.2-kb linear plasmid. We estimate that the nuclear genome encodes 8613 protein coding genes, and the mitochondrion encodes 15 genes and 26 tRNAs.

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.

Evaluation of commercially available RNA amplification kits for RNA sequencing using very low input amounts of total RNA

This article includes supplemental data. Please visit http://www.fasebj.org to obtain this information.Multiple recent publications on RNA sequencing (RNA-seq) have demonstrated the power of next-generation sequencing technologies in whole-transcriptome analysis. Vendor-specific protocols used for RNA library construction often require at least 100 ng total RNA. However, under certain conditions, much less RNA is available for library construction. In these cases, effective transcriptome profiling requires amplification of subnanogram amounts of RNA. Several commercial RNA amplification kits are available for amplification prior to library construction for next-generation sequencing, but these kits have not been comprehensively field evaluated for accuracy and performance of RNA-seq for picogram amounts of RNA. To address this, 4 types of amplification kits were tested with 3 different concentrations, from 5 ng to 50 pg, of a commercially available RNA. Kits were tested at multiple sites to assess reproducibility and ease of use. The human total reference RNA used was spiked with a control pool of RNA molecules in order to further evaluate quantitative recovery of input material. Additional control data sets were generated from libraries constructed following polyA selection or ribosomal depletion using established kits and protocols. cDNA was collected from the different sites, and libraries were synthesized at a single site using established protocols. Sequencing runs were carried out on the Illumina platform. Numerous metrics were compared among the kits and dilutions used. Overall, no single kit appeared to meet all the challenges of small input material. However, it is encouraging that excellent data can be recovered with even the 50 pg input total RNA.

Single-cell transcriptome analysis of avian neural crest migration reveals signatures of invasion and molecular transitions

Neural crest cells migrate throughout the embryo, but how cells move in a directed and collective manner has remained unclear. Here, we perform the first single-cell transcriptome analysis of cranial neural crest cell migration at three progressive stages in chick and identify and establish hierarchical relationships between cell position and time-specific transcriptional signatures. We determine a novel transcriptional signature of the most invasive neural crest Trailblazer cells that is consistent during migration and enriched for approximately 900 genes. Knockdown of several Trailblazer genes shows significant but modest changes to total distance migrated. However, in vivo expression analysis by RNAscope and immunohistochemistry reveals some salt and pepper patterns that include strong individual Trailblazer gene expression in cells within other subregions of the migratory stream. These data provide new insights into the molecular diversity and dynamics within a neural crest cell migratory stream that underlie complex directed and collective cell behaviors.