Cryptic species diversity in a widespread neotropical tree genus: The case of Cedrela odorata

PREMISE

Reconciling the use of taxonomy to partition morphological variation and describe genetic divergence within and among closely related species is a persistent challenge in phylogenetics. We reconstructed phylogenetic relationships among Cedrela odorata (Meliaceae) and five closely allied species to test the genetic basis for the current model of species delimitation in this economically valuable and threatened genus.

METHODS

We prepared a nuclear species tree with the program SNPhylo and 16,000 single-nucleotide polymorphisms from 168 Cedrela specimens. Based on clades present and ancestral patterns ADMIXTURE, we designed nine species delimitation models and compared each model to current taxonomy with Bayes factor delimitation. Timing of major lineage divergences was estimated with the program SNAPP.

RESULTS

The resulting analysis revealed that modern C. odorata evolved from two genetically distinct ancestral sources. All species delimitation models tested better fit the data than the model representing current taxonomic delimitation. Models with the greatest marginal likelihoods separated Mesoamerican C. odorata and South American C. odorata into two species and lumped C. angustifolia and C. montana as a single species. We estimated that Cedrela diversified in South America within the last 19 million years following one or more dispersal events from Mesoamerican lineages.

CONCLUSIONS

Our analyses show that the present taxonomic understanding within the genus obscures divergent lineages in C. odorata due in part to morphological differentiation and taxonomic distinctions that are not predictably associated with genetic divergence. A more accurate application of taxonomy to C. odorata and related species may aid in its conservation, management, and restoration efforts.

Plastome Structural Evolution and Homoplastic Inversions in Neo-Astragalus (Fabaceae)

The plastid genomes of photosynthetic green plants have largely maintained conserved gene content and order as well as structure over hundreds of millions of years of evolution. Several plant lineages, however, have departed from this conservation and contain many plastome structural rearrangements, which have been associated with an abundance of repeated sequences both overall and near rearrangement endpoints. We sequenced the plastomes of 25 taxa of Astragalus L. (Fabaceae), a large genus in the inverted repeat-lacking clade of legumes, to gain a greater understanding of the connection between repeats and plastome inversions. We found plastome repeat structure has a strong phylogenetic signal among these closely related taxa mostly in the New World clade of Astragalus called Neo-Astragalus. Taxa without inversions also do not differ substantially in their overall repeat structure from four taxa each with one large-scale inversion. For two taxa with inversion endpoints between the same pairs of genes, differences in their exact endpoints indicate the inversions occurred independently. Our proposed mechanism for inversion formation suggests the short inverted repeats now found near the endpoints of the four inversions may be there as a result of these inversions rather than their cause. The longer inverted repeats now near endpoints may have allowed the inversions first mediated by shorter microhomologous sequences to propagate, something that should be considered in explaining how any plastome rearrangement becomes fixed regardless of the mechanism of initial formation.

A draft genome and transcriptome of common milkweed (Asclepias syriaca) as resources for evolutionary, ecological, and molecular studies in milkweeds and Apocynaceae

Milkweeds (Asclepias) are used in wide-ranging studies including floral development, pollination biology, plant-insect interactions and co-evolution, secondary metabolite chemistry, and rapid diversification. We present a transcriptome and draft nuclear genome assembly of the common milkweed, Asclepias syriaca. This reconstruction of the nuclear genome is augmented by linkage group information, adding to existing chloroplast and mitochondrial genomic resources for this member of the Apocynaceae subfamily Asclepiadoideae. The genome was sequenced to 80.4× depth and the draft assembly contains 54,266 scaffolds ≥1 kbp, with N50 = 3,415 bp, representing 37% (156.6 Mbp) of the estimated 420 Mbp genome. A total of 14,474 protein-coding genes were identified based on transcript evidence, closely related proteins, and ab initio models, and 95% of genes were annotated. A large proportion of gene space is represented in the assembly, with 96.7% of Asclepias transcripts, 88.4% of transcripts from the related genus Calotropis, and 90.6% of proteins from Coffea mapping to the assembly. Scaffolds covering 75 Mbp of the Asclepias assembly formed 11 linkage groups. Comparisons of these groups with pseudochromosomes in Coffea found that six chromosomes show consistent stability in gene content, while one may have a long history of fragmentation and rearrangement. The progesterone 5β-reductase gene family, a key component of cardenolide production, is likely reduced in Asclepias relative to other Apocynaceae. The genome and transcriptome of common milkweed provide a rich resource for future studies of the ecology and evolution of a charismatic plant family.

eDNA as a tool for identifying freshwater species in sustainable forestry: A critical review and potential future applications

Environmental DNA (eDNA) is an emerging biological monitoring tool that can aid in assessing the effects of forestry and forest manufacturing activities on biota. Monitoring taxa across broad spatial and temporal scales is necessary to ensure forest management and forest manufacturing activities meet their environmental goals of maintaining biodiversity. Our objectives are to describe potential applications of eDNA across the wood products supply chain extending from regenerating forests, harvesting, and wood transport, to manufacturing facilities, and to review the current state of the science in this context. To meet our second objective, we summarize the taxa examined with targeted (PCR, qPCR or ddPCR) or metagenomic eDNA methods (eDNA metabarcoding), evaluate how estimated species richness compares between traditional field sampling and eDNA metabarcoding approaches, and compare the geographical representation of prior eDNA studies in freshwater ecosystems to global wood baskets. Potential applications of eDNA include evaluating the effects of forestry and forest manufacturing activities on aquatic biota, delineating fish-bearing versus non fish-bearing reaches, evaluating effectiveness of constructed road crossings for freshwater organism passage, and determining the presence of at-risk species. Studies using targeted eDNA approaches focused on fish, amphibians, and invertebrates, while metagenomic studies focused on fish, invertebrates, and microorganisms. Rare, threatened, or endangered species received the least attention in targeted eDNA research, but are arguably of greatest interest to sustainable forestry and forest manufacturing that seek to preserve freshwater biodiversity. Ultimately, using eDNA methods will enable forestry and forest manufacturing managers to have data-driven prioritization for conservation actions for all freshwater species.

Genomic resources for the Neotropical tree genus Cedrela (Meliaceae) and its relatives

BACKGROUND

Tree species in the genus Cedrela P. Browne are threatened by timber overexploitation across the Neotropics. Genetic identification of processed timber can be used to supplement wood anatomy to assist in the taxonomic and source validation of protected species and populations of Cedrela. However, few genetic resources exist that enable both species and source identification of Cedrela timber products. We developed several 'omic resources including a leaf transcriptome, organelle genome (cpDNA), and diagnostic single nucleotide polymorphisms (SNPs) that may assist the classification of Cedrela specimens to species and geographic origin and enable future research on this widespread Neotropical tree genus.

RESULTS

We designed hybridization capture probes to enrich for thousands of genes from both freshly preserved leaf tissue and from herbarium specimens across eight Meliaceae species. We first assembled a draft de novo transcriptome for C. odorata, and then identified putatively low-copy genes. Hybridization probes for 10,001 transcript models successfully enriched 9795 (98%) of these targets, and analysis of target capture efficiency showed that probes worked effectively for five Cedrela species, with each species showing similar mean on-target sequence yield and depth. The probes showed greater enrichment efficiency for Cedrela species relative to the other three distantly related Meliaceae species. We provide a set of candidate SNPs for species identification of four of the Cedrela species included in this analysis, and present draft chloroplast genomes for multiple individuals of eight species from four genera in the Meliaceae.

CONCLUSIONS

Deforestation and illegal logging threaten forest biodiversity globally, and wood screening tools offer enforcement agencies new approaches to identify illegally harvested timber. The genomic resources described here provide the foundation required to develop genetic screening methods for Cedrela species identification and source validation. Due to their transferability across the genus and family as well as demonstrated applicability for both fresh leaves and herbarium specimens, the genomic resources described here provide additional tools for studies examining the ecology and evolutionary history of Cedrela and related species in the Meliaceae.

Evolution at the tips: Asclepias phylogenomics and new perspectives on leaf surfaces

PREMISE OF THE STUDY

Leaf surface traits, such as trichome density and wax production, mediate important ecological processes such as anti-herbivory defense and water-use efficiency. We present a phylogenetic analysis of Asclepias plastomes as a framework for analyzing the evolution of trichome density and presence of epicuticular waxes.

METHODS

We produced a maximum-likelihood phylogeny using plastomes of 103 species of Asclepias. We reconstructed ancestral states and used model comparisons in a likelihood framework to analyze character evolution across Asclepias.

KEY RESULTS

We resolved the backbone of Asclepias, placing the Sonoran Desert clade and Incarnatae clade as successive sisters to the remaining species. We present novel findings about leaf surface evolution of Asclepias-the ancestor is reconstructed as waxless and sparsely hairy, a macroevolutionary optimal trichome density is supported, and the rate of evolution of trichome density has accelerated.

CONCLUSIONS

Increased sampling and selection of best-fitting models of evolution provide more resolved and robust estimates of phylogeny and character evolution than obtained in previous studies. Evolutionary inferences are more sensitive to character coding than model selection.

A time and a place for everything: phylogenetic history and geography as joint predictors of oak plastome phylogeny

Owing to high rates of introgressive hybridization, the plastid genome is poorly suited to fine-scale DNA barcoding and phylogenetic studies of the oak genus (Quercus, Fagaceae). At the tips of the oak plastome phylogeny, recent gene migration and reticulation generally cause topology to reflect geographic structure, while deeper branches reflect lineage divergence. In this study, we quantify the simple and partial effects of geographic proximity and nucleome-inferred phylogenetic history on oak plastome phylogeny at different evolutionary scales. Our study compares pairwise phylogenetic distances based on complete plastome sequences, pairwise phylogenetic distances from nuclear restriction site-associated DNA sequences (RADseq), and pairwise geographic distances for 34 individuals of the white oak clade representing 24 North American and Eurasian species. Within the North American white oak clade alone, phylogenetic history has essentially no effect on plastome variation, while geography explains 11%–21% of plastome phylogenetic variance. However, across multiple continents and clades, phylogeny predicts 30%–41% of plastome variation, geography 3%–41%. Tipwise attenuation of phylogenetic informativeness in the plastome means that in practical terms, plastome data has little use in solving phylogenetic questions, but can still be a useful barcoding or phylogenetic marker for resolving questions among major clades.

Hyb-Seq: Combining target enrichment and genome skimming for plant phylogenomics

PREMISE OF THE STUDY

Hyb-Seq, the combination of target enrichment and genome skimming, allows simultaneous data collection for low-copy nuclear genes and high-copy genomic targets for plant systematics and evolution studies. •

METHODS AND RESULTS

Genome and transcriptome assemblies for milkweed (Asclepias syriaca) were used to design enrichment probes for 3385 exons from 768 genes (>1.6 Mbp) followed by Illumina sequencing of enriched libraries. Hyb-Seq of 12 individuals (10 Asclepias species and two related genera) resulted in at least partial assembly of 92.6% of exons and 99.7% of genes and an average assembly length >2 Mbp. Importantly, complete plastomes and nuclear ribosomal DNA cistrons were assembled using off-target reads. Phylogenomic analyses demonstrated signal conflict between genomes. •

CONCLUSIONS

The Hyb-Seq approach enables targeted sequencing of thousands of low-copy nuclear exons and flanking regions, as well as genome skimming of high-copy repeats and organellar genomes, to efficiently produce genome-scale data sets for phylogenomics.

Horizontal transfer of DNA from the mitochondrial to the plastid genome and its subsequent evolution in milkweeds (apocynaceae)

Horizontal gene transfer (HGT) of DNA from the plastid to the nuclear and mitochondrial genomes of higher plants is a common phenomenon; however, plastid genomes (plastomes) are highly conserved and have generally been regarded as impervious to HGT. We sequenced the 158 kb plastome and the 690 kb mitochondrial genome of common milkweed (Asclepias syriaca [Apocynaceae]) and found evidence of intracellular HGT for a 2.4-kb segment of mitochondrial DNA to the rps2–rpoC2 intergenic spacer of the plastome. The transferred region contains an rpl2 pseudogene and is flanked by plastid sequence in the mitochondrial genome, including an rpoC2 pseudogene, which likely provided the mechanism for HGT back to the plastome through double-strand break repair involving homologous recombination. The plastome insertion is restricted to tribe Asclepiadeae of subfamily Asclepiadoideae, whereas the mitochondrial rpoC2 pseudogene is present throughout the subfamily, which confirms that the plastid to mitochondrial HGT event preceded the HGT to the plastome. Although the plastome insertion has been maintained in all lineages of Asclepiadoideae, it shows minimal evidence of transcription in A. syriaca and is likely nonfunctional. Furthermore, we found recent gene conversion of the mitochondrial rpoC2 pseudogene in Asclepias by the plastid gene, which reflects continued interaction of these genomes.

Navigating the tip of the genomic iceberg: Next-generation sequencing for plant systematics

PREMISE OF THE STUDY

Just as Sanger sequencing did more than 20 years ago, next-generation sequencing (NGS) is poised to revolutionize plant systematics. By combining multiplexing approaches with NGS throughput, systematists may no longer need to choose between more taxa or more characters. Here we describe a genome skimming (shallow sequencing) approach for plant systematics.

METHODS

Through simulations, we evaluated optimal sequencing depth and performance of single-end and paired-end short read sequences for assembly of nuclear ribosomal DNA (rDNA) and plastomes and addressed the effect of divergence on reference-guided plastome assembly. We also used simulations to identify potential phylogenetic markers from low-copy nuclear loci at different sequencing depths. We demonstrated the utility of genome skimming through phylogenetic analysis of the Sonoran Desert clade (SDC) of Asclepias (Apocynaceae).

KEY RESULTS

Paired-end reads performed better than single-end reads. Minimum sequencing depths for high quality rDNA and plastome assemblies were 40× and 30×, respectively. Divergence from the reference significantly affected plastome assembly, but relatively similar references are available for most seed plants. Deeper rDNA sequencing is necessary to characterize intragenomic polymorphism. The low-copy fraction of the nuclear genome was readily surveyed, even at low sequencing depths. Nearly 160000 bp of sequence from three organelles provided evidence of phylogenetic incongruence in the SDC.

CONCLUSIONS

Adoption of NGS will facilitate progress in plant systematics, as whole plastome and rDNA cistrons, partial mitochondrial genomes, and low-copy nuclear markers can now be efficiently obtained for molecular phylogenetics studies.

Transcriptome characterization and polymorphism detection between subspecies of big sagebrush (Artemisia tridentata)

Background

Big sagebrush (Artemisia tridentata) is one of the most widely distributed and ecologically important shrub species in western North America. This species serves as a critical habitat and food resource for many animals and invertebrates. Habitat loss due to a combination of disturbances followed by establishment of invasive plant species is a serious threat to big sagebrush ecosystem sustainability. Lack of genomic data has limited our understanding of the evolutionary history and ecological adaptation in this species. Here, we report on the sequencing of expressed sequence tags (ESTs) and detection of single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers in subspecies of big sagebrush.

Results

cDNA of A. tridentata sspp. tridentata and vaseyana were normalized and sequenced using the 454 GS FLX Titanium pyrosequencing technology. Assembly of the reads resulted in 20,357 contig consensus sequences in ssp. tridentata and 20,250 contigs in ssp. vaseyana. A BLASTx search against the non-redundant (NR) protein database using 29,541 consensus sequences obtained from a combined assembly resulted in 21,436 sequences with significant blast alignments (≤ 1e^-15). A total of 20,952 SNPs and 119 polymorphic SSRs were detected between the two subspecies. SNPs were validated through various methods including sequence capture. Validation of SNPs in different individuals uncovered a high level of nucleotide variation in EST sequences. EST sequences of a third, tetraploid subspecies (ssp. wyomingensis) obtained by Illumina sequencing were mapped to the consensus sequences of the combined 454 EST assembly. Approximately one-third of the SNPs between sspp. tridentata and vaseyana identified in the combined assembly were also polymorphic within the two geographically distant ssp. wyomingensis samples.

Conclusion

We have produced a large EST dataset for Artemisia tridentata, which contains a large sample of the big sagebrush leaf transcriptome. SNP mapping among the three subspecies suggest the origin of ssp. wyomingensis via mixed ancestry. A large number of SNP and SSR markers provide the foundation for future research to address questions in big sagebrush evolution, ecological genetics, and conservation using genomic approaches.

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

Building a model: developing genomic resources for common milkweed (Asclepias syriaca) with low coverage genome sequencing

Background

Milkweeds (Asclepias L.) have been extensively investigated in diverse areas of evolutionary biology and ecology; however, there are few genetic resources available to facilitate and compliment these studies. This study explored how low coverage genome sequencing of the common milkweed (Asclepias syriaca L.) could be useful in characterizing the genome of a plant without prior genomic information and for development of genomic resources as a step toward further developing A. syriaca as a model in ecology and evolution.

Results

A 0.5× genome of A. syriaca was produced using Illumina sequencing. A virtually complete chloroplast genome of 158,598 bp was assembled, revealing few repeats and loss of three genes: accD, clpP, and ycf1. A nearly complete rDNA cistron (18S-5.8S-26S; 7,541 bp) and 5S rDNA (120 bp) sequence were obtained. Assessment of polymorphism revealed that the rDNA cistron and 5S rDNA had 0.3% and 26.7% polymorphic sites, respectively. A partial mitochondrial genome sequence (130,764 bp), with identical gene content to tobacco, was also assembled. An initial characterization of repeat content indicated that Ty1/copia-like retroelements are the most common repeat type in the milkweed genome. At least one A. syriaca microread hit 88% of Catharanthus roseus (Apocynaceae) unigenes (median coverage of 0.29×) and 66% of single copy orthologs (COSII) in asterids (median coverage of 0.14×). From this partial characterization of the A. syriaca genome, markers for population genetics (microsatellites) and phylogenetics (low-copy nuclear genes) studies were developed.

Conclusions

The results highlight the promise of next generation sequencing for development of genomic resources for any organism. Low coverage genome sequencing allows characterization of the high copy fraction of the genome and exploration of the low copy fraction of the genome, which facilitate the development of molecular tools for further study of a target species and its relatives. This study represents a first step in the development of a community resource for further study of plant-insect co-evolution, anti-herbivore defense, floral developmental genetics, reproductive biology, chemical evolution, population genetics, and comparative genomics using milkweeds, and A. syriaca in particular, as ecological and evolutionary models.

Length polymorphism scanning is an efficient approach for revealing chloroplast DNA variation

Phylogeographic and population genetic screens of chloroplast DNA (cpDNA) provide insights into seed-based gene flow in angiosperms, yet studies are frequently hampered by the low mutation rate of this genome. Detection methods for intraspecific variation can be either direct (DNA sequencing) or indirect (PCR–RFLP), although no single method incorporates the best features of both approaches. We show that screening universal chloroplast amp li cons for length polymorphism provides an accurate and efficient method for identifying cpDNA variation. By sequencing 4500 bp of cpDNA from 17 accessions of Purshia tridentata (bitterbrush), we detected 9 haplotypes, 8 of which were identifiable by unique multilocus length combinations resolvable by automated fragment analysis. In silico estimates of PCR–RFLP for these loci show that 5 haplotypes would be resolved by agarose electrophoresis. A survey of 4 intraspecific data sets from diverse angiosperms revealed that length variation in cpDNA amplicons is nearly ubiquitous, and 61 of 67 haplotypes identified by direct sequencing could be identified by screening length variation. Combined with automated fluorescent detection, length polymorphism screening of universal cpDNA regions offers a simple screen for intraspecific variation that can be used across angiosperms with minimal optimization, providing detection limits that rival direct sequencing at a fraction of the cost.Key words: cpDNA, intraspecific polymorphism, population genetics, phylogeography, indels.

Quantitative analysis of transcript accumulation from genes duplicated by polyploidy using cDNA-SSCP

Repeated rounds of polyploidy have been commonplace in the lineages leading to modern eukaryotic genomes, giving rise to widespread gene duplication. Genes duplicated by polyploidy, or "homoeologs," may continue to be expressed at equal levels following polyploidization or their expression may be dramatically altered. In this report, we describe how SSCP analysis of RT-PCR products can be used to evaluate the expression status (presence and relative quantity) of highly similar homoelogous gene pairs from an allotetraploid genome. This cDNA-SSCP approach was used to evaluate transcript abundance in "synthetic tetraploid" mRNA pools (i.e., mixtures of diploid products) and three natural homoeologous gene pairs expressed in tetraploid cotton (Gossypium hirsutum) ovules. Results from replicated tests show that cDNA-SSCP reliably separates duplicated transcripts with 99% sequence identity. Most significantly, the method yields quantitative estimates of transcript ratios in template pools that range from equimolar to approximately 100:1.

Rate variation among nuclear genes and the age of polyploidy in Gossypium

Molecular evolutionary rate variation in Gossypium (cotton) was characterized using sequence data for 48 nuclear genes from both genomes of allotetraploid cotton, models of its diploid progenitors, and an outgroup. Substitution rates varied widely among the 48 genes, with silent and replacement substitution levels varying from 0.018 to 0.162 and from 0.000 to 0.073, respectively, in comparisons between orthologous Gossypium and outgroup sequences. However, about 90% of the genes had silent substitution rates spanning a more narrow threefold range. Because there was no evidence of rate heterogeneity among lineages for any gene and because rates were highly correlated in independent tests, evolutionary rate is inferred to be a property of each gene or its genetic milieu rather than the clade to which it belongs. Evidence from approximately 200,000 nucleotides (40,000 per genome) suggests that polyploidy in Gossypium led to a modest enhancement in rates of nucleotide substitution. Phylogenetic analysis for each gene yielded the topology expected from organismal history, indicating an absence of gene conversion or recombination among homoeologs subsequent to allopolyploid formation. Using the mean synonymous substitution rate calculated across the 48 genes, allopolyploid cotton is estimated to have formed circa 1.5 million years ago (MYA), after divergence of the diploid progenitors about 6.7 MYA.

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.

Feast and famine in plant genomes

Plant genomes vary over several orders of magnitude in size, even among closely related species, yet the origin, genesis and significance of this variation are not clear. Because DNA content varies over a sevenfold range among diploid species in the cotton genus (Gossypium) and its allies, this group offers opportunities for exploring patterns and mechanisms of genome size evolution. For example, the question has been raised whether plant genomes have a 'one-way ticket to genomic obesity', as a consequence of retroelement accumulation. Few empirical studies directly address this possibility, although it is consistent with recent insights gleaned from evolutionary genomic investigations. We used a phylogenetic approach to evaluate the directionality of genome size evolution among Gossypium species and their relatives in the cotton tribe (Gossypieae, Malvaceae). Our results suggest that both DNA content increase and decrease have occurred repeatedly during evolution. In contrast to a model of unidirectional genome size change, the frequency of inferred genome size contraction exceeded that of expansion. In conjunction with other evidence, this finding highlights the dynamic nature of plant genome size evolution, and suggests that poorly understood genomic contraction mechanisms operate on a more extensive scale that previously recognized. Moreover, the research sets the stage for fine-scale analysis of the evolutionary dynamics and directionality of change for the full spectrum of genomic constituents.

Rapid diversification of the cotton genus (Gossypium: Malvaceae) revealed by analysis of sixteen nuclear and chloroplast genes

Previous molecular phylogenetic studies have failed to resolve the branching order among the major cotton (Gossypium) lineages, and it has been unclear whether this reflects actual history (rapid radiation) or sampling properties of the genes evaluated. In this paper, we reconsider the phylogenetic relationships of diploid cotton genome groups using DNA sequences from 11 single-copy nuclear loci (10 293 base pairs [bp]), nuclear ribosomal DNA (695 bp), and four chloroplast loci (7370 bp). Results from individual loci and combined nuclear and chloroplast DNA partitions reveal that the cotton genome groups radiated in rapid succession following the formation of the genus. Maximum likelihood analysis of nuclear synonymous sites shows that this radiation occurred within a time span equivalent to 17% of the time since the separation of Gossypium from its nearest extant relatives in the genera Kokia and Gossypioides. Chloroplast and nuclear phylogenies differ significantly with respect to resolution of the basal divergence in the genus and to interrelationships among African cottons. This incongruence is due to limited character evolution in cpDNA and either previously unsuspected hybridization or unreliable phylogenetic performance of the cpDNA characters. This study highlights the necessity of using multiple, independent data sets for resolving phylogenetic relationships of rapidly diverged lineages.

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.

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.Key words: polyploidy, genome evolution, cotton, Gossypium, amplified fragment length polymorphism (AFLP).

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.