AftrRAD: a pipeline for accurate and efficient de novo assembly of RADseq data

An increase in studies using restriction site‐associated DNA sequencing (RADseq) methods has led to a need for both the development and assessment of novel bioinformatic tools that aid in the generation and analysis of these data. Here, we report the availability of AftrRAD, a bioinformatic pipeline that efficiently assembles and genotypes RADseq data, and outputs these data in various formats for downstream analyses. We use simulated and experimental data sets to evaluate AftrRAD's ability to perform accurate de novo assembly of loci, and we compare its performance with two other commonly used programs, stacks and pyrad. We demonstrate that AftrRAD is able to accurately assemble loci, while accounting for indel variation among alleles, in a more computationally efficient manner than currently available programs. AftrRAD run times are not strongly affected by the number of samples in the data set, making this program a useful tool when multicore systems are not available for parallel processing, or when data sets include large numbers of samples.     

Ray Meta: scalable de novo metagenome assembly and profiling

Voluminous parallel sequencing datasets, especially metagenomic experiments, require distributed computing for de novo assembly and taxonomic profiling. Ray Meta is a massively distributed metagenome assembler that is coupled with Ray Communities, which profiles microbiomes based on uniquely-colored k-mers. It can accurately assemble and profile a three billion read metagenomic experiment representing 1,000 bacterial genomes of uneven proportions in 15 hours with 1,024 processor cores, using only 1.5 GB per core. The software will facilitate the processing of large and complex datasets, and will help in generating biological insights for specific environments. Ray Meta is open source and available at http://denovoassembler.sf.net.     

Drosophila parthenogenesis: a tool to decipher centrosomal vs acentrosomal spindle assembly pathways

Development of unfertilized eggs in the parthenogenetic strain K23-O-im of Drosophila mercatorum requires the stochastic interactions of self-assembled centrosomes with the female chromatin. In a portion of the unfertilized eggs that do not assemble centrosomes, microtubules organize a bipolar anastral mitotic spindle around the chromatin like the one formed during the first female meiosis, suggesting that similar pathways may be operative. In the cytoplasm of eggs in which centrosomes do form, monastral and biastral spindles are found. Analysis by laser scanning confocal microscopy suggests that these spindles are derived from the stochastic interaction of astral microtubules directly with kinetochore regions or indirectly with kinetochore microtubules. Our findings are consistent with the idea that mitotic spindle assembly requires both acentrosomal and centrosomal pathways, strengthening the hypothesis that astral microtubules can dictate the organization of the spindle by capturing kinetochore microtubules.     

RADProc: A computationally efficient de novo locus assembler for population studies using RADseq data

Restriction site‐associated DNA sequencing (RADseq) is a powerful tool for genotyping of individuals, but the identification of loci and assignment of sequence reads is a crucial and often challenging step. The optimal parameter settings for a given de novo RADseq assembly vary between data sets and can be difficult and computationally expensive to determine. Here, we introduce RADProc, a software package that uses a graph data structure to represent all sequence reads and their similarity relationships. Storing sequence–comparison results in a graph eliminates unnecessary and redundant sequence similarity calculations. De novo locus formation for a given parameter set can be performed on the precomputed graph, making parameter sweeps far more efficient. RADProc also uses a clustering approach for faster nucleotide‐distance calculation. The performance of RADProc compares favourably with that of the widely used Stacks software. The run‐time comparisons between RADProc and Stacks for 32 different parameter settings using 20 green‐crab (Carcinus maenas) samples showed that RADProc took as little as 2 hr 40 min compared to 78 hr by Stacks, while 16 brown trout (Salmo trutta L.) samples were processed by RADProc and Stacks in 23 and 263 hr, respectively. Comparisons of the de novo loci formed, and catalog built using both the methods demonstrate that the improvement in processing speeds achieved by RADProc does not affect much the actual loci formed and the results of downstream analyses based on those loci.     

A genetic screen for mutations affecting gonad formation in Drosophila reveals a role for the slit/robo pathway

Organogenesis is a complex process requiring multiple cell types to associate with one another through correct cell contacts and in the correct location to achieve proper organ morphology and function. To better understand the mechanisms underlying gonad formation, we performed a mutagenesis screen in Drosophila and identified twenty-four genes required for gonadogenesis. These genes affect all different aspects of gonad formation and provide a framework for understanding the molecular mechanisms that control these processes. We find that gonad formation is regulated by multiple, independent pathways; some of these regulate the key cell adhesion molecule DE-cadherin, while others act through distinct mechanisms. In addition, we discover that the Slit/Roundabout pathway, best known for its role in regulating axonal guidance, is essential for proper gonad formation. Our findings shed light on the complexities of gonadogenesis and the genetic regulation required for proper organ formation.     

Evolutionary history inferred from the de novo assembly of a nonmodel organism,the blue‐eyed black lemur

Lemurs, the living primates most distantly related to humans, demonstrate incredible diversity in behaviour, life history patterns and adaptive traits. Although many lemur species are endangered within their native Madagascar, there is no high‐quality genome assembly from this taxon, limiting population and conservation genetic studies. One critically endangered lemur is the blue‐eyed black lemur Eulemur flavifrons. This species is fixed for blue irises, a convergent trait that evolved at least four times in primates and was subject to positive selection in humans, where 5′ regulatory variation of OCA2 explains most of the brown/blue eye colour differences. We built a de novo genome assembly for E. flavifrons, providing the most complete lemur genome to date, and a high confidence consensus sequence for close sister species E. macaco, the (brown‐eyed) black lemur. From diversity and divergence patterns across the genomes, we estimated a recent split time of the two species (160 Kya) and temporal fluctuations in effective population sizes that accord with known environmental changes. By looking for regions of unusually low diversity, we identified potential signals of directional selection in E. flavifrons at MITF, a melanocyte development gene that regulates OCA2 and has previously been associated with variation in human iris colour, as well as at several other genes involved in melanin biosynthesis in mammals. Our study thus illustrates how whole‐genome sequencing of a few individuals can illuminate the demographic and selection history of nonmodel species.     

The N-terminal prodomain of sV23 is essential for the assembly of a functional vitelline membrane network in Drosophila

The vitelline membrane is a specialized extracellular matrix that surrounds and protects the oocyte. Recent studies indicate that it also serves as a storage site for embryonic pattern determinants. sV23, a major vitelline membrane protein, is essential for the morphogenesis of the vitelline membrane as sV23 protein null mutants lay flaccid, infertile eggs. By analyzing a series of sV23 mutant transgenes in the sV23 protein null genetic background, we have shown that sV23 is secreted as a proprotein in functional excess and that C- and N-terminal prodomains are removed successively, following its deposition in the extracellular space. Although a target site for subtilisin-like convertases is essential for N-terminal processing, N-terminal processing is not necessary for the assembly of a functional vitelline membrane layer. While C-terminal truncations were tolerated, the removal of N-terminal sequences lead to the production of flaccid, infertile eggs with a soluble, rather than insoluble, vitelline membrane network. We propose that the hydrophobic N-terminal prodomain plays an early and essential role in aligning molecules within the vitelline membrane network, much like hydrophobic domains within elastin drive the assembly and alignment of molecules within elastin-based extracellular matrices.     

The Drosophila pericentrin-like protein is essential for cilia/flagella function, but appears to be dispensable for mitosis

Centrosomes consist of a pair of centrioles surrounded by an amorphous pericentriolar material (PCM). Proteins that contain a Pericentrin/AKAP450 centrosomal targeting (PACT) domain have been implicated in recruiting several proteins to the PCM. We show that the only PACT domain protein in Drosophila (the Drosophila pericentrin-like protein [D-PLP]) is associated with both the centrioles and the PCM, and is essential for the efficient centrosomal recruitment of all six PCM components that we tested. Surprisingly, however, all six PCM components are eventually recruited to centrosomes during mitosis in d-plp mutant cells, and mitosis is not dramatically perturbed. Although viable, d-plp mutant flies are severely uncoordinated, a phenotype usually associated with defects in mechanosensory neuron function. We show that the sensory cilia of these neurons are malformed and the neurons are nonfunctional in d-plp mutants. Moreover, the flagella in mutant sperm are nonmotile. Thus, D-PLP is essential for the formation of functional cilia and flagella in flies.     

Drosophila cells and ecdysterone: A model system for gene regulation

Summary WhenDrosophila cell lines are exposed to physiological doses of the steroid molting hormone, ecdysterone, they enter mitotic arrest and differentiate morphologically. These responses are accompanied by specific changes in gene expression. Several enzyme activities (acetylcholinesterase, β-galactosidase, dopa decarboxylase, and catalase) are induced and the synthesis of a cytoplasmic actin and the four small heat-shock proteins is initiated. Several of these ecdysterone inducible genes have been physically isolated and characterized, in several cases by DNA sequencing. Current studies focus on introducing cloned ecdysterone inducible genes into responsive cells by DNA mediated transfection. Once it is clear that these introduced genes acquire the normal pattern of hormone-regulated gene expression in the cell, in vitro mutagenesis can be used before transfection to modify their structure. Transient expression, then, can become a functional assay to define regions of DNA flanking the coding region of inducible genes that are needed for proper gene expression and regulation in cultured cells. This work has been supported by grants from the NIH (GM 22866, GM 33235, CA 23108) and the American Cancer Society (1N157).     

The Drosophila gene zfh2 is required to establish proximal-distal domains in the wing disc

Three main events characterize the development of the proximal-distal axis of the Drosophila wing disc: first, generation of nested circular domains defined by different combinations of gene expression; second, activation of wingless (wg) gene expression in a ring of cells; and third, an increase of cell number in each domain in response to Wg. The mechanisms by which these domains of gene expression are established and maintained are unknown. We have analyzed the role of the gene zinc finger homeodomain 2 (zfh2). We report that in discs lacking zfh2 the limits of the expression domains of the genes tsh, nub, rn, dve and nab coincide, and expression of wg in the wing hinge, is lost. We show that zfh2 expression is delimited distally by Vg, Nub and Dpp signalling, and proximally by Tsh and Dpp. Distal repression of zfh2 permits activation of nab in the wing blade and wg in the wing hinge. We suggest that the proximal-most wing fate, the hinge, is specified first and that later repression of zfh2 permits specification of the distal-most fate, the wing blade. We propose that proximal-distal axis development is achieved by a combination of two strategies: on one hand a process involving proximal to distal specification, with the wing hinge specified first followed later by the distal wing blade; on the other hand, early specification of the proximal-distal domains by different combinations of gene expression. The results we present here indicate that Zfh2 plays a critical role in both processes.     

Moleskin is essential for the formation of the myotendinous junction in Drosophila

It is the precise connectivity between skeletal muscles and their corresponding tendon cells to form a functional myotendinous junction (MTJ) that allows for the force generation required for muscle contraction and organismal movement. The Drosophila MTJ is composed of secreted extracellular matrix (ECM) proteins deposited between integrin-mediated hemi-adherens junctions on the surface of muscle and tendon cells. In this paper, we have identified a novel, cytoplasmic role for the canonical nuclear import protein Moleskin (Msk) in Drosophila embryonic somatic muscle attachment. Msk protein is enriched at muscle attachment sites in late embryogenesis and msk mutant embryos exhibit a failure in muscle–tendon cell attachment. Although the muscle–tendon attachment sites are reduced in size, components of the integrin complexes and ECM proteins are properly localized in msk mutant embryos. However, msk mutants fail to localize phosphorylated focal adhesion kinase (pFAK) to the sites of muscle–tendon cell junctions. In addition, the tendon cell specific proteins Stripe (Sr) and activated mitogen-activated protein kinase (MAPK) are reduced in msk mutant embryos. Our rescue experiments demonstrate that Msk is required in the muscle cell, but not in the tendon cells. Moreover, muscle attachment defects due to loss of Msk are rescued by an activated form of MAPK or the secreted epidermal growth factor receptor (Egfr) ligand Vein. Taken together, these findings provide strong evidence that Msk signals non-autonomously through the Vein-Egfr signaling pathway for late tendon cell late differentiation and/or maintenance.     

Implementation of a de novo genome-wide computational approach for updating Brachypodium miRNAs

Plant microRNAs (miRNAs) are single-stranded 20-22 nt small RNAs (sRNA) that are produced from their own genes. We have developed a de novo genome-wide approach for the computational identification of novel plant miRNAs based on the integration of the complete genome sequence with sRNA libraries. It comprises three modules - the clustering module identifies genomic regions that have two closely-located unidirectional sRNA clusters, the mirplan module explores the secondary structure of the genomic regions, and the duplex module predicts miRNA/miRNA* duplexes. We applied our approach to the Brachypodium genome and publicly available sRNA libraries and predicted 102 miRNAs. Our results extend the list of known miRNAs with 58 novel miRNAs and define the genomic loci of all predicted miRNAs. Because this approach considers specific features of plant miRNAs, it can be employed for the analysis of the genome and sRNA libraries generated for plant species to achieve systematic miRNA discovery.     

Restriction site‐associated DNA sequencing,genotyping error estimation and de novo assembly optimization for population genetic inference

Restriction site‐associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single‐nucleotide polymorphisms. As an empirical example, we use a double‐digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high‐altitude mountains in Mexico.     

spalt-induced specification of distinct dorsal and ventral domains is required for Drosophila tracheal patterning

Morphogenesis of the Drosophila tracheal system relies on different signalling pathways that have distinct roles in specifying both the migration of the tracheal cells and the particular morphological features of the primary branches. The current view is that the tracheal cells are initially specified as an equivalent group of cells whose diversification depends on signals from the surrounding cells. In this work, we show that the tracheal primordia are already specified as distinct dorsal and ventral cell populations. This subdivision depends on the activity of the spalt (sal) gene and occurs prior to the activity of the signalling pathways that dictate the development of the primary branches. Finally, we show that the specification of these two distinct cell populations, which are not defined by cell lineage, are critical for proper tracheal patterning. These results indicate that tracheal patterning depends not only on signalling from surrounding cells but also in the different response of the tracheal cells depending on their allocation to the dorsal or ventral domains.     

Genome assembly and annotation of a Drosophila simulans strain from Madagascar

Drosophila simulans is a close relative of the genetic model D. melanogaster. Its worldwide distribution in combination with the absence of segregating chromosomal inversions makes this species an increasingly attractive model to study the molecular signatures of adaptation in natural and experimental populations. In an effort to improve the genomic resources for D. simulans, we assembled and annotated the genome of a strain originating from Madagascar (M252), the ancestral range of D. simulans. The comparison of the M252 genome to other available D. simulans assemblies confirmed its high quality, but also highlighted genomic regions that are difficult to assemble with NGS data. The annotation of M252 provides a clear improvement with alternative splicing for 52% of the multiple‐exon genes, UTRs for 70% of the genes, 225 novel genes and 781 pseudogenes being reported. We anticipate that the M252 genome will be a valuable resource for many research questions.