Variable effects of captivity on microbiomes in populations of IUCN-endangered Blanding's turtles (Emydoidea blandingii).

AIMS: Microbiome composition is increasingly considered in species reintroduction efforts and may influence survival and reproductive success. Many turtle species are threatened by anthropogenic pressures and are frequently raised in captivity for reintroduction efforts, yet little is known about turtle microbiome composition in either wild or captive settings. Here, we investigated trends in microbiome composition of captive and wild IUCN-endangered Blanding's turtles (Emydoidea blandingii). METHODS AND RESULTS: We amplified and sequenced the V4 region of the 16S rDNA locus from plastron, cloaca, and water samples of wild E. blandingii adults and two populations of captive E. blandingii juveniles being raised for headstarting. Plastron, cloaca, and water-associated microbiomes differed strongly from each other and were highly variable among captive sites and between captive and wild sites. Across plastron, cloaca, and water-associated microbial communities, microbial diversity changed over time, but not in a predictable direction between captive sites. Plastron beta diversity correlated with growth rate in captive samples, indicating that external microbiomes may correlate with individual fitness. CONCLUSIONS: Our results indicate that external and internal microbiomes vary between captive and wild turtles and may reflect differences in fitness of captive-raised individuals.

Functional Dynamics within the Human Ribosome Regulate the Rate of Active Protein Synthesis.

The regulation of protein synthesis contributes to gene expression in both normal physiology and disease, yet kinetic investigations of the human translation mechanism are currently lacking. Using single-molecule fluorescence imaging methods, we have quantified the nature and timing of structural processes in human ribosomes during single-turnover and processive translation reactions. These measurements reveal that functional complexes exhibit dynamic behaviors and thermodynamic stabilities distinct from those observed for bacterial systems. Structurally defined sub-states of pre- and post-translocation complexes were sensitive to specific inhibitors of the eukaryotic ribosome, demonstrating the utility of this platform to probe drug mechanism. The application of three-color single-molecule fluorescence resonance energy transfer (smFRET) methods further revealed a long-distance allosteric coupling between distal tRNA binding sites within ribosomes bearing three tRNAs, which contributed to the rate of processive translation.

Implications of sequence variation on the evolution of rRNA.

The evolution of the multi-copy family of ribosomal RNA (rRNA) genes is unique in regard to its genetics and genome evolution. Paradoxically, rRNA genes are highly homogenized within and between individuals, yet they are globally distinct between species. Here, we discuss the implications for models of rRNA gene evolution in light of our recent discoveries that ribosomes bearing rRNA sequence variants can affect gene expression and physiology and that intra-individual rRNA alleles exhibit both context- and tissue-specific expression.

Correction to: Implications of sequence variation on the evolution of rRNA.

The article Implications of sequence variation on the evolution of rRNA, written by Matthew M. Parks, Chad M. Kurylo, Jake E. Batchelder, C. Theresa Vincent and Scott C. Blanchard, was originally published electronically on the publisher's internet portal (currently SpringerLink).

Signal, Uncertainty, and Conflict in Phylogenomic Data for a Diverse Lineage of Microbial Eukaryotes (Diatoms, Bacillariophyta).

Diatoms (Bacillariophyta) are a species-rich group of eukaryotic microbes diverse in morphology, ecology, and metabolism. Previous reconstructions of the diatom phylogeny based on one or a few genes have resulted in inconsistent resolution or low support for critical nodes. We applied phylogenetic paralog pruning techniques to a data set of 94 diatom genomes and transcriptomes to infer perennially difficult species relationships, using concatenation and summary-coalescent methods to reconstruct species trees from data sets spanning a wide range of thresholds for taxon and column occupancy in gene alignments. Conflicts between gene and species trees decreased with both increasing taxon occupancy and bootstrap cutoffs applied to gene trees. Concordance between gene and species trees was lowest for short internodes and increased logarithmically with increasing edge length, suggesting that incomplete lineage sorting disproportionately affects species tree inference at short internodes, which are a common feature of the diatom phylogeny. Although species tree topologies were largely consistent across many data treatments, concatenation methods appeared to outperform summary-coalescent methods for sparse alignments. Our results underscore that approaches to species-tree inference based on few loci are likely to be misled by unrepresentative sampling of gene histories, particularly in lineages that may have diversified rapidly. In addition, phylogenomic studies of diatoms, and potentially other hyperdiverse groups, should maximize the number of gene trees with high taxon occupancy, though there is clearly a limit to how many of these genes will be available.

ADAM10 controls the differentiation of the coronary arterial endothelium.

The coronary vasculature is crucial for normal heart function, yet much remains to be learned about its development, especially the maturation of coronary arterial endothelium. Here, we show that endothelial inactivation of ADAM10, a key regulator of Notch signaling, leads to defects in coronary arterial differentiation, as evidenced by dysregulated genes related to Notch signaling and arterial identity. Moreover, transcriptome analysis indicated reduced EGFR signaling in A10ΔEC coronary endothelium. Further analysis revealed that A10ΔEC mice have enlarged dysfunctional hearts with abnormal myocardial compaction, and increased expression of venous and immature endothelium markers. These findings provide the first evidence for a potential role for endothelial ADAM10 in cardioprotective homeostatic EGFR signaling and implicate ADAM10/Notch signaling in coronary arterial cell specification, which is vital for normal heart development and function. The ADAM10/Notch signaling pathway thus emerges as a potential therapeutic target for improving the regenerative capacity and maturation of the coronary vasculature.

An exact test for comparing a fixed quantitative property between gene sets.

Motivation: A significant difference in the distribution of a feature between two gene sets can provide insight into function or regulation. This statistical setting differs from much of hypothesis testing theory because the genome is often considered to be effectively fixed, finite and entirely known in commonly studied organisms, such as human. The Mann-Whitney U test is commonly employed in this scenario despite the assumptions of the test not being met, leading to unreliable and generally underpowered results. Permutation tests are also commonly employed for this purpose, but are computationally burdensome and are not tractable for obtaining small P values or for multiple comparisons. Results: We present an exact test for the null hypothesis that gene set membership is independent of the quantitative gene feature of interest. We derive an analytic expression for the randomization distribution of the median of the quantitative feature under the null hypothesis. Efficient implementation permits calculation of precise P values of arbitrary magnitude and makes thousands of simultaneous tests of transcriptome-sized gene sets computationally tractable. The flexibility of the hypothesis testing framework presented permits extension to a variety of related tests commonly found in genomics. The exact test is used to identify signatures of translation control and protein function in the human genome. Availability and implementation: The exact test presented here is implemented in R in the package kpmt available on CRAN. Contact: map2085@med.cornell.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Wnt5a Signals through DVL1 to Repress Ribosomal DNA Transcription by RNA Polymerase I.

Ribosome biogenesis is essential for cell growth and proliferation and is commonly elevated in cancer. Accordingly, numerous oncogene and tumor suppressor signaling pathways target rRNA synthesis. In breast cancer, non-canonical Wnt signaling by Wnt5a has been reported to antagonize tumor growth. Here, we show that Wnt5a rapidly represses rDNA gene transcription in breast cancer cells and generates a chromatin state with reduced transcription of rDNA by RNA polymerase I (Pol I). These effects were specifically dependent on Dishevelled1 (DVL1), which accumulates in nucleolar organizer regions (NORs) and binds to rDNA regions of the chromosome. Upon DVL1 binding, the Pol I transcription activator and deacetylase Sirtuin 7 (SIRT7) releases from rDNA loci, concomitant with disassembly of Pol I transcription machinery at the rDNA promoter. These findings reveal that Wnt5a signals through DVL1 to suppress rRNA transcription. This provides a novel mechanism for how Wnt5a exerts tumor suppressive effects and why disruption of Wnt5a signaling enhances mammary tumor growth in vivo.

Using controls to limit false discovery in the era of big data.

BACKGROUND: Procedures for controlling the false discovery rate (FDR) are widely applied as a solution to the multiple comparisons problem of high-dimensional statistics. Current FDR-controlling procedures require accurately calculated p-values and rely on extrapolation into the unknown and unobserved tails of the null distribution. Both of these intermediate steps are challenging and can compromise the reliability of the results. RESULTS: We present a general method for controlling the FDR that capitalizes on the large amount of control data often found in big data studies to avoid these frequently problematic intermediate steps. The method utilizes control data to empirically construct the distribution of the test statistic under the null hypothesis and directly compares this distribution to the empirical distribution of the test data. By not relying on p-values, our control data-based empirical FDR procedure more closely follows the foundational principles of the scientific method: that inference is drawn by comparing test data to control data. The method is demonstrated through application to a problem in structural genomics. CONCLUSIONS: The method described here provides a general statistical framework for controlling the FDR that is specifically tailored for the big data setting. By relying on empirically constructed distributions and control data, it forgoes potentially problematic modeling steps and extrapolation into the unknown tails of the null distribution. This procedure is broadly applicable insofar as controlled experiments or internal negative controls are available, as is increasingly common in the big data setting.

Phylogenomics reveals an extensive history of genome duplication in diatoms (Bacillariophyta).

PREMISE OF THE STUDY: Diatoms are one of the most species-rich lineages of microbial eukaryotes. Similarities in clade age, species richness, and primary productivity motivate comparisons to angiosperms, whose genomes have been inordinately shaped by whole-genome duplication (WGD). WGDs have been linked to speciation, increased rates of lineage diversification, and identified as a principal driver of angiosperm evolution. We synthesized a large but scattered body of evidence that suggests polyploidy may be common in diatoms as well. METHODS: We used gene counts, gene trees, and distributions of synonymous divergence to carry out a phylogenomic analysis of WGD across a diverse set of 37 diatom species. KEY RESULTS: Several methods identified WGDs of varying age across diatoms. Determining the occurrence, exact number, and placement of events was greatly impacted by uncertainty in gene trees. WGDs inferred from synonymous divergence of paralogs varied depending on how redundancy in transcriptomes was assessed, gene families were assembled, and synonymous distances (Ks) were calculated. Our results highlighted a need for systematic evaluation of key methodological aspects of Ks-based approaches to WGD inference. Gene tree reconciliations supported allopolyploidy as the predominant mode of polyploid formation, with strong evidence for ancient allopolyploid events in the thalassiosiroid and pennate diatom clades. CONCLUSIONS: Our results suggest that WGD has played a major role in the evolution of diatom genomes. We outline challenges in reconstructing paleopolyploid events in diatoms that, together with these results, offer a framework for understanding the impact of genome duplication in a group that likely harbors substantial genomic diversity.

Impacts of low coverage depths and post-mortem DNA damage on variant calling: a simulation study.

BACKGROUND: Massively parallel sequencing platforms, featuring high throughput and relatively short read lengths, are well suited to ancient DNA (aDNA) studies. Variant identification from short-read alignment could be hindered, however, by low DNA concentrations common to historic samples, which constrain sequencing depths, and post-mortem DNA damage patterns. RESULTS: We simulated pairs of sequences to act as reference and sample genomes at varied GC contents and divergence levels. Short-read sequence pools were generated from sample sequences, and subjected to varying levels of "post-mortem" damage by adjusting levels of fragmentation and fragmentation biases, transition rates at sequence ends, and sequencing depths. Mapping of sample read pools to reference sequences revealed several trends, including decreased alignment success with increased read length and decreased variant recovery with increased divergence. Variants were generally called with high accuracy, however identification of SNPs (single-nucleotide polymorphisms) was less accurate for high damage/low divergence samples. Modest increases in sequencing depth resulted in rapid gains in total variant recovery, and limited improvements to recovery of heterozygous variants. CONCLUSIONS: This in silico study suggests aDNA-associated damage patterns minimally impact variant call accuracy and recovery from short-read alignment, while modest increases in sequencing depth can greatly improve variant recovery.

Comparison of nuclear, plastid, and mitochondrial phylogenies and the origin of wild octoploid strawberry species.

PREMISE OF THE STUDY: Molecular phylogenies derived from all three plant genomes can provide insight into the evolutionary history of plant groups influenced by reticulation. We sought to reconstruct mitochondrial exome, chloroplast, and nuclear genome phylogenies for octoploid Fragaria and their diploid ancestors and to document patterns of incongruence between and within the cytoplasmic genomes and interpret these in the context of evolutionary origin of the octoploid strawberries. METHODS: Using a genome-skimming approach, we assembled chloroplast genomes and mitochondrial exomes, and we used the POLiMAPS method to assemble nuclear sequence for octoploid species and constructed phylogenies from all three genomes. We assessed incongruence between and within cytoplasmic genomes using topology-based phylogenetic incongruence tests. KEY RESULTS: The incongruent cytoplasmic genome phylogeny with respect to the placement of octoploids suggests potential breakage in linkage disequilibrium of cytoplasmic genomes during allopolyploid origin of the octoploids. Furthermore, a single mitochondrial chimeric gene with a putative role in cytoplasmic male sterility yields a phylogeny that is inconsistent with the rest of the mitochondrial genome but consistent with the chloroplast phylogeny, suggesting intracellular gene transfer between heteroplasmic mitochondria, possibly driven by selection to overcome the effects of mito-nuclear incompatibility in octoploid origins. CONCLUSIONS: This work expands on the current understanding of evolutionary history of the octoploid ancestors of cultivated strawberry. It demonstrates phylogenetic incongruence between cytoplasmic genomes in octoploids with respect to diploid ancestors, indicating breakage in linkage disequilibrium of cytoplasmic genomes. We discuss potential organism-level processes that may have contributed to the observed incongruence in Fragaria.

Turtle species and ecology drive carapace microbiome diversity in three seasonally interconnected wetland habitats.

Different species of freshwater turtles exhibit primary behaviours ranging from aerial basking to benthic bottom-walking, cycle between wet and dry conditions at different time intervals, and undertake short-distance overland movements between aquatic habitats. These behaviours in turn may impact the accumulation of microbes on external shell surfaces of turtles and provide novel niches for differentiation of microbial communities. We assessed microbial diversity using 16S and 18S rRNA metabarcoding on carapace surfaces of six species of freshwater turtles residing in three adjacent and seasonally interconnected wetland habitats in southeast Oklahoma (United States). Communities were highly diverse, with nearly 4200 prokaryotic and 500 micro-eukaryotic amplicon sequence variants recovered, and included taxa previously reported as common or differentially abundant on turtle shells. The 16S rRNA alpha diversity tended to be highest for two species of benthic turtles, while 18S rRNA alpha diversity was highest for two basking and one shallow-water benthic species. Beta diversity of communities was more strongly differentiated by turtle species than by collection site, and ordination patterns were largely reflective of turtle species' primary habits (i.e. benthic, basking, or benthic-basking). Our data support that freshwater turtles could play a role in microbial ecology and evolution in freshwater habitats and warrant additional exploration including in areas with high native turtle diversity and inter-habitat turtle movements.

Ancient and Modern Genomes Reveal Microsatellites Maintain a Dynamic Equilibrium Through Deep Time.

Microsatellites are widely used in population genetics, but their evolutionary dynamics remain poorly understood. It is unclear whether microsatellite loci drift in length over time. This is important because the mutation processes that underlie these important genetic markers are central to the evolutionary models that employ microsatellites. We identify more than 27 million microsatellites using a novel and unique dataset of modern and ancient Adélie penguin genomes along with data from 63 published chordate genomes. We investigate microsatellite evolutionary dynamics over 2 timescales: one based on Adélie penguin samples dating to ∼46.5 ka and the other dating to the diversification of chordates aged more than 500 Ma. We show that the process of microsatellite allele length evolution is at dynamic equilibrium; while there is length polymorphism among individuals, the length distribution for a given locus remains stable. Many microsatellites persist over very long timescales, particularly in exons and regulatory sequences. These often retain length variability, suggesting that they may play a role in maintaining phenotypic variation within populations.

Bioclimatic, ecological, and phenotypic intermediacy and high genetic admixture in a natural hybrid of octoploid strawberries.

PREMISE OF THE STUDY: Hybrid zones provide "natural laboratories" for understanding the processes of selection, reinforcement, and speciation. We sought to gain insight into the degree of introgression and the extent of ecological-phenotypic intermediacy in the natural hybrid strawberry, Fragaria × ananassa subsp. cuneifolia. • METHODS: We used whole-plastome sequencing to identify parental species-specific (Fragaria chiloensis and F. virginiana) chloroplast single-nucleotide polymorphisms and combined the use of these with nuclear microsatellite markers to genetically characterize the hybrid zone. We assessed the potential role of selection in the observed geographic patterns by bioclimatically characterizing the niche of the hybrid populations and phenotypically characterizing hybrid individuals of known genomic constitution. • KEY RESULTS: Significant admixture and little overall maternal bias in chloroplast or nuclear genomes suggest a high degree of interfertility among the parental and hybrid species and point to a long history of backcrossing and genetic mixing in the hybrid zone. Even though hybrids were phenotypically intermediate to the parental species, there was a discernible fingerprint of the parental genotype within hybrid individuals. Thus, although the pattern of introgression observed suggests geographic limitations to gene flow, it may be reinforced by selection for specific parental traits in the bioclimatically intermediate habitat occupied by the hybrid. • CONCLUSIONS: This work uncovered the genetic complexity underlying the hybrid zone of the wild relatives of the cultivated strawberry. It lays the foundation for experimental dissection of the causes of genomic introgression and nuclear-cytoplasmic disassociation, and for understanding other parts of Fragaria evolutionary history.

Separating the wheat from the chaff: mitigating the effects of noise in a plastome phylogenomic data set from Pinus L. (Pinaceae).

BACKGROUND: Through next-generation sequencing, the amount of sequence data potentially available for phylogenetic analyses has increased exponentially in recent years. Simultaneously, the risk of incorporating 'noisy' data with misleading phylogenetic signal has also increased, and may disproportionately influence the topology of weakly supported nodes and lineages featuring rapid radiations and/or elevated rates of evolution. RESULTS: We investigated the influence of phylogenetic noise in large data sets by applying two fundamental strategies, variable site removal and long-branch exclusion, to the phylogenetic analysis of a full plastome alignment of 107 species of Pinus and six Pinaceae outgroups. While high overall phylogenetic resolution resulted from inclusion of all data, three historically recalcitrant nodes remained conflicted with previous analyses. Close investigation of these nodes revealed dramatically different responses to data removal. Whereas topological resolution and bootstrap support for two clades peaked with removal of highly variable sites, the third clade resolved most strongly when all sites were included. Similar trends were observed using long-branch exclusion, but patterns were neither as strong nor as clear. When compared to previous phylogenetic analyses of nuclear loci and morphological data, the most highly supported topologies seen in Pinus plastome analysis are congruent for the two clades gaining support from variable site removal and long-branch exclusion, but in conflict for the clade with highest support from the full data set. CONCLUSIONS: These results suggest that removal of misleading signal in phylogenomic datasets can result not only in increased resolution for poorly supported nodes, but may serve as a tool for identifying erroneous yet highly supported topologies. For Pinus chloroplast genomes, removal of variable sites appears to be more effective than long-branch exclusion for clarifying phylogenetic hypotheses.

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.

Targeted enrichment strategies for next-generation plant biology.

PREMISE OF THE STUDY: The dramatic advances offered by modern DNA sequencers continue to redefine the limits of what can be accomplished in comparative plant biology. Even with recent achievements, however, plant genomes present obstacles that can make it difficult to execute large-scale population and phylogenetic studies on next-generation sequencing platforms. Factors like large genome size, extensive variation in the proportion of organellar DNA in total DNA, polyploidy, and gene number/redundancy contribute to these challenges, and they demand flexible targeted enrichment strategies to achieve the desired goals. METHODS: In this article, we summarize the many available targeted enrichment strategies that can be used to target partial-to-complete organellar genomes, as well as known and anonymous nuclear targets. These methods fall under four categories: PCR-based enrichment, hybridization-based enrichment, restriction enzyme-based enrichment, and enrichment of expressed gene sequences. KEY RESULTS: Examples of plant-specific applications exist for nearly all methods described. While some methods are well established (e.g., transcriptome sequencing), other promising methods are in their infancy (hybridization enrichment). A direct comparison of methods shows that PCR-based enrichment may be a reasonable strategy for accessing small genomic targets (e.g., ≤50 kbp), but that hybridization and transcriptome sequencing scale more efficiently if larger targets are desired. CONCLUSIONS: While the benefits of targeted sequencing are greatest in plants with large genomes, nearly all comparative projects can benefit from the improved throughput offered by targeted multiplex DNA sequencing, particularly as the amount of data produced from a single instrument approaches a trillion bases per run.

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.

Newly developed primers for complete ycf1 amplification in Pinus (Pinaceae) chloroplasts with possible family-wide utility.

PREMISE OF THE STUDY: Primers were designed to amplify the highly variable locus ycf1 from all 11 subsections of Pinus to facilitate plastome assemblies based on short sequence reads as well as future phylogenetic and population genetic analyses. METHODS AND RESULTS: Primer design was based on alignment of 33 Pinus and four Pinaceae plastomes with mostly incomplete ycf1 sequences. Sanger sequencing of 12 Pinus accessions resulted in open reading frames ranging in size from 5.2 to 6.1 kbp. Highest sequence diversity was identified in two regions totaling 5.5 kbp aligned length, which can be targeted in all pine subsections with three primer combinations. Preliminary results suggest the primers described also amplify homologous targets in the broader Pinaceae. CONCLUSIONS: The successful design and implementation of PCR primers spanning the large, variable locus ycf1 in Pinus represents the development of a valuable tool in pine genetic studies, and should facilitate studies throughout Pinaceae.