The impact of age,biogenesis, and genomic clustering on Drosophila microRNA evolution

The molecular evolutionary signatures of miRNAs inform our understanding of their emergence, biogenesis, and function. The known signatures of miRNA evolution have derived mostly from the analysis of deeply conserved, canonical loci. In this study, we examine the impact of age, biogenesis pathway, and genomic arrangement on the evolutionary properties of Drosophila miRNAs. Crucial to the accuracy of our results was our curation of high-quality miRNA alignments, which included nearly 150 corrections to ortholog calls and nucleotide sequences of the global 12-way Drosophilid alignments currently available. Using these data, we studied primary sequence conservation, normalized free-energy values, and types of structure-preserving substitutions. We expand upon common miRNA evolutionary patterns that reflect fundamental features of miRNAs that are under functional selection. We observe that melanogaster-subgroup-specific miRNAs, although recently emerged and rapidly evolving, nonetheless exhibit evolutionary signatures that are similar to well-conserved miRNAs and distinct from other structured noncoding RNAs and bulk conserved non-miRNA hairpins. This provides evidence that even young miRNAs may be selected for regulatory activities. More strikingly, we observe that mirtrons and clustered miRNAs both exhibit distinct evolutionary properties relative to solo, well-conserved miRNAs, even after controlling for sequence depth. These studies highlight the previously unappreciated impact of biogenesis strategy and genomic location on the evolutionary dynamics of miRNAs, and affirm that miRNAs do not evolve as a unitary class.     

Genome Engineering of Drosophila with the CRISPR RNA-Guided Cas9 Nuclease

We have adapted a bacterial CRISPR RNA/Cas9 system to precisely engineer the Drosophila genome and report that Cas9-mediated genomic modifications are efficiently transmitted through the germline. This RNA-guided Cas9 system can be rapidly programmed to generate targeted alleles for probing gene function in Drosophila.     

Organellar genomes of the four-toothed moss,Tetraphis pellucida

Background

Mosses are the largest of the three extant clades of gametophyte-dominant land plants and remain poorly studied using comparative genomic methods. Major monophyletic moss lineages are characterised by different types of a spore dehiscence apparatus called the peristome, and the most important unsolved problem in higher-level moss systematics is the branching order of these peristomate clades. Organellar genome sequencing offers the potential to resolve this issue through the provision of both genomic structural characters and a greatly increased quantity of nucleotide substitution characters, as well as to elucidate organellar evolution in mosses. We publish and describe the chloroplast and mitochondrial genomes of Tetraphis pellucida, representative of the most phylogenetically intractable and morphologically isolated peristomate lineage.

Results

Assembly of reads from Illumina SBS and Pacific Biosciences RS sequencing reveals that the Tetraphis chloroplast genome comprises 127,489 bp and the mitochondrial genome 107,730 bp. Although genomic structures are similar to those of the small number of other known moss organellar genomes, the chloroplast lacks the petN gene (in common with Tortula ruralis) and the mitochondrion has only a non-functional pseudogenised remnant of nad7 (uniquely amongst known moss chondromes).

Conclusions

Structural genomic features exist with the potential to be informative for phylogenetic relationships amongst the peristomate moss lineages, and thus organellar genome sequences are urgently required for exemplars from other clades. The unique genomic and morphological features of Tetraphis confirm its importance for resolving one of the major questions in land plant phylogeny and for understanding the evolution of the peristome, a likely key innovation underlying the diversity of mosses. The functional loss of nad7 from the chondrome is now shown to have occurred independently in all three bryophyte clades as well as in the early-diverging tracheophyte Huperzia squarrosa.     

Ack kinase regulates CTP synthase filaments during Drosophila oogenesis

The enzyme CTP synthase (CTPS) dynamically assembles into macromolecular filaments in bacteria, yeast, Drosophila, and mammalian cells, but the role of this morphological reorganization in regulating CTPS activity is controversial. During Drosophila oogenesis, CTPS filaments are transiently apparent in ovarian germline cells during a period of intense genomic endoreplication and stockpiling of ribosomal RNA. Here, we demonstrate that CTPS filaments are catalytically active and that their assembly is regulated by the non-receptor tyrosine kinase DAck, the Drosophila homologue of mammalian Ack1 (activated cdc42-associated kinase 1), which we find also localizes to CTPS filaments. Egg chambers from flies deficient in DAck or lacking DAck catalytic activity exhibit disrupted CTPS filament architecture and morphological defects that correlate with reduced fertility. Furthermore, ovaries from these flies exhibit reduced levels of total RNA, suggesting that DAck may regulate CTP synthase activity. These findings highlight an unexpected function for DAck and provide insight into a novel pathway for the developmental control of an essential metabolic pathway governing nucleotide biosynthesis.     

The adaptive role of transposable elements in the Drosophila genome

Transposable elements (TEs) are short DNA sequences with the capacity to move between different sites in the genome. This ability provides them with the capacity to mutate the genome in many different ways, from subtle regulatory mutations to gross genomic rearrangements. The potential adaptive significance of TEs was recognized by those involved in their initial discovery although it was hotly debated afterwards. For more than two decades, TEs were considered to be intragenomic parasites leading to almost exclusively detrimental effects to the host genome. The sequencing of the Drosophila melanogaster genome provided an unprecedented opportunity to study TEs and led to the identification of the first TE-induced adaptations in this species. These studies were followed by a systematic genome-wide search for adaptive insertions that allowed for the first time to infer that TEs contribute substantially to adaptive evolution. This study also revealed that there are at least twice as many TE-induced adaptations that remain to be identified. To gain a better understanding of the adaptive role of TEs in the genome we clearly need to (i) identify as many adaptive TEs as possible in a range of Drosophila species as well as (ii) carry out in-depth investigations of the effects of adaptive TEs on as many phenotypes as possible.     

Cytological Analysis of Spermatogenesis: Live and Fixed Preparations of Drosophila Testes

Drosophila melanogaster is a powerful model system that has been widely used to elucidate a variety of biological processes. For example, studies of both the female and male germ lines of Drosophila have contributed greatly to the current understanding of meiosis as well as stem cell biology. Excellent protocols are available in the literature for the isolation and imaging of Drosophila ovaries and testes^3-12. Herein, methods for the dissection and preparation of Drosophila testes for microscopic analysis are described with an accompanying video demonstration. A protocol for isolating testes from the abdomen of adult males and preparing slides of live tissue for analysis by phase-contrast microscopy as well as a protocol for fixing and immunostaining testes for analysis by fluorescence microscopy are presented. These techniques can be applied in the characterization of Drosophila mutants that exhibit defects in spermatogenesis as well as in the visualization of subcellular localizations of proteins.     

2D and 3D Chromosome Painting in Malaria Mosquitoes

Fluorescent in situ hybridization (FISH) of whole arm chromosome probes is a robust technique for mapping genomic regions of interest, detecting chromosomal rearrangements, and studying three-dimensional (3D) organization of chromosomes in the cell nucleus. The advent of laser capture microdissection (LCM) and whole genome amplification (WGA) allows obtaining large quantities of DNA from single cells. The increased sensitivity of WGA kits prompted us to develop chromosome paints and to use them for exploring chromosome organization and evolution in non-model organisms. Here, we present a simple method for isolating and amplifying the euchromatic segments of single polytene chromosome arms from ovarian nurse cells of the African malaria mosquito Anopheles gambiae. This procedure provides an efficient platform for obtaining chromosome paints, while reducing the overall risk of introducing foreign DNA to the sample. The use of WGA allows for several rounds of re-amplification, resulting in high quantities of DNA that can be utilized for multiple experiments, including 2D and 3D FISH. We demonstrated that the developed chromosome paints can be successfully used to establish the correspondence between euchromatic portions of polytene and mitotic chromosome arms in An. gambiae. Overall, the union of LCM and single-chromosome WGA provides an efficient tool for creating significant amounts of target DNA for future cytogenetic and genomic studies.     

Fine-Tuning of the Actin Cytoskeleton and Cell Adhesion During Drosophila Development by the Unconventional Guanine Nucleotide Exchange Factors Myoblast City and Sponge

The evolutionarily conserved Dock proteins function as unconventional guanine nucleotide exchange factors (GEFs). Upon binding to engulfment and cell motility (ELMO) proteins, Dock–ELMO complexes activate the Rho family of small GTPases to mediate a diverse array of biological processes, including cell motility, apoptotic cell clearance, and axon guidance. Overlapping expression patterns and functional redundancy among the 11 vertebrate Dock family members, which are subdivided into four families (Dock A, B, C, and D), complicate genetic analysis. In both vertebrate and invertebrate systems, the actin dynamics regulator, Rac, is the target GTPase of the Dock-A subfamily. However, it remains unclear whether Rac or Rap1 are the in vivo downstream GTPases of the Dock-B subfamily. Drosophila melanogaster is an excellent genetic model organism for understanding Dock protein function as its genome encodes one ortholog per subfamily: Myoblast city (Mbc; Dock A) and Sponge (Spg; Dock B). Here we show that the roles of Spg and Mbc are not redundant in the Drosophila somatic muscle or the dorsal vessel. Moreover, we confirm the in vivo role of Mbc upstream of Rac and provide evidence that Spg functions in concert with Rap1, possibly to regulate aspects of cell adhesion. Together these data show that Mbc and Spg can have different downstream GTPase targets. Our findings predict that the ability to regulate downstream GTPases is dependent on cellular context and allows for the fine-tuning of actin cytoskeletal or cell adhesion events in biological processes that undergo cell morphogenesis.     

Genome bias influences amino acid choices: analysis of amino acid substitution and re-compilation of substitution matrices exclusive to an AT-biased genome

The genomic era has seen a remarkable increase in the number of genomes being sequenced and annotated. Nonetheless, annotation remains a serious challenge for compositionally biased genomes. For the preliminary annotation, popular nucleotide and protein comparison methods such as BLAST are widely employed. These methods make use of matrices to score alignments such as the amino acid substitution matrices. Since a nucleotide bias leads to an overall bias in the amino acid composition of proteins, it is possible that a genome with nucleotide bias may have introduced atypical amino acid substitutions in its proteome. Consequently, standard matrices fail to perform well in sequence analysis of these genomes. To address this issue, we examined the amino acid substitution in the AT-rich genome of Plasmodium falciparum, chosen as a reference and reconstituted a substitution matrix in the genome's context. The matrix was used to generate protein sequence alignments for the parasite proteins that improved across the functional regions. We attribute this to the consistency that may have been achieved amid the target and background frequencies calculated exclusively in our study. This study has important implications on annotation of proteins that are of experimental interest but give poor sequence alignments with standard conventional matrices.     

Drosophila genome project: One-hit coverage in yeast artificial chromosomes

We present a strategy for assembling a physical map of the genome of Drosophila melanogaster based on yeast artificial chromosomes (YACs). In this paper we report 500 YACs containing inserts of Drosophila DNA averaging 200 kb that have been assigned positions on the physical map by means of in situ hybridization with salivary gland chromosomes. The cloned DNA fragments have randomly sheared ends (DY clones) or ends generated by partial digestion with either NotI (N clones) or EcoRI (E clones). Relative to the euchromatic portion of the genome, the size distribution and genomic positions of the clones reveal no significant bias in the completeness or randomness of genome coverage. The 500 mapped euchromatic clones contain an aggregate of approximately 100 million base pairs of DNA, which is approximately one genome equivalent of Drosophila euchromatin.by W. Hennig     

When peptides fly: Advances in Drosophila proteomics

In the past decade, improvements in genome annotation, protein fractionation methods and mass spectrometry instrumentation resulted in rapid growth of Drosophila proteomics. This review presents the current status of proteomics research in the fly. Areas that have seen major advances in recent years include efforts to map and catalog the Drosophila proteome and high-throughput as well as targeted studies to analyze protein–protein interactions and post-translational modifications. Stable isotope labeling of flies and other applications of quantitative proteomics have opened up new possibilities for functional analyses. It is clear that proteomics is becoming an indispensable tool in Drosophila systems biology research that adds a unique dimension to studying gene function.     

The transposable elements of the Drosophila melanogaster euchromatin: a genomics perspective

Background

Transposable elements are found in the genomes of nearly all eukaryotes. The recent completion of the Release 3 euchromatic genomic sequence of Drosophila melanogaster by the Berkeley Drosophila Genome Project has provided precise sequence for the repetitive elements in the Drosophila euchromatin. We have used this genomic sequence to describe the euchromatic transposable elements in the sequenced strain of this species.

Results

We identified 85 known and eight novel families of transposable element varying in copy number from one to 146. A total of 1,572 full and partial transposable elements were identified, comprising 3.86% of the sequence. More than two-thirds of the transposable elements are partial. The density of transposable elements increases an average of 4.7 times in the centromere-proximal regions of each of the major chromosome arms. We found that transposable elements are preferentially found outside genes; only 436 of 1,572 transposable elements are contained within the 61.4 Mb of sequence that is annotated as being transcribed. A large proportion of transposable elements is found nested within other elements of the same or different classes. Lastly, an analysis of structural variation from different families reveals distinct patterns of deletion for elements belonging to different classes.

Conclusions

This analysis represents an initial characterization of the transposable elements in the Release 3 euchromatic genomic sequence of D. melanogaster for which comparison to the transposable elements of other organisms can begin to be made. These data have been made available on the Berkeley Drosophila Genome Project website for future analyses.     

Partial venom gland transcriptome of a Drosophila parasitoid wasp, Leptopilina heterotoma, reveals novel and shared bioactive profiles with stinging Hymenoptera

Analysis of natural host-parasite relationships reveals the evolutionary forces that shape the delicate and unique specificity characteristic of such interactions. The accessory long gland-reservoir complex of the wasp Leptopilina heterotoma (Figitidae) produces venom with virus-like particles. Upon delivery, venom components delay host larval development and completely block host immune responses. The host range of this Drosophila endoparasitoid notably includes the highly-studied model organism, Drosophila melanogaster. Categorization of 827 unigenes, using similarity as an indicator of putative homology, reveals that approximately 25% are novel or classified as hypothetical proteins. Most of the remaining unigenes are related to processes involved in signaling, cell cycle, and cell physiology including detoxification, protein biogenesis, and hormone production. Analysis of L. heterotoma's predicted venom gland proteins demonstrates conservation among endo- and ectoparasitoids within the Apocrita (e.g., this wasp and the jewel wasp Nasonia vitripennis) and stinging aculeates (e.g., the honey bee and ants). Enzyme and KEGG pathway profiling predicts that kinases, esterases, and hydrolases may contribute to venom activity in this unique wasp. To our knowledge, this investigation is among the first functional genomic studies for a natural parasitic wasp of Drosophila. Our findings will help explain how L. heterotoma shuts down its hosts' immunity and shed light on the molecular basis of a natural arms race between these insects.     

The generation of chromosomal deletions to provide extensive coverage and subdivision of the Drosophila melanogaster genome

Background

Chromosomal deletions are used extensively in Drosophila melanogaster genetics research. Deletion mapping is the primary method used for fine-scale gene localization. Effective and efficient deletion mapping requires both extensive genomic coverage and a high density of molecularly defined breakpoints across the genome.

Results

A large-scale resource development project at the Bloomington Drosophila Stock Center has improved the choice of deletions beyond that provided by previous projects. FLP-mediated recombination between FRT-bearing transposon insertions was used to generate deletions, because it is efficient and provides single-nucleotide resolution in planning deletion screens. The 793 deletions generated pushed coverage of the euchromatic genome to 98.4%. Gaps in coverage contain haplolethal and haplosterile genes, but the sizes of these gaps were minimized by flanking these genes as closely as possible with deletions. In improving coverage, a complete inventory of haplolethal and haplosterile genes was generated and extensive information on other haploinsufficient genes was compiled. To aid mapping experiments, a subset of deletions was organized into a Deficiency Kit to provide maximal coverage efficiently. To improve the resolution of deletion mapping, screens were planned to distribute deletion breakpoints evenly across the genome. The median chromosomal interval between breakpoints now contains only nine genes and 377 intervals contain only single genes.

Conclusions

Drosophila melanogaster now has the most extensive genomic deletion coverage and breakpoint subdivision as well as the most comprehensive inventory of haploinsufficient genes of any multicellular organism. The improved selection of chromosomal deletion strains will be useful to nearly all Drosophila researchers.     

Discovery of 17 conserved structural RNAs in fungi

Many non-coding RNAs with known functions are structurally conserved: their intramolecular secondary and tertiary interactions are maintained across evolutionary time. Consequently, the presence of conserved structure in multiple sequence alignments can be used to identify candidate functional non-coding RNAs. Here, we present a bioinformatics method that couples iterative homology search with covariation analysis to assess whether a genomic region has evidence of conserved RNA structure. We used this method to examine all unannotated regions of five well-studied fungal genomes (Saccharomyces cerevisiae, Candida albicans, Neurospora crassa, Aspergillus fumigatus, and Schizosaccharomyces pombe). We identified 17 novel structurally conserved non-coding RNA candidates, which include four H/ACA box small nucleolar RNAs, four intergenic RNAs and nine RNA structures located within the introns and untranslated regions (UTRs) of mRNAs. For the two structures in the 3′ UTRs of the metabolic genes GLY1 and MET13, we performed experiments that provide evidence against them being eukaryotic riboswitches.     

Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-540) contains supplementary material, which is available to authorized users.     

Analysis of covariation in an SH3 domain sequence alignment: applications in tertiary contact prediction and the design of compensating hydrophobic core substitutions

We have analyzed sequence covariation in an alignment of 266 non-redundant SH3 domain sequences using chi-squared statistical methods. Artifactual covariations arising from close evolutionary relationships among certain sequence subgroups were eliminated using empirically derived sequence diversity thresholds. This covariation detection method was able to predict residue-residue contacts (side-chain centres of mass within 8 A) in the structure of the SH3 domain with an accuracy of 85 %, which is greater than that achieved in many previous covariation studies. In examining the positions involved most frequently in covariations, we discovered a dramatic over-representation of a subset of five hydrophobic core positions. This covariation information was used to design second and third site substitutions that could compensate for highly destabilizing hydrophobic core substitutions in the Fyn SH3 domain, thus providing experimental data to validate the covariation analysis. The testing of our covariation detection method on 15 other alignments showed that the accuracy of contact prediction is highly variable depending on which sequence alignment is used, and useful levels of prediction accuracy were obtained with only approximately one-third of alignments. The results presented here provide insight into the difficulties inherent in covariation analysis, and suggest that it may have limited usefulness in tertiary structure prediction. On the other hand, our ability to use covariation analysis to design stabilizing combinations of hydrophobic core substitutions attests to its potential utility for gaining deeper insight into the stability determinants and functional mechanisms of proteins with known three-dimensional structures.     

Heterochromatic sequences in a Drosophila whole-genome shotgun assembly

Background

Most eukaryotic genomes include a substantial repeat-rich fraction termed heterochromatin, which is concentrated in centric and telomeric regions. The repetitive nature of heterochromatic sequence makes it difficult to assemble and analyze. To better understand the heterochromatic component of the Drosophila melanogaster genome, we characterized and annotated portions of a whole-genome shotgun sequence assembly.

Results

WGS3, an improved whole-genome shotgun assembly, includes 20.7 Mb of draft-quality sequence not represented in the Release 3 sequence spanning the euchromatin. We annotated this sequence using the methods employed in the re-annotation of the Release 3 euchromatic sequence. This analysis predicted 297 protein-coding genes and six non-protein-coding genes, including known heterochromatic genes, and regions of similarity to known transposable elements. Bacterial artificial chromosome (BAC)-based fluorescence in situ hybridization analysis was used to correlate the genomic sequence with the cytogenetic map in order to refine the genomic definition of the centric heterochromatin; on the basis of our cytological definition, the annotated Release 3 euchromatic sequence extends into the centric heterochromatin on each chromosome arm.

Conclusions

Whole-genome shotgun assembly produced a reliable draft-quality sequence of a significant part of the Drosophila heterochromatin. Annotation of this sequence defined the intron-exon structures of 30 known protein-coding genes and 267 protein-coding gene models. The cytogenetic mapping suggests that an additional 150 predicted genes are located in heterochromatin at the base of the Release 3 euchromatic sequence. Our analysis suggests strategies for improving the sequence and annotation of the heterochromatic portions of the Drosophila and other complex genomes.