The Effects of Geography, Environment and Phylogeny on Community Assembly and Gene Flow Dynamics

It is increasingly evident that evolutionary processes play a role in how ecological communities are assembled. However the extend to which evolution influences how plants respond to spatial and environmental gradients and interact with each other is less clear. In this dissertation I leverage evolutionary tools and thinking to understand how space and environment affect community composition and patterns of gene flow in a unique system of Atlantic rainforest and restinga (sandy coastal plains) habitats in Southeastern Brazil.

In chapter one I investigate how space and environment affect the population genetic structure and gene flow of Aechmea nudicaulis, a bromeliad species that co-occurs in forest and restinga habitats. I genotyped seven microsatellite loci and sequenced one chloroplast DNA region for individuals collected in 7 pairs of forest / restinga sites. Bayesian genetic clustering analyses show that populations of A. nudicaulis are geographically structured in northern and southern populations, a pattern consistent with broader scale phylogeographic dynamics of the Atlantic rainforest. On the other hand, explicit migration models based on the coalescent estimate that inter-habitat gene flow is less common than gene flow between populations in the same habitat type, despite their geographic discontinuity. I conclude that there is evidence for repeated colonization of the restingas from forest populations even though the steep environmental gradient between habitats is a stronger barrier to gene flow than geographic distance.

In chapter two I use data on 2800 individual plants finely mapped in a restinga plot and on first-year survival of 500 seedlings to understand the roles of phylogeny, functional traits and abiotic conditions in the spatial structuring of that community. I demonstrate that phylogeny is a poor predictor of functional traits in and that convergence in these traits is pervasive. In general, the community is not phylogenetically structured, with at best 14% of the plots deviating significantly from the null model. The functional traits SLA, leaf dry matter content (LDMC), and maximum height also showed no clear pattern of spatial structuring. On the other hand, leaf area is strongly overdispersed across all spatial scales. Although leaf area overdispersion would be generally taken as evidence of competition, I argue that interpretation is probably misleading. Finally, I show that seedling survival is dramatically increased when they grow shaded by an adult individual, suggesting that seedlings are being facilitated. Phylogenetic distance to their adult neighbor has no influence on rates of survival though. Taken together, these results indicate that phylogeny has very limited influence on the fine scale assembly of restinga communities.

Using Bird Distributions to Assess Extinction Risk and Identify Conservation Priorities in Biodiversity Hotspots

Habitat loss, fragmentation, and degradation threaten the World’s ecosystems and species. These, and other threats, will likely be exacerbated by climate change. Due to a limited budget for conservation, we are forced to prioritize a few areas over others. These places are selected based on their uniqueness and vulnerability. One of the most famous examples is the biodiversity hotspots: areas where large quantities of endemic species meet alarming rates of habitat loss. Most of these places are in the tropics, where species have smaller ranges, diversity is higher, and ecosystems are most threatened.

Species distributions are useful to understand ecological theory and evaluate extinction risk. Small-ranged species, or those endemic to one place, are more vulnerable to extinction than widely distributed species. However, current range maps often overestimate the distribution of species, including areas that are not within the suitable elevation or habitat for a species. Consequently, assessment of extinction risk using these maps could underestimate vulnerability.

In order to be effective in our quest to conserve the World’s most important places we must: 1) Translate global and national priorities into practical local actions, 2) Find synergies between biodiversity conservation and human welfare, 3) Evaluate the different dimensions of threats, in order to design effective conservation measures and prepare for future threats, and 4) Improve the methods used to evaluate species’ extinction risk and prioritize areas for conservation. The purpose of this dissertation is to address these points in Colombia and other global biodiversity hotspots.

In Chapter 2, I identified the global, strategic conservation priorities and then downscaled to practical local actions within the selected priorities in Colombia. I used existing range maps of 171 bird species to identify priority conservation areas that would protect the greatest number of species at risk in Colombia (endemic and small-ranged species). The Western Andes had the highest concentrations of such species—100 in total—but the lowest densities of national parks. I then adjusted the priorities for this region by refining these species ranges by selecting only areas of suitable elevation and remaining habitat. The estimated ranges of these species shrank by 18–100% after accounting for habitat and suitable elevation. Setting conservation priorities on the basis of currently available range maps excluded priority areas in the Western Andes and, by extension, likely elsewhere and for other taxa. By incorporating detailed maps of remaining natural habitats, I made practical recommendations for conservation actions. One recommendation was to restore forest connections to a patch of cloud forest about to become isolated from the main Andes.

For Chapter 3, I identified areas where bird conservation met ecosystem service protection in the Central Andes of Colombia. Inspired by the November 11th (2011) landslide event near Manizales, and the current poor results of Colombia’s Article 111 of Law 99 of 1993 as a conservation measure in this country, I set out to prioritize conservation and restoration areas where landslide prevention would complement bird conservation in the Central Andes. This area is one of the most biodiverse places on Earth, but also one of the most threatened. Using the case of the Rio Blanco Reserve, near Manizales, I identified areas for conservation where endemic and small-range bird diversity was high, and where landslide risk was also high. I further prioritized restoration areas by overlapping these conservation priorities with a forest cover map. Restoring forests in bare areas of high landslide risk and important bird diversity yields benefits for both biodiversity and people. I developed a simple landslide susceptibility model using slope, forest cover, aspect, and stream proximity. Using publicly available bird range maps, refined by elevation, I mapped concentrations of endemic and small-range bird species. I identified 1.54 km2 of potential restoration areas in the Rio Blanco Reserve, and 886 km2 in the Central Andes region. By prioritizing these areas, I facilitate the application of Article 111 which requires local and regional governments to invest in land purchases for the conservation of watersheds.

Chapter 4 dealt with elevational ranges of montane birds and the impact of lowland deforestation on their ranges in the Western Andes of Colombia, an important biodiversity hotspot. Using point counts and mist-nets, I surveyed six altitudinal transects spanning 2200 to 2800m. Three transects were forested from 2200 to 2800m, and three were partially deforested with forest cover only above 2400m. I compared abundance-weighted mean elevation, minimum elevation, and elevational range width. In addition to analyzing the effect of deforestation on 134 species, I tested its impact within trophic guilds and habitat preference groups. Abundance-weighted mean and minimum elevations were not significantly different between forested and partially deforested transects. Range width was marginally different: as expected, ranges were larger in forested transects. Species in different trophic guilds and habitat preference categories showed different trends. These results suggest that deforestation may affect species’ elevational ranges, even within the forest that remains. Climate change will likely exacerbate harmful impacts of deforestation on species’ elevational distributions. Future conservation strategies need to account for this by protecting connected forest tracts across a wide range of elevations.

In Chapter 5, I refine the ranges of 726 species from six biodiversity hotspots by suitable elevation and habitat. This set of 172 bird species for the Atlantic Forest, 138 for Central America, 100 for the Western Andes of Colombia, 57 for Madagascar, 102 for Sumatra, and 157 for Southeast Asia met the criteria for range size, endemism, threat, and forest use. Of these 586 species, the Red List deems 108 to be threatened: 15 critically endangered, 29 endangered, and 64 vulnerable. When ranges are refined by elevational limits and remaining forest cover, 10 of those critically endangered species have ranges < 100km2, but then so do 2 endangered species, seven vulnerable, and eight non-threatened ones. Similarly, 4 critically endangered species, 20 endangered, and 12 vulnerable species have refined ranges < 5000km2, but so do 66 non-threatened species. A striking 89% of these species I have classified in higher threat categories have <50% of their refined ranges inside protected areas. I find that for 43% of the species I assessed, refined range sizes fall within thresholds that typically have higher threat categories than their current assignments. I recommend these species for closer inspection by those who assess risk. These assessments are not only important on a species-by-species basis, but by combining distributions of threatened species, I create maps of conservation priorities. They differ significantly from those created from unrefined ranges.

Microbial inactivation of Pseudomonas putida and Pichia pastoris using gene silencing.

Antisense deoxyoligonucleotide (ASO) gene silencing was investigated as a potential disinfection tool for industrial and drinking water treatment application. ASOs bind with their reverse complementary mRNA transcripts thereby blocking protein translation. While ASO silencing has mainly been studied in medicine, it may be useful for modulating gene expression and inactivating microorganisms in environmental applications. In this proof of concept work, gene targets were sh ble (zeocin resistance) and todE (catechol-2,3-dioxygenase) in Pichia pastoris and npt (kanamycin resistance) in Pseudomonas putida. A maximum 0.5-fold decrease in P. pastoris cell numbers was obtained following a 120 min incubation with single-stranded DNA (ssDNA) concentrations ranging from 0.2 to 200 nM as compared to the no ssDNA control. In P. putida, a maximum 5.2-fold decrease was obtained after 90 min with 400 nM ssDNA. While the silencing efficiencies varied for the 25 targets tested, these results suggest that protein activity as well as microbial growth can be altered using ASO gene silencing-based tools. If successful, this technology has the potential to eliminate some of the environmental and health issues associated with the use of strong chemical biocides. However, prior to its dissemination, more research is needed to increase silencing efficiency and develop effective delivery methods.

Morphologically cryptic biological species within the liverwort Frullania asagrayana.

UNLABELLED: PREMISE OF THE STUDY: The Frullania tamarisci complex includes eight Holarctic liverwort species. One of these, F. asagrayana, is distributed broadly throughout eastern North America from Canada to the Gulf Coast. Preliminary genetic data suggested that the species includes two groups of populations. This study was designed to test whether the two groups are reproductively isolated biological species. • METHODS: Eighty-eight samples from across the range of F. asagrayana, plus 73 samples from one population, were genotyped for 13 microsatellite loci. Sequences for two plastid loci and nrITS were obtained from 13 accessions. Genetic data were analyzed using coalescent models and Bayesian inference. • KEY RESULTS: Frullania asagrayana is sequence-invariant at the two plastid loci and ITS2, but two clear groups were resolved by microsatellites. The two groups are largely reproductively isolated, but there is a low level of gene flow from the southern to the northern group. No gene flow was detected in the other direction. A local population was heterogeneous but displayed strong genetic structure. • CONCLUSIONS: The genetic structure of F. asagrayana in eastern North America reflects morphologically cryptic differentiation between reproductively isolated groups of populations, near-panmixis within groups, and clonal propagation at local scales. Reproductive isolation between groups that are invariant at the level of nucleotide sequences shows that caution must be exercised in making taxonomic and evolutionary inferences from reciprocal monophyly (or lack thereof) between putative species.

The scale of population structure in Arabidopsis thaliana.

The population structure of an organism reflects its evolutionary history and influences its evolutionary trajectory. It constrains the combination of genetic diversity and reveals patterns of past gene flow. Understanding it is a prerequisite for detecting genomic regions under selection, predicting the effect of population disturbances, or modeling gene flow. This paper examines the detailed global population structure of Arabidopsis thaliana. Using a set of 5,707 plants collected from around the globe and genotyped at 149 SNPs, we show that while A. thaliana as a species self-fertilizes 97% of the time, there is considerable variation among local groups. This level of outcrossing greatly limits observed heterozygosity but is sufficient to generate considerable local haplotypic diversity. We also find that in its native Eurasian range A. thaliana exhibits continuous isolation by distance at every geographic scale without natural breaks corresponding to classical notions of populations. By contrast, in North America, where it exists as an exotic species, A. thaliana exhibits little or no population structure at a continental scale but local isolation by distance that extends hundreds of km. This suggests a pattern for the development of isolation by distance that can establish itself shortly after an organism fills a new habitat range. It also raises questions about the general applicability of many standard population genetics models. Any model based on discrete clusters of interchangeable individuals will be an uneasy fit to organisms like A. thaliana which exhibit continuous isolation by distance on many scales.

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

Delimiting species without nuclear monophyly in Madagascar's mouse lemurs.

BACKGROUND: Speciation begins when populations become genetically separated through a substantial reduction in gene flow, and it is at this point that a genetically cohesive set of populations attain the sole property of species: the independent evolution of a population-level lineage. The comprehensive delimitation of species within biodiversity hotspots, regardless of their level of divergence, is important for understanding the factors that drive the diversification of biota and for identifying them as targets for conservation. However, delimiting recently diverged species is challenging due to insufficient time for the differential evolution of characters--including morphological differences, reproductive isolation, and gene tree monophyly--that are typically used as evidence for separately evolving lineages. METHODOLOGY: In this study, we assembled multiple lines of evidence from the analysis of mtDNA and nDNA sequence data for the delimitation of a high diversity of cryptically diverged population-level mouse lemur lineages across the island of Madagascar. Our study uses a multi-faceted approach that applies phylogenetic, population genetic, and genealogical analysis for recognizing lineage diversity and presents the most thoroughly sampled species delimitation of mouse lemur ever performed. CONCLUSIONS: The resolution of a large number of geographically defined clades in the mtDNA gene tree provides strong initial evidence for recognizing a high diversity of population-level lineages in mouse lemurs. We find additional support for lineage recognition in the striking concordance between mtDNA clades and patterns of nuclear population structure. Lineages identified using these two sources of evidence also exhibit patterns of population divergence according to genealogical exclusivity estimates. Mouse lemur lineage diversity is reflected in both a geographically fine-scaled pattern of population divergence within established and geographically widespread taxa, as well as newly resolved patterns of micro-endemism revealed through expanded field sampling into previously poorly and well-sampled regions.

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.

Delivery and Scavenging of Nucleic Acids by Polycationic Polymers

Electrostatic interaction is a strong force that attracts positively and negatively charged molecules to each other. Such an interaction is formed between positively charged polycationic polymers and negatively charged nucleic acids. In this dissertation, the electrostatic attraction between polycationic polymers and nucleic acids is exploited for applications in oral gene delivery and nucleic acid scavenging. An enhanced nanoparticle for oral gene delivery of a human Factor IX (hFIX) plasmid is developed using the polycationic polysaccharide, chitosan (Ch), in combination with protamine sulfate (PS) to treat hemophilia B. For nucleic acid scavenging purposes, the development of an effective nucleic acid scavenging nanofiber platform is described for dampening hyper-inflammation and reducing the formation of biofilms.

Non-viral gene therapy may be an attractive alternative to chronic protein replacement therapy. Orally administered non-viral gene vectors have been investigated for more than one decade with little progress made beyond the initial studies. Oral administration has many benefits over intravenous injection including patient compliance and overall cost; however, effective oral gene delivery systems remain elusive. To date, only chitosan carriers have demonstrated successful oral gene delivery due to chitosan’s stability via the oral route. In this study, we increase the transfection efficiency of the chitosan gene carrier by adding protamine sulfate to the nanoparticle formulation. The addition of protamine sulfate to the chitosan nanoparticles results in up to 42x higher in vitro transfection efficiency than chitosan nanoparticles without protamine sulfate. Therapeutic levels of hFIX protein are detected after oral delivery of Ch/PS/phFIX nanoparticles in 5/12 mice in vivo, ranging from 3 -132 ng/mL, as compared to levels below 4 ng/mL in 1/12 mice given Ch/phFIX nanoparticles. These results indicate the protamine sulfate enhances the transfection efficiency of chitosan and should be considered as an effective ternary component for applications in oral gene delivery.

Dying cells release nucleic acids (NA) and NA-complexes that activate the inflammatory pathways of immune cells. Sustained activation of these pathways contributes to chronic inflammation related to autoimmune diseases including systemic lupus erythematosus, rheumatoid arthritis, and inflammatory bowel disease. Studies have shown that certain soluble, cationic polymers can scavenge extracellular nucleic acids and inhibit RNA-and DNA-mediated activation of Toll-like receptors (TLRs) and inflammation. In this study, the cationic polymers are incorporated onto insoluble nanofibers, enabling local scavenging of negatively charged pro-inflammatory species such as damage-associated molecular pattern (DAMP) molecules in the extracellular space, reducing cytotoxicity related to unwanted internalization of soluble cationic polymers. In vitro data show that electrospun nanofibers grafted with cationic polymers, termed nucleic acid scavenging nanofibers (NASFs), can scavenge nucleic acid-based agonists of TLR 3 and TLR 9 directly from serum and prevent the production of NF-ĸB, an immune system activating transcription factor while also demonstrating low cytotoxicity. NASFs formed from poly (styrene-alt-maleic anhydride) conjugated with 1.8 kDa branched polyethylenimine (bPEI) resulted in randomly aligned fibers with diameters of 486±9 nm. NASFs effectively eliminate the immune stimulating response of NA based agonists CpG (TLR 9) and poly (I:C) (TLR 3) while not affecting the activation caused by the non-nucleic acid TLR agonist pam3CSK4. Results in a more biologically relevant context of doxorubicin-induced cell death in RAW cells demonstrates that NASFs block ~25-40% of NF-ĸβ response in Ramos-Blue cells treated with RAW extracellular debris, ie DAMPs, following doxorubicin treatment. Together, these data demonstrate that the formation of cationic NASFs by a simple, replicable, modular technique is effective and that such NASFs are capable of modulating localized inflammatory responses.

An understandable way to clinically apply the NASF is as a wound bandage. Chronic wounds are a serious clinical problem that is attributed to an extended period of inflammation as well as the presence of biofilms. An NASF bandage can potentially have two benefits in the treatment of chronic wounds by reducing the inflammation and preventing biofilm formation. NASF can prevent biofilm formation by reducing the NA present in the wound bed, therefore removing large components of what the bacteria use to develop their biofilm matrix, the extracellular polymeric substance, without which the biofilm cannot develop. The NASF described above is used to show the effect of the nucleic acid scavenging technology on in vitro and in vivo biofilm formation of P. aeruginosa, S. aureus, and S. epidermidis biofilms. The in vitro studies demonstrated that the NASFs were able to significantly reduce the biofilm formation in all three bacterial strains. In vivo studies of the NASF on mouse wounds infected with biofilm show that the NASF retain their functionality and are able to scavenge DNA, RNA, and protein from the wound bed. The NASF remove DNA that are maintaining the inflammatory state of the open wound and contributing to the extracellular polymeric substance (EPS), such as mtDNA, and also removing proteins that are required for bacteria/biofilm formation and maintenance such as chaperonin, ribosomal proteins, succinyl CoA-ligase, and polymerases. However, the NASF are not successful at decreasing the wound healing time because their repeated application and removal disrupts the wound bed and removes proteins required for wound healing such as fibronectin, vibronectin, keratin, and plasminogen. Further optimization of NASF treatment duration and potential combination treatments should be tested to reduce the unwanted side effects of increased wound healing time.