单细胞自由生物基因组全序列的测定和比较基因组研究

近几年来五种单细胞生物的基因组计划得以完成。本文介绍了从五种生物的全基因组序列获得的一些成果,包括全基因组鸟枪法测序、基因组分析和新比较基因组等三个方面,并对生物基因组计划的研究方法作一些探讨。     

Phylogeny and sex chromosome evolution of Palaeognathae

Many paleognaths (ratites and tinamous) have a pair of homomorphic ZW sex chromosomes in contrast to the highly differentiated sex chromosomes of most other birds. To understand the evolutionary causes for the different tempos of sex chromosome evolution, we produced female genomes of 12 paleognathous species and reconstructed the phylogeny and the evolutionary history of paleognathous sex chromosomes. We uncovered that Palaeognathae sex chromosomes had undergone stepwise recombination suppression and formed a pattern of “evolutionary strata”. Nine of the 15 studied species' sex chromosomes have maintained homologous recombination in their long pseudoautosomal regions extending more than half of the entire chromosome length. We found that in the older strata, the W chromosome suffered more serious functional gene loss. Their homologous Z-linked regions, compared with other genomic regions, have produced an excess of species-specific autosomal duplicated genes that evolved female-specific expression, in contrast to their broadly expressed progenitors. We speculate such “defeminization” of Z chromosome with underrepresentation of female-biased genes and slow divergence of sex chromosomes of paleognaths might be related to their distinctive mode of sexual selection targeting females rather than males, which evolved in their common ancestors.     

Sequencing chromosome 5D of Aegilops tauschii and comparison with its allopolyploid descendant bread wheat (Triticum aestivum)

Flow cytometric sorting of individual chromosomes and chromosome‐based sequencing reduces the complexity of large, repetitive Triticeae genomes. We flow‐sorted chromosome 5D of Aegilops tauschii, the D genome donor of bread wheat and sequenced it by Roche 454 GS FLX platform to approximately 2.2x coverage. Repetitive sequences represent 81.09% of the survey sequences of this chromosome, and Class I retroelements are the prominent type, with a particular abundance of LTR/Gypsy superfamily. Nonrepetitive sequences were assembled to cover 17.76% of the total chromosome regions. Up to 6188 nonrepetitive gene loci were predicted to be encoded by the 5D chromosome. The numbers and chromosomal distribution patterns of tRNA genes suggest abundance in tRNA^Lys and tRNA^Met species, while the nonrepetitive assembly reveals tRNA^Ala species as the most abundant type. A comparative analysis of the genomic sequences of bread wheat and Aegilops chromosome 5D indicates conservation of gene content. Orthologous unique genes, matching Aegilops 5D sequences, numbered 3730 in barley, 5063 in Brachypodium, 4872 in sorghum and 4209 in rice. In this study, we provide a chromosome‐specific view into the structure and organization of the 5D chromosome of Ae. tauschii, the D genome ancestor of bread wheat. This study contributes to our understanding of the chromosome‐level evolution of the wheat genome and presents a valuable resource in wheat genomics due to the recent hybridization of Ae. tauschii genome with its tetraploid ancestor.     

Construction of a dairy microbial genome catalog opens new perspectives for the metagenomic analysis of dairy fermented products

Background

Microbial communities of traditional cheeses are complex and insufficiently characterized. The origin, safety and functional role in cheese making of these microbial communities are still not well understood. Metagenomic analysis of these communities by high throughput shotgun sequencing is a promising approach to characterize their genomic and functional profiles. Such analyses, however, critically depend on the availability of appropriate reference genome databases against which the sequencing reads can be aligned.

Results

We built a reference genome catalog suitable for short read metagenomic analysis using a low-cost sequencing strategy. We selected 142 bacteria isolated from dairy products belonging to 137 different species and 67 genera, and succeeded to reconstruct the draft genome of 117 of them at a standard or high quality level, including isolates from the genera Kluyvera, Luteococcus and Marinilactibacillus, still missing from public database. To demonstrate the potential of this catalog, we analysed the microbial composition of the surface of two smear cheeses and one blue-veined cheese, and showed that a significant part of the microbiota of these traditional cheeses was composed of microorganisms newly sequenced in our study.

Conclusions

Our study provides data, which combined with publicly available genome references, represents the most expansive catalog to date of cheese-associated bacteria. Using this extended dairy catalog, we revealed the presence in traditional cheese of dominant microorganisms not deliberately inoculated, mainly Gram-negative genera such as Pseudoalteromonas haloplanktis or Psychrobacter immobilis, that may contribute to the characteristics of cheese produced through traditional methods.

Electronic supplementary material

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

Sequencing and assembly of low copy and genic regions of isolated Triticum aestivum chromosome arm 7DS

The genome of bread wheat (Triticum aestivum) is predicted to be greater than 16 Gbp in size and consist predominantly of repetitive elements, making the sequencing and assembly of this genome a major challenge. We have reduced genome sequence complexity by isolating chromosome arm 7DS and applied second‐generation technology and appropriate algorithmic analysis to sequence and assemble low copy and genic regions of this chromosome arm. The assembly represents approximately 40% of the chromosome arm and all known 7DS genes. Comparison of the 7DS assembly with the sequenced genomes of rice (Oryza sativa) and Brachypodium distachyon identified large regions of conservation. The syntenic relationship between wheat, B. distachyon and O. sativa, along with available genetic mapping data, has been used to produce an annotated draft 7DS syntenic build, which is publicly available at http://www.wheatgenome.info . Our results suggest that the sequencing of isolated chromosome arms can provide valuable information of the gene content of wheat and is a step towards whole‐genome sequencing and variation discovery in this important crop.     

Divergent microsatellite evolution in the human and chimpanzee lineages

Comparison of the complete human genome sequence to one of its closest relatives, the chimpanzee genome, provides a unique opportunity for exploring recent evolutionary events affecting the microsatellites in these species. A simple assumption on microsatellite distribution is that the total length of perfect repeats is constant compared to that of imperfect ones regardless of the repeat sequence. In this paper, we show that this is valid for most of the chimpanzee genome but not for a number of human chromosomes. Our results suggest accelerated evolution of microsatellites in the human genome relative to the chimpanzee lineage.     

The structure of the Aegilops tauschii genepool and the evolution of hexaploid wheat

 Polymorphism in the lengths of restriction fragments at 53 single-copy loci, the rRNA locus Nor3, and the high-molecular-weight glutenin locus Glu1 was investigated in the D genome of hexaploid Triticum aestivum and that of Aegilops tauschii, the source of the T. aestivum D genome. The distribution of genetic variation in Ae. tauschii suggests gene flow between Ae. tauschii ssp. strangulata and ssp. tauschii in Iran but less in Transcaucasia. The “strangulata” genepool is wider than it appears on the basis of morphology and includes ssp. strangulata in Transcaucasia and southeastern (SE) Caspian Iran and ssp. tauschii in north-central Iran and southwestern (SW) Caspian Iran. In the latter region, Ae. tauschii morphological varieties ‘meyeri’ and ‘typica’ are equidistant to ssp. strangulata in Transcaucasia, and both belong to the “strangulata” genepool. A model of the evolution of Ae. tauschii is presented. On the geographic region basis, the D genomes of all investigated forms of T. aestivum are most closely related to the “strangulata” genepool in Transcaucasia, Armenia in particular, and SW Caspian Iran. It is suggested that the principal area of the origin of T. aestivum is Armenia, but the SW coastal area of the Caspian Sea and a corridor between the two areas may have played a role as well. Little genetic differentiation was found among the D genomes of all investigated free-threshing and hulled forms of T. aestivum, and all appear to share a single D-genome genepool, in spite of the fact that several Ae. tauschii parents were involved in the evolution of T. aestivum. Received: 17 November 1997 / Accepted: 17 March 1998     

植物多倍体基因组的形成与进化

多倍化是植物进化变异的自然现象,也是促进植物发生进化改变的重要力量。在被子植物中,约 70％的种类在进化史中曾发生过一次或多次多倍化的过程。目前的研究结果表明,自然界绝大多数多倍体是通过未减数配子的融合而形成的,并且很多多倍体种是通过多次独立的多倍化过程而重复发生的。由多倍化所导致的重复基因在多倍体基因组中可能有三种不同的命运,即：保持原有的功能、基因沉默或分化并执行新的功能。多倍化以后,重复基因组的进化动态则主要表现在染色体重排和“染色体二倍化”、不同基因组之间的相互渗透、以及核-质之间的相互作用等方面。     

Characterizing the composition and evolution of homoeologous genomes in hexaploid wheat through BAC-end sequencing on chromosome 3B

Bread wheat (Triticum aestivum) is one of the most important crops worldwide. However, because of its large, hexaploid, highly repetitive genome it is a challenge to develop efficient means for molecular analysis and genetic improvement in wheat. To better understand the composition and molecular evolution of the hexaploid wheat homoeologous genomes and to evaluate the potential of BAC-end sequences (BES) for marker development, we have followed a chromosome-specific strategy and generated 11 Mb of random BES from chromosome 3B, the largest chromosome of bread wheat. The sequence consisted of about 86% of repetitive elements, 1.2% of coding regions, and 13% remained unknown. With 1.2% of the sequence length corresponding to coding sequences, 6000 genes were estimated for chromosome 3B. New repetitive sequences were identified, including a Triticineae-specific tandem repeat (Fat) that represents 0.6% of the B-genome and has been differentially amplified in the homoeologous genomes before polyploidization. About 10% of the BES contained junctions between nested transposable elements that were used to develop chromosome-specific markers for physical and genetic mapping. Finally, sequence comparison with 2.9 Mb of random sequences from the D-genome of Aegilops tauschii suggested that the larger size of the B-genome is due to a higher content in repetitive elements. It also indicated which families of transposable elements are mostly responsible for differential expansion of the homoeologous wheat genomes during evolution. Our data demonstrate that BAC-end sequencing from flow-sorted chromosomes is a powerful tool for analysing the structure and evolution of polyploid and highly repetitive genomes.     

Comparative analysis of methylthioalkylmalate synthase (MAM) gene family and flanking DNA sequences in Brassica oleracea and Arabidopsis thaliana

Gene BoGSL-PRO is associated with presence of 3-carbon side-chain glucosinolates (GSL). This gene is a member of the methylthioalkylmalate synthase (MAM) gene family. A BAC clone of Brassica oleracea, B21F5, containing this gene, was sequenced, annotated and compared to its corresponding region in Arabidopsis thaliana. Twelve protein-coding genes and 10 transposable elements were found in this clone. The corresponding region in A. thaliana chromosome I has 14 genes and no transposable elements. Analysis of MAM gene family in both species, which also include genes controlling 4-carbon side-chain GSL, separated the genes in two groups based on exon numbers and function. Phylogenetic analysis of the amino acid sequences encoded by these genes suggest that these two groups were produced by a duplication that must have occurred before the divergence of the Rosid and Asterid lineages of angiosperms. Comparison with putative orthologs from several prokaryotes further suggest that the members of the gene family with 10 exons, which encode proteins involved in 4-carbon side-chain GSL biosynthesis, were derived via truncation of the 3′ end from ancestral genes more similar in length to those with 12 exons, which encode proteins involved in 3-carbon side-chain GSL biosynthesis. Lower gene density in B. oleracea compared to A. thaliana is due in part to presence of transposable elements (TE) mostly in inter-genic regions.     

Comparative RFLP mapping of Hordeum vulgare and Triticum tauschii

Hordeum vulgare (barley) and Triticum tauschii are related, but sexually incompatible, species. This study was conducted to determine the extent of homology between the genomes of barley and T. tauschii using a common set of restriction fragment length polymorphism (RFLP) markers. Results showed that >95% of low-copy sequences are shared, but 42% of the conserved sequences showed copy-number differences. Sixty-three loci were mapped in T. tauschii using RFLP markers previously mapped in barley. A comparison of RFLP marker order showed that, in general, barley and T. tauschii have conserved linkage groups, with markers in the same linear orders. However, six of the seven linkage groups of T. tauschii contained markers which mapped to unrelated (i.e., non-homoeologous) barley chromosomes. Additionally, four of the T. tauschii linkage groups contained markers that were switched in order with respect to barley. All the chromosome segments differing between T. tauschii and barley contained markers that were detected by multi-copy probes. The results suggest that the observed differences between the T. tauschii and barley genomes were brought about by duplications or deletions of segments in one or both species. The implications of these findings for genetic mapping, breeding, and plant genome evolution are discussed.Published with the approval of the Director of the Colorado State University/Agricultural Experiment Station     

Sequencing,annotation and comparative analysis of nine BACs of the giant panda(Ailuropoda melanoleuca)

A 10-fold BAC library for the giant panda was constructed and nine BACs were selected to generate finish sequences.These BACs could be used as a validation resource for the de novo assembly accuracy of the whole genome shotgun sequencing reads of the giant panda newly generated by Illumina GA sequencing technology.Complete Sanger sequencing,assembly,annotation and comparative analysis were carried out on the selected BACs of a joint length 878 kb.Homologue search and de novo prediction methods were used to ...     

Variation and evolution of biomolecular systems: searching for functional relevance

The availability of genome sequences and functional genomics data from multiple species enables us to compare the composition of biomolecular systems like biochemical pathways and protein complexes between species. Here, we review small- and large-scale, "genomics-based" approaches to biomolecular systems variation. In general, caution is required when comparing the results of bioinformatics analyses of genomes or of functional genomics data between species. Limitations to the sensitivity of sequence analysis tools and the noisy nature of genomics data tend to lead to systematic overestimates of the amount of variation. Nevertheless, the results from detailed manual analyses, and of large-scale analyses that filter out systematic biases, point to a large amount of variation in the composition of biomolecular systems. Such observations challenge our understanding of the function of the systems and their individual components and can potentially facilitate the identification and functional characterization of sub-systems within a system. Mapping the inter-species variation of complex biomolecular systems on a phylogenetic species tree allows one to reconstruct their evolution.     

Genome sequence and transcriptome of Sorbus pohuashanensis provide insights into population evolution and leaf sunburn response

Sorbus pohuashanensis (Hance) Hedl. is a potential horticulture and medicinal plant, but its genomic and genetic backgrounds remain unknown. Here, we sequence and assemble the S. pohuashanensis reference genome using PacBio long reads. Based on the new reference genome, we resequence a core collection of 22 Sorbus spp. samples, which are divided into 2 groups (G1 and G2) based on phylogenetic and PCA analyses. These phylogenetic clusters are highly consistent with their classification based on leaf shape. Natural hybridization between the G1 and G2 groups is evidenced by a sample (R21) with a highly heterozygous genotype. Nucleotide diversity (π) analysis shows that G1 has a higher diversity than G2 and that G2 originated from G1. During the evolution process, the gene families involved in photosynthesis pathways expanded and the gene families involved in energy consumption contracted. RNA-seq data suggests that flavonoid biosynthesis and heat-shock protein (HSP)-heat-shock factor (HSF) pathways play important roles in protection against sunburn. This study provides new insights into the evolution of Sorbus spp. genomes. In addition, the genomic resources, and the identified genetic variations, especially those related to stress resistance, will help future efforts to produce and breed Sorbus spp.     

Systematic Cross-biospecimen Evaluation of DNA Extraction Kits for Long- and Short-read Multi-metagenomic Sequencing Studies

High-quality DNA extraction is a crucial step in metagenomic studies. Bias by different isolation kits impairs the comparison across datasets. A trending topic is, however, the analysis of multiple metagenomes from the same patients to draw a holistic picture of microbiota associated with diseases. We thus collected bile, stool, saliva, plaque, sputum, and conjunctival swab samples and performed DNA extraction with three commercial kits. For each combination of the specimen type and DNA extraction kit, 20-gigabase (Gb) metagenomic data were generated using short-read sequencing. While profiles of the specimen types showed close proximity to each other, we observed notable differences in the alpha diversity and composition of the microbiota depending on the DNA extraction kits. No kit outperformed all selected kits on every specimen. We reached consistently good results using the Qiagen QiAamp DNA Microbiome Kit. Depending on the specimen, our data indicate that over 10 Gb of sequencing data are required to achieve sufficient resolution, but DNA-based identification is superior to identification by mass spectrometry. Finally, long-read nanopore sequencing confirmed the results (correlation coefficient > 0.98). Our results thus suggest using a strategy with only one kit for studies aiming for a direct comparison of multiple microbiotas from the same patients.     

Directed sequencing and annotation of three Dicentrarchus labrax L. chromosomes by applying Sanger- and pyrosequencing technologies on pooled DNA of comparatively mapped BAC clones

Dicentrarchus labrax is one of the major marine aquaculture species in the European Union. In this study, we have developed a directed-sequencing strategy to sequence three sea bass chromosomes and compared results with other teleosts.Three BAC DNA pools were created from sea bass BAC clones that mapped to stickleback chromosomes/groups V, XVII and XXI. The pools were sequenced to 17-39x coverage by pyrosequencing. Data assembly was supported by Sanger reads and mate pair data and resulted in superscaffolds of 13.2 Mb, 17.5 Mb and 13.7 Mb respectively. Annotation features of the superscaffolds include 1477 genes. We analyzed size change of exon, intron and intergenic sequence between teleost species and deduced a simple model for the evolution of genome composition in teleost lineage.Combination of second generation sequencing technologies, Sanger sequencing and genome partitioning strategies allows “high-quality draft assemblies” of chromosome-sized superscaffolds, which are crucial for the prediction and annotation of complete genes.     

一株新的产吡咯喹啉醌甲基营养菌全基因组测序和比较基因组学分析

【背景】由于甲基营养菌被发现的时间较短,而且可以生产吡咯喹啉醌(pyrroloquinoline quinone,PQQ)的甲基杆菌属细菌只有少数菌株的全基因组序列被公布,增加了该类细菌基因组学和生物代谢途径研究的难度。【目的】将本实验室筛选的PQQ生产菌经多种诱变方式处理,用于提高PQQ的发酵产量。对高产突变菌株进行全基因组解析,以探究甲基杆菌PQQ合成的分子机制,为后续分子育种提供序列背景信息。【方法】将野生型PQQ生产菌株进行紫外诱变、亚硝基胍诱变、甲基磺酸乙酯诱变、硫酸二乙酯诱变和紫外-氯化锂复合诱变。将突变菌株利用PromethION三代测序平台和MGISEQ-2000二代测序平台测序,然后进行组装和功能注释。组装得到的全基因组序列与模式菌株扭脱甲基杆菌AM1 (Methylobacterium extorquens AM1)进行比较基因组学分析。【结果】经11轮诱变获得一株突变菌株NI91,其PQQ产量为19.49mg/L,相较原始菌株提高44.91%。突变菌株NI91的基因组由一个5 409 262 bp的染色体组成,共编码4 957个蛋白,与模式菌株M. extorqu...     

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.     

Genomic characterization of Helicobacter hepaticus: ordered cosmid library and comparative sequence analysis

Helicobacter hepaticus is an important pathogen in laboratory mice and induces the development of liver tumors and gastrointestinal disease in susceptible strains of mice. In this study, a miniset of 36 cosmid clones from a genomic library of H. hepaticus was ordered and grouped into four large contigs representing approximately 1 Mb of the H. hepaticus genome using PCR, DNA sequencing, Southern and dot-blot hybridization and pulsed-field gel electrophoresis. From the 200-300 terminal nucleotide sequences of 38 cosmid clones, 56 coding regions were predicted, of which 51 were found to have orthologs in the public databases and five appeared to be unique to H. hepaticus. Of these 51 genes, 36 have orthologs in Helicobacter pylori and 25 display the highest sequence similarity to H. pylori. However, chromosomal positions of these genes are not conserved between these two helicobacters. In addition, 10 H. hepaticus genes had the highest sequence similarity to orthologs in Campylobacter jejuni. The GC content in a randomly selected 21-kb H. hepaticus genomic sequence was 35.8%, which approximates the average between H. pylori (39%) and C. jejuni (30.6%). These results demonstrate that: (1) H. hepaticus is more closely related to H. pylori than C. jejuni; (2) significant genomic alterations exist between H. hepaticus and H. pylori, including gene organization, protein sequences and GC content, probably in part due to specific adaptation to distinct ecological niches.