Genes and genomes: Sequencing 5-methylcytosine residues in genomic DNA

To analyse the biological role of 5-methylation of cytosine residues in DNA requires precise and efficient methods for detecting individual 5-methylcytosines (5-MeCs) in genomic DNA. The methods developed over the past decade rely on either differential enzymatic or chemical cleavage of DNA, or more recently on differential sensitivity to chemical conversion of one base to another. The most commonly used methods for studying the methylation profile of DNA, including the bisulphite base-conversion method, are reviewed.     

Sequencing of Australian wild rice genomes reveals ancestral relationships with domesticated rice

The related A genome species of the Oryza genus are the effective gene pool for rice. Here, we report draft genomes for two Australian wild A genome taxa: O. rufipogon‐like population, referred to as Taxon A, and O. meridionalis‐like population, referred to as Taxon B. These two taxa were sequenced and assembled by integration of short‐ and long‐read next‐generation sequencing (NGS) data to create a genomic platform for a wider rice gene pool. Here, we report that, despite the distinct chloroplast genome, the nuclear genome of the Australian Taxon A has a sequence that is much closer to that of domesticated rice (O. sativa) than to the other Australian wild populations. Analysis of 4643 genes in the A genome clade showed that the Australian annual, O. meridionalis, and related perennial taxa have the most divergent (around 3 million years) genome sequences relative to domesticated rice. A test for admixture showed possible introgression into the Australian Taxon A (diverged around 1.6 million years ago) especially from the wild indica/O. nivara clade in Asia. These results demonstrate that northern Australia may be the centre of diversity of the A genome Oryza and suggest the possibility that this might also be the centre of origin of this group and represent an important resource for rice improvement.     

Sequencing of 21 varicella-zoster virus genomes reveals two novel genotypes and evidence of recombination

Genotyping of 21 varicella-zoster virus (VZV) strains using a scattered single nucleotide polymorphism (SNP) method revealed ambiguous SNPs and two nontypeable isolates. For a further genetic characterization, the genomes of all strains were sequenced using the 454 technology. Almost-complete genome sequences were assembled, and most remaining gaps were closed with Sanger sequencing. Phylogenetic analysis of 42 genomes revealed five established and two novel VZV genotypes, provisionally termed VIII and IX. Genotypes VIII and IX are distinct from the previously reported provisional genotypes VI and VII as judged from the SNP pattern. The alignments showed evidence of ancient recombination events in the phylogeny of clade 4 and recent recombinations within single strains: 3/2005 (clade 1), 11 and 405/2007 (clade 3), 8 and DR (clade 4), CA123 and 413/2000 (clade 5), and strains of the novel genotypes VIII and IX. Bayesian tree inference of the thymidine kinase and the polymerase genes of the VZV clades and other varicelloviruses revealed that VZV radiation began some 110,000 years ago, which correlates with the out-of-Africa dispersal of modern humans. The split of ancestral clades 2/4 and 1/3/5/VIII/IX shows the greatest node height.     

Large-scale genomic 2D visualization reveals extensive CG-AT skew correlation in bird genomes

Background  

Bird genomes have very different compositional structure compared with other warm-blooded animals. The variation in the base skew rules in the vertebrate genomes remains puzzling, but it must relate somehow to large-scale genome evolution. Current research is inclined to relate base skew with mutations and their fixation. Here we wish to explore base skew correlations in bird genomes, to develop methods for displaying and quantifying such correlations at different scales, and to discuss possible explanations for the peculiarities of the bird genomes in skew correlation.     

Novel insights into the diversity of catabolic metabolism from ten haloarchaeal genomes

Background

The extremely halophilic archaea are present worldwide in saline environments and have important biotechnological applications. Ten complete genomes of haloarchaea are now available, providing an opportunity for comparative analysis.

Methodology/Principal Findings

We report here the comparative analysis of five newly sequenced haloarchaeal genomes with five previously published ones. Whole genome trees based on protein sequences provide strong support for deep relationships between the ten organisms. Using a soft clustering approach, we identified 887 protein clusters present in all halophiles. Of these core clusters, 112 are not found in any other archaea and therefore constitute the haloarchaeal signature. Four of the halophiles were isolated from water, and four were isolated from soil or sediment. Although there are few habitat-specific clusters, the soil/sediment halophiles tend to have greater capacity for polysaccharide degradation, siderophore synthesis, and cell wall modification. Halorhabdus utahensis and Haloterrigena turkmenica encode over forty glycosyl hydrolases each, and may be capable of breaking down naturally occurring complex carbohydrates. H. utahensis is specialized for growth on carbohydrates and has few amino acid degradation pathways. It uses the non-oxidative pentose phosphate pathway instead of the oxidative pathway, giving it more flexibility in the metabolism of pentoses.

Conclusions

These new genomes expand our understanding of haloarchaeal catabolic pathways, providing a basis for further experimental analysis, especially with regard to carbohydrate metabolism. Halophilic glycosyl hydrolases for use in biofuel production are more likely to be found in halophiles isolated from soil or sediment.     

A haplotype map of allohexaploid wheat reveals distinct patterns of selection on homoeologous genomes

BackgroundBread wheat is an allopolyploid species with a large, highly repetitive genome. To investigate the impact of selection on variants distributed among homoeologous wheat genomes and to build a foundation for understanding genotype-phenotype relationships, we performed population-scale re-sequencing of a diverse panel of wheat lines.ResultsA sample of 62 diverse lines was re-sequenced using the whole exome capture and genotyping-by-sequencing approaches. We describe the allele frequency, functional significance, and chromosomal distribution of 1.57 million single nucleotide polymorphisms and 161,719 small indels. Our results suggest that duplicated homoeologous genes are under purifying selection. We find contrasting patterns of variation and inter-variant associations among wheat genomes; this, in addition to demographic factors, could be explained by differences in the effect of directional selection on duplicated homoeologs. Only a small fraction of the homoeologous regions harboring selected variants overlapped among the wheat genomes in any given wheat line. These selected regions are enriched for loci associated with agronomic traits detected in genome-wide association studies.ConclusionsEvidence suggests that directional selection in allopolyploids rarely acted on multiple parallel advantageous mutations across homoeologous regions, likely indicating that a fitness benefit could be obtained by a mutation at any one of the homoeologs. Additional advantageous variants in other homoelogs probably either contributed little benefit, or were unavailable in populations subjected to directional selection. We hypothesize that allopolyploidy may have increased the likelihood of beneficial allele recovery by broadening the set of possible selection targets.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0606-4) contains supplementary material, which is available to authorized users.     

Sequencing complete mitochondrial and plastid genomes

Organelle genomics has become an increasingly important research field, with applications in molecular modeling, phylogeny, taxonomy, population genetics and biodiversity. Typically, research projects involve the determination and comparative analysis of complete mitochondrial and plastid genome sequences, either from closely related species or from a taxonomically broad range of organisms. Here, we describe two alternative organelle genome sequencing protocols. The "random genome sequencing" protocol is suited for the large majority of organelle genomes irrespective of their size. It involves DNA fragmentation by shearing (nebulization) and blunt-end cloning of the resulting fragments into pUC or BlueScript-type vectors. This protocol excels in randomness of clone libraries as well as in time and cost-effectiveness. The "long-PCR-based genome sequencing" protocol is specifically adapted for DNAs of low purity and quantity, and is particularly effective for small organelle genomes. Library construction by either protocol can be completed within 1 week.     

Comparative genomic analysis reveals evolutionary characteristics and patterns of microRNA clusters in vertebrates

MicroRNAs (miRNAs) are a class of small non-coding RNAs that can play important regulatory roles in many important biological processes. Although clustering patterns of miRNA clusters have been uncovered in animals, the origin and evolution of miRNA clusters in vertebrates are still poorly understood. Here, we performed comparative genomic analyses to construct 51 sets of orthologous miRNA clusters (SOMCs) across seven test vertebrate species, a collection of miRNA clusters from two or more species that are likely to have evolved from a common ancestral miRNA cluster, and used these to systematically examine the evolutionary characteristics and patterns of miRNA clusters in vertebrates. We found that miRNA clusters are continuously generated, and most of them tend to be conserved and maintained in vertebrate genomes, although some adaptive gains and losses of miRNA cluster have occurred during evolution. Furthermore, miRNA clusters appeared relatively early in the evolutionary history might suffer from more complicated adaptive gain-and-loss than those young miRNA clusters. Detailed analysis showed that genomic duplication events of ancestral miRNAs or miRNA clusters are likely to be major driving force and apparently contribute to origin and evolution of miRNA clusters. Comparison of conserved with lineage-specific miRNA clusters revealed that the contribution of duplication events for the formation of miRNA cluster appears to be more important for conserved miRNA clusters than lineage-specific. Our study provides novel insights for further exploring the origins and evolution of miRNA clusters in vertebrates at a genome scale.     

Sequencing intractable DNA to close microbial genomes

Advancement in high throughput DNA sequencing technologies has supported a rapid proliferation of microbial genome sequencing projects, providing the genetic blueprint for in-depth studies. Oftentimes, difficult to sequence regions in microbial genomes are ruled "intractable" resulting in a growing number of genomes with sequence gaps deposited in databases. A procedure was developed to sequence such problematic regions in the "non-contiguous finished" Desulfovibrio desulfuricans ND132 genome (6 intractable gaps) and the Desulfovibrio africanus genome (1 intractable gap). The polynucleotides surrounding each gap formed GC rich secondary structures making the regions refractory to amplification and sequencing. Strand-displacing DNA polymerases used in concert with a novel ramped PCR extension cycle supported amplification and closure of all gap regions in both genomes. The developed procedures support accurate gene annotation, and provide a step-wise method that reduces the effort required for genome finishing.     

Sequencing crop genomes: approaches and applications

Many challenges face plant scientists, in particular those working on crop production, such as a projected increase in population, decrease in water and arable land, changes in weather patterns and predictability. Advances in genome sequencing and resequencing can and should play a role in our response to meeting these challenges. However, several barriers prevent rapid and effective deployment of these tools to a wide variety of crops. Because of the complexity of crop genomes, de novo sequencing with next-generation sequencing technologies is a process fraught with difficulties that then create roadblocks to the utilization of these genome sequences for crop improvement. Collecting rapid and accurate phenotypes in crop plants is a hindrance to integrating genomics with crop improvement, and advances in informatics are needed to put these tools in the hands of the scientists on the ground.     

Sequencing genomes from single cells by polymerase cloning

Genome sequencing currently requires DNA from pools of numerous nearly identical cells (clones), leaving the genome sequences of many difficult-to-culture microorganisms unattainable. We report a sequencing strategy that eliminates culturing of microorganisms by using real-time isothermal amplification to form polymerase clones (plones) from the DNA of single cells. Two Escherichia coli plones, analyzed by Affymetrix chip hybridization, demonstrate that plonal amplification is specific and the bias is randomly distributed. Whole-genome shotgun sequencing of Prochlorococcus MIT9312 plones showed 62% coverage of the genome from one plone at a sequencing depth of 3.5x, and 66% coverage from a second plone at a depth of 4.7x. Genomic regions not revealed in the initial round of sequencing are recovered by sequencing PCR amplicons derived from plonal DNA. The mutation rate in single-cell amplification is <2 x 10(5), better than that of current genome sequencing standards. Polymerase cloning should provide a critical tool for systematic characterization of genome diversity in the biosphere.     

Compositional patterns in reptilian genomes

Sauropsids form a complex group of vertebrates including squamates (lizards and snakes), turtles, crocodiles, sphenodon and birds (which are often considered as a separate class). Although avian genomes have been relatively well studied, the genomes of the other groups have remained only sparsely characterized. Moreover, the nuclear sequences available in databanks are still very limited. In the present study, we have analysed the compositional patterns, i.e. the GC (molar fraction of guanine and cytosine in DNA) distributions, of 31 reptilian (particularly snake) genomes by analytical ultracentrifugation of DNAs in CsCl gradients. The profiles were characterized by their modal buoyant density rho(o), mean buoyant density < rho>, asymmetry < rho>- rho(o), and heterogeneity H. The modal buoyant density distribution of reptilian DNAs clearly distinguishes two groups. The snakes fall in the same range of modal densities as most mammals, whereas crocodiles, turtles and lizards show higher values (>1.700 g/cm(3)). As far as the more important compositional properties of asymmetry and heterogeneity are concerned, previous studies showed that amphibians and fishes share relatively low values, whereas birds and mammals are characterized by highly heterogeneous and asymmetric patterns (with the exception of Muridae, which have a lower heterogeneity). The present results show that the snake genomes cover a broad range of asymmetry and heterogeneity values, whereas the genomes of crocodiles and turtles cover a narrow range that is intermediate between those of fishes/amphibians and those of mammals/birds.     

Evolutionary patterns in prokaryotic genomes

Prokaryotic genomics is shifting towards comparative approaches to unravel how and why genomes change over time. Both phylogenetic and population genetics approaches are required to dissect the relative roles of selection and drift under these conditions. Lineages evolve adaptively by selection of changes in extant genomes and the way this occurs is being explored from a systemic and evolutionary perspective to understand how mutations relate with gene repertoire changes and how both are contextualized in cellular networks. Through an increased appreciation of genome dynamics in given ecological contexts, a more detailed picture of the genetic basis of prokaryotic evolution is emerging.     

Sequencing of avian influenza virus genomes following random amplification

Increasing surveillance for the avian influenza virus (AIV) has underscored the need for quickly and precisely characterizing isolates of this highly variable target. Random amplification, sequencing, and assembly of total RNA from nonpurified virus overcomes the need for specific primers for DNA microarray or PCR protocols.