Fc γ receptor IIIa/CD16a processing correlates with the expression of glycan-related genes in human natural killer cells.

Many therapeutic monoclonal antibodies require binding to Fc γ receptors (FcγRs) for full effect and increasing the binding affinity increases efficacy. Preeminent among the five activating human FcγRs is FcγRIIIa/CD16a expressed by natural killer (NK) cells. CD16a is heavily processed, and recent reports indicate that the composition of the five CD16a asparagine(N)-linked carbohydrates (glycans) impacts affinity. These observations indicate that specific manipulation of CD16a N-glycan composition in CD16a-expressing effector cells including NK cells may improve treatment efficacy. However, it is unclear if modifying the expression of select genes that encode processing enzymes in CD16a-expressing effector cells is sufficient to affect N-glycan composition. We identified substantial processing differences using a glycoproteomics approach by comparing CD16a isolated from two NK cell lines, NK92 and YTS, with CD16a expressed by HEK293F cells and previous reports of CD16a from primary NK cells. Gene expression profiling by RNA-Seq and qRT-PCR revealed expression levels for glycan-modifying genes that correlated with CD16a glycan composition. These results identified a high degree of variability between the processing of the same human protein by different human cell types. N-glycan processing correlated with the expression of glycan-modifying genes and thus explained the substantial differences in CD16a processing by NK cells of different origins.

New data on Henneguya postexilis Minchew, 1977, a parasite of channel catfish Ictalurus punctatus, with notes on resolution of molecular markers for myxozoan phylogeny.

Previous morphological and histological data are supplemented with molecular and ultrastructural data for a Henneguya sp. isolated from farm-raised channel catfish Ictalurus punctatus in Mississippi, USA. Myxospores were cryptic, encapsulated within a thin layer of epithelium in the gill lamellae with spore measurements consistent with the original description of Henneguya postexilis Minchew, 1977. Myxospores were 42.7-49.1 µm in total length with spore bodies 12.1-17.2 × 3.6-4.8 × 2.9-3 µm. Polar capsules were of unequal length, with the longer capsule being 4.4-6.7 × 1.1-1.6 µm and the shorter capsule being 4.4-6.4 × 1.1-1.6 µm. Polar tubules had 6-8 turns. Caudal processes were 25.7-38.1 µm in length. Spores were encapsulated in a thin layer of epithelium in the gill lamellae. Molecular data from the most commonly used markers for myxozoan identification and phylogeny, partial 18S small subunit ribosomal gene (SSU), partial 28S large subunit ribosomal gene (LSU), and elongation factor 2 (EF2) were generated for H. postexilis. Additionally, novel data for LSU and EF2 were generated for archived myxozoan specimens from farm-raised catfish (H. mississippiensis, H. ictaluri, H. exilis, H. adiposa, H. sutherlandi, H. bulbosus, Unicauda fimbrethilae), as well as archived specimens from wild fish (H. laseeae [from Pylodictis olivaris], Hennegoides flockae [from Aphredoderus sayanus], Myxobolus cloutmani [from Cycleptus elongatus]. These include the first EF2 sequence data for the genera Hennegoides and Unicauda. Phylogenetic analyses using these data placed H. postexilis in well supported clades with other ictalurid-infecting Henneguya species. Phylogenetic signal assessments on these analyses suggest that while SSU provided the greatest phylogenetic signal, LSU yielded comparable signal, supporting previous work implying this region may be underutilised in reconstructing myxobolid phylogenies.

A unique insight into the MiRNA profile during genital chlamydial infection.

BACKGROUND: Genital C. trachomatis infection may cause pelvic inflammatory disease (PID) that can lead to tubal factor infertility (TFI). Understanding the pathogenesis of chlamydial complications including the pathophysiological processes within the female host genital tract is important in preventing adverse pathology. MicroRNAs regulate several pathophysiological processes of infectious and non-infectious etiologies. In this study, we tested the hypothesis that the miRNA profile of single and repeat genital chlamydial infections will be different and that these differences will be time dependent. Thus, we analyzed and compared differentially expressed mice genital tract miRNAs after single and repeat chlamydia infections using a C. muridarum mouse model. Mice were sacrificed and their genital tract tissues were collected at 1, 2, 4, and 8 weeks after a single and repeat chlamydia infections. Histopathology, and miRNA sequencing were performed. RESULTS: Histopathology presentation showed that the oviduct and uterus of reinfected mice were more inflamed, distended and dilated compared to mice infected once. The miRNAs expression profile was different in the reproductive tissues after a reinfection, with a greater number of miRNAs expressed after reinfection. Also, the number of miRNAs expressed each week after chlamydia infection and reinfection varied, with weeks eight and one having the highest number of differentially expressed miRNAs for chlamydia infection and reinfection respectively. Ten miRNAs; mmu-miR-378b, mmu-miR-204-5p, mmu-miR-151-5p, mmu-miR-142-3p, mmu-miR-128-3p, mmu-miR-335-3p, mmu-miR-195a-3p, mmu-miR-142-5p, mmu-miR-106a-5p and mmu-miR-92a-3p were common in both primary chlamydia infection and reinfection. Pathway analysis showed that, amongst other functions, the differentially regulated miRNAs control pathways involved in cellular and tissue development, disease conditions and toxicity. CONCLUSIONS: This study provides insights into the changes in miRNA expression over time after chlamydia infection and reinfection, as well as the pathways they regulate to determine pathological outcomes. The miRNAs networks generated in our study shows that there are differences in the focus molecules involved in significant biological functions in chlamydia infection and reinfection, implying that chlamydial pathogenesis occurs differently for each type of infection and that this could be important when determining treatments regime and disease outcome. The study underscores the crucial role of host factors in chlamydia pathogenesis.

Activation of defence pathways in Scots pine bark after feeding by pine weevil (Hylobius abietis).

BACKGROUND: During their lifetime, conifer trees are exposed to numerous herbivorous insects. To protect themselves against pests, trees have developed a broad repertoire of protective mechanisms. Many of the plant's defence reactions are activated upon an insect attack, and the underlying regulatory mechanisms are not entirely understood yet, in particular in conifer trees. Here, we present the results of our studies on the transcriptional response and the volatile compounds production of Scots pine (Pinus sylvestris) upon the large pine weevil (Hylobius abietis) feeding. RESULTS: Transcriptional response of Scots pine to the weevil attack was investigated using a novel customised 36.4 K Pinus taeda microarray. The weevil feeding caused large-scale changes in the pine transcriptome. In total, 774 genes were significantly up-regulated more than 4-fold (p≤0.05), whereas 64 genes were significantly down-regulated more than 4-fold. Among the up-regulated genes, we could identify genes involved in signal perception, signalling pathways, transcriptional regulation, plant hormone homeostasis, secondary metabolism and defence responses. The weevil feeding on stem bark of pine significantly increased the total emission of volatile organic compounds from the undamaged stem bark area. The emission levels of monoterpenes and sesquiterpenes were also increased. Interestingly, we could not observe any correlation between the increased production of the terpenoid compounds and expression levels of the terpene synthase-encoding genes. CONCLUSIONS: The obtained data provide an important insight into the transcriptional response of conifer trees to insect herbivory and illustrate the massive changes in the host transcriptome upon insect attacks. Moreover, many of the induced pathways are common between conifers and angiosperms. The presented results are the first ones obtained by the use of a microarray platform with an extended coverage of pine transcriptome (36.4 K cDNA elements). The platform will further facilitate the identification of resistance markers with the direct relevance for conifer tree breeding.

Noctilisin, a Venom Glycopeptide of Sirex noctilio (Hymenoptera: Siricidae), Causes Needle Wilt and Defense Gene Responses in Pines.

During oviposition, female Sirex noctilio (F.) (Siricidae) woodwasps inject their conifer hosts with a venom gland secretion. The secretion induces a variety of host physiological changes that facilitate subsequent lethal infection by a symbiotic fungus. A heat-stable factor that can migrate from the site of oviposition in the trunk through the xylem to needles in the crown of attacked pines was purified by size-fractionation and reversed-phase-high-performance liquid chromatography using activity assays based on defense gene induction as well as the needle wilt response in pine shoot explants. An 11-amino acid, posttranslationally modified peptide (SEGPROGTKRP) encoded by the most abundant transcript recovered from S. noctilio venom gland tissue comprised the backbone of the 1,850 Da active factor. Posttranslational modifications included hydroxylation of a Pro residue at position 6 as well as O-glycosylation of Ser and Thr residues at positions 1 and 8, respectively. The O-linked sugars were identical α-linked N-acetylgalactosamine residues modified at the C6 position by addition of phosphoethanolamine. In contrast to the native peptide, a synthetic version of the hydroxylated peptide backbone lacking the glycosyl side chains failed to induce pine defense genes or cause needle wilt in excised shoots. This peptide, hereafter called noctilisin, is related to the O-glycosylated short-chain proline-rich antimicrobial peptides exemplified by drosocin. The noctilisin structure contains motifs which may explain how it avoids detection by pine defense systems.

A SNP resource for Douglas-fir: de novo transcriptome assembly and SNP detection and validation.

BACKGROUND: Douglas-fir (Pseudotsuga menziesii), one of the most economically and ecologically important tree species in the world, also has one of the largest tree breeding programs. Although the coastal and interior varieties of Douglas-fir (vars. menziesii and glauca) are native to North America, the coastal variety is also widely planted for timber production in Europe, New Zealand, Australia, and Chile. Our main goal was to develop a SNP resource large enough to facilitate genomic selection in Douglas-fir breeding programs. To accomplish this, we developed a 454-based reference transcriptome for coastal Douglas-fir, annotated and evaluated the quality of the reference, identified putative SNPs, and then validated a sample of those SNPs using the Illumina Infinium genotyping platform. RESULTS: We assembled a reference transcriptome consisting of 25,002 isogroups (unique gene models) and 102,623 singletons from 2.76 million 454 and Sanger cDNA sequences from coastal Douglas-fir. We identified 278,979 unique SNPs by mapping the 454 and Sanger sequences to the reference, and by mapping four datasets of Illumina cDNA sequences from multiple seed sources, genotypes, and tissues. The Illumina datasets represented coastal Douglas-fir (64.00 and 13.41 million reads), interior Douglas-fir (80.45 million reads), and a Yakima population similar to interior Douglas-fir (8.99 million reads). We assayed 8067 SNPs on 260 trees using an Illumina Infinium SNP genotyping array. Of these SNPs, 5847 (72.5%) were called successfully and were polymorphic. CONCLUSIONS: Based on our validation efficiency, our SNP database may contain as many as ~200,000 true SNPs, and as many as ~69,000 SNPs that could be genotyped at ~20,000 gene loci using an Infinium II array-more SNPs than are needed to use genomic selection in tree breeding programs. Ultimately, these genomic resources will enhance Douglas-fir breeding and allow us to better understand landscape-scale patterns of genetic variation and potential responses to climate change.

Transcriptomic analysis highlights epigenetic and transcriptional regulation during zygotic embryo development of Pinus pinaster.

BACKGROUND: It is during embryogenesis that the plant body plan is established and the meristems responsible for all post-embryonic growth are specified. The molecular mechanisms governing conifer embryogenesis are still largely unknown. Their elucidation may contribute valuable information to clarify if the distinct features of embryo development in angiosperms and gymnosperms result from differential gene regulation. To address this issue, we have performed the first transcriptomic analysis of zygotic embryo development in a conifer species (Pinus pinaster) focusing our study in particular on regulatory genes playing important roles during plant embryo development, namely epigenetic regulators and transcription factors. RESULTS: Microarray analysis of P. pinaster zygotic embryogenesis was performed at five periods of embryo development from early developing to mature embryos. Our results show that most changes in transcript levels occurred in the first and the last embryo stage-to-stage transitions, namely early to pre-cotyledonary embryo and cotyledonary to mature embryo. An analysis of functional categories for genes that were differentially expressed through embryogenesis highlighted several epigenetic regulation mechanisms. While putative orthologs of transcripts associated with mechanisms that target transposable elements and repetitive sequences were strongly expressed in early embryogenesis, PRC2-mediated repression of genes seemed more relevant during late embryogenesis. On the other hand, functions related to sRNA pathways appeared differentially regulated across all stages of embryo development with a prevalence of miRNA functions in mid to late embryogenesis. Identification of putative transcription factor genes differentially regulated between consecutive embryo stages was strongly suggestive of the relevance of auxin responses and regulation of auxin carriers during early embryogenesis. Such responses could be involved in establishing embryo patterning. Later in development, transcripts with homology to genes acting on modulation of auxin flow and determination of adaxial-abaxial polarity were up-regulated, as were putative orthologs of genes required for meristem formation and function as well as establishment of organ boundaries. Comparative analysis with A. thaliana embryogenesis also highlighted genes involved in auxin-mediated responses, as well as epigenetic regulation, indicating highly correlated transcript profiles between the two species. CONCLUSIONS: This is the first report of a time-course transcriptomic analysis of zygotic embryogenesis in a conifer. Taken together our results show that epigenetic regulation and transcriptional control related to auxin transport and response are critical during early to mid stages of pine embryogenesis and that important events during embryogenesis seem to be coordinated by putative orthologs of major developmental regulators in angiosperms.

Linoleic acid-induced ANGPTL4 inhibits C2C12 skeletal muscle differentiation by suppressing Wnt/ß-catenin.

Skeletal muscle differentiation is an essential process in embryonic development as well as regeneration and repair throughout the lifespan. It is well-known that dietary fat intake impacts biological and physiological function in skeletal muscle, however, understanding of the contribution of nutritional factors in skeletal muscle differentiation is limited. Therefore, the objective of the current study was to evaluate the effects of free fatty acids (FFAs) on skeletal muscle differentiation in vitro. We used C2C12 murine myoblasts and treated them with various FFAs, which revealed a unique response of angiopoietin-like protein-4 (ANGPTL4) with linoleic acid (LA) treatment that was associated with reduced differentiation. LA significantly inhibited myotube formation and lowered the protein expression of myogenic regulatory factors, including MyoD and MyoG and increased Pax7 during cell differentiation. Next, recombinant ANGPTL4 protein or siRNA knockdown of ANGPTL4 was employed to examine its role in skeletal muscle differentiation, and we confirmed that ANGPTL4 knockdown at day two and six of differentiation restored myotube formation in the presence of LA. RNA-sequencing analysis revealed that ANGPTL4-mediated inhibition of skeletal muscle differentiation at day two as well as LA at day two or -6 led to a reduction in Wnt/ß-catenin signaling pathways. We confirmed that LA reduced Wnt11 and Axin2 while increasing expression of the Wnt inhibitor, Dkk2. ANGPTL4 knockdown increased ß-catenin protein in the nucleus in response to LA and increased Axin2 and Wnt11 expression. Taken together, these results demonstrate that LA induced ANGPTL4 inhibits C2C12 differentiation by suppressing Wnt/ß-catenin signaling.

The phenylalanine ammonia lyase (PAL) gene family shows a gymnosperm-specific lineage.

BACKGROUND: Phenylalanine ammonia lyase (PAL) is a key enzyme of the phenylpropanoid pathway that catalyzes the deamination of phenylalanine to trans-cinnamic acid, a precursor for the lignin and flavonoid biosynthetic pathways. To date, PAL genes have been less extensively studied in gymnosperms than in angiosperms. Our interest in PAL genes stems from their potential role in the defense responses of Pinus taeda, especially with respect to lignification and production of low molecular weight phenolic compounds under various biotic and abiotic stimuli. In contrast to all angiosperms for which reference genome sequences are available, P. taeda has previously been characterized as having only a single PAL gene. Our objective was to re-evaluate this finding, assess the evolutionary history of PAL genes across major angiosperm and gymnosperm lineages, and characterize PAL gene expression patterns in Pinus taeda. METHODS: We compiled a large set of PAL genes from the largest transcript dataset available for P. taeda and other conifers. The transcript assemblies for P. taeda were validated through sequencing of PCR products amplified using gene-specific primers based on the putative PAL gene assemblies. Verified PAL gene sequences were aligned and a gene tree was estimated. The resulting gene tree was reconciled with a known species tree and the time points for gene duplication events were inferred relative to the divergence of major plant lineages. RESULTS: In contrast to angiosperms, gymnosperms have retained a diverse set of PAL genes distributed among three major clades that arose from gene duplication events predating the divergence of these two seed plant lineages. Whereas multiple PAL genes have been identified in sequenced angiosperm genomes, all characterized angiosperm PAL genes form a single clade in the gene PAL tree, suggesting they are derived from a single gene in an ancestral angiosperm genome. The five distinct PAL genes detected and verified in P. taeda were derived from a combination of duplication events predating and postdating the divergence of angiosperms and gymnosperms. CONCLUSIONS: Gymnosperms have a more phylogenetically diverse set of PAL genes than angiosperms. This inference has contrasting implications for the evolution of PAL gene function in gymnosperms and angiosperms.

Microarray analysis and scale-free gene networks identify candidate regulators in drought-stressed roots of loblolly pine (P. taeda L.).

BACKGROUND: Global transcriptional analysis of loblolly pine (Pinus taeda L.) is challenging due to limited molecular tools. PtGen2, a 26,496 feature cDNA microarray, was fabricated and used to assess drought-induced gene expression in loblolly pine propagule roots. Statistical analysis of differential expression and weighted gene correlation network analysis were used to identify drought-responsive genes and further characterize the molecular basis of drought tolerance in loblolly pine. RESULTS: Microarrays were used to interrogate root cDNA populations obtained from 12 genotype × treatment combinations (four genotypes, three watering regimes). Comparison of drought-stressed roots with roots from the control treatment identified 2445 genes displaying at least a 1.5-fold expression difference (false discovery rate = 0.01). Genes commonly associated with drought response in pine and other plant species, as well as a number of abiotic and biotic stress-related genes, were up-regulated in drought-stressed roots. Only 76 genes were identified as differentially expressed in drought-recovered roots, indicating that the transcript population can return to the pre-drought state within 48 hours. Gene correlation analysis predicts a scale-free network topology and identifies eleven co-expression modules that ranged in size from 34 to 938 members. Network topological parameters identified a number of central nodes (hubs) including those with significant homology (E-values ≤ 2 × 10-30) to 9-cis-epoxycarotenoid dioxygenase, zeatin O-glucosyltransferase, and ABA-responsive protein. Identified hubs also include genes that have been associated previously with osmotic stress, phytohormones, enzymes that detoxify reactive oxygen species, and several genes of unknown function. CONCLUSION: PtGen2 was used to evaluate transcriptome responses in loblolly pine and was leveraged to identify 2445 differentially expressed genes responding to severe drought stress in roots. Many of the genes identified are known to be up-regulated in response to osmotic stress in pine and other plant species and encode proteins involved in both signal transduction and stress tolerance. Gene expression levels returned to control values within a 48-hour recovery period in all but 76 transcripts. Correlation network analysis indicates a scale-free network topology for the pine root transcriptome and identifies central nodes that may serve as drivers of drought-responsive transcriptome dynamics in the roots of loblolly pine.

Corticosterone and testosterone treatment influence expression of gene pathways linked to meiotic segregation in preovulatory follicles of the domestic hen.

Decades of work indicate that female birds can control their offspring sex ratios in response to environmental and social cues. In laying hens, hormones administered immediately prior to sex chromosome segregation can exert sex ratio skews, indicating that these hormones may act directly on the germinal disc to influence which sex chromosome is retained in the oocyte and which is discarded into an unfertilizable polar body. We aimed to uncover the gene pathways involved in this process by testing whether treatments with testosterone or corticosterone that were previously shown to influence sex ratios elicit changes in the expression of genes and/or gene pathways involved in the process of meiotic segregation. We injected laying hens with testosterone, corticosterone, or control oil 5h prior to ovulation and collected germinal discs from the F1 preovulatory follicle in each hen 1.5h after injection. We used RNA-sequencing (RNA-seq) followed by DESeq2 and gene set enrichment analyses to identify genes and gene pathways that were differentially expressed between germinal discs of control and hormone-treated hens. Corticosterone treatment triggered downregulation of 13 individual genes, as well as enrichment of gene sets related to meiotic spindle organization and chromosome segregation, and additional gene sets that function in ion transport. Testosterone treatment triggered upregulation of one gene, and enrichment of one gene set that functions in nuclear chromosome segregation. This work indicates that corticosterone can be a potent regulator of meiotic processes and provides potential gene targets on which corticosterone and/or testosterone may act to influence offspring sex ratios in birds.

Effects of diethylcarbamazine and ivermectin treatment on Brugia malayi gene expression in infected gerbils (Meriones unguiculatus).

Lymphatic filariasis (LF) threatens nearly 20% of the world's population and has handicapped one-third of the 120 million people currently infected. Current control and elimination programs for LF rely on mass drug administration of albendazole plus diethylcarbamazine (DEC) or ivermectin. Only the mechanism of action of albendazole is well understood. To gain a better insight into antifilarial drug action in vivo, we treated gerbils harbouring patent Brugia malayi infections with 6 mg kg-1 DEC, 0.15 mg kg-1 ivermectin or 1 mg kg-1 albendazole. Treatments had no effect on the numbers of worms present in the peritoneal cavity of treated animals, so effects on gene expression were a direct result of the drug and not complicated by dying parasites. Adults and microfilariae were collected 1 and 7 days post-treatment and RNA isolated for transcriptomic analysis. The experiment was repeated three times. Ivermectin treatment produced the most differentially expressed genes (DEGs), 113. DEC treatment yielded 61 DEGs. Albendazole treatment resulted in little change in gene expression, with only 6 genes affected. In total, nearly 200 DEGs were identified with little overlap between treatment groups, suggesting that these drugs may interfere in different ways with processes important for parasite survival, development, and reproduction.

Variation among human, veterinary and environmental Mycobacterium chelonae-abscessus complex isolates observed using core genome phylogenomic analysis, targeted gene comparison, and anti-microbial susceptibility patterns.

Mycobacterium chelonae is a member of the Mycobacterium chelonae-abscessus complex and a cause of opportunistic disease in fish, reptiles, birds, and mammals including humans. Isolates in the complex are often difficult to identify and have differing antimicrobial susceptibilities. Thirty-one previously identified rapidly-growing, non-tuberculous Mycobacterium sp. isolates cultured from biofilms, fish, reptiles, mammals, including humans, and three ATCC reference strains were evaluated with nine M. chelonae-abscessus complex whole genome sequences from GenBank by phylogenomic analysis, targeted gene comparisons, and in-vitro antimicrobial susceptibility patterns to assess strain variation among isolates from different sources. Results revealed minimal genetic variation among the M. chelonae strains. However, the core genomic alignment and SNP pattern of the complete 16S rRNA sequence clearly separated the turtle type strain ATCC 35752T from the clinical isolates and human reference strain "M. chelonae chemovar niacinogenes" ATCC 19237, providing evidence of two distinct subspecies. Concatenation of the partial rpoB (752 bp) and complete hsp65 (1,626 bp) sequence produced the same species/subspecies delineations as the core phylogeny. Partial rpoB and hsp65 sequences identified all the clinical isolates to the appropriate species level when respective cut-offs of 98% and 98.4% identity to the M. chelonae type strain ATCC 35752T were employed. The human strain, ATCC19237, was the most representative strain for the evaluated human, veterinary, and environmental strains. Additionally, two isolates were identified as Mycobacterium saopaulense, its first identification in a non-fish or non-human host.

Recognition and killing of Brugia malayi microfilariae by human immune cells is dependent on the parasite sample and is not altered by ivermectin treatment.

Mass administration of macrocyclic lactones targets the transmission of the causative agents of lymphatic filariasis to their insect vectors by rapidly clearing microfilariae (Mf) from the circulation. It has been proposed that the anti-filarial action of these drugs may be mediated through the host immune system. We recently developed an in vitro assay for monitoring the attachment to and killing of B. malayi Mf by human neutrophils (PMNs) and monocytes (PBMCs), however, the levels of both cell to worm attachment and leukocyte mediated Mf killing varied greatly between individual experiments. To determine whether differences in an individual's immune cells or the Mf themselves might account for the variability in survival, PMNs and PBMCs were isolated from 12 donors every week for 4 weeks and the cells used for survival assays with a different batch of Mf, thereby keeping donors constant but varying the Mf sample. Results from these experiments indicate that, overall, killing is Mf-rather than donor-dependent. To assess whether ivermectin (IVM) or diethylcarbamazine (DEC) increase killing, Mf were incubated either alone or with immune cells in the presence of IVM or DEC. Neither drug induced a significant difference in the survival of Mf whether cultured with or without cells, with the exception of DEC at 2 h post incubation. In addition, human PBMCs and PMNs were incubated with IVM or DEC for 1 h or 16 h prior to RNA extraction and Illumina sequencing. Although donor-to-donor variation may mask subtle differences in gene expression, principle component analysis of the RNASeq data indicates that there is no significant change in the expression of any genes from the treated cells versus controls. Together these data suggest that IVM and DEC have little direct effect on immune cells involved in the rapid clearance of Mf from the circulation.

Water stress-responsive genes in loblolly pine (Pinus taeda) roots identified by analyses of expressed sequence tag libraries.

Drought stress is the principal cause of seedling mortality in pine forests of the southeastern United States and in many other forested regions around the globe. As part of a larger effort to discover loblolly pine genes, this study subjected rooted cuttings of three unrelated pine genotypes to three watering regimens. Expressed sequence tags (ESTs) were obtained from both the 3' and 5' ends of 12,918 randomly selected cDNAs generated from root tissues. These ESTs were clustered to identify 6,765 unique transcripts (UniScripts) derived from 6,202 putative unique genes (UniGenes-S). Tentative annotations were assigned on the basis of BLASTX comparisons to the Protein Information Resource Nonredundant Reference (PIR-NREF) database. Expression levels of 42 UniScripts varied with high statistical significance with respect to treatment. Many of them resembled gene products shown to be important for drought tolerance in other species, including dehydrins, endochitinases, cytochrome P450 enzymes, pathogenesis-related proteins and various late-embryogenesis abundant (LEA) gene products. Similarly, expression levels of 110 UniScripts varied with high statistical significance among genotypes, indicating that gene expression patterns in this species are much more dependent on genotype than on treatment. Most of the water stress-induced pine UniScripts that appeared to encode products resembling drought tolerance factors in other species were most highly induced in a single genotype, suggesting that particularly useful adaptive alleles for drought tolerance might exist within the collection of cDNAs characterized from this genotype. Mining and visualizing the complete data set, as well as downloading of both EST and UniScript contig sequences, are possible using MAGIC Gene Discovery at http://fungen.org/genediscovery/.

SAGE profiling and demonstration of differential gene expression along the axial developmental gradient of lignifying xylem in loblolly pine (Pinus taeda).

Wood formation has been studied extensively at the cellular and biochemical levels, but remains poorly understood with respect to gene expression and regulation. As a first step toward identifying genes specifically involved in wood formation and characterizing their roles in determining wood quality, serial analysis of gene expression (SAGE) was used to quantify gene expression in lignifying xylem from a single, 10-year-old loblolly pine (Pinus taeda L.). Two SAGE libraries were generated based on lignifying xylem isolated from either the upper (crown) or lower (base) portions of the trunk. Over 85,000 tags representing a maximum of 27,398 expressed genes were analyzed from the crown wood library, and more than 65,000 tags, representing a maximum of 25,983 expressed genes, were analyzed in the base wood library. Combining these data sets to reflect the sum of genes expressed in lignifying xylem, 150,855 tags were cataloged, representing a maximum of 42,641 different genes. Currently, this study represents the most extensive analysis of its kind in a higher plant and provides a quantitative description of the transcriptome representing the lignifying xylem of a 10-year-old loblolly pine.

Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies.

BACKGROUND: The size and complexity of conifer genomes has, until now, prevented full genome sequencing and assembly. The large research community and economic importance of loblolly pine, Pinus taeda L., made it an early candidate for reference sequence determination. RESULTS: We develop a novel strategy to sequence the genome of loblolly pine that combines unique aspects of pine reproductive biology and genome assembly methodology. We use a whole genome shotgun approach relying primarily on next generation sequence generated from a single haploid seed megagametophyte from a loblolly pine tree, 20-1010, that has been used in industrial forest tree breeding. The resulting sequence and assembly was used to generate a draft genome spanning 23.2 Gbp and containing 20.1 Gbp with an N50 scaffold size of 66.9 kbp, making it a significant improvement over available conifer genomes. The long scaffold lengths allow the annotation of 50,172 gene models with intron lengths averaging over 2.7 kbp and sometimes exceeding 100 kbp in length. Analysis of orthologous gene sets identifies gene families that may be unique to conifers. We further characterize and expand the existing repeat library based on the de novo analysis of the repetitive content, estimated to encompass 82% of the genome. CONCLUSIONS: In addition to its value as a resource for researchers and breeders, the loblolly pine genome sequence and assembly reported here demonstrates a novel approach to sequencing the large and complex genomes of this important group of plants that can now be widely applied.

Biomarker genes highlight intraspecific and interspecific variations in the responses of Pinus taeda L. and Pinus radiata D. Don to Sirex noctilio F. acid gland secretions.

Sirex noctilio F., a Eurasian horntail woodwasp recently introduced into North America, oviposits in pines and other conifers and in the process spreads a phytopathogenic fungus that serves as a food source for its larvae. During oviposition the woodwasp also deposits mucus produced in its acid (venom) gland that alters pine defense responses and facilitates infection by the fungus. A 26,496-feature loblolly pine cDNA microarray was used to survey gene expression of pine tissue responding to S. noctilio venom. Six genes were selected for further assessment by quantitative real-time polymerase chain reaction (qRT-PCR), including one that encoded an apparent PR-4 protein and another that encoded a thaumatin-like protein. Expression of both was strongly induced in response to venom, while expression of an apparent actin gene (ACT1) was stable in response to the venom. The pattern of gene response was similar in Pinus taeda L. and Pinus radiata D. Don, but the magnitude of response in P. radiata was significantly stronger for each of the induced genes. The magnitude of the biomarker gene response to venom also varied according to genotype within these two species. The qRT-PCR assay was used to demonstrate that the primary bioactive component in S. noctilio venom is a polypeptide.

An improved method of RNA isolation from loblolly pine (P. taeda L.) and other conifer species.

Tissues isolated from conifer species, particularly those belonging to the Pinaceae family, such as loblolly pine (Pinus taeda L.), contain high concentrations of phenolic compounds and polysaccharides that interfere with RNA purification. Isolation of high-quality RNA from these species requires rigorous tissue collection procedures in the field and the employment of an RNA isolation protocol comprised of multiple organic extraction steps in order to isolate RNA of sufficient quality for microarray and other genomic analyses. The isolation of high-quality RNA from field-collected loblolly pine samples can be challenging, but several modifications to standard tissue and RNA isolation procedures greatly improve results. The extent of general RNA degradation increases if samples are not properly collected and transported from the field, especially during large-scale harvests. Total RNA yields can be increased significantly by pulverizing samples in a liquid nitrogen freezer mill prior to RNA isolation, especially when samples come from woody tissues. This is primarily due to the presence of oxidizing agents, such as phenolic compounds, and polysaccharides that are both present at high levels in extracts from the woody tissues of most conifer species. If not removed, these contaminants can carry over leading to problems, such as RNA degradation, that result in low yields and a poor quality RNA sample. Carryover of phenolic compounds, as well as polysaccharides, can also reduce or even completely eliminate the activity of reverse transcriptase or other polymerases commonly used for cDNA synthesis. In particular, RNA destined to be used as template for double-stranded cDNA synthesis in the generation of cDNA libraries, single-stranded cDNA synthesis for PCR or qPCR's, or for the synthesis of microarray target materials must be of the highest quality if researchers expect to obtain optimal results. RNA isolation techniques commonly employed for many other plant species are often insufficient in their ability to remove these contaminants from conifer samples and thus do not yield total RNA samples suitable for downstream manipulations. In this video we demonstrate methods for field collection of conifer tissues, beginning with the felling of a forty year-old tree, to the harvesting of phloem, secondary xylem, and reaction wood xylem. We also demonstrate an RNA isolation protocol that has consistently yielded high-quality RNA for subsequent enzymatic manipulations.

Processing the loblolly pine PtGen2 cDNA microarray.

PtGen2 is a 26,496 feature cDNA microarray containing amplified loblolly pine ESTs. The array is produced in our laboratory for use by researchers studying gene expression in pine and other conifer species. PtGen2 was developed as a result of our gene discovery efforts in loblolly pine, and is comprised of sequences identified primarily from root tissues, but also from needle and stem. PtGen2 has been tested by hybridizing different Cy-dye labeled conifer target cDNAs, using both amplified and non-amplified indirect labeling methods, and also tested with a number of hybridization and washing conditions. This video focuses on the handling and processing of slides before and after pre-hybridization, as well as after hybridization, using some modifications to procedures developed previously. Also included, in text form only, are the protocols used for the generation, labeling and clean up of target cDNA s, as well as information on software used for downstream data processing. PtGen2 is printed with a proprietary print buffer that contains high concentrations of salt that can be difficult to remove completely. The slides are washed first in a warm SDS solution prior to pre-hybridization. After pre-hybridization, the slides are washed vigorously in several changes of water to complete removal of remaining salts. LifterSlips are then cleaned and positioned on the slides and labeled cDNA is carefully loaded onto the microarray by way of capillary action which provides for even distribution of the sample across the slide, and reduces the chance of bubble incorporation. Hybridization of targets to the array is done at 48 degrees C in high humidity conditions. After hybridization, a series of standard washes are done at 53 degrees C and room temperature for extended times. Processing PtGen2 slides using this technique reduces salt and SDS-derived artifacts often seen when the array is processed less rigorously. Hybridizing targets derived from several different conifer RNA sources, this processing protocol yielded fewer artifacts, reduced background, and provided better consistency among different experimental groups of arrays.