The Roles of Phenotypic Plasticity and Plant-Microbe Interactions in the Evolution of Complex Traits in Boechera stricta

All organisms live in complex habitats that shape the course of their evolution by altering the phenotype expressed by a given genotype (a phenomenon known as phenotypic plasticity) and simultaneously by determining the evolutionary fitness of that phenotype. In some cases, phenotypic evolution may alter the environment experienced by future generations. This dissertation describes how genetic and environmental variation act synergistically to affect the evolution of glucosinolate defensive chemistry and flowering time in Boechera stricta, a wild perennial herb. I focus particularly on plant-associated microbes as a part of the plant’s environment that may alter trait evolution and in turn be affected by the evolution of those traits. In the first chapter I measure glucosinolate production and reproductive fitness of over 1,500 plants grown in common gardens in four diverse natural habitats, to describe how patterns of plasticity and natural selection intersect and may influence glucosinolate evolution. I detected extensive genetic variation for glucosinolate plasticity and determined that plasticity may aid colonization of new habitats by moving phenotypes in the same direction as natural selection. In the second chapter I conduct a greenhouse experiment to test whether naturally-occurring soil microbial communities contributed to the differences in phenotype and selection that I observed in the field experiment. I found that soil microbes cause plasticity of flowering time but not glucosinolate production, and that they may contribute to natural selection on both traits; thus, non-pathogenic plant-associated microbes are an environmental feature that could shape plant evolution. In the third chapter, I combine a multi-year, multi-habitat field experiment with high-throughput amplicon sequencing to determine whether B. stricta-associated microbial communities are shaped by plant genetic variation. I found that plant genotype predicts the diversity and composition of leaf-dwelling bacterial communities, but not root-associated bacterial communities. Furthermore, patterns of host genetic control over associated bacteria were largely site-dependent, indicating an important role for genotype-by-environment interactions in microbiome assembly. Together, my results suggest that soil microbes influence the evolution of plant functional traits and, because they are sensitive to plant genetic variation, this trait evolution may alter the microbial neighborhood of future B. stricta generations. Complex patterns of plasticity, selection, and symbiosis in natural habitats may impact the evolution of glucosinolate profiles in Boechera stricta.

Distinct routes to metastasis: plasticity-dependent and plasticity-independent pathways.

The cascade that culminates in macrometastases is thought to be mediated by phenotypic plasticity, including epithelial-mesenchymal and mesenchymal-epithelial transitions (EMT and MET). Although there is substantial support for the role of EMT in driving cancer cell invasion and dissemination, much less is known about the importance of MET in the later steps of metastatic colonization. We created novel reporters, which integrate transcriptional and post-transcriptional regulation, to test whether MET is required for metastasis in multiple in vivo cancer models. In a model of carcinosarcoma, metastasis occurred via an MET-dependent pathway; however, in two prostate carcinoma models, metastatic colonization was MET independent. Our results provide evidence for both MET-dependent and MET-independent metastatic pathways.

Behavior Genetics and Post Genomics

The science of genetics is undergoing a paradigm shift. Recent discoveries, including the activity of retrotransposons, the extent of copy number variations, somatic and chromosomal mosaicism, and the nature of the epigenome as a regulator of DNA expressivity, are challenging a series of dogmas concerning the nature of the genome and the relationship between genotype and phenotype. DNA, once held to be the unchanging template of heredity, now appears subject to a good deal of environmental change; considered to be identical in all cells and tissues of the body, there is growing evidence that somatic mosaicism is the normal human condition; and treated as the sole biological agent of heritability, we now know that the epigenome, which regulates gene expressivity, can be inherited via the germline. These developments are particularly significant for behavior genetics for at least three reasons: First, these phenomena appear to be particularly prevalent in the human brain, and likely are involved in much of human behavior; second, they have important implications for the validity of heritability and gene association studies, the methodologies that largely define the discipline of behavior genetics; and third, they appear to play a critical role in development during the perinatal period, and in enabling phenotypic plasticity in offspring in particular. I examine one of the central claims to emerge from the use of heritability studies in the behavioral sciences, the principle of “minimal shared maternal effects,” in light of the growing awareness that the maternal perinatal environment is a critical venue for the exercise of adaptive phenotypic plasticity. This consideration has important implications for both developmental and evolutionary biology

Multiple phenotypic changes in mice after knockout of the B3gnt5 gene, encoding Lc3 synthase--a key enzyme in lacto-neolacto ganglioside synthesis.

BACKGROUND: Ganglioside biosynthesis occurs through a multi-enzymatic pathway which at the lactosylceramide step is branched into several biosynthetic series. Lc3 synthase utilizes a variety of galactose-terminated glycolipids as acceptors by establishing a glycosidic bond in the beta-1,3-linkage to GlcNaAc to extend the lacto- and neolacto-series gangliosides. In order to examine the lacto-series ganglioside functions in mice, we used gene knockout technology to generate Lc3 synthase gene B3gnt5-deficient mice by two different strategies and compared the phenotypes of the two null mouse groups with each other and with their wild-type counterparts. RESULTS: B3gnt5 gene knockout mutant mice appeared normal in the embryonic stage and, if they survived delivery, remained normal during early life. However, about 9% developed early-stage growth retardation, 11% died postnatally in less than 2 months, and adults tended to die in 5-15 months, demonstrating splenomegaly and notably enlarged lymph nodes. Without lacto-neolacto series gangliosides, both homozygous and heterozygous mice gradually displayed fur loss or obesity, and breeding mice demonstrated reproductive defects. Furthermore, B3gnt5 gene knockout disrupted the functional integrity of B cells, as manifested by a decrease in B-cell numbers in the spleen, germinal center disappearance, and less efficiency to proliferate in hybridoma fusion. CONCLUSIONS: These novel results demonstrate unequivocally that lacto-neolacto series gangliosides are essential to multiple physiological functions, especially the control of reproductive output, and spleen B-cell abnormality. We also report the generation of anti-IgG response against the lacto-series gangliosides 3'-isoLM1 and 3',6'-isoLD1.

Independent evolutionary origins of landlocked alewife populations and rapid parallel evolution of phenotypic traits.

Alewife, Alosa pseudoharengus, populations occur in two discrete life-history variants, an anadromous form and a landlocked (freshwater resident) form. Landlocked populations display a consistent pattern of life-history divergence from anadromous populations, including earlier age at maturity, smaller adult body size, and reduced fecundity. In Connecticut (USA), dams constructed on coastal streams separate anadromous spawning runs from lake-resident landlocked populations. Here, we used sequence data from the mtDNA control region and allele frequency data from five microsatellite loci to ask whether coastal Connecticut landlocked alewife populations are independently evolved from anadromous populations or whether they share a common freshwater ancestor. We then used microsatellite data to estimate the timing of the divergence between anadromous and landlocked populations. Finally, we examined anadromous and landlocked populations for divergence in foraging morphology and used divergence time estimates to calculate the rate of evolution for foraging traits. Our results indicate that landlocked populations have evolved multiple times independently. Tests of population divergence and estimates of gene flow show that landlocked populations are genetically isolated, whereas anadromous populations exchange genes. These results support a 'phylogenetic raceme' model of landlocked alewife divergence, with anadromous populations forming an ancestral core from which landlocked populations independently diverged. Divergence time estimates suggest that landlocked populations diverged from a common anadromous ancestor no longer than 5000 years ago and perhaps as recently as 300 years ago, depending on the microsatellite mutation rate assumed. Examination of foraging traits reveals landlocked populations to have significantly narrower gapes and smaller gill raker spacings than anadromous populations, suggesting that they are adapted to foraging on smaller prey items. Estimates of evolutionary rates (in haldanes) indicate rapid evolution of foraging traits, possibly in response to changes in available resources.

Experimental evidence that phenotypic divergence in predators drives community divergence in prey.

Studies of adaptive divergence have traditionally focused on the ecological causes of trait diversification, while the ecological consequences of phenotypic divergence remain relatively unexplored. Divergence in predator foraging traits, in particular, has the potential to impact the structure and dynamics of ecological communities. To examine the effects of predator trait divergence on prey communities, we exposed zooplankton communities in lake mesocosms to predation from either anadromous or landlocked (freshwater resident) alewives, which have undergone recent and rapid phenotypic differentiation in foraging traits (gape width, gill raker spacing, and prey size-selectivity). Anadromous alewives, which exploit large prey items, significantly reduced the mean body size, total biomass, species richness, and diversity of crustacean zooplankton relative to landlocked alewives, which exploit smaller prey. The zooplankton responses observed in this experiment are consistent with patterns observed in lakes. This study provides direct evidence that phenotypic divergence in predators, even in its early stages, can play a critical role in determining prey community structure.

Phenex: ontological annotation of phenotypic diversity.

BACKGROUND: Phenotypic differences among species have long been systematically itemized and described by biologists in the process of investigating phylogenetic relationships and trait evolution. Traditionally, these descriptions have been expressed in natural language within the context of individual journal publications or monographs. As such, this rich store of phenotype data has been largely unavailable for statistical and computational comparisons across studies or integration with other biological knowledge. METHODOLOGY/PRINCIPAL FINDINGS: Here we describe Phenex, a platform-independent desktop application designed to facilitate efficient and consistent annotation of phenotypic similarities and differences using Entity-Quality syntax, drawing on terms from community ontologies for anatomical entities, phenotypic qualities, and taxonomic names. Phenex can be configured to load only those ontologies pertinent to a taxonomic group of interest. The graphical user interface was optimized for evolutionary biologists accustomed to working with lists of taxa, characters, character states, and character-by-taxon matrices. CONCLUSIONS/SIGNIFICANCE: Annotation of phenotypic data using ontologies and globally unique taxonomic identifiers will allow biologists to integrate phenotypic data from different organisms and studies, leveraging decades of work in systematics and comparative morphology.

To plasticity and back again.

Both the gain and the loss of flexibility in the development of phenotypes have led to an increased diversity of physical forms in nematode worms.

Evaluation of an epithelial plasticity biomarker panel in men with localized prostate cancer.

BACKGROUND: Given the potential importance of epithelial plasticity (EP) to cancer metastasis, we sought to investigate biomarkers related to EP in men with localized prostate cancer (PC) for the association with time to PSA recurrence and other clinical outcomes after surgery. METHODS: Men with localized PC treated with radical prostatectomy at the Durham VA Medical Center and whose prostatectomy tissues were included in a tissue microarray (TMA) linked to long-term outcomes. We performed immunohistochemical studies using validated antibodies against E-cadherin and Ki-67 and mesenchymal biomarkers including N-cadherin, vimentin, SNAIL, ZEB1 and TWIST. Association studies were conducted for each biomarker with baseline clinical/pathologic characteristics an risk of PSA recurrence over time. RESULTS: Two hundred and five men contributed TMA tissue and had long-term follow-up (median 11 years). Forty-three percent had PSA recurrence; three died of PC. The majority had high E-cadherin expression (86%); 14% had low/absent E-cadherin expression. N-cadherin was rarely expressed (<4%) and we were unable to identify an E-to-N-cadherin switch as independently prognostic. No associations with clinical risk group, PSA recurrence or Gleason sum were noted for SNAIL, ZEB1, vimentin or TWIST, despite heterogeneous expression between patients. We observed an association of higher Ki-67 expression with Gleason sum (P=0.043), National Comprehensive Cancer Network risk (P=0.013) and PSA recurrence (hazard ratio 1.07, P=0.016). CONCLUSIONS: The expression of EP biomarkers in this cohort of men with a low risk of PC-specific mortality was not associated with aggressive features or PSA relapse after surgery.

BDNF-TrkB Signaling in Single-Spine Structural Plasticity

Multiple lines of evidence reveal that activation of the tropomyosin related kinase B (TrkB) receptor is a critical molecular mechanism underlying status epilepticus (SE) induced epilepsy development. However, the cellular consequences of such signaling remain unknown. To this point, localization of SE-induced TrkB activation to CA1 apical dendritic spines provides an anatomic clue pointing to Schaffer collateral-CA1 synaptic plasticity as one potential cellular consequence of TrkB activation. Here, we combine two-photon glutamate uncaging with two photon fluorescence lifetime imaging microscopy (2pFLIM) of fluorescence resonance energy transfer (FRET)-based sensors to specifically investigate the roles of TrkB and its canonical ligand brain derived neurotrophic factor (BDNF) in dendritic spine structural plasticity (sLTP) of CA1 pyramidal neurons in cultured hippocampal slices of rodents. To begin, we demonstrate a critical role for post-synaptic TrkB and post-synaptic BDNF in sLTP. Building on these findings, we develop a novel FRET-based sensor for TrkB activation that can report both BDNF and non-BDNF activation in a specific and reversible manner. Using this sensor, we monitor the spatiotemporal dynamics of TrkB activity during single-spine sLTP. In response to glutamate uncaging, we report a rapid (onset less than 1 minute) and sustained (lasting at least 20 minutes) activation of TrkB in the stimulated spine that depends on N-methyl-D-aspartate receptor (NMDAR)-Ca2+/Calmodulin dependent kinase II (CaMKII) signaling as well as post-synaptically synthesized BDNF. Consistent with these findings, we also demonstrate rapid, glutamate uncaging-evoked, time-locked release of BDNF from single dendritic spines using BDNF fused to superecliptic pHluorin (SEP). Finally, to elucidate the molecular mechanisms by which TrkB activation leads to sLTP, we examined the dependence of Rho GTPase activity - known mediators of sLTP - on BDNF-TrkB signaling. Through the use of previously described FRET-based sensors, we find that the activities of ras-related C3 botulinum toxin substrate 1 (Rac1) and cell division control protein 42 (Cdc42) require BDNF-TrkB signaling. Taken together, these findings reveal a spine-autonomous, autocrine signaling mechanism involving NMDAR-CaMKII dependent BDNF release from stimulated dendritic spines leading to TrkB activation and subsequent activation of the downstream molecules Rac1 and Cdc42 in these same spines that proves critical for sLTP. In conclusion, these results highlight structural plasticity as one cellular consequence of CA1 dendritic spine TrkB activation that may potentially contribute to larger, circuit-level changes underlying SE-induced epilepsy.

Calcium/Calmodulin-Dependent Protein Kinase II Serves as a Biochemical Integrator of Calcium Signals for the Induction of Synaptic Plasticity

Repetitive Ca2+ transients in dendritic spines induce various forms of synaptic plasticity by transmitting information encoded in their frequency and amplitude. CaMKII plays a critical role in decoding these Ca2+ signals to initiate long-lasting synaptic plasticity. However, the properties of CaMKII that mediate Ca2+ decoding in spines remain elusive. Here, I measured CaMKII activity in spines using fast-framing two-photon fluorescence lifetime imaging. Following each repetitive Ca2+ elevations, CaMKII activity increased in a stepwise manner. This signal integration, at the time scale of seconds, critically depended on Thr286 phosphorylation. In the absence of Thr286 phosphorylation, only by increasing the frequency of repetitive Ca2+ elevations could high peak CaMKII activity or plasticity be induced. In addition, I measured the association between CaMKII and Ca2+/CaM during spine plasticity induction. Unlike CaMKII activity, association of Ca2+/CaM to CaMKII plateaued at the first Ca2+ elevation event. This result indicated that integration of Ca2+ signals was initiated by the binding of Ca2+/CaM and amplified by the subsequent increases in Thr286-phosphorylated form of CaMKII. Together, these findings demonstrate that CaMKII functions as a leaky integrator of repetitive Ca2+ signals during the induction of synaptic plasticity, and that Thr286 phosphorylation is critical for defining the frequencies of such integration.

Spatiotemporal Dynamics of CaMKI During Structural Plasticity of Single Dendritic Spines

Multifunctional calcium/calmodulin dependent protein kinases (CaMKs) are key regulators of spine structural plasticity and long-term potentiation (LTP) in neurons. CaMKs have promiscuous and overlapping substrate recognition motifs, and are distinguished in their regulatory role based on differences in the spatiotemporal dynamics of activity. While the function and activity of CaMKII in synaptic plasticity has been extensively studied, that of CaMKI, another major class of CaMK required for LTP, still remain elusive.

Here, we develop a Förster’s Resonance Energy Transfer (FRET) based sensor to measure the spatiotemporal activity dynamics of CaMK1. We monitored CaMKI activity using 2-photon fluorescence lifetime imaging, while inducing LTP in single dendritic spines of rat (Rattus Norvegicus, strain Sprague Dawley) hippocampal CA1 pyramidal neurons using 2-photon glutamate uncaging. Using RNA-interference and pharmacological means, we also characterize the role of CaMKI during spine structural plasticity.

We found that CaMKI was rapidly and transiently activated with a rise time of ~0.3 s and decay time of ~1 s in response to each uncaging pulse. Activity of CaMKI spread out of the spine. Phosphorylation of CaMKI by CaMKK was required for this spreading and for the initial phase of structural LTP. Combined with previous data showing that CaMKII is restricted to the stimulated spine and required for long-term maintenance of structural LTP, these results suggest that CaMK diversity allows the same incoming signal – calcium – to independently regulate distinct phases of LTP by activating different CaMKs with distinct spatiotemporal dynamics.