Mapping Genetic Diversity of Cherimoya (Annona cherimola Mill.): Application of Spatial Analysis for Conservation and Use of Plant Genetic Resources.

There is a growing call for inventories that evaluate geographic patterns in diversity of plant genetic resources maintained on farm and in species' natural populations in order to enhance their use and conservation. Such evaluations are relevant for useful tropical and subtropical tree species, as many of these species are still undomesticated, or in incipient stages of domestication and local populations can offer yet-unknown traits of high value to further domestication. For many outcrossing species, such as most trees, inbreeding depression can be an issue, and genetic diversity is important to sustain local production. Diversity is also crucial for species to adapt to environmental changes. This paper explores the possibilities of incorporating molecular marker data into Geographic Information Systems (GIS) to allow visualization and better understanding of spatial patterns of genetic diversity as a key input to optimize conservation and use of plant genetic resources, based on a case study of cherimoya (Annona cherimola Mill.), a Neotropical fruit tree species. We present spatial analyses to (1) improve the understanding of spatial distribution of genetic diversity of cherimoya natural stands and cultivated trees in Ecuador, Bolivia and Peru based on microsatellite molecular markers (SSRs); and (2) formulate optimal conservation strategies by revealing priority areas for in situ conservation, and identifying existing diversity gaps in ex situ collections. We found high levels of allelic richness, locally common alleles and expected heterozygosity in cherimoya's putative centre of origin, southern Ecuador and northern Peru, whereas levels of diversity in southern Peru and especially in Bolivia were significantly lower. The application of GIS on a large microsatellite dataset allows a more detailed prioritization of areas for in situ conservation and targeted collection across the Andean distribution range of cherimoya than previous studies could do, i.e. at province and department level in Ecuador and Peru, respectively.

Investigation of the bottleneck leading to the domestication of maize

Maize (Zea mays ssp. mays) is genetically diverse, yet it is also morphologically distinct from its wild relatives. These two observations are somewhat contradictory: the first observation is consistent with a large historical population size for maize, but the latter observation is consistent with strong, diversity-limiting selection during maize domestication. In this study, we sampled sequence diversity, coupled with simulations of the coalescent process, to study the dynamics of a population bottleneck during the domestication of maize. To do this, we determined the DNA sequence of a 1,400-bp region of the Adh1 locus from 19 individuals representing maize, its presumed progenitor (Z. mays ssp. parviglumis), and a more distant relative (Zea luxurians). The sequence data were used to guide coalescent simulations of population bottlenecks associated with domestication. Our study confirms high genetic diversity in maize—maize contains 75% of the variation found in its progenitor and is more diverse than its wild relative, Z. luxurians—but it also suggests that sequence diversity in maize can be explained by a bottleneck of short duration and very small size. For example, the breadth of genetic diversity in maize is consistent with a founding population of only 20 individuals when the domestication event is 10 generations in length.

Pro-phenol oxidase activating proteinase from an insect, Manduca sexta: A bacteria-inducible protein similar to Drosophila easter

Activation of pro-phenol oxidase (proPO) in insects and crustaceans is important in defense against wounding and infection. The proPO zymogen is activated by a specific proteolytic cleavage. PO oxidizes phenolic compounds to produce quinones, which may help to kill pathogens and can also be used for synthesis of melanin to seal wounds and encapsulate parasites. We have isolated from the tobacco hornworm, Manduca sexta, a serine proteinase that activates proPO, and have cloned its cDNA. The isolated proPO activating proteinase (PAP) hydrolyzed artificial substrates but required other protein factors for proPO activation, suggesting that proPO-activating enzyme may exist as a protein complex, one component of which is PAP. PAP (44 kDa) is composed of two disulfide-linked polypeptide chains (31 kDa and 13 kDa). A cDNA for PAP was isolated from a hemocyte library, by using a PCR-generated probe based on the amino-terminal amino acid sequence of the 31-kDa catalytic domain. PAP belongs to a family of arthropod serine proteinases containing a carboxyl-terminal proteinase domain and an amino-terminal “clip” domain. The member of this family most similar in sequence to PAP is the product of the easter gene from Drosophila melanogaster. PAP mRNA was present at a low level in larval hemocytes and fat body, but became much more abundant in fat body after insects were injected with Escherichia coli. Sequence data and 3H-diisopropyl fluorphosphate labeling results suggest that the same PAP exists in hemolymph and cuticle.

Calibrating bacterial evolution

Attempts to calibrate bacterial evolution have relied on the assumption that rates of molecular sequence divergence in bacteria are similar to those of higher eukaryotes, or to those of the few bacterial taxa for which ancestors can be reliably dated from ecological or geological evidence. Despite similarities in the substitution rates estimated for some lineages, comparisons of the relative rates of evolution at different classes of nucleotide sites indicate no basis for their universal application to all bacteria. However, there is evidence that bacteria have a constant genome-wide mutation rate on an evolutionary time scale but that this rate differs dramatically from the rate estimated by experimental methods.

Respiratory microbiota resistance and resilience to pulmonary exacerbation and subsequent antimicrobial intervention

This study was supported by the UK Natural Environment Research Council (NE/H019456/1) to CJvdG, by the Wellcome Trust (WT 098051) to AWW and JP for sequencing costs, and by The Anna Trust (KB2008) to KDB. AWW and The Rowett Institute of Nutrition and Health, University of Aberdeen, receive core funding support from the Scottish Government Rural and Environmental Science and Analysis Service (RESAS). We thank Paul Scott, Richard Rance and the Wellcome Trust Sanger Institute’s sequencing team for generating 16S rRNA gene sequence data.

Bootstrap confidence levels for phylogenetic trees

Evolutionary trees are often estimated from DNA or RNA sequence data. How much confidence should we have in the estimated trees? In 1985, Felsenstein [Felsenstein, J. (1985) Evolution 39, 783–791] suggested the use of the bootstrap to answer this question. Felsenstein’s method, which in concept is a straightforward application of the bootstrap, is widely used, but has been criticized as biased in the genetics literature. This paper concerns the use of the bootstrap in the tree problem. We show that Felsenstein’s method is not biased, but that it can be corrected to better agree with standard ideas of confidence levels and hypothesis testing. These corrections can be made by using the more elaborate bootstrap method presented here, at the expense of considerably more computation.

Duplicated genes evolve independently after polyploid formation in cotton

Of the many processes that generate gene duplications, polyploidy is unique in that entire genomes are duplicated. This process has been important in the evolution of many eukaryotic groups, and it occurs with high frequency in plants. Recent evidence suggests that polyploidization may be accompanied by rapid genomic changes, but the evolutionary fate of discrete loci recently doubled by polyploidy (homoeologues) has not been studied. Here we use locus-specific isolation techniques with comparative mapping to characterize the evolution of homoeologous loci in allopolyploid cotton (Gossypium hirsutum) and in species representing its diploid progenitors. We isolated and sequenced 16 loci from both genomes of the allopolyploid, from both progenitor diploid genomes and appropriate outgroups. Phylogenetic analysis of the resulting 73.5 kb of sequence data demonstrated that for all 16 loci (14.7 kb/genome), the topology expected from organismal history was recovered. In contrast to observations involving repetitive DNAs in cotton, there was no evidence of interaction among duplicated genes in the allopolyploid. Polyploidy was not accompanied by an obvious increase in mutations indicative of pseudogene formation. Additionally, differences in rates of divergence among homoeologues in polyploids and orthologues in diploids were indistinguishable across loci, with significant rate deviation restricted to two putative pseudogenes. Our results indicate that most duplicated genes in allopolyploid cotton evolve independently of each other and at the same rate as those of their diploid progenitors. These indications of genic stasis accompanying polyploidization provide a sharp contrast to recent examples of rapid genomic evolution in allopolyploids.

Ribosomes can slide over and beyond “hungry” codons, resuming protein chain elongation many nucleotides downstream

In cells subjected to moderate aminoacyl-tRNA limitation, the peptidyl-tRNA–ribosome complex stalled at the “hungry” codon can slide well beyond it on the messenger RNA and resume translation further downstream. This behavior is proved by unequivocal amino acid sequence data, showing a protein that lacks the bypassed sequence encoded between the hungry codon and specific landing sites. The landing sites are codons cognate to the anticodon of the peptidyl-tRNA. The efficiency of this behavior can be as high as 10–20% but declines with the length of the slide. Interposition of “trap” sites (nonproductive landing sites) in the bypassed region reduces the frequency of successful slides, confirming that the ribosome–peptidyl-tRNA complex passes through the untranslated region of the message. This behavior appears to be quite general: it can occur at the two kinds of hungry codons tested, AUA and AAG; the sliding peptidyl-tRNA can be any of three species tested, phenylalanine, tyrosine, or leucine tRNA; the peptidyl component can be either of two very different peptide sequences; and translation can resume at any of the three codons tested.

Coalescent estimates of HIV-1 generation time in vivo

The generation time of HIV Type 1 (HIV-1) in vivo has previously been estimated using a mathematical model of viral dynamics and was found to be on the order of one to two days per generation. Here, we describe a new method based on coalescence theory that allows the estimate of generation times to be derived by using nucleotide sequence data and a reconstructed genealogy of sequences obtained over time. The method is applied to sequences obtained from a long-term nonprogressing individual at five sampling occasions. The estimate of viral generation time using the coalescent method is 1.2 days per generation and is close to that obtained by mathematical modeling (1.8 days per generation), thus strengthening confidence in estimates of a short viral generation time. Apart from the estimation of relevant parameters relating to viral dynamics, coalescent modeling also allows us to simulate the evolutionary behavior of samples of sequences obtained over time.

Hypoxia-induced gene expression profiling in the euryoxic fish Gillichthys mirabilis

Hypoxia is important in both biomedical and environmental contexts and necessitates rapid adaptive changes in metabolic organization. Mammals, as air breathers, have a limited capacity to withstand sustained exposure to hypoxia. By contrast, some aquatic animals, such as certain fishes, are routinely exposed and resistant to severe environmental hypoxia. Understanding the changes in gene expression in fishes exposed to hypoxic stress could reveal novel mechanisms of tolerance that may shed new light on hypoxia and ischemia in higher vertebrates. Using cDNA microarrays, we have studied gene expression in a hypoxia-tolerant burrow-dwelling goby fish, Gillichthys mirabilis. We show that a coherent picture of a complex transcriptional response can be generated for a nonmodel organism for which sequence data were unavailable. We demonstrate that: (i) although certain shifts in gene expression mirror changes in mammals, novel genes are differentially expressed in fish; and (ii) tissue-specific patterns of expression reflect the different metabolic roles of tissues during hypoxia.

RISSC: a novel database for ribosomal 16S–23S RNA genes spacer regions

A novel database, under the acronym RISSC (Ribosomal Intergenic Spacer Sequence Collection), has been created. It compiles more than 1600 entries of edited DNA sequence data from the 16S–23S ribosomal spacers present in most prokaryotes and organelles (e.g. mitochondria and chloroplasts) and is accessible through the Internet (http://ulises.umh.es/RISSC), where systematic searches for specific words can be conducted, as well as BLAST-type sequence searches. Additionally, a characteristic feature of this region, the presence/absence and nature of tRNA genes within the spacer, is included in all the entries, even when not previously indicated in the original database. All these combined features could provide a useful documen­tation tool for studies on evolution, identification, typing and strain characterization, among others.

PseudoBase: structural information on RNA pseudoknots

PseudoBase is a database containing structural, functional and sequence data related to RNA pseudo­knots. It can be reached at http://wwwbio.LeidenUniv.nl/∼Batenburg/PKB.html. For each pseudoknot, thirteen items are stored, for example the relevant sequence, the stem positions of the pseudoknot, the EMBL accession number of the sequence and the support that can be given regarding the reliability of the pseudo­knot. Since the last publication, information on sizes of the stems and the loops in the pseudoknots has been added. Also added are alternative entries that produce surveys of where the pseudoknots are, sorted according to stem size or loop size.

The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species

While genome sequencing projects are advancing rapidly, EST sequencing and analysis remains a primary research tool for the identification and categorization of gene sequences in a wide variety of species and an important resource for annotation of genomic sequence. The TIGR Gene Indices (http://www.tigr.org/tdb/tgi.shtml) are a collection of species-specific databases that use a highly refined protocol to analyze EST sequences in an attempt to identify the genes represented by that data and to provide additional information regarding those genes. Gene Indices are constructed by first clustering, then assembling EST and annotated gene sequences from GenBank for the targeted species. This process produces a set of unique, high-fidelity virtual transcripts, or Tentative Consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to mapping and genomic sequence data, to provide links between orthologous and paralogous genes and as a resource for comparative sequence analysis.