Multiplexed microsatellite recovery using massively parallel sequencing.

Conservation and management of natural populations requires accurate and inexpensive genotyping methods. Traditional microsatellite, or simple sequence repeat (SSR), marker analysis remains a popular genotyping method because of the comparatively low cost of marker development, ease of analysis and high power of genotype discrimination. With the availability of massively parallel sequencing (MPS), it is now possible to sequence microsatellite-enriched genomic libraries in multiplex pools. To test this approach, we prepared seven microsatellite-enriched, barcoded genomic libraries from diverse taxa (two conifer trees, five birds) and sequenced these on one lane of the Illumina Genome Analyzer using paired-end 80-bp reads. In this experiment, we screened 6.1 million sequences and identified 356,958 unique microreads that contained di- or trinucleotide microsatellites. Examination of four species shows that our conversion rate from raw sequences to polymorphic markers compares favourably to Sanger- and 454-based methods. The advantage of multiplexed MPS is that the staggering capacity of modern microread sequencing is spread across many libraries; this reduces sample preparation and sequencing costs to less than $400 (USD) per species. This price is sufficiently low that microsatellite libraries could be prepared and sequenced for all 1373 organisms listed as 'threatened' and 'endangered' in the United States for under $0.5 M (USD).

Finding a (pine) needle in a haystack: chloroplast genome sequence divergence in rare and widespread pines.

Critical to conservation efforts and other investigations at low taxonomic levels, DNA sequence data offer important insights into the distinctiveness, biogeographic partitioning and evolutionary histories of species. The resolving power of DNA sequences is often limited by insufficient variability at the intraspecific level. This is particularly true of studies involving plant organelles, as the conservative mutation rate of chloroplasts and mitochondria makes it difficult to detect polymorphisms necessary to track genealogical relationships among individuals, populations and closely related taxa, through space and time. Massively parallel sequencing (MPS) makes it possible to acquire entire organelle genome sequences to identify cryptic variation that would be difficult to detect otherwise. We are using MPS to evaluate intraspecific chloroplast-level divergence across biogeographic boundaries in narrowly endemic and widespread species of Pinus. We focus on one of the world's rarest pines - Torrey pine (Pinus torreyana) - due to its conservation interest and because it provides a marked contrast to more widespread pine species. Detailed analysis of nearly 90% ( approximately 105 000 bp each) of these chloroplast genomes shows that mainland and island populations of Torrey pine differ at five sites in their plastome, with the differences fixed between populations. This is an exceptionally low level of divergence (1 polymorphism/ approximately 21 kb), yet it is comparable to intraspecific divergence present in widespread pine species and species complexes. Population-level organelle genome sequencing offers new vistas into the timing and magnitude of divergence within species, and is certain to provide greater insight into pollen dispersal, migration patterns and evolutionary dynamics in plants.

Evolution and expression of MYB genes in diploid and polyploid cotton.

R2R3-MYB transcription factors have been implicated in a diversity of plant-specific processes. Among the functions attributed to myb factors is the determination of cell shape, including regulation of trichome length and density. Because myb transcription factors are likely to play a role in cotton fiber development, the molecular evolutionary properties of six MYB genes previously shown to be expressed in cotton fiber initiation were examined. In accordance with their presumed central role, each of the genes display conservative substitution patterns and limited sequence divergence in diploid members of the genus Gossypium, and this pattern is conserved in allotetraploid cottons. In contrast to highly reiterated rDNA repeats, GhMYB homologues (duplicated gene pairs) exhibit no evidence of concerted evolution, but instead appear to evolve independently in the allopolyploid nucleus. Expression patterns for the MYB genes were examined in several organs to determine if there have been changes in expression patterns between the diploids (G. raimondii and G. arboreum) and the tetraploid (G. hirsutum) or between the duplicated copies in the tetraploid. Spatial and temporal expression patterns appear to have been evolutionarily conserved, both during divergence of the diploid parents of allopolyploid cotton and following polyploid formation. However, the duplicated copies of MYB1 in the tetraploid are not expressed at equal levels or equivalently in all organs, suggesting possible functional differentiation.

PCR-mediated recombination in amplification products derived from polyploid cotton.

PCR recombination describes a process of in vitro chimera formation from non-identical templates. The key requirement of this process is the inclusion of two partially homologous templates in one reaction, a condition met when amplifying any locus from polyploid organisms and members of multigene families from diploid organisms. Because polyploids possess two or more divergent genomes ("homoeologues") in a common nucleus, intergenic chimeras can form during the PCR amplification of any gene. Here we report a high frequency of PCR-induced recombination for four low-copy genes from allotetraploid cotton ( Gossypium hirsutum). Amplification products from these genes ( Myb3, Myb5, G1262 and CesA1) range in length from 860 to 4,050 bp. Intergenomic recombinants were formed frequently, accounting for 23 of the 74 (31.1%) amplicons evaluated, with the frequency of recombination in individual reactions ranging from 0% to approximately 89%. Inspection of the putative recombination zones failed to reveal sequence-specific attributes that promote recombination. The high levels of observed in vitro recombination indicate that the tacit assumption of exclusive amplification of target templates may often be violated, particularly from polyploid genomes. This conclusion has profound implications for population and evolutionary genetic studies, where unrecognized artifactually recombinant molecules may bias results or alter interpretations.

Polyploid formation in cotton is not accompanied by rapid genomic changes.

Recent work has demonstrated that allopolyploid speciation in plants may be associated with non-Mendelian genomic changes in the early generations following polyploid synthesis. To address the question of whether rapid genomic changes also occur in allopolyploid cotton (Gossypium) species, amplified fragment length polymorphism (AFLP) analysis was performed to evaluate nine sets of newly synthesized allotetraploid and allohexaploid plants, their parents, and the selfed progeny from colchicine-doubled synthetics. Using both methylation-sensitive and methylation-insensitive enzymes, the extent of fragment additivity in newly combined genomes was ascertained for a total of approximately 22,000 genomic loci. Fragment additivity was observed in nearly all cases, with the few exceptions most likely reflecting parental heterozygosity or experimental error. In addition, genomic Southern analysis on six sets of synthetic allopolyploids probed with five retrotransposons also revealed complete additivity. Because no alterations were observed using methylation-sensitive isoschizomers, epigenetic changes following polyploid synthesis were also minimal. These indications of genomic additivity and epigenetic stasis during allopolyploid formation provide a contrast to recent evidence from several model plant allopolyploids, most notably wheat and Brassica, where rapid and unexplained genomic changes have been reported. In addition, the data contrast with evidence from repetitive DNAs in Gossypium, some of which are subject to non-Mendelian molecular evolutionary phenomena in extant polyploids. These contrasts indicate polyploid speciation in plants is accompanied by a diverse array of molecular evolutionary phenomena, which will vary among both genomic constituents and taxa.

Comparative development of fiber in wild and cultivated cotton.

One of the most striking examples of plant hairs is the single-celled epidermal seed trichome of cultivated cotton. The developmental morphology of these commercial "fibers" has been well-characterized in Gossypium hirsutum, but little is known about the pattern and tempo of fiber development in wild Gossypium species, all of which have short, agronomically inferior fiber. To identify developmental differences that account for variation in fiber length, and to place these differences in a phylogenetic context, we conducted SEM studies of ovules at and near the time of flowering, and generated growth curves for cultivated and wild diploid and tetraploid species. Trichome initiation was found to be similar in all taxa, with few notable differences in trichome density or early growth. Developmental profiles of the fibers of most wild species are similar, with fiber elongation terminating at about two weeks post-anthesis. In contrast, growth is extended to three weeks in the A- and F-genome diploids. This prolonged elongation period is diagnosed as a key evolutionary event in the origin of long fiber. A second evolutionary innovation is that absolute growth rate is higher in species with long fibers. Domestication of species is associated with a further prolongation of elongation at both the diploid and allopolyploid levels, suggesting the effects of parallel artificial selection. Comparative analysis of fiber growth curves lends developmental support to previous quantitative genetic suggestions that genes for fiber "improvement" in tetraploid cotton were contributed by the agronomically inferior D-genome diploid parent.

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.

The tortoise and the hare: choosing between noncoding plastome and nuclear Adh sequences for phylogeny reconstruction in a recently diverged plant group.

Phylogenetic resolution is often low within groups of recently diverged taxa due to a paucity of phylogenetically informative characters. We tested the relative utility of seven noncoding cpDNA regions and a pair of homoeologous nuclear genes for resolving recent divergences, using tetraploid cottons (Gossypium) as a model system. The five tetraploid species of Gossypium are a monophyletic assemblage derived from an allopolyploidization event that probably occurred within the last 0.5-2 million years. Previous analysis of cpDNA restriction site data provided only partial resolution within this clade despite a large number of enzymes employed. We sequenced three cpDNA introns (rpl16, rpoC1, ndhA) and four cpDNA spacers (accD-psaI, trnL-trnF, trnT-trnL, atpB-rbcL) for a total of over 7 kb of sequence per taxon, yet obtained only four informative nucleotide substitutions (0.05%) resulting in incomplete phylogenetic resolution. In addition, we sequenced a 1.65-kb region of a homoeologous pair of nuclear-encoded alcohol dehydrogenase (Adh) genes. In contrast with the cpDNA sequence data, the Adh homoeologues yielded 25 informative characters (0.76%) and provided a robust and completely resolved topology that is concordant with previous cladistic and phenetic analyses. The enhanced resolution obtained using the nuclear genes reflects an approximately three- to sixfold increase in nucleotide substitution rate relative to the plastome spacers and introns.

Polymorphism and concerted evolution in a tandemly repeated gene family: 5S ribosomal DNA in diploid and allopolyploid cottons.

5S RNA genes and their nontranscribed spacers are tandemly repeated in plant genomes at one or more chromosomal loci. To facilitate an understanding of the forces that govern 5S rDNA evolution, copy-number estimation and DNA sequencing were conducted for a phylogenetically well-characterized set of 16 diploid species of cotton (Gossypium) and 4 species representing allopolyploid derivatives of the diploids. Copy number varies over twentyfold in the genus, from approximately 1,000 to 20,000 copies/2C genome. When superimposed on the organismal phylogeny, these data reveal examples of both array expansion and contraction. Across species, a mean of 12% of nucleotide positions are polymorphic within individual arrays, for both gene and spacer sequences. This shows, in conjunction with phylogenetic evidence for ancestral polymorphisms that survive speciation events, that intralocus concerted evolutionary forces are relatively weak and that the rate of interrepeat homogenization is approximately equal to the rate of speciation. Evidence presented also shows that duplicated 5S rDNA arrays in allopolyploids have retained their subgenomic identity since polyploid formation, thereby indicating that interlocus concerted evolution has not been an important factor in the evolution of these arrays. A descriptive model, one which incorporates the opposing forces of mutation and homogenization within a selective framework, is outlined to account for the empirical data presented. Weak homogenizing forces allow equivalent levels of sequence polymorphism to accumulate in the 5S gene and spacer sequences, but fixation of mutations is nearly prohibited in the 5S gene. As a consequence, fixed interspecific differences are statistically underrepresented for 5S genes. This result explains the apparent paradox that despite similar levels of gene and spacer diversity, phylogenetic analysis of spacer sequences yields highly resolved trees, whereas analyses based on 5S gene sequences do not.

Inhibition of reverse transcriptase from feline immunodeficiency virus by analogs of 2'-deoxyadenosine-5'-triphosphate.

The replication of feline immunodeficiency virus (FIV) in cultured cells was inhibited by 2',3'-dideoxyadenosine (ddA) and by 9-(2-phosphonylmethoxyethyl)adenine (PMEA) with IC50 values of 0.98 and 0.95 microM, respectively. The effects of the presumed active forms of these inhibitors, ddATP and PMEA-diphosphate (PMEApp), upon the FIV reverse transcriptase (RT) were examined with two different template-primer systems. Both of these compounds were potent inhibitors of the FIV RT in reactions with primed phi X-174 DNA, yielding Ki values of 8.8 nM for ddATP and 5.0 nM for PMEApp. However, they were both poor inhibitors of the reaction with poly(rU)-oligo(dA); concentrations of ddATP or PMEApp greater than 10 microM were required to inhibit this reaction by 50%. Further analysis of the reaction with poly(rU)-oligo(dA) revealed that even in the absence of inhibitors the primers were extended by less than 20 nucleotides. In contrast, high molecular weight products were obtained in reactions with phi X-174 DNA. These results suggest that the reaction of FIV RT with poly(rU)-oligo(dA) is not highly processive. The high degree of termination encountered during this reaction with poly(rU)-oligo(dA) may be responsible for the low inhibitory potential of ddATP and PMEApp.

RNase H activity associated with reverse transcriptase from feline immunodeficiency virus.

Reverse transcription of retroviral genomes requires the action of an RNase H for template switching and primer generation. In this report, we compare enzymatic properties of the RNase H associated with the reverse transcriptase (RT) from feline immunodeficiency virus (FIV) and that from human immunodeficiency virus (HIV). Both enzymes displayed substrate preference for poly[3H](rG) . poly(dC) hybird over poly[3H](rA) . poly(dT) and cation preference for Mg2+ over Mn2+. Activity of the FIV RNase H upon poly(rG) . poly(dC) produced hydrolysis products from 1 to 6 nucleotides in length, similar to that reported for HIV. Dextran sulfates were effective inhibitors of both the FIV and HIV RNase H and RT activities. Nearly identical inhibition constants (0.12 nM) were obtained for all enzyme activities with dextran sulfate 500,000, while different inhibition constants were observed with dextran sulfate 8,000. Our results suggest that FIV and HIV RTs contain a conserved region that is sensitive to the larger dextran sulfate and that dextran sulfate 8,000 may interact at a different site or by a different mechanism.

Direct comparisons of inhibitor sensitivities of reverse transcriptases from feline and human immunodeficiency viruses.

The sensitivities of reverse transcriptases (RTs) from feline immunodeficiency virus (FIV) and human immunodeficiency virus type 1 (HIV) were directly compared. The two enzymes had similar sensitivities to three analogs of dTTP, namely, 3'-azido-3'-deoxythymidine 5'-triphosphate, 2',3'-dideoxythymidine 5'-triphosphate, and 2',3'-dideoxy-2',3'-didehydrothymidine 5'-triphosphate. Each of these analogs demonstrated competitive inhibition of both enzymes. Ki values for inhibition of FIV RT by these three inhibitors were 3.3, 6.7, and 1.8 nM, respectively; Ki values for inhibition of the HIV enzyme were 6.5, 5.9, and 8.3 nM, respectively. Ratios of the Ki for the inhibitor to the Km for the substrate were also determined for each inhibitor, and no differences between the two enzymes greater than threefold were observed. Inhibition constants for 3'-amino-3'-deoxythymidine 5'-triphosphate and 3'-fluoro-3'-deoxythymidine 5'-triphosphate were determined for FIV RT, and these were similar to published values for HIV RT. The activities of three dideoxynucleoside 5'-triphosphates against FIV RT were determined; ddGTP was slightly more potent than ddTTP, whereas both were much more effective than ddCTP. The activity of a noncompetitive inhibitor, phosphonoformate, was also examined with the FIV enzyme; it was much more active with poly(rA)-oligo(dT) as the template-primer than with poly(rC)-oligo(dG) or poly(rI)-oligo(dC).

Characterization of reverse transcriptase from feline immunodeficiency virus.

Reverse transcriptase has been purified from feline immunodeficiency virus (FIV) by DEAE-cellulose and phosphocellulose chromatography. The purified enzyme consists of a single protein with a Mr of 67,000. When proteolysis is not minimized during purification, a fragment of Mr 54,000 is also observed. This is similar to the reverse transcriptase from human immunodeficiency virus type 1 (HIV), which consists of a polypeptide of Mr 66,000; when proteolysis is not minimized during purification, a fragment of Mr 51,000 is also observed. In direct comparisons, the FIV reverse transcriptase is very similar to the HIV reverse transcriptase in template specificity and requirements for Mg2+. In contrast to these similarities, the FIV and HIV reverse transcriptases are substantially different in primary sequence, as determined by peptide mapping.