Spatial isolation and genetic differentiation in naturally fragmented plant populations of the Swiss Alps

Aims The effect Of anthropogenic landscape fragmentation on the genetic diversity and adaptive potential of plant populations is a major issue in conservation biology. However, little is known about the partitioning of genetic diversity in alpine species, which occur in naturally fragmented habitats. Here, we, investigate molecular patterns of three alpine plants (Epilobium fleischeri, Geum reptans and Campanula thyrsoides) across Switzerland and ask whether Spatial isolation has led to high levels of populations differentiation, increasing over distance, and a decrease of within-population variability. We further hypothesize that file contrasting potential for long-distance dispersal (LDD) of Seed in these Species will considerably influence and explain diversity partitioning. Methods For each study species, we Sampled 20-23 individuals from each of 20-32 populations across entire Switzerland. We applied Random Amplified Polymorphic Dimorphism markers to assess genetic diversity within (Nei's expected heterozygosity, H-e; percentage of polymorphic hands, P-P) and among (analysis of molecular variance, Phi(st)) populations and correlated population size and altitude with within-populalion diversity. Spatial patterns of genetic relatedness were investigated using Mantel tests and standardized major axis regression as well as unweighted pair group method with arithmetic mean cluster analyses and Monmonier's algorithm. To avoid known biases, We standardized the numbers of populations, individuals and markers using multiple random reductions. We modelled LDD with a high alpine wind data set using the terminal velocity and height of seed release as key parameters. Additionally, we assessed a number of important life-history traits and factors that potentially influence genetic diversity partitioning (e.g. breeding system, longevity and population size). Important findings For all three species, We found a significant isolation-by-distance relationship but only a moderately high differentiation among populations (Phi(st): 22.7, 48 and 16.8%, for E. fleischeri, G. reptans and C. thyrsoides, respectively). Within-population diversity (H-c: 0.19-0.21, P-p: 62-75%) was not reduced in comparison to known results from lowland species and even small populations with < 50 reproductive individuals contained high levels of genetic diversity. We further found no indication that a high long-distance seed dispersal potential enhances genetic connectivity among populations. Gene flow seems to have a strong stochastic component causing large dissimilarity between population pairs irrespective of the spatial distance. Our results suggest that other life-history traits, especially the breeding System, may play an important role in genetic diversity partitioning. We conclude that spatial isolation in the alpine environment has a strong influence on population relatedness but that a number of factors can considerably influence the strength of this relationship.

DOUBLE-STRAND BREAK REPAIR PATHWAYS IN DNA STRUCTURE-INDUCED GENETIC INSTABILITY

Genetic instability in mammalian cells can occur by many different mechanisms. In the absence of exogenous sources of DNA damage, the DNA structure itself has been implicated in genetic instability. When the canonical B-DNA helix is naturally altered to form a non-canonical DNA structure such as a Z-DNA or H-DNA, this can lead to genetic instability in the form of DNA double-strand breaks (DSBs) (1, 2). Our laboratory found that the stability of these non-B DNA structures was different in mammals versus Escherichia coli (E.coli) bacteria (1, 2). One explanation for the difference between these species may be a result of how DSBs are repaired within each species. Non-homologous end-joining (NHEJ) is primed to repair DSBs in mammalian cells, while bacteria that lack NHEJ (such as E.coli), utilize homologous recombination (HR) to repair DSBs. To investigate the role of the error-prone NHEJ repair pathway in DNA structure-induced genetic instability, E.coli cells were modified to express genes to allow for a functional NHEJ system under different HR backgrounds. The Mycobacterium tuberculosis NHEJ sufficient system is composed of Ku and Ligase D (LigD) (3). These inducible NHEJ components were expressed individually and together in E.coli cells, with or without functional HR (RecA/RecB), and the Z-DNA and H-DNA-induced mutations were characterized. The Z-DNA structure gave rise to higher mutation frequencies compared to the controls, regardless of the DSB repair pathway(s) available; however, the type of mutants produced after repair was greatly dictated on the available DSB repair system, indicated by the shift from 2% large-scale deletions in the total mutant population to 24% large-scale deletions when NHEJ was present (4). This suggests that NHEJ has a role in the large deletions induced by Z-DNA-forming sequences. H-DNA structure, however, did not exhibit an increase in mutagenesis in the newly engineered E.coli environment, suggesting the involvement of other factors in regulating H-DNA formation/stability in bacterial cells. Accurate repair by established DNA DSB repair pathways is essential to maintain the stability of eukaryotic and prokaryotic genomes and our results suggest that an error-prone NHEJ pathway was involved in non-B DNA structure-induced mutagenesis in both prokaryotes and eukaryotes.

Using spatial linear models with SAR and CAR structure to examine Texas lung cancer incidence rates

Scholars have found that socioeconomic status was one of the key factors that influenced early-stage lung cancer incidence rates in a variety of regions. This thesis examined the association between median household income and lung cancer incidence rates in Texas counties. A total of 254 individual counties in Texas with corresponding lung cancer incidence rates from 2004 to 2008 and median household incomes in 2006 were collected from the National Cancer Institute Surveillance System. A simple linear model and spatial linear models with two structures, Simultaneous Autoregressive Structure (SAR) and Conditional Autoregressive Structure (CAR), were used to link median household income and lung cancer incidence rates in Texas. The residuals of the spatial linear models were analyzed with Moran's I and Geary's C statistics, and the statistical results were used to detect similar lung cancer incidence rate clusters and disease patterns in Texas.^

Making genetic biodiversity measurable : a review of statistical multivariate methods to study variability at gene level

Measures of agro-ecosystems genetic variability are essential to sustain scientific-based actions and policies tending to protect the ecosystem services they provide. To build the genetic variability datum it is necessary to deal with a large number and different types of variables. Molecular marker data is highly dimensional by nature, and frequently additional types of information are obtained, as morphological and physiological traits. This way, genetic variability studies are usually associated with the measurement of several traits on each entity. Multivariate methods are aimed at finding proximities between entities characterized by multiple traits by summarizing information in few synthetic variables. In this work we discuss and illustrate several multivariate methods used for different purposes to build the datum of genetic variability. We include methods applied in studies for exploring the spatial structure of genetic variability and the association of genetic data to other sources of information. Multivariate techniques allow the pursuit of the genetic variability datum, as a unifying notion that merges concepts of type, abundance and distribution of variability at gene level.

Mapping Genetic Diversity of Cherimoya (Annona cherimola Mill.): Application of Spatial Analysis for Conservation and Use of Plant Genetic Resources.

There is a growing call for inventories that evaluate geographic patterns in diversity of plant genetic resources maintained on farm and in species' natural populations in order to enhance their use and conservation. Such evaluations are relevant for useful tropical and subtropical tree species, as many of these species are still undomesticated, or in incipient stages of domestication and local populations can offer yet-unknown traits of high value to further domestication. For many outcrossing species, such as most trees, inbreeding depression can be an issue, and genetic diversity is important to sustain local production. Diversity is also crucial for species to adapt to environmental changes. This paper explores the possibilities of incorporating molecular marker data into Geographic Information Systems (GIS) to allow visualization and better understanding of spatial patterns of genetic diversity as a key input to optimize conservation and use of plant genetic resources, based on a case study of cherimoya (Annona cherimola Mill.), a Neotropical fruit tree species. We present spatial analyses to (1) improve the understanding of spatial distribution of genetic diversity of cherimoya natural stands and cultivated trees in Ecuador, Bolivia and Peru based on microsatellite molecular markers (SSRs); and (2) formulate optimal conservation strategies by revealing priority areas for in situ conservation, and identifying existing diversity gaps in ex situ collections. We found high levels of allelic richness, locally common alleles and expected heterozygosity in cherimoya's putative centre of origin, southern Ecuador and northern Peru, whereas levels of diversity in southern Peru and especially in Bolivia were significantly lower. The application of GIS on a large microsatellite dataset allows a more detailed prioritization of areas for in situ conservation and targeted collection across the Andean distribution range of cherimoya than previous studies could do, i.e. at province and department level in Ecuador and Peru, respectively.

Genetic characterization of chestnut (Castanea sativa Mill.) orchards and traditional nut varieties in El Bierzo, a glacial refuge and major cultivation site in Northwestern Spain

This paper presents a detailed genetic study of Castanea sativa in El Bierzo, a major nut production region with interesting features. It is located within a glacial refuge at one extreme of the distribution area (northwest Spain); it has a centenary tradition of chestnut management; and more importantly, it shows an unusual degree of genetic isolation. Seven nuclear microsatellite markers were selected to analyze the genetic variability and structure of 169 local trees grafted for nut production. We analyzed in the same manner 62 local nuts. The selected loci were highly discriminant for the genotypes studied, giving a combined probability of identity of 6.1 × 10−6. An unprecedented density of trees was sampled for this project over the entire region, and nuts were collected representing 18 cultivars marketed by local producers. Several instances of misclassification by local growers were detected. Fixation index estimates and analysis of molecular variance (AMOVA) data are supportive of an unexpectedly high level of genetic differentiation in El Bierzo, larger than that estimated in a previous study with broader geographical scope but based on limited local sampling (Pereira-Lorenzo et al., Tree Genet Genomes 6: 701–715, 2010a). Likewise, we have determined that clonality due to grafting had been previously overestimated. In line with these observations, no significant spatial structure was found using both a model-based Bayesian procedure and Mantel’s tests. Taken together, our results evidence the need for more fine-scale genetic studies if conservation strategies are to be efficiently improved.

Mycorrhizal preferences and fine spatial structure of the epiphytic orchid Epidendrum rhopalostele

• Premise of the study: The presence of compatible fungi is necessary for epiphytic orchid recruitment. Thus, identifying associated mycorrhizal fungi at the population level is essential for orchid conservation. Recruitment patterns may also be conditioned by factors such as seed dispersal range and specific environmental characteristics. • Methods: In a forest plot, all trees with a diameter at breast height >1 cm and all individuals of the epiphytic orchid Epidendrum rhopalostele were identified and mapped. Additionally, one flowering individual of E. rhopalostele per each host tree was randomly selected for root sampling and DNA extraction. • Key results: A total of 239 E. rhopalostele individuals were located in 25 of the 714 potential host trees. Light microscopy of sampled roots showed mycorrhizal fungi in 22 of the 25 sampled orchids. Phylogenetic analysis of ITS1-5.8S-ITS2 sequences yielded two Tulasnella clades. In four cases, plants were found to be associated with both clades. The difference between univariate and bivariate K functions was consistent with the random labeling null model at all spatial scales, indicating that trees hosting clades A and B of Tulasnella are not spatially segregated. The analysis of the inhomogenous K function showed that host trees are not clustered, suggesting no limitations to population-scale dispersal. χ2 analysis of contingency tables showed that E. rhopalostele is more frequent on dead trees than expected. • Conclusions: Epidendrum rhopalostele establishes mycorrhizal associations with at least two different Tulasnella species. The analysis of the distribution patterns of this orchid suggests a microsite preference for dead trees and no seed dispersal limitation.

Spatial network surrogates for disentangling complex system structure from spatial embedding of nodes

ACKNOWLEDGMENTS MW and RVD have been supported by the German Federal Ministry for Education and Research (BMBF) via the Young Investigators Group CoSy-CC2 (grant no. 01LN1306A). JFD thanks the Stordalen Foundation and BMBF (project GLUES) for financial support. JK acknowledges the IRTG 1740 funded by DFG and FAPESP. MT Gastner is acknowledged for providing his data on the airline, interstate, and Internet network. P Menck thankfully provided his data on the Scandinavian power grid. We thank S Willner on behalf of the entire zeean team for providing the data on the world trade network. All computations have been performed using the Python package pyunicorn [41] that is available at https://github.com/pik-copan/pyunicorn.

Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

Biodiversity of Costa Rican salamanders: Implications of high levels of genetic differentiation and phylogeographic structure for species formation

Although salamanders are characteristic amphibians in Holarctic temperate habitats, in tropical regions they have diversified evolutionarily only in tropical America. An adaptive radiation centered in Middle America occurred late in the history of a single clade, the supergenus Bolitoglossa (Plethodontidae), and large numbers of species now occur in diverse habitats. Sublineages within this clade decrease in number from the northern to southern parts of Middle America, and in Costa Rica, there are but three. Despite this phylogenetic constraint, Costa Rica has many species; the number of salamander species on one local elevational transect in the Cordillera de Talamanca may be the largest for any such transect in the world. Extraordinary variation in sequences of the mitochondrial gene cytochrome b within a clade of the genus Bolitoglossa in Costa Rica reveals strong phylogeographic structure within a single species, Bolitoglossa pesrubra. Allozymic variation in 19 proteins reveals a pattern largely concordant with the mitochondrial DNA phylogeography. More species exist than are currently recognized. Diversification occurs in restricted geographic areas and involves sharp geographic and elevational differentiation and zonation. In their degree of genetic differentiation at a local scale, these species of the deep tropics exceed the known variation of extratropical salamanders, which also differ in being less restricted in elevational range. Salamanders display “tropicality” in that although speciose, they are usually local in distribution and rare. They display strong ecological and physiological differentiation that may contribute importantly to morphological divergence and species formation.

Homologous genetic recombination as an intrinsic dynamic property of a DNA structure induced by RecA/Rad51-family proteins: A possible advantage of DNA over RNA as genomic material

Heteroduplex joints are general intermediates of homologous genetic recombination in DNA genomes. A heteroduplex joint is formed between a single-stranded region (or tail), derived from a cleaved parental double-stranded DNA, and homologous regions in another parental double-stranded DNA, in a reaction mediated by the RecA/Rad51-family of proteins. In this reaction, a RecA/Rad51-family protein first forms a filamentous complex with the single-stranded DNA, and then interacts with the double-stranded DNA in a search for homology. Studies of the three-dimensional structures of single-stranded DNA bound either to Escherichia coli RecA or Saccharomyces cerevisiae Rad51 have revealed a novel extended DNA structure. This structure contains a hydrophobic interaction between the 2′ methylene moiety of each deoxyribose and the aromatic ring of the following base, which allows bases to rotate horizontally through the interconversion of sugar puckers. This base rotation explains the mechanism of the homology search and base-pair switch between double-stranded and single-stranded DNA during the formation of heteroduplex joints. The pivotal role of the 2′ methylene-base interaction in the heteroduplex joint formation is supported by comparing the recombination of RNA genomes with that of DNA genomes. Some simple organisms with DNA genomes induce homologous recombination when they encounter conditions that are unfavorable for their survival. The extended DNA structure confers a dynamic property on the otherwise chemically and genetically stable double-stranded DNA, enabling gene segment rearrangements without disturbing the coding frame (i.e., protein-segment shuffling). These properties may give an extensive evolutionary advantage to DNA.

Genetic and comparative analyses reveal an alternative secondary structure in the region of nt 912 of Escherichia coli 16S rRNA.

Mutations at position 912 of Escherichia coli 16S rRNA result in two notable phenotypes. The C-->U transition confers resistance to streptomycin, a translational-error-inducing antibiotic, while a C-->G transversion causes marked retardation of cell growth rate. Starting with the slow-growing G912 mutant, random mutagenesis was used to isolate a second site mutation that restored growth nearly to the wild-type rate. The second site mutation was identified as a G-->C transversion at position 885 in 16S rRNA. Cells containing the G912 mutation had an increased doubling time, abnormal sucrose gradient ribosome/subunit profile, increased sensitivity to spectinomycin, dependence upon streptomycin for growth in the presence of spectinomycin, and slower translation rate, whereas cells with the G912/C885 double mutation were similar to wild type in these assays. Comparative analysis showed there was significant covariation between positions 912 and 885. Thus the second-site suppressor analysis, the functional assays, and the comparative data suggest that the interaction between nt 912 and nt 885 is conserved and necessary for normal ribosome function. Furthermore, the comparative data suggest that the interaction extends to include G885-G886-G887 pairing with C912-U911-C910. An alternative secondary structure element for the central domain of 16S rRNA is proposed.

Primary structure of hepatocyte nuclear factor/forkhead homologue 4 and characterization of gene expression in the developing respiratory and reproductive epithelium.

Members of the winged helix/forkhead family of transcription factors are believed to play a role in cell-specific gene expression. A cDNA encoding a member of this family of proteins, termed hepatocyte nuclear factor/forkhead homologue 4 (HFH-4), has been isolated from rat lung and rat testis cDNA libraries. This cDNA contains an open reading frame of 421 amino acids with a conserved DNA binding domain and several potential transactivating regions. During murine lung development, a single species of HFH-4-specific transcript (2.4 kb long) is first detected precisely at the start of the late pseudoglandular stage (embryonic day 14.5) and, by in situ hybridization, is specifically localized to the proximal pulmonary epithelium. The unique temporal and spatial pattern of HFH-4 gene expression in the developing lung defines this protein as a marker for the initiation of bronchial epithelial cell differentiation and suggests that it may play an important role in cell fate determination during lung development. In addition to expression in the pulmonary epithelium, RNA blot analysis reveals 2.4-kb HFH-4 transcripts in the testis and oviduct. By using mice with genetic defects in spermatogenesis, HFH-4 expression in the testis is found to be associated with the appearance of haploid germ cells and in situ hybridization studies indicate that HFH-4 expression is confined to stages I-VII of spermatogenesis. This pattern of HFH-4 gene expression during the early stages of differentiation of haploid germ cells suggests that HFH-4 may play a role in regulating stage-specific gene expression and cell-fate determination during lung development and in spermatogenesis.

Grouping genetic operators for the delineation of functional areas based on spatial interaction

The delineation of functional economic areas, or market areas, is a problem of high practical relevance, since the delineation of functional sets such as economic areas in the US, Travel-to-Work Areas in the United Kingdom, and their counterparts in other OECD countries are the basis of many statistical operations and policy making decisions at local level. This is a combinatorial optimisation problem defined as the partition of a given set of indivisible spatial units (covering a territory) into regions characterised by being (a) self-contained and (b) cohesive, in terms of spatial interaction data (flows, relationships). Usually, each region must reach a minimum size and self-containment level, and must be continuous. Although these optimisation problems have been typically solved through greedy methods, a recent strand of the literature in this field has been concerned with the use of evolutionary algorithms with ad hoc operators. Although these algorithms have proved to be successful in improving the results of some of the more widely applied official procedures, they are so time consuming that cannot be applied directly to solve real-world problems. In this paper we propose a new set of group-based mutation operators, featuring general operations over disjoint groups, tailored to ensure that all the constraints are respected during the operation to improve efficiency. A comparative analysis of our results with those from previous approaches shows that the proposed algorithm systematically improves them in terms of both quality and processing time, something of crucial relevance since it allows dealing with most large, real-world problems in reasonable time.

Genetic variation and inference of demographic histories in non-model species

Both long-term environmental changes such as those driven by the glacial cycles and more recent anthropogenic impacts have had major effects on the past demography in wild organisms. Within species, these changes are reflected in the amount and distribution of neutral genetic variation. In this thesis, mitochondrial and microsatellite DNA was analysed to investigate how environmental and anthropogenic factors have affected genetic diversity and structure in four ecologically different animal species. Paper I describes the post-glacial recolonisation history of the speckled-wood butterfly (Pararge aegeria) in Northern Europe. A decrease in genetic diversity with latitude and a marked population structure were uncovered, consistent with a hypothesis of repeated founder events during the postglacial recolonisation. Moreover, Approximate Bayesian Computation analyses indicate that the univoltine populations in Scandinavia and Finland originate from recolonisations along two routes, one on each side of the Baltic. Paper II aimed to investigate how past sea-level rises affected the population history of the convict surgeonfish (Acanthurus triostegus) in the Indo-Pacific. Assessment of the species’ demographic history suggested a population expansion that occurred approximately at the end of the last glaciation. Moreover, the results demonstrated an overall lack of phylogeographic structure, probably due to the high dispersal rates associated with the species’ pelagic larval stage. Populations at the species’ eastern range margin were significantly differentiated from other populations, which likely is a consequence of their geographic isolation. In Paper III, we assessed the effect of human impact on the genetic variation of European moose (Alces alces) in Sweden. Genetic analyses revealed a spatial structure with two genetic clusters, one in northern and one in southern Sweden, which were separated by a narrow transition zone. Moreover, demographic inference suggested a recent population bottleneck. The inferred timing of this bottleneck coincided with a known reduction in population size in the 19th and early 20th century due to high hunting pressure. In Paper IV, we examined the effect of an indirect but well-described human impact, via environmental toxic chemicals (PCBs), on the genetic variation of Eurasian otters (Lutra lutra) in Sweden. Genetic clustering assignment revealed differentiation between otters in northern and southern Sweden, but also in the Stockholm region. ABC analyses indicated a decrease in effective population size in both northern and southern Sweden. Moreover, comparative analyses of historical and contemporary samples demonstrated a more severe decline in genetic diversity in southern Sweden compared to northern Sweden, in agreement with the levels of PCBs found.