Accurate recognition of cis-regulatory motifs with the correct lengths in prokaryotic genomes

We present a new computational method for solving a classical problem, the identification problem of cis-regulatory motifs in a given set of promoter sequences, based on one key new idea. Instead of scoring candidate motifs individually like in all the existing motif-finding programs, our method scores groups of candidate motifs with similar sequences, called motif closures, using a P-value, which has substantially improved the prediction reliability over the existing methods. Our new P-value scoring scheme is sequence length independent, hence allowing direct comparisons among predicted motifs with different lengths on the same footing. We have implemented this method as a Motif Recognition Computer (MREC) program, and have extensively tested MREC on both simulated and biological data from prokaryotic genomes. Our test results indicate that MREC can accurately pick out the actual motif with the correct length as the best scoring candidate for the vast majority of the cases in our test set. We compared our prediction results with two motif-finding programs Cosmo and MEME, and found that MREC outperforms both programs across all the test cases by a large margin. The MREC program is available at http://csbl.bmb.uga.edu/~bingqiang/MREC1/.     

Common physical map of four Brucella bacteriophage genomes

Abstract DNAs isolated from four strains of Brucella bacteriophages were studied by restriction endonuclease mapping and Southern blot analysis. In all strains the genome was composed of a 38 kb (25.1 × 10⁶ dalton) double-stranded circular DNA. The physical map was the same for the four genomes and Southern blot hybridization of restriction endonuclease fragments with the Tbilissi strain DNA as a probe showed complete homology between the four DNAs. Thus, the four phage strains appear to be identical, the specific host range of each originating from minor changes in phage or Brucella receptors or both.     

A high-density consensus map of A and B wheat genomes

A durum wheat consensus linkage map was developed by combining segregation data from six mapping populations. All of the crosses were derived from durum wheat cultivars, except for one accession of T. ssp. dicoccoides. The consensus map was composed of 1,898 loci arranged into 27 linkage groups covering all 14 chromosomes. The length of the integrated map and the average marker distance were 3,058.6 and 1.6?cM, respectively. The order of the loci was generally in agreement with respect to the individual maps and with previously published maps. When the consensus map was aligned to the deletion bin map, 493 markers were assigned to specific bins. Segregation distortion was found across many durum wheat chromosomes, with a higher frequency for the B genome. This high-density consensus map allowed the scanning of the genome for chromosomal rearrangements occurring during the wheat evolution. Translocations and inversions that were already known in literature were confirmed, and new putative rearrangements are proposed. The consensus map herein described provides a more complete coverage of the durum wheat genome compared with previously developed maps. It also represents a step forward in durum wheat genomics and an essential tool for further research and studies on evolution of the wheat genome.     

Reconstruction of a composite comparative map composed of ten legume genomes

The Fabaceae (legume family) is the third largest and the second of agricultural importance among flowering plant groups. In this study, we report the reconstruction of a composite comparative map composed of ten legume genomes, including seven species from the galegoid clade (Medicago truncatula, Medicago sativa, Lens culinaris, Pisum sativum, Lotus japonicus, Cicer arietinum, Vicia faba) and three species from the phaseoloid clade (Vigna radiata, Phaseolus vulgaris, Glycine max). To accomplish this comparison, a total of 209 cross-species gene-derived markers were employed. The comparative analysis resulted in a single extensive genetic/genomic network composed of 93 chromosomes or linkage groups, from which 110 synteny blocks and other evolutionary events (e.g., 13 inversions) were identified. This comparative map also allowed us to deduce several large scale evolutionary events, such as chromosome fusion/fission, with which might explain differences in chromosome numbers among compared species or between the two clades. As a result, useful properties of cross-species genic markers were re-verified as an efficient tool for cross-species translation of genomic information, and similar approaches, combined with a high throughput bioinformatic marker design program, should be effective for applying the knowledge of trait-associated genes to other important crop species for breeding purposes. Here, we provide a basic comparative framework for the ten legume species, and expect to be usefully applied towards the crop improvement in legume breeding.     

Data compression and genomes: a two-dimensional life domain map

We define the complexity of DNA sequences as the information content per nucleotide, calculated by means of some Lempel-Ziv data compression algorithm. It is possible to use the statistics of the complexity values of the functional regions of different complete genomes to distinguish among genomes of different domains of life (Archaea, Bacteria and Eukarya). We shall focus on the distribution function of the complexity of non-coding regions. We show that the three domains may be plotted in separate regions within the two-dimensional space where the axes are the skewness coefficient and the curtosis coefficient of the aforementioned distribution. Preliminary results on 15 genomes are introduced.     

A high-resolution map of synteny disruptions in gibbon and human genomes

Gibbons are part of the same superfamily (Hominoidea) as humans and great apes, but their karyotype has diverged faster from the common hominoid ancestor. At least 24 major chromosome rearrangements are required to convert the presumed ancestral karyotype of gibbons into that of the hominoid ancestor. Up to 28 additional rearrangements distinguish the various living species from the common gibbon ancestor. Using the northern white-cheeked gibbon (2n = 52) (Nomascus leucogenys leucogenys) as a model, we created a high-resolution map of the homologous regions between the gibbon and human. The positions of 100 synteny breakpoints relative to the assembled human genome were determined at a resolution of about 200 kb. Interestingly, 46% of the gibbon–human synteny breakpoints occur in regions that correspond to segmental duplications in the human lineage, indicating a common source of plasticity leading to a different outcome in the two species. Additionally, the full sequences of 11 gibbon BACs spanning evolutionary breakpoints reveal either segmental duplications or interspersed repeats at the exact breakpoint locations. No specific sequence element appears to be common among independent rearrangements. We speculate that the extraordinarily high level of rearrangements seen in gibbons may be due to factors that increase the incidence of chromosome breakage or fixation of the derivative chromosomes in a homozygous state.     

The evolutionary significance of cis-regulatory mutations

For decades, evolutionary biologists have argued that changes in cis-regulatory sequences constitute an important part of the genetic basis for adaptation. Although originally based on first principles, this claim is now empirically well supported: numerous studies have identified cis-regulatory mutations with functionally significant consequences for morphology, physiology and behaviour. The focus has now shifted to considering whether cis-regulatory and coding mutations make qualitatively different contributions to phenotypic evolution. Cases in which parallel mutations have produced parallel trait modifications in particular suggest that some phenotypic changes are more likely to result from cis-regulatory mutations than from coding mutations.     

Mapping unexplored genomes: a genetic linkage map of the Hawaiian cricket Laupala

As with many organisms of evolutionary interest, the Hawaiian cricket Laupala genome is not well characterized genetically. Mapping such an unexplored genome therefore presents challenges not often faced in model genetic organisms and not well covered in the literature. We discuss the evolutionary merits of Laupala as a model for speciation studies involving prezygotic change, our choice of marker system for detecting genetic variation, and the initial genetic expectations pertaining to the construction of any unknown genomic map in general and to the Laupala linkage map construction in particular. We used the technique of amplified fragment length polymorphism (AFLP) to develop a linkage map of Laupala. We utilized both EcoRI/MseI- and EcoRI/PstI-digested genomic DNA to generate AFLP bands and identified 309 markers that segregated among F(2) interspecific hybrid individuals. The map is composed of 231 markers distributed over 11 and 7 species-specific autosomal groups together with a number of putative X chromosome linkage groups. The integration of codominant markers enabled the identification of five homologous linkage groups corresponding to five of the seven autosomal chromosomal pairs found in Laupala.     

On the possibility of obtaining a physical map of genomes by photoelectron imaging.

Photoelectron imaging provides the possibility of a new method of mapping chromosomes. The basic concept is to cause DNA to emit electrons under the action of UV light. The criteria which must be met to map genomes by photoelectron imaging are set forth and discussed. Forming an image of the DNA by accelerating and focusing the electrons is a necessary but not sufficient condition for genome mapping. Equally important is to identify wavelengths of UV light which will cause selective emission from the base pairs, adenine-thymine and guanine-cytosine. The resulting image would then contain a modulation in the image brightness along the DNA duplex. By examining the photoelectron current from uniform films of homopolymers, a wavelength region is identified where marked differences in emission from base pairs is observed. At 160 nm, for example, the relative electron emission from a film of poly(dGdC) is approximately 5 times greater than for an equivalent film of poly(dAdT). Using the experimental data and known sequences, photoelectron gene maps are calculated for the bacteriophage lambda and for a short interspersed repetitive DNA sequence (an Alu repeat) of the human genome. The results suggest that a 5-nm physical map of chromosomes generated by photoelectron imaging would be informative and useful in mapping human and other large genomes.     

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.     

A deer (subfamily Cervinae) genetic linkage map and the evolution of ruminant genomes

Comparative maps between ruminant species and humans are increasingly important tools for the discovery of genes underlying economically important traits. In this article we present a primary linkage map of the deer genome derived from an interspecies hybrid between red deer (Cervus elaphus) and Père David's deer (Elaphurus davidianus). The map is approximately 2500 cM long and contains >600 markers including both evolutionary conserved type I markers and highly polymorphic type II markers (microsatellites). Comparative mapping by annotation and sequence similarity (COMPASS) was demonstrated to be a useful tool for mapping bovine and ovine ESTs in deer. Using marker order as a phylogenetic character and comparative map information from human, mouse, deer, cattle, and sheep, we reconstructed the karyotype of the ancestral Pecoran mammal and identified the chromosome rearrangements that have occurred in the sheep, cattle, and deer lineages. The deer map and interspecies hybrid pedigrees described here are a valuable resource for (1) predicting the location of orthologs to human genes in ruminants, (2) mapping QTL in farmed and wild deer populations, and (3) ruminant phylogenetic studies.     

Bacterial cis-regulatory RNA structures

The review considers mechanism of bacterial gene expression based on formation of alternative RNA structures, such as riboswitches, attenuators, T-boxes, etc. These structures are classified by mechanism of action. Evolution and interaction of regulatory systems are discussed.     

A cis-regulatory logic simulator

Background  

A major goal of computational studies of gene regulation is to accurately predict the expression of genes based on the cis-regulatory content of their promoters. The development of computational methods to decode the interactions among cis-regulatory elements has been slow, in part, because it is difficult to know, without extensive experimental validation, whether a particular method identifies the correct cis-regulatory interactions that underlie a given set of expression data. There is an urgent need for test expression data in which the interactions among cis-regulatory sites that produce the data are known. The ability to rapidly generate such data sets would facilitate the development and comparison of computational methods that predict gene expression patterns from promoter sequence.