Strain-specific genomic regions of Ruminococcus flavefaciens FD-1 as revealed by combinatorial random-phase genome sequencing and suppressive subtractive hybridization

Two closely related strains of the Gram-positive, cellulolytic ruminal bacterium Ruminococcus flavefaciens were compared at the genomic level by suppressive subtractive hybridization. The two strains investigated in this study differ by 1.94% in their respective 16S rDNA genes. Three hundred and eighty-four PCR-amplified products were cloned and then screened for their strain identity by dot blot hybridization. Based on redundancy percentages of the clones sequenced, 9.5% of the genome of the R. flavefaciens FD-1 strain is not present in the JM1 strain. The majority of identities of individual cloned subtracted products (642 bp average length) bore no relation to deposited sequences in GenBank (42% of the subtracted library), whereas of those with putative assigned functions 7% are loosely associated with fibre-degradation, 6% with insertion elements, transposons and phage-like ORFs, 5% with cell membrane associated proteins and 3% with signal transduction. Subtracted sequences were then supplemented with the draft (2 x coverage) genome sequence of R. flavefaciens FD-1 to indicate potential regions of rearrangement within the genome, including a novel insertion sequence.     

Structural dynamics of cereal mitochondrial genomes as revealed by complete nucleotide sequencing of the wheat mitochondrial genome

The application of a new gene-based strategy for sequencing the wheat mitochondrial genome shows its structure to be a 452528 bp circular molecule, and provides nucleotide-level evidence of intra-molecular recombination. Single, reciprocal and double recombinant products, and the nucleotide sequences of the repeats that mediate their formation have been identified. The genome has 55 genes with exons, including 35 protein-coding, 3 rRNA and 17 tRNA genes. Nucleotide sequences of seven wheat genes have been determined here for the first time. Nine genes have an exon–intron structure. Gene amplification responsible for the production of multicopy mitochondrial genes, in general, is species-specific, suggesting the recent origin of these genes. About 16, 17, 15, 3.0 and 0.2% of wheat mitochondrial DNA (mtDNA) may be of genic (including introns), open reading frame, repetitive sequence, chloroplast and retro-element origin, respectively. The gene order of the wheat mitochondrial gene map shows little synteny to the rice and maize maps, indicative that thorough gene shuffling occurred during speciation. Almost all unique mtDNA sequences of wheat, as compared with rice and maize mtDNAs, are redundant DNA. Features of the gene-based strategy are discussed, and a mechanistic model of mitochondrial gene amplification is proposed.     

Development of genomic resources in support of sequencing, assembly, and annotation of the catfish genome

Major progress has been made in catfish genomics including construction of high-density genetic linkage maps, BAC-based physical maps, and integration of genetic linkage and physical maps. Large numbers of ESTs have been generated from both channel catfish and blue catfish. Microarray platforms have been developed for the analysis of genome expression. Genome repeat structures are studied, laying grounds for whole genome sequencing. USDA recently approved funding of the whole genome sequencing project of catfish using the next generation sequencing technologies. Generation of the whole genome sequence is a historical landmark of catfish research as it opens the real first step of the long march toward genetic enhancement. The research community needs to be focused on aquaculture performance and production traits, take advantage of the unprecedented genome information and technology, and make real progress toward genetic improvements of aquaculture brood stocks.     

Preferred analysis methods for single genomic regions in RNA sequencing revealed by processing the shape of coverage

The informational content of RNA sequencing is currently far from being completely explored. Most of the analyses focus on processing tables of counts or finding isoform deconvolution via exon junctions. This article presents a comparison of several techniques that can be used to estimate differential expression of exons or small genomic regions of expression, based on their coverage function shapes. The problem is defined as finding the differentially expressed exons between two samples using local expression profile normalization and statistical measures to spot the differences between two profile shapes. Initial experiments have been done using synthetic data, and real data modified with synthetically created differential patterns. Then, 160 pipelines (5 types of generator × 4 normalizations × 8 difference measures) are compared. As a result, the best analysis pipelines are selected based on linearity of the differential expression estimation and the area under the ROC curve. These platform-independent techniques have been implemented in the Bioconductor package rnaSeqMap. They point out the exons with differential expression or internal splicing, even if the counts of reads may not show this. The areas of application include significant difference searches, splicing identification algorithms and finding suitable regions for QPCR primers.     

The evolution of the avian genome as revealed by comparative molecular cytogenetics

Birds are characterised by feathers, flight, a small genome and a very distinctive karyotype. Despite the large numbers of chromosomes, the diploid count of 2n approximately 80 has remained remarkably constant with 63% of birds where 2n = 74-86, 24% with 2n = 66-74 and extremes of 2n = 40 and 2n = 142. Of these, the most studied is the chicken (2n = 78), and molecular cytogenetic probes generated from this species have been used to further understand the evolution of the avian genome. The ancestral karyotype is, it appears, very similar to that of the chicken, with chicken chromosomes 1, 2, 3, 4q, 5, 6, 7, 8, 9, 4p and Z representing the ancestral avian chromosomes 1-10 + Z; chromosome 4 being the most ancient. Avian evolution occurred primarily in three stages: the divergence of the group represented by extant ratites (emu, ostrich etc.) from the rest; divergence of the Galloanserae (chicken, turkey, duck, goose etc.)--the most studied group; and divergence of the 'land' and 'water' higher birds. Other than sex chromosome differentiation in the first divergence there are no specific changes associated with any of these evolutionary milestones although certain families and orders have undergone multiple fusions (and some fissions), which has reduced their chromosome number; the Falconiformes are the best described. Most changes, overall, seem to involve chromosomes 1, 2, 4, 10 and Z where the Z changes are intrachromosomal; there are also some recurring (convergent) events. Of these, the most puzzling involves chromosomes 4 and 10, which appear to have undergone multiple fissions and/or fusions throughout evolution - three possible hypotheses are presented to explain the findings. We conclude by speculating as to the reasons for the strange behaviour of these chromosomes as well as the role of telomeres and nuclear organisation in avian evolution.     

The origin and evolution of sex chromosomes,revealed by sequencing of the Silene latifolia female genome

1. Download : Download high-res image (180KB)
2. Download : Download full-size image

     

Flexible genomic islands as drivers of genome evolution

1. Download : Download high-res image (244KB)
2. Download : Download full-size image

     

Interactions among genomic structure, function, and evolution revealed by comprehensive analysis of the Arabidopsis thaliana genome

The genome in a higher organism consists of a number of types of nucleotide sequence-specialized components, with each having tens of thousands of members or elements. It is crucial for our understanding of how a genome as an entity is organized, functions, and evolves to determine how these components are organized in the genome and how they relate with each other; however, no such knowledge is available. Here, we report a comprehensive analysis of the organization and interaction of all 40 components constituting the genome of the plant model species, Arabidopsis thaliana, at the whole-genome and chromosome levels. The 40 components include (i) 6 genome structural components consisting of GC%, genes, retrotransposons, DNA transposons, simple repeats, and low complex repeats; (ii) 3 evolutionarily critical features consisting of recombination rate, nucleotide substitutions, and nucleotide insertions/deletions; and (iii) 31 categories of genes with different functions and numbers of functions. We show that the distributions of 39 of the 40 components of the genome (excepting GC%) deviate significantly from the random distribution model and different types of the genome components are significantly correlated. These results remained to be true even when the genomic regions, such as centromeric regions, where transposable and repeat elements are abundant were excluded from the analyses. These findings suggest that DNA molecules contained in the Arabidopsis genome are each organized and structured from their constituting components in an unambiguous manner and that different types of the components that constitute or characterize the genome interact. The analysis also showed that each chromosome consists of a similar set of the components at similar densities, suggesting that the unique organization and interaction pattern of the components in each chromosome may represent, at least in part, the identity of a chromosome or a genome at the genome level, thus partly accounting for the phenotypic variation among different species. The data also provide comprehensive and new insights into many phenomena significant in genome biology, with which we particularly discuss the variation of genetic recombination. The variation of genetic recombination rate along a chromosomal arm is shaped, not only by the distribution of simple repeats, retrotransposons, DNA transposons, and nucleotide substitutions, but also by the functions of genes contained, especially those with multiple functions, suggesting that variation of genetic recombination along a chromosomal arm is the result of interactions among the components constituting local genome structure, function, and evolution.     

Whole genome sequencing and genome annotation of the wild edible mushroom,Russula griseocarnosa

Russula griseocarnosa is a species of edible ectomycorrhizal fungi with medicinal properties that grows in southern China. Total DNA was isolated from a fresh fruiting body of R. griseocarnosa and subjected to sequencing using Illumina Hiseq with the PacBio RS sequencing platform. Here, we present the 64.81 Mb draft genome map of R. griseocarnosa based on 471 scaffolds and 16,128 coding protein genes. The gene annotation of protein coding genes was used to obtain corresponding annotations by blastp. Phylogenetic analysis revealed a close evolutionary relationship of R. griseocarnosa to Heterobasidion irregulare and Stereum hirsutum in the core Russulales clade. The R. griseocarnosa genome encodes a repertoire of enzymes engaged in carbohydrate and polysaccharide metabolism, along with cytochrome P450s and secondary metabolite biosynthesis. The genome content of R. griseocarnosa provides insights into the genetic basis of its reported medicinal properties and serves as a reference for comparative genomics of fungi.     

Rapid genome evolution revealed by comparative sequence analysis of orthologous regions from four triticeae genomes

Bread wheat (Triticum aestivum) is an allohexaploid species, consisting of three subgenomes (A, B, and D). To study the molecular evolution of these closely related genomes, we compared the sequence of a 307-kb physical contig covering the high molecular weight (HMW)-glutenin locus from the A genome of durum wheat (Triticum turgidum, AABB) with the orthologous regions from the B genome of the same wheat and the D genome of the diploid wheat Aegilops tauschii (Anderson et al., 2003; Kong et al., 2004). Although gene colinearity appears to be retained, four out of six genes including the two paralogous HMW-glutenin genes are disrupted in the orthologous region of the A genome. Mechanisms involved in gene disruption in the A genome include retroelement insertions, sequence deletions, and mutations causing in-frame stop codons in the coding sequences. Comparative sequence analysis also revealed that sequences in the colinear intergenic regions of these different genomes were generally not conserved. The rapid genome evolution in these regions is attributable mainly to the large number of retrotransposon insertions that occurred after the divergence of the three wheat genomes. Our comparative studies indicate that the B genome diverged prior to the separation of the A and D genomes. Furthermore, sequence comparison of two distinct types of allelic variations at the HMW-glutenin loci in the A genomes of different hexaploid wheat cultivars with the A genome locus of durum wheat indicates that hexaploid wheat may have more than one tetraploid ancestor.     

Next-generation sequencing and genome evolution in allopolyploids

?Premise of the study: Hybridization and polyploidization (allopolyploidy) are ubiquitous in the evolution of plants, but tracing the origins and subsequent evolution of the constituent genomes of allopolyploids has been challenging. Genome doubling greatly complicates genetic analyses, and this has long hindered investigation in that most allopolyploid species are "nonmodel" organisms. However, recent advances in sequencing and genomics technologies now provide unprecedented opportunities to analyze numerous genetic markers in multiple individuals in any organism. ?Methods: Here we review the application of next-generation sequencing technologies to the study of three aspects of allopolyploid genome evolution: duplicated gene loss and expression in two recently formed Tragopogon allopolyploids, intergenomic interactions and chromosomal evolution in Tragopogon miscellus, and repetitive DNA evolution in Nicotiana allopolyploids. ?Key results: For the first time, we can explore on a genomic scale the evolutionary processes that are ongoing in natural allopolyploids and not be restricted to well-studied crops and genetic models. ?Conclusions: These approaches can be easily and inexpensively applied to many other plant species-making any evolutionarily provocative system a new "model" system.     

Efficient identification of genomic insertions and flanking regions through whole-genome sequencing in three transgenic soybean events

Transgenic Research - Genomic insertions and flanking regions of transgenes in host genomes constitute a critical component of precise molecular characterization and event-specific detection, which...     

Presence of three rice sucrose synthase genes as revealed by cloning and sequencing of cDNA

By sequencing cDNA clones, we have concluded that three distinct sucrose genes are expressed in rice (Oryza sativa cv. Tainong 67). When the amino acid sequences deduced from these cDNAs as well as those of known sucrose synthase are compared, the highest divergence is found in the C-termini. The most suitable DNA sequences for use as specific for the mRNA derived from these genes have been suggested.     

Rapid evolution and complex structural organization in genomic regions harboring multiple prolamin genes in the polyploid wheat genome

Genes encoding wheat prolamins belong to complicated multi-gene families in the wheat genome. To understand the structural complexity of storage protein loci, we sequenced and analyzed orthologous regions containing both gliadin and LMW-glutenin genes from the A and B genomes of a tetraploid wheat species, Triticum turgidum ssp. durum. Despite their physical proximity to one another, the gliadin genes and LMW-glutenin genes are organized quite differently. The gliadin genes are found to be more clustered than the LMW-glutenin genes which are separated from each other by much larger distances. The separation of the LMW-glutenin genes is the result of both the insertion of large blocks of repetitive DNA owing to the rapid amplification of retrotransposons and the presence of genetic loci interspersed between them. Sequence comparisons of the orthologous regions reveal that gene movement could be one of the major factors contributing to the violation of microcolinearity between the homoeologous A and B genomes in wheat. The rapid sequence rearrangements and differential insertion of repetitive DNA has caused the gene islands to be not conserved in compared regions. In addition, we demonstrated that the i-type LMW-glutenin originated from a deletion of 33-bps in the 5′ coding region of the m-type gene. Our results show that multiple rounds of segmental duplication of prolamin genes have driven the amplification of the ω-gliadin genes in the region; such segmental duplication could greatly increase the repetitive DNA content in the genome depending on the amount of repetitive DNA present in the original duplicate region. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Shuangcheng Gao and Yong Qiang Gu contributed equally to the work.     

Commonly conserved genetic fragments revealed by genome profiling can serve as tracers of evolution

We developed a method to produce, identify and analyze DNA fragments for the purpose of taxonomic classification. Genome profiling (GP) is a strategy that identifies genomic DNA fragments common to closely related species without prior knowledge of the DNA sequence. Random PCR, one of the key technologies of GP, is used to produce fragments and may be used even when there are mutations at the priming site. These fragments can then be distinguished based on the information of mobility and melting pattern when subjected to temperature gradient gel electrophoresis (TGGE). Corresponding fragments among several species, designated as commonly conserved genetic fragments (CCGFs), likely have the same genetic origin or correspond to the same gene. The criteria for identification of CCGFs has been defined and presented here. To assess this prediction, some of the fragments were sequenced and were confirmed to be CCGFs. We show that genome profiles bearing evolutionarily conserved CCGFs can be used to classify organisms and trace evolutionary pathways, among other profound applications.     

Genome evolution in diploid and tetraploid Coffea species as revealed by comparative analysis of orthologous genome segments

Sequence comparison of orthologous regions enables estimation of the divergence between genomes, analysis of their evolution and detection of particular features of the genomes, such as sequence rearrangements and transposable elements. Despite the economic importance of Coffea species, little genomic information is currently available. Coffea is a relatively young genus that includes more than one hundred diploid species and a single tetraploid species. Three Coffea orthologous regions of 470-900 kb were analyzed and compared: both subgenomes of allotetraploid Coffea arabica (contributed by the diploid species Coffea eugenioides and Coffea canephora) and the genome of diploid C. canephora. Sequence divergence was calculated on global alignments or on coding and non-coding sequences separately. A search for transposable elements detected 43 retrotransposons and 198 transposons in the sequences analyzed. Comparative insertion analysis made it possible to locate 165 TE insertions in the phylogenetic tree of the three genomes/subgenomes. In the tetraploid C. arabica, a homoeologous non-reciprocal transposition (HNRT) was detected and characterized: a 50 kb region of the C. eugenioides derived subgenome replaced the C. canephora derived counterpart. Comparative sequence analysis on three Coffea genomes/subgenomes revealed almost perfect gene synteny, low sequence divergence and a high number of shared transposable elements. Compared to the results of similar analysis in other genera (Aegilops/Triticum and Oryza), Coffea genomes/subgenomes appeared to be dramatically less diverged, which is consistent with the relatively recent radiation of the Coffea genus. Based on nucleotide substitution frequency, the HNRT was dated at 10,000-50,000 years BP, which is also the most recent estimation of the origin of C. arabica.