Phylogenetic taxonomy in Drosophila: Problems and prospects

The genus Drosophila is one of the best-studied model systems in modern biology, with twelve fully sequenced genomes available. In spite of the large number of genetic and genomic resources, little is known concerning the phylogenetic relationships, ecology, and evolutionary history of all but a few species. Recent molecular systematic studies have shown that this genus is comprised of at least three independent lineages and that several other genera are actually embedded within Drosophila. This genus accounts for over 2000 described, and many more undescribed, species. While some Drosophila researchers are advocating dividing this genus into three or more separate genera, others favor maintaining Drosophila as a single large genus. With the recent sequencing of the genomes of multiple Drosophila species and their expanding use in comparative biology, it is critical that the Drosophila research community understands the taxonomic framework underlying the naming and relationships of these species. The subdivision of this genus has significant biological implications, ranging from the accurate annotation of single genes to understanding how ecological adaptations have occurred over the history of the group.     

A physical map of the bovine genome

Background

Cattle are important agriculturally and relevant as a model organism. Previously described genetic and radiation hybrid (RH) maps of the bovine genome have been used to identify genomic regions and genes affecting specific traits. Application of these maps to identify influential genetic polymorphisms will be enhanced by integration with each other and with bacterial artificial chromosome (BAC) libraries. The BAC libraries and clone maps are essential for the hybrid clone-by-clone/whole-genome shotgun sequencing approach taken by the bovine genome sequencing project.

Results

A bovine BAC map was constructed with HindIII restriction digest fragments of 290,797 BAC clones from animals of three different breeds. Comparative mapping of 422,522 BAC end sequences assisted with BAC map ordering and assembly. Genotypes and pedigree from two genetic maps and marker scores from three whole-genome RH panels were consolidated on a 17,254-marker composite map. Sequence similarity allowed integrating the BAC and composite maps with the bovine draft assembly (Btau3.1), establishing a comprehensive resource describing the bovine genome. Agreement between the marker and BAC maps and the draft assembly is high, although discrepancies exist. The composite and BAC maps are more similar than either is to the draft assembly.

Conclusion

Further refinement of the maps and greater integration into the genome assembly process may contribute to a high quality assembly. The maps provide resources to associate phenotypic variation with underlying genomic variation, and are crucial resources for understanding the biology underpinning this important ruminant species so closely associated with humans.     

Construction and analysis of bacterial artificial chromosome libraries from a marine microbial assemblage

Cultivation-independent surveys of ribosomal RNA genes have revealed the existence of novel microbial lineages, many with no known cultivated representatives. Ribosomal RNA-based analyses, however, often do not provide significant information beyond phylogenetic affiliation. Analysis of large genome fragments recovered directly from microbial communities represents one promising approach for characterizing uncultivated microbial species better. To assess further the utility of this approach, we constructed large-insert bacterial artificial chromosome (BAC) libraries from the genomic DNA of planktonic marine microbial assemblages. The BAC libraries we prepared had average insert sizes of 80 kb, with maximal insert sizes > 150 kb. A rapid screening method assessing the phylogenetic diversity and representation in the library was developed and applied. In general, representation in the libraries agreed well with previous culture-independent surveys based on polymerase chain reaction (PCR)amplified rRNA fragments. A significant fraction of the genome fragments in the BAC libraries originated from as yet uncultivated microbial species, thought to be abundant and widely distributed in the marine environment. One entire BAC insert, derived from an uncultivated, surface-dwelling euryarchaeote, was sequenced completely. The planktonic euryarchaeal genome fragment contained some typical archaeal genes, as well as unique open reading frames (ORFs) suggesting novel function. In total, our results verify the utility of BAC libraries for providing access to the genomes of as yet uncultivated microbial species. Further analysis of these BAC libraries has the potential to provide significant insight into the genomic potential and ecological roles of many indigenous microbial species, cultivated or not.     

Construction, characterization and FISH mapping of a bacterial artificial chromosome library of Chinese pangolin (Manis pentadactyla)

Chinese pangolins as a representative species in the order Pholidota have highly specified morphological characters and occupy an important place in the mammalian phylogenetic tree. To obtain genomic data for this species, we have constructed a bacterial artificial chromosome (BAC) library of Chinese pangolin. The library contains 208,272 clones with an average insert size of 122.1 kb and represents approximately eight times the Chinese pangolin haploid genome (if we assume that the Chinese pangolins have a genome size similar to human). One hundred and twenty randomly-selected BAC clones were mapped onto Chinese pangolin chromosomes by fluorescence in situ hybridization (FISH), showing a largely unbiased chromosomal distribution. Several clones containing repetitive DNA and ribosomal DNA genes were also found. The BAC library and FISH mapped BAC clones are useful resources for comparative genomics and cytogenetics of mammals and in particular, the ongoing genome sequencing project of Chinese pangolins.     

Using comparative genomics to reorder the human genome sequence into a virtual sheep genome

Background

Is it possible to construct an accurate and detailed subgene-level map of a genome using bacterial artificial chromosome (BAC) end sequences, a sparse marker map, and the sequences of other genomes?

Results

A sheep BAC library, CHORI-243, was constructed and the BAC end sequences were determined and mapped with high sensitivity and low specificity onto the frameworks of the human, dog, and cow genomes. To maximize genome coverage, the coordinates of all BAC end sequence hits to the cow and dog genomes were also converted to the equivalent human genome coordinates. The 84,624 sheep BACs (about 5.4-fold genome coverage) with paired ends in the correct orientation (tail-to-tail) and spacing, combined with information from sheep BAC comparative genome contigs (CGCs) built separately on the dog and cow genomes, were used to construct 1,172 sheep BAC-CGCs, covering 91.2% of the human genome. Clustered non-tail-to-tail and outsize BACs located close to the ends of many BAC-CGCs linked BAC-CGCs covering about 70% of the genome to at least one other BAC-CGC on the same chromosome. Using the BAC-CGCs, the intrachromosomal and interchromosomal BAC-CGC linkage information, human/cow and vertebrate synteny, and the sheep marker map, a virtual sheep genome was constructed. To identify BACs potentially located in gaps between BAC-CGCs, an additional set of 55,668 sheep BACs were positioned on the sheep genome with lower confidence. A coordinate conversion process allowed us to transfer human genes and other genome features to the virtual sheep genome to display on a sheep genome browser.

Conclusion

We demonstrate that limited sequencing of BACs combined with positioning on a well assembled genome and integrating locations from other less well assembled genomes can yield extensive, detailed subgene-level maps of mammalian genomes, for which genomic resources are currently limited.     

Dissecting plant genomes with the PLAZA comparative genomics platform

With the arrival of low-cost, next-generation sequencing, a multitude of new plant genomes are being publicly released, providing unseen opportunities and challenges for comparative genomics studies. Here, we present PLAZA 2.5, a user-friendly online research environment to explore genomic information from different plants. This new release features updates to previous genome annotations and a substantial number of newly available plant genomes as well as various new interactive tools and visualizations. Currently, PLAZA hosts 25 organisms covering a broad taxonomic range, including 13 eudicots, five monocots, one lycopod, one moss, and five algae. The available data consist of structural and functional gene annotations, homologous gene families, multiple sequence alignments, phylogenetic trees, and colinear regions within and between species. A new Integrative Orthology Viewer, combining information from different orthology prediction methodologies, was developed to efficiently investigate complex orthology relationships. Cross-species expression analysis revealed that the integration of complementary data types extended the scope of complex orthology relationships, especially between more distantly related species. Finally, based on phylogenetic profiling, we propose a set of core gene families within the green plant lineage that will be instrumental to assess the gene space of draft or newly sequenced plant genomes during the assembly or annotation phase.     

Construction of a BAC library and its application to the identification of simple sequence repeats in peach [ Prunus persica (L.) Batsch

The Rosaceae contains many economically important crop species, but their genomes are not well characterized, and comparative genetic mapping lags well behind that of other families. To facilitate genome comparisons and gene discovery in the Rosaceae, we have begun the development of genomic resources for peach as the model genome for this family. First, we developed a simplified, cost-effective method for constructing BAC libraries, particularly appropriate for plant species of relatively minor economic importance. Second, we used the library to investigate the abundance and local distribution of simple sequence repeats (SSRs) in peach. Our results indicate that microsatellite loci are locally much more highly abundant than previously estimated, and BAC sequencing results suggest that microsatellite repeats are not randomly distributed within gene-containing regions of the peach genome. This makes it relatively easy to identify SSRs in peach by hybridization to BAC clones, and even by random sequencing of BAC clones, not known a priori to contain SSRs.     

BAC library construction and BAC end sequencing of five Drosophila species: the comparative map with the D. melanogaster genome

We constructed and characterized arrayed bacterial artificial chromosome (BAC) libraries of five Drosophila species (D. melanogaster, D. simulans, D. sechellia, D. auraria, and D. ananassae), which are genetically well characterized in the studies of meiosis, evolution, population genetics, and developmental biology. The BAC libraries comprise 8,000 to 12,500 clones for each species, estimated to cover the most of the genomes. We sequenced both ends of most of these BAC clones with a success rate of 91%. Of these, 53,701 clones consisting of non-repetitive BAC end sequences (BESs) were mapped with reference of the public D. melanogaster genome sequences. The BES mapping estimated that the BAC libraries of D. auraria and D. ananassae covered 47% and 57% of the D. melanogaster genome, respectively, and those of D. melanogaster, D. sechellia, and D. simulans covered 94-97%. The low coverage by BESs of D. auraria and D. ananassae may be due to the high sequence divergence with D. melanogaster. From the comparative BES mapping, 111 possible breakpoints of chromosomal rearrangements were identified in these four species. The breakpoints of the major chromosome rearrangement between D. simulans and D. melanogaster on the third chromosome were determined within 20 kb in 84E and 30 kb in 93E/F. Corresponding breakpoints were also identified in D. sechellia. The BAC clones described here will be an important addition to the Drosophila genomic resources.     

Construction and characterization of bacterial artificial chromosome library of black-handed spider monkey (Ateles geoffroyi).

The large-insert genomic DNA library is a critical resource for genome-wide genetic dissection of target species. We constructed a high-redundancy bacterial artificial chromosome (BAC) library of a New World monkey species, the black-handed spider monkey (Ateles geoffroyi). A total of 193 152 BAC clones were generated in this library. The average insert size of the BAC clones was estimated to be 184.6 kb with the small inserts (50-100 kb) accounting for less than 3% and the non-recombinant clones only 1.2%. Assuming a similar genome size with humans, the spider monkey BAC library has about 11x genome coverage. In addition, by end sequencing of randomly selected BAC clones, we generated 367 sequence tags for the library. When blasted against human genome, they showed a good correlation between the number of hit clones and the size of the chromosomes, an indication of unbiased chromosomal distribution of the library. This black-handed spider monkey BAC library would serve as a valuable resource in comparative genomic study and large-scale genome sequencing of nonhuman primates.     

Construction of a bacterial artificial chromosome library from the spikemoss Selaginella moellendorffii: a new resource for plant comparative genomics

Background

The lycophytes are an ancient lineage of vascular plants that diverged from the seed plant lineage about 400 Myr ago. Although the lycophytes occupy an important phylogenetic position for understanding the evolution of plants and their genomes, no genomic resources exist for this group of plants.

Results

Here we describe the construction of a large-insert bacterial artificial chromosome (BAC) library from the lycophyte Selaginella moellendorffii. Based on cell flow cytometry, this species has the smallest genome size among the different lycophytes tested, including Huperzia lucidula, Diphaiastrum digita, Isoetes engelmanii and S. kraussiana. The arrayed BAC library consists of 9126 clones; the average insert size is estimated to be 122 kb. Inserts of chloroplast origin account for 2.3% of the clones. The BAC library contains an estimated ten genome-equivalents based on DNA hybridizations using five single-copy and two duplicated S. moellendorffii genes as probes.

Conclusion

The S. moellenforffii BAC library, the first to be constructed from a lycophyte, will be useful to the scientific community as a resource for comparative plant genomics and evolution.     

Genome size diversity in the family Drosophilidae

Flies in the genus Drosophila have been the dominant model organisms in genetics for over a century and, with a dozen complete sequences now available, continue as such in modern comparative genomics. Surprisingly, estimates of genome size for this genus have been relatively sparse, covering less than 2% of species. Here, best practice flow cytometric genome size estimates are reported for both male and female flies from 67 species from six genera in the family Drosophilidae, including 55 species from the genus Drosophila. Direct and phylogenetically corrected correlation analyses indicate that genome size is positively correlated with temperature-controlled duration of development in Drosophila, and there is indication that genome size may be positively related to body size and sperm length in this genus. These findings may provide some explanation for the streamlined genomes found in these insects, and complement recent work demonstrating possible selective constraints on further deletion of noncoding DNA.     

Generation of Physical Map Contig-Specific Sequences Useful for Whole Genome Sequence Scaffolding

Along with the rapid advances of the nextgen sequencing technologies, more and more species are added to the list of organisms whose whole genomes are sequenced. However, the assembled draft genome of many organisms consists of numerous small contigs, due to the short length of the reads generated by nextgen sequencing platforms. In order to improve the assembly and bring the genome contigs together, more genome resources are needed. In this study, we developed a strategy to generate a valuable genome resource, physical map contig-specific sequences, which are randomly distributed genome sequences in each physical contig. Two-dimensional tagging method was used to create specific tags for 1,824 physical contigs, in which the cost was dramatically reduced. A total of 94,111,841 100-bp reads and 315,277 assembled contigs are identified containing physical map contig-specific tags. The physical map contig-specific sequences along with the currently available BAC end sequences were then used to anchor the catfish draft genome contigs. A total of 156,457 genome contigs (~79% of whole genome sequencing assembly) were anchored and grouped into 1,824 pools, in which 16,680 unique genes were annotated. The physical map contig-specific sequences are valuable resources to link physical map, genetic linkage map and draft whole genome sequences, consequently have the capability to improve the whole genome sequences assembly and scaffolding, and improve the genome-wide comparative analysis as well. The strategy developed in this study could also be adopted in other species whose whole genome assembly is still facing a challenge.     

Comparative Genomic Analysis of Sulfurospirillum cavolei MES Reconstructed from the Metagenome of an Electrosynthetic Microbiome

Sulfurospirillum spp. play an important role in sulfur and nitrogen cycling, and contain metabolic versatility that enables reduction of a wide range of electron acceptors, including thiosulfate, tetrathionate, polysulfide, nitrate, and nitrite. Here we describe the assembly of a Sulfurospirillum genome obtained from the metagenome of an electrosynthetic microbiome. The ubiquity and persistence of this organism in microbial electrosynthesis systems suggest it plays an important role in reactor stability and performance. Understanding why this organism is present and elucidating its genetic repertoire provide a genomic and ecological foundation for future studies where Sulfurospirillum are found, especially in electrode-associated communities. Metabolic comparisons and in-depth analysis of unique genes revealed potential ecological niche-specific capabilities within the Sulfurospirillum genus. The functional similarities common to all genomes, i.e., core genome, and unique gene clusters found only in a single genome were identified. Based upon 16S rRNA gene phylogenetic analysis and average nucleotide identity, the Sulfurospirillum draft genome was found to be most closely related to Sulfurospirillum cavolei. Characterization of the draft genome described herein provides pathway-specific details of the metabolic significance of the newly described Sulfurospirillum cavolei MES and, importantly, yields insight to the ecology of the genus as a whole. Comparison of eleven sequenced Sulfurospirillum genomes revealed a total of 6246 gene clusters in the pan-genome. Of the total gene clusters, 18.5% were shared among all eleven genomes and 50% were unique to a single genome. While most Sulfurospirillum spp. reduce nitrate to ammonium, five of the eleven Sulfurospirillum strains encode for a nitrous oxide reductase (nos) cluster with an atypical nitrous-oxide reductase, suggesting a utility for this genus in reduction of the nitrous oxide, and as a potential sink for this potent greenhouse gas.     

Genome evolution in Reptilia: <Emphasis Type="Italic">in silico</Emphasis> chicken mapping of 12,000 BAC-end sequences from two reptiles and a basal bird

Background

With the publication of the draft chicken genome and the recent production of several BAC clone libraries from non-avian reptiles and birds, it is now possible to undertake more detailed comparative genomic studies in Reptilia. Of interest in particular are the genomic events that transformed the large, repeat-rich genomes of mammals and non-avian reptiles into the minimalist chicken genome. We have used paired BAC end sequences (BESs) from the American alligator (Alligator mississippiensis), painted turtle (Chrysemys picta) and emu (Dromaius novaehollandiae) to investigate patterns of sequence divergence, gene and retroelement content, and microsynteny between these species and chicken.

Results

From a total of 11,967 curated BESs, we successfully mapped 725, 773 and 2597 sequences in alligator, turtle, and emu, respectively, to sites in the draft chicken genome using a stringent BLAST protocol. Most commonly, sequences mapped to a single site in the chicken genome. Of 1675, 1828 and 2936 paired BESs obtained for alligator, turtle, and emu, respectively, a total of 34 (alligator, 2%), 24 (turtle, 1.3%) and 479 (emu, 16.3%) pairs were found to map with high confidence and in the correct orientation and with BAC-sized intermarker distances to single chicken chromosomes, including 25 such paired hits in emu mapping to the chicken Z chromosome. By determining the insert sizes of a subset of BAC clones from these three species, we also found a significant correlation between the intermarker distance in alligator and turtle and in chicken, with slopes as expected on the basis of the ratio of the genome sizes.

Conclusion

Our results suggest that a large number of small-scale chromosomal rearrangements and deletions in the lineage leading to chicken have drastically reduced the number of detected syntenies observed between the chicken and alligator, turtle, and emu genomes and imply that small deletions occurring widely throughout the genomes of reptilian and avian ancestors led to the ~50% reduction in genome size observed in birds compared to reptiles. We have also mapped and identified likely gene regions in hundreds of new BAC clones from these species.

     

Streamlining universal single‐copy orthologue and ultraconserved element design: A case study in Collembola

Genomic data sets are increasingly central to ecological and evolutionary biology, but far fewer resources are available for invertebrates. Powerful new computational tools and the rapidly decreasing cost of Illumina sequencing are beginning to change this, enabling rapid genome assembly and reference marker extraction. We have developed and tested a practical workflow for developing genomic resources in nonmodel groups with real‐world data on Collembola (springtails), one of the most dominant soil animals on Earth. We designed universal molecular marker sets, single‐copy orthologues (BUSCO s) and ultraconserved elements (UCEs), using three existing and 11 newly generated genomes. Both marker types were tested in silico via marker capture success and phylogenetic performance. The new genomes were assembled with Illumina short reads and 9,585?14,743 protein‐coding genes were predicted with ab initio and protein homology evidence. We identified 1,997 benchmarking universal single‐copy orthologues (BUSCO s) across 14 genomes and created and assessed a custom BUSCO data set for extracting single‐copy genes. We also developed a new UCE probe set containing 46,087 baits targeting 1,885 loci. We successfully captured 1,437?1,865 BUSCO s and 975?1,186 UCEs across 14 genomes. Phylogenomic reconstructions using these markers proved robust, giving new insight on deep‐time collembolan relationships. Our study demonstrates the feasibility of generating thousands of universal markers from highly efficient whole‐genome sequencing, providing a valuable resource for genome‐scale investigations in evolutionary biology and ecology.     

Genome-informed Bradyrhizobium taxonomy: where to from here?

Bradyrhizobium is thought to be the largest and most diverse rhizobial genus, but this is not reflected in the number of described species. Although it was one of the first rhizobial genera recognised, its taxonomy remains complex. Various contemporary studies are showing that genome sequence information may simplify taxonomic decisions. Therefore, the growing availability of genomes for Bradyrhizobium will likely aid in the delineation and characterization of new species. In this study, we addressed two aims: first, we reviewed the availability and quality of available genomic resources for Bradyrhizobium. This was achieved by comparing genome sequences in terms of sequencing technologies used and estimated level of completeness for inclusion in genome-based phylogenetic analyses. Secondly, we utilized these genomes to investigate the taxonomic standing of Bradyrhizobium in light of its diverse lifestyles. Although genome sequences differed in terms of their quality and completeness, our data indicate that the use of these genome sequences is adequate for taxonomic purposes. By using these resources, we inferred a fully resolved, well-supported phylogeny. It separated Bradyrhizobium into seven lineages, three of which corresponded to the so-called supergroups known for the genus. Wide distribution of key lifestyle traits such as nodulation, nitrogen fixation and photosynthesis revealed that these traits have complicated evolutionary histories. We present the first robust Bradyrhizobium species phylogeny based on genome sequence information for investigating the evolution of this important assemblage of bacteria. Furthermore, this study provides the basis for using genome sequence information as a resource to make important taxonomic decisions, particularly at the species and genus levels.     

稻属植物的基因组进化

水稻所在的稻属(Oryza)共有24个左右的物种。由于野生稻含有大量的优良农艺性状基因, 在水稻遗传学研究中日益受到重视。随着国际稻属基因组计划的开展, 越来越多的稻属基因组序列被测定, 稻属成为进行比较、功能和进化基因组学研究的模式系统。近期开展的一系列研究对稻属不同基因组区段以及全基因组序列的比较分析, 揭示了稻属在基因组大小、基因移动、多倍体进化、常染色质到异染色质的转化以及着丝粒区域的进化等方面的分子机制。转座子的活性以及转座子因非均等重组或非法重组而造成的删除, 对稻属基因组的扩增和收缩具有重要作用。DNA双链断裂修复介导的基因移动, 特别是非同源末端连接, 是稻属基因组非共线性基因形成的主要来源。稻属基因组从常染色质到异染色质的转换过程, 伴随着转座子的大量扩增、基因片段的区段性和串联重复以及从基因组其他位置不断捕获异染色质基因。对稻属不同物种间基因拷贝数、特异基因和重要农艺性状基因的进化等研究, 可揭示稻属不同物种间表型和适应性差异的分子基础, 将加速水稻的育种和改良。     

Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish

ABSTRACT: BACKGROUND: The availability of a large number of recently sequenced vertebrate genomes opens new avenues to integrate cytogenetics and genomics in comparative and evolutionary studies. Cytogenetic mapping can offer alternative means to identify conserved synteny shared by distinct genomes and also to define genome regions that are still not fine characterized even after wide-ranging nucleotide sequence efforts. An efficient way to perform comparative cytogenetic mapping is based on BAC clones mapping by fluorescence in situ hybridization. In this report, to address the knowledge gap on the genome evolution in cichlid fishes, BAC clones of an Oreochromis niloticus library covering the linkage groups (LG) 1, 3, 5, and 7 were mapped onto the chromosomes of 9 African cichlid species. The cytogenetic mapping data were also integrated with BAC-end sequences information of O. niloticus and comparatively analyzed against the genome of other fish species and vertebrates. RESULTS: The location of BACs from LG1, 3, 5, and 7 revealed a strong chromosomal conservation among the analyzed cichlid species genomes, which evidenced a synteny of the markers of each LG. Comparative in silico analysis also identified large genomic blocks that were conserved in distantly related fish groups and also in other vertebrates. CONCLUSIONS: Although it has been suggested that fishes contain plastic genomes with high rates of chromosomal rearrangements and probably low rates of synteny conservation, our results evidence that large syntenic chromosome segments have been maintained conserved during evolution, at least for the considered markers. Additionally, our current cytogenetic mapping efforts integrated with genomic approaches conduct to a new perspective to address important questions involving chromosome evolution in fishes.