ALBERT

Description: Genomic structural variation (SV), a common hallmark of cancer, has important predictive and therapeutic implications. However, accurately detecting SV using high-throughput sequencing data remains challenging, especially for ‘targeted’ resequencing efforts. This is critically important in the clinical setting where targeted resequencing is frequently being applied to rapidly assess clinically actionable mutations in tumor biopsies in a cost-effective manner. We present BreaKmer, a novel approach that uses a ‘kmer’ strategy to assemble misaligned sequence reads for predicting insertions, deletions, inversions, tandem duplications and translocations at base-pair resolution in targeted resequencing data. Variants are predicted by realigning an assembled consensus sequence created from sequence reads that were abnormally aligned to the reference genome. Using targeted resequencing data from tumor specimens with orthogonally validated SV, non-tumor samples and whole-genome sequencing data, BreaKmer had a 97.4% overall sensitivity for known events and predicted 17 positively validated, novel variants. Relative to four publically available algorithms, BreaKmer detected SV with increased sensitivity and limited calls in non-tumor samples, key features for variant analysis of tumor specimens in both the clinical and research settings.

Description: Recent studies on microRNA (miRNA) evolution focused mainly on the comparison of miRNA complements between animal clades. However, evolution of miRNAs within such groups is poorly explored despite the availability of comparable data that in some cases lack only a few key taxa. For flatworms (Platyhelminthes), miRNA complements are available for some free-living flatworms and all major parasitic lineages, except for the Monogenea. We present the miRNA complement of the monogenean flatworm Gyrodactylus salaris that facilitates a comprehensive analysis of miRNA evolution in Platyhelminthes. Using the newly designed bioinformatics pipeline miRCandRef, the miRNA complement was disentangled from next-generation sequencing of small RNAs and genomic DNA without a priori genome assembly. It consists of 39 miRNA hairpin loci of conserved miRNA families, and 22 novel miRNAs. A comparison with the miRNA complements of Schmidtea mediterranea (Turbellaria), Schistosoma japonicum (Trematoda), and Echinococcus granulosus (Cestoda) reveals a substantial loss of conserved bilaterian, protostomian, and lophotrochozoan miRNAs. Eight of the 46 expected conserved miRNAs were lost in all flatworms, 16 in Neodermata and 24 conserved miRNAs could not be detected in the cestode and the trematode. Such a gradual loss of miRNAs has not been reported before for other animal phyla. Currently, little is known about miRNAs in Platyhelminthes, and for the majority of the lost miRNAs there is no prediction of function. As suggested earlier they might be related to morphological simplifications. The presence and absence of 153 conserved miRNAs was compared for platyhelminths and 32 other metazoan taxa. Phylogenetic analyses support the monophyly of Platyhelminthes (Turbellaria + Neodermata [Monogenea {Trematoda + Cestoda}]).

Description: We perform a suite of multimass cosmological zoom simulations of individual dark matter haloes and explore how to best select Lagrangian regions for resimulation without contaminating the halo of interest with low-resolution particles. Such contamination can lead to significant errors in the gas distribution of hydrodynamical simulations, as we show. For a fixed Lagrange volume, we find that the chance of contamination increases systematically with the level of zoom. In order to avoid contamination, the Lagrangian volume selected for resimulation must increase monotonically with the resolution difference between parent box and the zoom region. We provide a simple formula for selecting Lagrangian regions (in units of the halo virial volume) as a function of the level of zoom required. We also explore the degree to which a halo's Lagrangian volume correlates with other halo properties (concentration, spin, formation time, shape, etc.) and find no significant correlation. There is a mild correlation between Lagrange volume and environment, such that haloes living in the most clustered regions have larger Lagrangian volumes. Nevertheless, selecting haloes to be isolated is not the best way to ensure inexpensive zoom simulations. We explain how one can safely choose haloes with the smallest Lagrangian volumes, which are the least expensive to resimulate, without biasing one's sample.

Description: Convergent evolution occurs when the same trait arises independently in multiple lineages. In most cases of phenotypic convergence such transitions are adaptive, so finding the underlying molecular causes of convergence can provide insight into the process of adaptation. Convergent evolution at the genomic level also lends itself to study by comparative methods, although molecular convergence can also occur by chance, adding noise to this process. Parker et al. studied convergence across the genomes of several mammals, including echolocating bats and dolphins (Parker J, Tsagkogeorga G, Cotton JA, Liu Y, Provero P, Stupka E, Rossiter SJ. 2013. Genome-wide signatures of convergent evolution in echolocating mammals. Nature 502:228–231). On the basis of a null distribution of site-specific likelihood support (SSLS) generated using simulated topologies, they concluded that there was evidence for genome-wide adaptive convergence between echolocating taxa. Here, we demonstrate that methods based on SSLS do not adequately measure convergence, and reiterate the use of an empirical null model that directly compares convergent substitutions between all pairs of species. We find that when the proper comparisons are made there is no surprising excess of convergence between echolocating mammals, even in sensory genes.

Description: Gene duplication is a major driver of organismal adaptation and evolution and plays an important role in multiple human diseases. Whole-genome analyses have shown similar and high rates of gene duplication across a variety of eukaryotic species. Most of these studies, however, did not address the possible impact of interlocus gene conversion (IGC) on the evolution of gene duplicates. Because IGC homogenizes pairs of duplicates, widespread conversion would cause gene duplication events that happened long ago to appear more recent, resulting in artificially high estimates of duplication rates. Although the majority of genome-wide studies (including in the budding yeast Saccharomyces cerevisiae [Scer]) point to levels of IGC between paralogs ranging from 2% to 18%, Gao and Innan (Gao LZ, Innan H. 2004. Very low gene duplication rate in the yeast genome. Science 306:1367–1370.) found that gene conversion in yeast affected 〉80% of paralog pairs. If conversion rates really are this high, it would imply that the rate of gene duplication in eukaryotes is much lower than previously reported. In this work, we apply four different methodologies—including one approach that closely mirrors Gao and Innan’s method—to estimate the level of IGC in Scer. Our analyses point to a maximum conversion level of 13% between paralogs in this species, in close agreement with most estimates of IGC in eukaryotes. We also show that the exceedingly high levels of conversion found previously derive from application of an accurate method to an inappropriate data set. In conclusion, our work provides the most striking evidence to date supporting the reduced incidence of IGC among Scer paralogs and sets up a framework for future analyses in other eukaryotes.

Description: Deleterious mutations contribute to polymorphism even when selection effectively prevents their fixation. The efficacy of selection in removing deleterious mitochondrial mutations from populations depends on the effective population size ( N e ) of the mitochondrial DNA and the degree to which a lack of recombination magnifies the effects of linked selection. Using complete mitochondrial genomes from Drosophila melanogaster and nuclear data available from the same samples, we reexamine the hypothesis that nonrecombining animal mitochondrial DNA harbor an excess of deleterious polymorphisms relative to the nuclear genome. We find no evidence of recombination in the mitochondrial genome, and the much-reduced level of mitochondrial synonymous polymorphism relative to nuclear genes is consistent with a reduction in N e . Nevertheless, we find that the neutrality index, a measure of the excess of nonsynonymous polymorphism relative to the neutral expectation, is only weakly significantly different between mitochondrial and nuclear loci. This difference is likely the result of the larger proportion of beneficial mutations in X-linked relative to autosomal loci, and we find little to no difference between mitochondrial and autosomal neutrality indices. Reanalysis of published data from Homo sapiens reveals a similar lack of a difference between the two genomes, although previous studies have suggested a strong difference in both species. Thus, despite a smaller N e , mitochondrial loci of both flies and humans appear to experience similar efficacies of purifying selection as do loci in the recombining nuclear genome.

Description: The evolution of a pair of chromosomes that differ in appearance between males and females (heteromorphic sex chromosomes) has occurred repeatedly across plants and animals. Recent work has shown that the male heterogametic (XY) and female heterogametic (ZW) sex chromosomes evolved independently from different pairs of homomorphic autosomes in the common ancestor of birds and mammals but also that X and Z chromosomes share many convergent molecular features. However, little is known about how often heteromorphic sex chromosomes have either evolved convergently from different autosomes or in parallel from the same pair of autosomes and how universal patterns of molecular evolution on sex chromosomes really are. Among winged insects with sequenced genomes, there are male heterogametic species in both the Diptera (e.g., Drosophila melanogaster ) and the Coleoptera ( Tribolium castaneum ), female heterogametic species in the Lepidoptera ( Bombyx mori ), and haplodiploid species in the Hymenoptera (e.g., Nasonia vitripennis ). By determining orthologous relationships among genes on the X and Z chromosomes of insects with sequenced genomes, we are able to show that these chromosomes are not homologous to one another but are homologous to autosomes in each of the other species. These results strongly imply that heteromorphic sex chromosomes have evolved independently from different pairs of ancestral chromosomes in each of the insect orders studied. We also find that the convergently evolved X chromosomes of Diptera and Coleoptera share genomic features with each other and with vertebrate X chromosomes, including excess gene movement from the X to the autosomes. However, other patterns of molecular evolution—such as increased codon bias, decreased gene density, and the paucity of male-biased genes on the X—differ among the insect X and Z chromosomes. Our results provide evidence for both differences and nearly universal similarities in patterns of evolution among independently derived sex chromosomes.

Description: Substitution rates vary between species, and many explanations regarding the causes of this variation have been proposed. Here we consider how new genomic data on the per-generation mutation rate impinge on proposed hypotheses for substitution rate variation in primates. We propose that the generation-time effect as it is usually understood cannot explain the observed rate variation, but instead that selection for decreased somatic mutation rates can. By considering the disparate causes underlying mutation rate changes in recent human history, we also show that the per-generation mutation rate is increasing even as the per-cell-division rate is decreasing.

Description: Current de novo whole-genome sequencing approaches often are inadequate for organisms lacking substantial preexisting genetic data. Problems with these methods are manifest as: large numbers of scaffolds that are not ordered within chromosomes or assigned to individual chromosomes, misassembly of allelic sequences as separate loci when the individual(s) being sequenced are heterozygous, and the collapse of recently duplicated sequences into a single locus, regardless of levels of heterozygosity. Here we propose a new approach for producing de novo whole-genome sequences—which we call recombinant population genome construction—that solves many of the problems encountered in standard genome assembly and that can be applied in model and nonmodel organisms. Our approach takes advantage of next-generation sequencing technologies to simultaneously barcode and sequence a large number of individuals from a recombinant population. The sequences of all recombinants can be combined to create an initial de novo assembly, followed by the use of individual recombinant genotypes to correct assembly splitting/collapsing and to order and orient scaffolds within linkage groups. Recombinant population genome construction can rapidly accelerate the transformation of nonmodel species into genome-enabled systems by simultaneously producing a high-quality genome assembly and providing genomic tools ( e.g. , high-confidence single-nucleotide polymorphisms) for immediate applications. In populations segregating for important functional traits, this approach also enables simultaneous mapping of quantitative trait loci. We demonstrate our method using simulated Illumina data from a recombinant population of Caenorhabditis elegans and show that the method can produce a high-fidelity, high-quality genome assembly for both parents of the cross.

Description: Genetic differentiation across populations that is maintained in the presence of gene flow is a hallmark of spatially varying selection. In Drosophila melanogaster , the latitudinal clines across the eastern coasts of Australia and North America appear to be examples of this type of selection, with recent studies showing that a substantial portion of the D. melanogaster genome exhibits allele frequency differentiation with respect to latitude on both continents. As of yet there has been no genome-wide examination of differentiated copy-number variants (CNVs) in these geographic regions, despite their potential importance for phenotypic variation in Drosophila and other taxa. Here, we present an analysis of geographic variation in CNVs in D. melanogaster. We also present the first genomic analysis of geographic variation for copy-number variation in the sister species, D. simulans , in order to investigate patterns of parallel evolution in these close relatives. In D. melanogaster we find hundreds of CNVs, many of which show parallel patterns of geographic variation on both continents, lending support to the idea that they are influenced by spatially varying selection. These findings support the idea that polymorphic CNVs contribute to local adaptation in D. melanogaster . In contrast, we find very few CNVs in D. simulans that are geographically differentiated in parallel on both continents, consistent with earlier work suggesting that clinal patterns are weaker in this species.