Transcriptome features of stone cell development in weevil-resistant and susceptible Sitka spruce.

Stone cells are a specialized, highly lignified cell type found in both angiosperms and gymnosperms. In conifers, abundance of stone cells in the cortex provides a robust constitutive physical defense against stem feeding insects. Stone cells are a major insect-resistance trait in Sitka spruce (Picea sitchensis), occurring in dense clusters in apical shoots of trees resistant (R) to spruce weevil (Pissodes strobi) but being rare in susceptible (S) trees. To learn more about molecular mechanisms of stone cell formation in conifers, we used laser microdissection and RNA sequencing to develop cell-type-specific transcriptomes of developing stone cells from R and S trees. Using light, immunohistochemical, and fluorescence microscopy, we also visualized the deposition of cellulose, xylan, and lignin associated with stone cell development. A total of 1293 genes were differentially expressed at higher levels in developing stone cells relative to cortical parenchyma. Genes with potential roles in stone cell secondary cell wall formation (SCW) were identified and their expression evaluated over a time course of stone cell formation in R and S trees. The expression of several transcriptional regulators was associated with stone cell formation, including a NAC family transcription factor and several genes annotated as MYB transcription factors with known roles in SCW formation.

Organellar Genomes of White Spruce (Picea glauca): Assembly and Annotation.

The genome sequences of the plastid and mitochondrion of white spruce (Picea glauca) were assembled from whole-genome shotgun sequencing data using ABySS. The sequencing data contained reads from both the nuclear and organellar genomes, and reads of the organellar genomes were abundant in the data as each cell harbors hundreds of mitochondria and plastids. Hence, assembly of the 123-kb plastid and 5.9-Mb mitochondrial genomes were accomplished by analyzing data sets primarily representing low coverage of the nuclear genome. The assembled organellar genomes were annotated for their coding genes, ribosomal RNA, and transfer RNA. Transcript abundances of the mitochondrial genes were quantified in three developmental tissues and five mature tissues using data from RNA-seq experiments. C-to-U RNA editing was observed in the majority of mitochondrial genes, and in four genes, editing events were noted to modify ACG codons to create cryptic AUG start codons. The informatics methodology presented in this study should prove useful to assemble organellar genomes of other plant species using whole-genome shotgun sequencing data.

Reduced Genetic Diversity and Increased Dispersal in Guigna (Leopardus guigna) in Chilean Fragmented Landscapes.

Landscape fragmentation is often a major cause of species extinction as it can affect a wide variety of ecological processes. The impact of fragmentation varies among species depending on many factors, including their life-history traits and dispersal abilities. Felids are one of the groups most threatened by fragmented landscapes because of their large home ranges, territorial behavior, and low population densities. Here, we model the impacts of habitat fragmentation on patterns of genetic diversity in the guigna (Leopardus guigna), a small felid that is closely associated with the heavily human-impacted temperate rainforests of southern South America. We assessed genetic variation in 1798 base pairs of mitochondrial DNA sequences, 15 microsatellite loci, and 2 sex chromosome genes and estimated genetic diversity, kinship, inbreeding, and dispersal in 38 individuals from landscapes with differing degrees of fragmentation on Chiloé Island in southern Chile. Increased fragmentation was associated with reduced genetic diversity, but not with increased kinship or inbreeding. However, in fragmented landscapes, there was a weaker negative correlation between pairwise kinship and geographic distance, suggesting increased dispersal distances. These results highlight the importance of biological corridors to maximize connectivity in fragmented landscapes and contribute to our understanding of the broader genetic consequences of habitat fragmentation, especially for forest-specialist carnivores.

Multiple mating and family structure of the western tent caterpillar, Malacosoma californicum pluviale: impact on disease resistance.

BACKGROUND: Levels of genetic diversity can strongly influence the dynamics and evolutionary changes of natural populations. Survival and disease resistance have been linked to levels of genetic diversity in eusocial insects, yet these relationships remain untested in gregarious insects where disease transmission can be high and selection for resistance is likely to be strong. METHODOLOGY/PRINCIPAL FINDINGS: Here we use 8 microsatellite loci to examine genetic variation in 12 families of western tent caterpillars, Malacosoma californicum pluviale from four different island populations to determine the relationship of genetic variability to survival and disease resistance. In addition these genetic markers were used to elucidate the population structure of western tent caterpillars. Multiple paternity was revealed by microsatellite markers, with the number of sires estimated to range from one to three per family (mean ± SE  = 1.92±0.23). Observed heterozygosity (H(O)) of families was not associated to the resistance of families to a nucleopolyhedrovirus (NPV) (r = 0.161, F(1,12)  = 0.271, P = 0.614), a major cause of mortality in high-density populations, but was positively associated with larval survival (r = 0.635, F(1,10)  = 5.412, P = 0.048). Genetic differentiation among the families was high (F(ST) = 0.269, P<0.0001), and families from the same island were as differentiated as were families from other islands. CONCLUSION/SIGNIFICANCE: We have been able to describe and characterize 8 microsatellite loci, which demonstrate patterns of variation within and between families of western tent caterpillars. We have discovered an association between larval survival and family-level heterozygosity that may be relevant to the population dynamics of this cyclic forest lepidopteran, and this will be the topic of future work.

Genetic analysis of cabbage loopers, Trichoplusia ni (Lepidoptera: Noctuidae), a seasonal migrant in western North America.

Long-range migrations of many wind-borne noctuid moths will have been influenced by the expansion of agriculture that provides greater availability of food plants along the migratory route. The migratory, agricultural pest, Trichoplusia ni (cabbage looper) over-winters in southern California and each summer migrates as far north as British Columbia. We explored the degree of genetic connectivity of populations over this migratory range. Preliminary investigation of seven mitochondrial gene regions found little to no variation among 13 populations, while partial regions of the NADH dehydrogenase subunits 1 and 4 in 42 individuals revealed eight and six haplotypes, respectively. The pattern of haplotype distribution indicated genetic homogeneity of persistent populations in California but weak differentiation among populations further north. Four highly variable amplified fragment length polymorphism primer combinations generated 167 polymorphic bands, with heterozygosity levels ranging from 0.250 to 0.302. Pairwise F ST values and clustering analyses also showed similarty among populations in California with some differentiation among populations initiated from the annual migration. Overall, some differentiation occurs among temporary, annual migratory populations but no pattern occurs with distance from the source population. Population subdivision in British Columbia associated with greenhouses has the greatest impact on genetic differentiation.

Spatial and temporal changes in genetic structure of greenhouse and field populations of cabbage looper, Trichoplusia ni.

Trichoplusia ni is a subtropical moth that migrates annually from southern California to southern British Columbia, Canada where it invades vegetable greenhouses and field crops. The heated greenhouse environment has altered the natural extinction-recolonization dynamics of T. ni populations, and allows year-round persistence in some locations. In addition, the extensive use of the biopesticide, Bacillus thuringiensis subspecies kurstaki (Bt) in some greenhouses has selected for resistance. Here we investigated the genetic structure of T. ni populations in British Columbia greenhouses and in field populations in California and British Columbia using amplified fragment length polymorphisms (AFLP) as related to patterns of Bt resistance. The majority of British Columbia field populations were similar to the California field populations, the potential source of migrants. However populations in two geographic areas with high concentrations of greenhouses showed local genetic differentiation. Some of these populations experienced severe bottlenecks over-winter and following Bt sprays. Greenhouse populations showed a pattern of isolation by distance and a strong positive relationship between genetic differentiation and levels of Bt resistance. These patterns indicate that greenhouses that sometimes support year-round populations of T. ni and the ensuing strong bottlenecking effects following winter cleanups and Bt application cause genetic differentiation of T. ni populations. Long distance migrants to field populations contribute to genetic homogeneity of these.

Development of 10 polymorphic microsatellite loci isolated from the mountain beaver, Aplodontia rufa rufa (Rafinesque).

We developed 10 microsatellite markers for the mountain beaver, Aplodontia rufa rufa. In three populations of A. r. rufa, the number of alleles for these loci ranged from monomorphic to nine. Average observed heterozygosities in these populations ranged from 0.29 to 0.60. We also tested previously published markers from the endangered subspecies A. r. nigra in A. r. rufa populations.

A conifer genomics resource of 200,000 spruce (Picea spp.) ESTs and 6,464 high-quality, sequence-finished full-length cDNAs for Sitka spruce (Picea sitchensis).

BACKGROUND: Members of the pine family (Pinaceae), especially species of spruce (Picea spp.) and pine (Pinus spp.), dominate many of the world's temperate and boreal forests. These conifer forests are of critical importance for global ecosystem stability and biodiversity. They also provide the majority of the world's wood and fiber supply and serve as a renewable resource for other industrial biomaterials. In contrast to angiosperms, functional and comparative genomics research on conifers, or other gymnosperms, is limited by the lack of a relevant reference genome sequence. Sequence-finished full-length (FL)cDNAs and large collections of expressed sequence tags (ESTs) are essential for gene discovery, functional genomics, and for future efforts of conifer genome annotation. RESULTS: As part of a conifer genomics program to characterize defense against insects and adaptation to local environments, and to discover genes for the production of biomaterials, we developed 20 standard, normalized or full-length enriched cDNA libraries from Sitka spruce (P. sitchensis), white spruce (P. glauca), and interior spruce (P. glauca-engelmannii complex). We sequenced and analyzed 206,875 3'- or 5'-end ESTs from these libraries, and developed a resource of 6,464 high-quality sequence-finished FLcDNAs from Sitka spruce. Clustering and assembly of 147,146 3'-end ESTs resulted in 19,941 contigs and 26,804 singletons, representing 46,745 putative unique transcripts (PUTs). The 6,464 FLcDNAs were all obtained from a single Sitka spruce genotype and represent 5,718 PUTs. CONCLUSION: This paper provides detailed annotation and quality assessment of a large EST and FLcDNA resource for spruce. The 6,464 Sitka spruce FLcDNAs represent the third largest sequence-verified FLcDNA resource for any plant species, behind only rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), and the only substantial FLcDNA resource for a gymnosperm. Our emphasis on capturing FLcDNAs and ESTs from cDNA libraries representing herbivore-, wound- or elicitor-treated induced spruce tissues, along with incorporating normalization to capture rare transcripts, resulted in a rich resource for functional genomics and proteomics studies. Sequence comparisons against five plant genomes and the non-redundant GenBank protein database revealed that a substantial number of spruce transcripts have no obvious similarity to known angiosperm gene sequences. Opportunities for future applications of the sequence and clone resources for comparative and functional genomics are discussed.

Genomics of hybrid poplar (Populus trichocarpax deltoides) interacting with forest tent caterpillars (Malacosoma disstria): normalized and full-length cDNA libraries, expressed sequence tags, and a cDNA microarray for the study of insect-induced defences in poplar.

As part of a genomics strategy to characterize inducible defences against insect herbivory in poplar, we developed a comprehensive suite of functional genomics resources including cDNA libraries, expressed sequence tags (ESTs) and a cDNA microarray platform. These resources are designed to complement the existing poplar genome sequence and poplar (Populus spp.) ESTs by focusing on herbivore- and elicitor-treated tissues and incorporating normalization methods to capture rare transcripts. From a set of 15 standard, normalized or full-length cDNA libraries, we generated 139,007 3'- or 5'-end sequenced ESTs, representing more than one-third of the c. 385,000 publicly available Populus ESTs. Clustering and assembly of 107,519 3'-end ESTs resulted in 14,451 contigs and 20,560 singletons, altogether representing 35,011 putative unique transcripts, or potentially more than three-quarters of the predicted c. 45,000 genes in the poplar genome. Using this EST resource, we developed a cDNA microarray containing 15,496 unique genes, which was utilized to monitor gene expression in poplar leaves in response to herbivory by forest tent caterpillars (Malacosoma disstria). After 24 h of feeding, 1191 genes were classified as up-regulated, compared to only 537 down-regulated. Functional classification of this induced gene set revealed genes with roles in plant defence (e.g. endochitinases, Kunitz protease inhibitors), octadecanoid and ethylene signalling (e.g. lipoxygenase, allene oxide synthase, 1-aminocyclopropane-1-carboxylate oxidase), transport (e.g. ABC proteins, calreticulin), secondary metabolism [e.g. polyphenol oxidase, isoflavone reductase, (-)-germacrene D synthase] and transcriptional regulation [e.g. leucine-rich repeat transmembrane kinase, several transcription factor classes (zinc finger C3H type, AP2/EREBP, WRKY, bHLH)]. This study provides the first genome-scale approach to characterize insect-induced defences in a woody perennial providing a solid platform for functional investigation of plant-insect interactions in poplar.

Robust simple sequence repeat markers for spruce (Picea spp.) from expressed sequence tags.

Traditionally, simple sequence repeat (SSR)markers have been developed from libraries of genomic DNA. However, the large, repetitive nature of conifer genomes makes development of robust, single-copy SSR markers from genomic DNA difficult. Expressed sequence tags (ESTs), or sequences of messenger RNA, offer the opportunity to exploit single, low-copy, conserved sequence motifs for SSR development. From a 20,275-unigene spruce EST set, we identified 44 candidate EST-SSR markers. Of these, 25 amplified and were polymorphic in white, Sitka, and black spruce; 20 amplified in all 3 spruce species tested; the remaining five amplified in all except one species. In addition, 101 previously described spruce SSRs (mostly developed from genomic DNA), were tested. Of these, 17 amplified across white,Sitka, and black spruce. The 25 EST-SSRs had approximately 9% less heterozygosity than the 17 genomic-derived SSRs (mean H=0.65 vs 0.72), but appeared to have less null alleles, as evidenced by much lower apparent inbreeding (mean F=0.046 vs 0.126). These robust SSRs are of particular use in comparative studies,and as the EST-SSRs are within the expressed portion of the genome, they are more likely to be associated with a particular gene of interest, improving their utility for quantitative trait loci mapping and allowing detection of selective sweeps at specific genes.

Single-copy, species-transferable microsatellite markers developed from loblolly pine ESTs.

Microsatellites, or simple sequence repeats (SSRs), are usually regarded as the "markers of choice" in population genetics research because they exhibit high variability. The development cost of these markers is usually high. In addition, microsatellite primers developed for one species often do not cross-amplify in related species, requiring separate development for each species. However, microsatellites found in expressed sequence tags (ESTs) might better cross-amplify as they reside in or near conserved coding DNA. In this study, we identified 14 Pinus taeda (loblolly pine) EST-SSRs from public EST databases and tested for their cross-species transferability to P. contorta ssp. latifolia, P. ponderosa, and P. sylvestris. As part of our development of a P. contorta microsatellite set, we also compared their transferability to that of 99 traditional microsatellite markers developed in P. taeda and tested on P. contorta ssp. latifolia. Compared to traditional microsatellites, EST-SSRs had higher transfer rates across pine species; however, the level of polymorphism of microsatellites derived from ESTs was lower. Sequence analyses revealed that the frequencies of insertions/deletions and base substitutions were lower in EST-SSRs than in other types of microsatellites, confirming that EST-SSRs are more conserved than traditional SSRs. Our results also provide a battery of 23 polymorphic, robust microsatellite primer pairs for lodgepole pine.