Gene content evolution in the arthropods

Arthropods comprise the largest and most diverse phylum on Earth and play vital roles in nearly every ecosystem. Their diversity stems in part from variations on a conserved body plan, resulting from and recorded in adaptive changes in the genome. Dissection of the genomic record of sequence change enables broad questions regarding genome evolution to be addressed, even across hyper-diverse taxa within arthropods. Using 76 whole genome sequences representing 21 orders spanning more than 500 million years of arthropod evolution, we document changes in gene and protein domain content and provide temporal and phylogenetic context for interpreting these innovations. We identify many novel gene families that arose early in the evolution of arthropods and during the diversification of insects into modern orders. We reveal unexpected variation in patterns of DNA methylation across arthropods and examples of gene family and protein domain evolution coincident with the appearance of notable phenotypic and physiological adaptations such as flight, metamorphosis, sociality, and chemoperception. These analyses demonstrate how large-scale comparative genomics can provide broad new insights into the genotype to phenotype map and generate testable hypotheses about the evolution of animal diversity

Study of the rate and spectrum of spontaneous mutations

Mutations are the initial force responsible for all aspects of genetic variation, and are a central part to evolution in all organisms. Yet despite its importance, the previously high cost that is associated with surveying mutations at a genome-wide scale has limited the understanding of the mutation process in eukaryotes. However, recent high-throughput sequencing technology has greatly reduced the cost of surveying mutations. By applying high-throughput sequencing to mutation accumulation experiments, we have begun to characterize the genome-wide mutation spectrum of eukaryotes. Across all eukaryotes, we observe a biased rate of G/C-\u3e A/T mutations that exceeds the number of A/T-\u3eG/C mutations. This finding is consistent with spontaneous deamination of cytosine or methylated cytosine. Alternate forces such as selection or G/C biased gene conversion must be driving eukaryotic genomes toward a higher G/C composition than expected from mutation bias. In Paramecium tetraurelia, we observe a nuclear mutation rate ∼75 fold lower than previously expected. When the base substitution rate per generation is extrapolated to the rate per expressed sexual cycle, it is equivalent to that observed in multicellular species with comparable genome sizes. This suggests that natural selection operates at germline expression, and favors a minimum rate that opposes random genetic drift. Using a natural population of Daphnia pulex we catalogue simple sequence repeats (SSR) and determine average heterozygosity of each SSR type. We find that SSR heterozygosity is motif specific, and positively correlated with repeat number as well as motif length. We identify a motif-dependent end-nucleotide polymorphism bias. Our observations also confirm the high frequency of multiple unit variation at large microsatellite loci. We observe that structural variants in C. elegans and S. cerevisiae, which are ∼1000 fold larger than base substitution rates on a per nucleotide basis, occur on the same order of magnitude as base substitutions. The rate and direction of structural gains and losses differ between yeast and C. elegans, and we hypothesize that the rate of structural variants corresponds with the coding portion of the genome. We also confirm a high rate of gene inversion and gene loss in the life history of C. elegans

Genome of the Avirulent Human-Infective Trypanosome—Trypanosoma rangeli

Background: Trypanosoma rangeli is a hemoflagellate protozoan parasite infecting humans and other wild and domestic mammals across Central and South America. It does not cause human disease, but it can be mistaken for the etiologic agent of Chagas disease, Trypanosoma cruzi. We have sequenced the T. rangeli genome to provide new tools for elucidating the distinct and intriguing biology of this species and the key pathways related to interaction with its arthropod and mammalian hosts.  Methodology/Principal Findings: The T. rangeli haploid genome is ,24 Mb in length, and is the smallest and least repetitive trypanosomatid genome sequenced thus far. This parasite genome has shorter subtelomeric sequences compared to those of T. cruzi and T. brucei; displays intraspecific karyotype variability and lacks minichromosomes. Of the predicted 7,613 protein coding sequences, functional annotations could be determined for 2,415, while 5,043 are hypothetical proteins, some with evidence of protein expression. 7,101 genes (93%) are shared with other trypanosomatids that infect humans. An ortholog of the dcl2 gene involved in the T. brucei RNAi pathway was found in T. rangeli, but the RNAi machinery is non-functional since the other genes in this pathway are pseudogenized. T. rangeli is highly susceptible to oxidative stress, a phenotype that may be explained by a smaller number of anti-oxidant defense enzymes and heatshock proteins.  Conclusions/Significance: Phylogenetic comparison of nuclear and mitochondrial genes indicates that T. rangeli and T. cruzi are equidistant from T. brucei. In addition to revealing new aspects of trypanosome co-evolution within the vertebrate and invertebrate hosts, comparative genomic analysis with pathogenic trypanosomatids provides valuable new information that can be further explored with the aim of developing better diagnostic tools and/or therapeutic targets

De novo assembly of the olive fruit fly (Bactrocera oleae) genome with linked-reads and long-read technologies minimizes gaps and provides exceptional Y chromosome assembly

Background: The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly’s biology and proposing alternative control methods to pesticide use. Results: Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR. Conclusions: The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome’s organization, function and evolution and is poised to provide avenues for sterile insect technique approaches

An integrated 4249 marker FISH/RH map of the canine genome

BACKGROUND: The 156 breeds of dog recognized by the American Kennel Club offer a unique opportunity to map genes important in genetic variation. Each breed features a defining constellation of morphological and behavioral traits, often generated by deliberate crossing of closely related individuals, leading to a high rate of genetic disease in many breeds. Understanding the genetic basis of both phenotypic variation and disease susceptibility in the dog provides new ways in which to dissect the genetics of human health and biology. RESULTS: To facilitate both genetic mapping and cloning efforts, we have constructed an integrated canine genome map that is both dense and accurate. The resulting resource encompasses 4249 markers, and was constructed using the RHDF5000-2 whole genome radiation hybrid panel. The radiation hybrid (RH) map features a density of one marker every 900 Kb and contains 1760 bacterial artificial chromosome clones (BACs) localized to 1423 unique positions, 851 of which have also been mapped by fluorescence in situ hybridization (FISH). The two data sets show excellent concordance. Excluding the Y chromosome, the map features an RH/FISH mapped BAC every 3.5 Mb and an RH mapped BAC-end, on average, every 2 Mb. For 2233 markers, the orthologous human genes have been established, allowing the identification of 79 conserved segments (CS) between the dog and human genomes, dramatically extending the length of most previously described CS. CONCLUSIONS: These results provide a necessary resource for the canine genome mapping community to undertake positional cloning experiments and provide new insights into the comparative canine-human genome maps

High-Throughput Sequencing of Three Lemnoideae (Duckweeds) Chloroplast Genomes from Total DNA

BACKGROUND: Chloroplast genomes provide a wealth of information for evolutionary and population genetic studies. Chloroplasts play a particularly important role in the adaption for aquatic plants because they float on water and their major surface is exposed continuously to sunlight. The subfamily of Lemnoideae represents such a collection of aquatic species that because of photosynthesis represents one of the fastest growing plant species on earth. METHODS: We sequenced the chloroplast genomes from three different genera of Lemnoideae, Spirodela polyrhiza, Wolffiella lingulata and Wolffia australiana by high-throughput DNA sequencing of genomic DNA using the SOLiD platform. Unfractionated total DNA contains high copies of plastid DNA so that sequences from the nucleus and mitochondria can easily be filtered computationally. Remaining sequence reads were assembled into contiguous sequences (contigs) using SOLiD software tools. Contigs were mapped to a reference genome of Lemna minor and gaps, selected by PCR, were sequenced on the ABI3730xl platform. CONCLUSIONS: This combinatorial approach yielded whole genomic contiguous sequences in a cost-effective manner. Over 1,000-time coverage of chloroplast from total DNA were reached by the SOLiD platform in a single spot on a quadrant slide without purification. Comparative analysis indicated that the chloroplast genome was conserved in gene number and organization with respect to the reference genome of L. minor. However, higher nucleotide substitution, abundant deletions and insertions occurred in non-coding regions of these genomes, indicating a greater genomic dynamics than expected from the comparison of other related species in the Pooideae. Noticeably, there was no transition bias over transversion in Lemnoideae. The data should have immediate applications in evolutionary biology and plant taxonomy with increased resolution and statistical power

De novo draft assembly of the Botrylloides leachii genome provides further insight into tunicate evolution

Tunicates are marine invertebrates that compose the closest phylogenetic group to the vertebrates. These chordates present a particularly diverse range of regenerative abilities and life-history strategies. Consequently, tunicates provide an extraordinary perspective into the emergence and diversity of these traits. Here we describe the genome sequencing, annotation and analysis of the Stolidobranchian Botrylloides leachii. We have produced a high-quality 159 Mb assembly, 82% of the predicted 194  Mb genome. Analysing genome size, gene number, repetitive elements, orthologs clustering and gene ontology terms show that B. leachii has a genomic architecture similar to that of most solitary tunicates, while other recently sequenced colonial ascidians have undergone genome expansion. In addition, ortholog clustering has identified groups of candidate genes for the study of colonialism and whole-body regeneration. By analysing the structure and composition of conserved gene linkages, we observed examples of cluster breaks and gene dispersions, suggesting that several lineage-specific genome rearrangements occurred during tunicate evolution. We also found lineage-specific gene gain and loss within conserved cell-signalling pathways. Such examples of genetic changes within conserved cell-signalling pathways commonly associated with regeneration and development that may underlie some of the diverse regenerative abilities observed in tunicates. Overall, these results provide a novel resource for the study of tunicates and of colonial ascidians