The effect of a startle-eliciting device on the foraging success of individual harbor seals (Phoca vitulina).

Pinniped predation on commercially and ecologically important prey has been a source of conflict for centuries. In the Salish Sea, harbor seals (Phoca vitulina) are suspected of impeding the recovery of culturally and ecologically critical Pacific salmon (Oncorhynchus spp.). In Fall 2020, a novel deterrent called Targeted Acoustic Startle Technology (TAST) was deployed at Whatcom Creek to deter harbor seals from preying on fall runs of hatchery chum (O. keta) and Chinook (O. tshawytscha) salmon in Bellingham, Washington, USA. Field observations were conducted in 2020 to compare the presence and foraging success of individual harbor seals across sound exposure (TAST-on) and control (TAST-off) conditions. Observations conducted the previous (2019) and following (2021) years were used to compare the effects observed in 2020 to two control years. Using photo-identification, individual seals were associated with foraging successes across all 3 years of the study. Generalized linear mixed models showed a significant 45.6% reduction in the duration (min) individuals remained at the creek with TAST on, and a significant 43.8% reduction in the overall foraging success of individuals. However, the observed effect of TAST varied across individual seals. Seals that were observed regularly within one season were more likely to return the year after, regardless of TAST treatment. Generalized linear models showed interannual variation in the number of seals present and salmon consumed. However, the effect of TAST in 2020 was greater than the observed variation across years. Our analyses suggest TAST can be an effective tool for managing pinniped predation, although alternate strategies such as deploying TAST longer-term and using multi-unit setups to increase coverage could help strengthen its effects. Future studies should further examine the individual variability found in this study.

Using Drosophila to identify naturally occurring genetic modifiers of amyloid beta 42- and tau-induced toxicity.

Alzheimer's disease is characterized by 2 pathological proteins, amyloid beta 42 and tau. The majority of Alzheimer's disease cases in the population are sporadic and late-onset Alzheimer's disease, which exhibits high levels of heritability. While several genetic risk factors for late-onset Alzheimer's disease have been identified and replicated in independent studies, including the ApoE ε4 allele, the great majority of the heritability of late-onset Alzheimer's disease remains unexplained, likely due to the aggregate effects of a very large number of genes with small effect size, as well as to biases in sample collection and statistical approaches. Here, we present an unbiased forward genetic screen in Drosophila looking for naturally occurring modifiers of amyloid beta 42- and tau-induced ommatidial degeneration. Our results identify 14 significant SNPs, which map to 12 potential genes in 8 unique genomic regions. Our hits that are significant after genome-wide correction identify genes involved in neuronal development, signal transduction, and organismal development. Looking more broadly at suggestive hits (P < 10-5), we see significant enrichment in genes associated with neurogenesis, development, and growth as well as significant enrichment in genes whose orthologs have been identified as significantly or suggestively associated with Alzheimer's disease in human GWAS studies. These latter genes include ones whose orthologs are in close proximity to regions in the human genome that are associated with Alzheimer's disease, but where a causal gene has not been identified. Together, our results illustrate the potential for complementary and convergent evidence provided through multitrait GWAS in Drosophila to supplement and inform human studies, helping to identify the remaining heritability and novel modifiers of complex diseases.

Disruption of a ∼23-24 nucleotide small RNA pathway elevates DNA damage responses in Tetrahymena thermophila.

Endogenous RNA interference (RNAi) pathways regulate a wide range of cellular processes in diverse eukaryotes, yet in the ciliated eukaryote, Tetrahymena thermophila, the cellular purpose of RNAi pathways that generate ∼23-24 nucleotide (nt) small (s)RNAs has remained unknown. Here, we investigated the phenotypic and gene expression impacts on vegetatively growing cells when genes involved in ∼23-24 nt sRNA biogenesis are disrupted. We observed slower proliferation and increased expression of genes involved in DNA metabolism and chromosome organization and maintenance in sRNA biogenesis mutants RSP1Δ, RDN2Δ, and RDF2Δ. In addition, RSP1Δ and RDN2Δ cells frequently exhibited enlarged chromatin extrusion bodies, which are nonnuclear, DNA-containing structures that may be akin to mammalian micronuclei. Expression of homologous recombination factor Rad51 was specifically elevated in RSP1Δ and RDN2Δ strains, with Rad51 and double-stranded DNA break marker γ-H2A.X localized to discrete macronuclear foci. In addition, an increase in Rad51 and γ-H2A.X foci was also found in knockouts of TWI8, a macronucleus-localized PIWI protein. Together, our findings suggest that an evolutionarily conserved role for RNAi pathways in maintaining genome integrity may be extended even to the early branching eukaryotic lineage that gave rise to Tetrahymena thermophila.

Evolutionary network genomics of wood formation in a phylogenetic survey of angiosperm forest trees.

Wood formation was present in early angiosperms, but has been highly modified through evolution to generate the anatomical diversity seen in extant angiosperm lineages. In this project, we modeled changes in gene coexpression relationships associated with the evolution of wood formation in a phylogenetic survey of 13 angiosperm tree species. Gravitropic stimulation was used as an experimental treatment to alter wood formation and also perturb gene expression. Gene transcript abundances were determined using RNA sequencing of developing wood tissues from upright trees, and from the top (tension wood) and bottom (opposite wood) tissues of gravistimulated trees. A network-based approach was employed to align gene coexpression networks across species based on orthologous relationships. A large-scale, multilayer network was modeled that identified both lineage-specific gene coexpression modules and modules conserved across multiple species. Functional annotation and analysis of modules identified specific regulatory processes associated with conserved modules, including regulation of hormones, protein phosphorylation, meristem development and epigenetic processes. Our results provide novel insights into the evolution and development of wood formation, and demonstrate the ability to identify biological processes and genes important for the evolution of a foundational trait in nonmodel, undomesticated forest trees.

A comprehensive genomic scan reveals gene dosage balance impacts on quantitative traits in Populus trees.

Gene dosage variation and the associated changes in gene expression influence a wide variety of traits, ranging from cancer in humans to yield in plants. It is also expected to affect important traits of ecological and agronomic importance in forest trees, but this variation has not been systematically characterized or exploited. Here we performed a comprehensive scan of the Populus genome for dosage-sensitive loci affecting quantitative trait variation for spring and fall phenology and biomass production. The study population was a large collection of clonally propagated F1 hybrid lines of Populus that saturate the genome 10-fold with deletions and insertions (indels) of known sizes and positions. As a group, the phenotypic means of the indel lines consistently differed from control nonindel lines, with an overall negative effect of both insertions and deletions on all biomass-related traits but more diverse effects and an overall wider phenotypic distribution of the indel lines for the phenology-related traits. We also investigated the correlation between gene dosage at specific chromosomal locations and phenotype, to identify dosage quantitative trait loci (dQTL). Such dQTL were detected for most phenotypes examined, but stronger effect dQTL were identified for the phenology-related traits than for the biomass traits. Our genome-wide screen for dosage sensitivity in a higher eukaryote demonstrates the importance of global genomic balance and the impact of dosage on life history traits.

Heritable Phytohormone Profiles of Poplar Genotypes Vary in Resistance to a Galling Aphid.

Insect galls are highly specialized structures arising from atypical development of plant tissue induced by insects. Galls provide the insect enhanced nutrition and protection against natural enemies and environmental stresses. Galls are essentially plant organs formed by an intimate biochemical interaction between the gall-inducing insect and its host plant. Because galls are plant organs, their development is likely to be governed by phytohormones involved in normal organogenesis. We characterized concentrations of both growth and defensive phytohormones in ungalled control leaves and galls induced by the aphid Pemphigus betae on narrowleaf cottonwood Populus angustifolia that differ genotypically in resistance to this insect. We found that susceptible trees differed from resistant trees in constitutive concentrations of both growth and defense phytohormones. Susceptible trees were characterized by significantly higher constitutive cytokinin concentrations in leaves, significantly greater ability of aphids to elicit cytokinin increases, and significantly lower constitutive defense phytohormone concentrations than observed in resistant trees. Phytohormone concentrations in both constitutive and induced responses in galled leaves exhibited high broad-sense heritability that, respectively, ranged from 0.39 to 0.93 and from 0.28 to 0.66, suggesting that selection can act upon these traits and that they might vary across the landscape. Increased cytokinin concentrations may facilitate forming strong photosynthate sinks in the galls, a requirement for galling insect success. By characterizing for the first time the changes in 15 phytohormones belonging to five different classes, this study offers a better overview of the signaling alteration occurring in galls that has likely been important for their ecology and evolution. Copyright © 2019 The Author(s). This is an open-access article distributed under the CC BY-NC-ND 4.0 International license .

Large effect quantitative trait loci for salicinoid phenolic glycosides in Populus: Implications for gene discovery.

Genomic studies have been used to identify genes underlying many important plant secondary metabolic pathways. However, genes for salicinoid phenolic glycosides (SPGs)-ecologically important compounds with significant commercial, cultural, and medicinal applications-remain largely undescribed. We used a linkage map derived from a full-sib population of hybrid cottonwoods (Populus spp.) to search for quantitative trait loci (QTL) for the SPGs salicortin and HCH-salicortin. SSR markers and primer sequences were used to anchor the map to the V3.0 P. trichocarpa genome. We discovered 21 QTL for the two traits, including a major QTL for HCH-salicortin (R2 = .52) that colocated with a QTL for salicortin on chr12. Using the V3.0 Populus genome sequence, we identified 2,983 annotated genes and 1,480 genes of unknown function within our QTL intervals. We note ten candidate genes of interest, including a BAHD-type acyltransferase that has been potentially linked to Populus SPGs. Our results complement other recent studies in Populus with implications for gene discovery and the evolution of defensive chemistry in a model genus. To our knowledge, this is the first study to use a full-sib mapping population to identify QTL intervals and gene lists associated with SPGs.

Transcriptional and temporal response of Populus stems to gravi-stimulation.

Plants modify development in response to external stimuli, to produce new growth that is appropriate for environmental conditions. For example, gravi-stimulation of leaning branches in angiosperm trees results in modifications of wood development, to produce tension wood that pulls leaning stems upright. Here, we use gravi-stimulation and tension wood response to dissect the temporal changes in gene expression underlying wood formation in Populus stems. Using time-series analysis of seven time points over a 14-d experiment, we identified 8,919 genes that were differentially expressed between tension wood (upper) and opposite wood (lower) sides of leaning stems. Clustering of differentially expressed genes showed four major transcriptional responses, including gene clusters whose transcript levels were associated with two types of tissue-specific impulse responses that peaked at about 24-48 h, and gene clusters with sustained changes in transcript levels that persisted until the end of the 14-d experiment. Functional enrichment analysis of those clusters suggests they reflect temporal changes in pathways associated with hormone regulation, protein localization, cell wall biosynthesis and epigenetic processes. Time-series analysis of gene expression is an underutilized approach for dissecting complex developmental responses in plants, and can reveal gene clusters and mechanisms influencing development.

Identifying gene coexpression networks underlying the dynamic regulation of wood-forming tissues in Populus under diverse environmental conditions.

Trees modify wood formation through integration of environmental and developmental signals in complex but poorly defined transcriptional networks, allowing trees to produce woody tissues appropriate to diverse environmental conditions. In order to identify relationships among genes expressed during wood formation, we integrated data from new and publically available datasets in Populus. These datasets were generated from woody tissue and include transcriptome profiling, transcription factor binding, DNA accessibility and genome-wide association mapping experiments. Coexpression modules were calculated, each of which contains genes showing similar expression patterns across experimental conditions, genotypes and treatments. Conserved gene coexpression modules (four modules totaling 8398 genes) were identified that were highly preserved across diverse environmental conditions and genetic backgrounds. Functional annotations as well as correlations with specific experimental treatments associated individual conserved modules with distinct biological processes underlying wood formation, such as cell-wall biosynthesis, meristem development and epigenetic pathways. Module genes were also enriched for DNase I hypersensitivity footprints and binding from four transcription factors associated with wood formation. The conserved modules are excellent candidates for modeling core developmental pathways common to wood formation in diverse environments and genotypes, and serve as testbeds for hypothesis generation and testing for future studies.

Transcript profiling of a novel plant meristem, the monocot cambium.

While monocots lack the ability to produce a vascular cambium or woody growth, some monocot lineages evolved a novel lateral meristem, the monocot cambium, which supports secondary radial growth of stems. In contrast to the vascular cambium found in woody angiosperm and gymnosperm species, the monocot cambium produces secondary vascular bundles, which have an amphivasal organization of tracheids encircling a central strand of phloem. Currently there is no information concerning the molecular genetic basis of the development or evolution of the monocot cambium. Here we report high-quality transcriptomes for monocot cambium and early derivative tissues in two monocot genera, Yucca and Cordyline. Monocot cambium transcript profiles were compared to those of vascular cambia and secondary xylem tissues of two forest tree species, Populus trichocarpa and Eucalyptus grandis. Monocot cambium transcript levels showed that there are extensive overlaps between the regulation of monocot cambia and vascular cambia. Candidate regulatory genes that vary between the monocot and vascular cambia were also identified, and included members of the KANADI and CLE families involved in polarity and cell-cell signaling, respectively. We suggest that the monocot cambium may have evolved in part through reactivation of genetic mechanisms involved in vascular cambium regulation.

Genetic variation in NIN1 and C/VIF1 genes is significantly associated with Populus angustifolia resistance to a galling herbivore, Pemphigus betae.

The identification of genes associated with ecologically important traits provides information on the potential genetic mechanisms underlying the responses of an organism to its natural environment. In this study, we investigated the genetic basis of host plant resistance to the gall-inducing aphid, Pemphigus betae, in a natural population of 154 narrowleaf cottonwoods (Populus angustifolia). We surveyed genetic variation in two genes putatively involved in sink-source relations and a phenology gene that co-located in a previously identified quantitative trait locus for resistance to galling. Using a candidate gene approach, three major findings emerged. First, natural variation in tree resistance to galling was repeatable. Sampling of the same tree genotypes 20 years after the initial survey in 1986 show that 80% of the variation in resistance was due to genetic differences among individuals. Second, we identified significant associations at the single nucleotide polymorphism and haplotype levels between the plant neutral invertase gene NIN1 and tree resistance. Invertases are a class of sucrose hydrolyzing enzymes and play an important role in plant responses to biotic stress, including the establishment of nutrient sinks. These associations with NIN1 were driven by a single nucleotide polymorphism (NIN1_664) located in the second intron of the gene and in an orthologous sequence to two known regulatory elements. Third, haplotypes from an inhibitor of invertase (C/VIF1) were significantly associated with tree resistance. The identification of genetic variation in these two genes provides a starting point to understand the possible genetic mechanisms that contribute to tree resistance to gall formation. We also build on previous work demonstrating that genetic differences in sink-source relationships of the host influence the ability of P. betae to manipulate the flow of nutrients and induce a nutrient sink.

Creation and Genomic Analysis of Irradiation Hybrids in Populus.

Establishing efficient functional genomic systems for creating and characterizing genetic variation in forest trees is challenging. Here we describe protocols for creating novel gene-dosage variation in Populus through gamma-irradiation of pollen, followed by genomic analysis to identify chromosomal regions that have been deleted or inserted in each F1 individual. Irradiated pollen is used in a controlled, interspecific cross to create F1 progeny that carry deletions and insertions of chromosomal regions. These insertions and deletions result in novel changes in gene dosage that in turn affect both qualitative and quantitative phenotypic traits. The protocols described here outline the processes involved in optimizing irradiation levels and performing controlled crosses, sowing seed and propagating seedlings, and genome-wide sequence-based analysis of deletions and insertions in the F1 progeny. The same approach could be applied to other vegetatively propagated species. © 2016 by John Wiley & Sons, Inc.

Transcriptional and Hormonal Regulation of Gravitropism of Woody Stems in Populus.

Angiosperm trees reorient their woody stems by asymmetrically producing a specialized xylem tissue, tension wood, which exerts a strong contractile force resulting in negative gravitropism of the stem. Here, we show, in Populus trees, that initial gravity perception and response occurs in specialized cells through sedimentation of starch-filled amyloplasts and relocalization of the auxin transport protein, PIN3. Gibberellic acid treatment stimulates the rate of tension wood formation and gravibending and enhances tissue-specific expression of an auxin-responsive reporter. Gravibending, maturation of contractile fibers, and gibberellic acid (GA) stimulation of tension wood formation are all sensitive to transcript levels of the Class I KNOX homeodomain transcription factor-encoding gene ARBORKNOX2 (ARK2). We generated genome-wide transcriptomes for trees in which gene expression was perturbed by gravistimulation, GA treatment, and modulation of ARK2 expression. These data were employed in computational analyses to model the transcriptional networks underlying wood formation, including identification and dissection of gene coexpression modules associated with wood phenotypes, GA response, and ARK2 binding to genes within modules. We propose a model for gravitropism in the woody stem in which the peripheral location of PIN3-expressing cells relative to the cambium results in auxin transport toward the cambium in the top of the stem, triggering tension wood formation, while transport away from the cambium in the bottom of the stem triggers opposite wood formation.

A System for Dosage-Based Functional Genomics in Poplar.

Altering gene dosage through variation in gene copy number is a powerful approach to addressing questions regarding gene regulation, quantitative trait loci, and heterosis, but one that is not easily applied to sexually transmitted species. Elite poplar (Populus spp) varieties are created through interspecific hybridization, followed by clonal propagation. Altered gene dosage relationships are believed to contribute to hybrid performance. Clonal propagation allows for replication and maintenance of meiotically unstable ploidy or structural variants and provides an alternative approach to investigating gene dosage effects not possible in sexually propagated species. Here, we built a genome-wide structural variation system for dosage-based functional genomics and breeding of poplar. We pollinated Populus deltoides with gamma-irradiated Populus nigra pollen to produce >500 F1 seedlings containing dosage lesions in the form of deletions and insertions of chromosomal segments (indel mutations). Using high-precision dosage analysis, we detected indel mutations in ∼55% of the progeny. These indels varied in length, position, and number per individual, cumulatively tiling >99% of the genome, with an average of 10 indels per gene. Combined with future phenotype and transcriptome data, this population will provide an excellent resource for creating and characterizing dosage-based variation in poplar, including the contribution of dosage to quantitative traits and heterosis.

A resource for characterizing genome-wide binding and putative target genes of transcription factors expressed during secondary growth and wood formation in Populus.

Identifying transcription factor target genes is essential for modeling the transcriptional networks underlying developmental processes. Here we report a chromatin immunoprecipitation sequencing (ChIP-seq) resource consisting of genome-wide binding regions and associated putative target genes for four Populus homeodomain transcription factors expressed during secondary growth and wood formation. Software code (programs and scripts) for processing the Populus ChIP-seq data are provided within a publically available iPlant image, including tools for ChIP-seq data quality control and evaluation adapted from the human Encyclopedia of DNA Elements (ENCODE) project. Basic information for each transcription factor (including members of Class I KNOX, Class III HD ZIP, BEL1-like families) binding are summarized, including the number and location of binding regions, distribution of binding regions relative to gene features, associated putative target genes, and enriched functional categories of putative target genes. These ChIP-seq data have been integrated within the Populus Genome Integrative Explorer (PopGenIE) where they can be analyzed using a variety of web-based tools. We present an example analysis that shows preferential binding of transcription factor ARBORKNOX1 to the nearest neighbor genes in a pre-calculated co-expression network module, and enrichment for meristem-related genes within this module including multiple orthologs of Arabidopsis KNOTTED-like Arabidopsis 2/6.

The Populus ARBORKNOX1 homeodomain transcription factor regulates woody growth through binding to evolutionarily conserved target genes of diverse function.

The class I KNOX homeodomain transcription factor ARBORKNOX1 (ARK1) is a key regulator of vascular cambium maintenance and cell differentiation in Populus. Currently, basic information is lacking concerning the distribution, functional characteristics, and evolution of ARK1 binding in the Populus genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) technology to identify ARK1 binding loci genome-wide in Populus. Computational analyses evaluated the distribution of ARK1 binding loci, the function of genes associated with bound loci, the effect of ARK1 binding on transcript levels, and evolutionary conservation of ARK1 binding loci. ARK1 binds to thousands of loci which are highly enriched proximal to the transcriptional start sites of genes of diverse functions. ARK1 target genes are significantly enriched in paralogs derived from the whole-genome salicoid duplication event. Both ARK1 and a maize (Zea mays) homolog, KNOTTED1, preferentially target evolutionarily conserved genes. However, only a small portion of ARK1 target genes are significantly differentially expressed in an ARK1 over-expression mutant. This study describes the functional characteristics and evolution of DNA binding by a transcription factor in an undomesticated tree, revealing complexities similar to those shown for transcription factors in model animal species.

Evaluation of experimental design and computational parameter choices affecting analyses of ChIP-seq and RNA-seq data in undomesticated poplar trees.

BACKGROUND: One of the great advantages of next generation sequencing is the ability to generate large genomic datasets for virtually all species, including non-model organisms. It should be possible, in turn, to apply advanced computational approaches to these datasets to develop models of biological processes. In a practical sense, working with non-model organisms presents unique challenges. In this paper we discuss some of these challenges for ChIP-seq and RNA-seq experiments using the undomesticated tree species of the genus Populus. RESULTS: We describe specific challenges associated with experimental design in Populus, including selection of optimal genotypes for different technical approaches and development of antibodies against Populus transcription factors. Execution of the experimental design included the generation and analysis of Chromatin immunoprecipitation-sequencing (ChIP-seq) data for RNA polymerase II and transcription factors involved in wood formation. We discuss criteria for analyzing the resulting datasets, determination of appropriate control sequencing libraries, evaluation of sequencing coverage needs, and optimization of parameters. We also describe the evaluation of ChIP-seq data from Populus, and discuss the comparison between ChIP-seq and RNA-seq data and biological interpretations of these comparisons. CONCLUSIONS: These and other "lessons learned" highlight the challenges but also the potential insights to be gained from extending next generation sequencing-supported network analyses to undomesticated non-model species.

Hybridization among dominant tree species correlates positively with understory plant diversity.

PREMISE OF THE STUDY: Elucidating the factors that determine the abundance and distribution of species remains a central goal of ecology. It is well recognized that genetic differences among individual species can affect the distribution and species interactions of dependent taxa, but the ecological effects of genetic differences on taxa of the same trophic level remain much less understood. Our goal was to test the hypothesis that differences between related overstory tree species and their hybrids can influence the understory plant community in wild settings. METHODS: We conducted vegetation surveys in a riparian community with the overstory dominated by Populus fremontii, P. angustifolia, and their natural hybrids (referred to as cross types) along the Weber River in north central Utah, USA. Understory diversity and community composition, as well as edaphic properties, were compared under individual trees. KEY RESULTS: Diversity metrics differ under the three different tree cross types such that a greater species richness, diversity, and cover of understory plants exist under the hybrids compared with either of the parental taxa (30-54%, 40-48%, and 35-74% greater, respectively). The community composition of the understory also varied by cross type, whereby additional understory plant species cluster with hybrids, not with parental species. CONCLUSIONS: Genetic composition dictated by hybridization in the overstory can play a role in structuring the associated understory plants in natural communities-where a hybridized overstory correlates with a species-rich understory-and thus can have cascading effects on community members of the same trophic level. The underlying mechanism requires further investigation.

A genetic basis for the manipulation of sink-source relationships by the galling aphid Pemphigus batae.

We examined how the galling aphid Pemphigus batae manipulates resource translocation patterns of resistant and susceptible narrowleaf cottonwood Populus angustifolia. Using carbon-14 ((14)C)-labeling experiments in common garden trials, five patterns emerged. First, although aphid galls on resistant and susceptible genotypes did not differ in their capacity to intercept assimilates exported from the leaf they occupied, aphids sequestered 5.8-fold more assimilates from surrounding leaves on susceptible tree genotypes compared to resistant genotypes. Second, gall sinks on the same side of a shoot as a labeled leaf were 3.4-fold stronger than gall sinks on the opposite side of a shoot, which agrees with patterns of vascular connections among leaves of the same shoot (orthostichy). Third, plant genetic-based traits accounted for 26% of the variation in sink strength of gall sinks and 41% of the variation in sink strength of a plant's own bud sinks. Fourth, tree susceptibility to aphid gall formation accounted for 63% of the variation in (14)C import, suggesting strong genetic control of sink-source relationships. Fifth, competition between two galls was observed on a susceptible but not a resistant tree. On the susceptible tree distal aphids intercepted 1.5-fold more (14)C from the occupied leaf than did basal aphids, but basal aphids compensated for the presence of a distal competitor by almost doubling import to the gall from surrounding leaves. These findings and others, aimed at identifying candidate genes for resistance, argue the importance of including plant genetics in future studies of the manipulation of translocation patterns by phytophageous insects.