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).

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.

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.

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.