Genetic effects of chronic habitat fragmentation on tree species: the case of Sorbus aucuparia in a deforested Scottish landscape

Sustainable forest restoration and management practices require a thorough understanding of the influence that habitat fragmentation has on the processes shaping genetic variation and its distribution in tree populations. We quantified genetic variation at isozyme markers and chloroplast DNA (cpDNA), analysed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) in severely fragmented populations of Sorbus aucuparia (Rosaceae) in a single catchment (Moffat) in southern Scotland. Remnants maintain surprisingly high levels of gene diversity (H-E) for isozymes (H-E = 0.195) and cpDNA markers (H-E = 0.490). Estimates are very similar to those from non-fragmented populations in continental Europe, even though the latter were sampled over a much larger spatial scale. Overall, no genetic bottleneck or departures from random mating were detected in the Moffat fragments. However, genetic differentiation among remnants was detected for both types of marker (isozymes Theta(n) = 0.043, cpDNA Theta(c) = 0.131; G-test, P-value < 0.001). In this self-incompatible, insect-pollinated, bird-dispersed tree species, the estimated ratio of pollen flow to seed flow between fragments is close to 1 (r = 1.36). Reduced pollen-mediated gene flow is a likely consequence of habitat fragmentation, but effective seed dispersal by birds is probably helping to maintain high levels of genetic diversity within remnants and reduce genetic differentiation between them.

Limited effect of anthropogenic habitat fragmentation on molecular diversity in a rain forest skink, Gnypetoscincus queenslandiae

To examine the effects of recent habitat fragmentation, we assayed genetic diversity in a rain forest endemic lizard, the prickly forest skink (Gnypetoscincus queenslandiae), from seven forest fragments and five sites in continuous forest on the Atherton tableland of northeastern Queensland, Australia. The rain forest in this region was fragmented by logging and clearing for dairy farms in the early 1900s and most forest fragments studied have been isolated for 50-80 years or nine to 12 skink generations. We genotyped 411 individuals at nine microsatellite DNA loci and found fewer alleles per locus in prickly forest skinks from small rain forest fragments and a lower ratio of allele number to allele size range in forest fragments than in continuous forest, indicative of a decrease in effective population size. In contrast, and as expected for populations with small neighbourhood sizes, neither heterozygosity nor variance in allele size differed between fragments and sites in continuous forests. Considering measures of among population differentiation, there was no increase in F-ST among fragments and a significant isolation by distance pattern was identified across all 12 sites. However, the relationship between genetic (F-ST) and geographical distance was significantly stronger for continuous forest sites than for fragments, consistent with disruption of gene flow among the latter. The observed changes in genetic diversity within and among populations are small, but in the direction predicted by the theory of genetic erosion in recently fragmented populations. The results also illustrate the inherent difficulty in detecting genetic consequences of recent habitat fragmentation, even in genetically variable species, and especially when effective population size and dispersal rates are low.

The effects of habitat fragmentation due to forestry plantation establishment on the demography and genetic variation of a marsupial carnivore, Antechinus agilis

We conducted a demographic and genetic study to investigate the effects of fragmentation due to the establishment of an exotic softwood plantation on populations of a small marsupial carnivore, the agile antechinus (Antechinus agilis), and the factors influencing the persistence of those populations in the fragmented habitat. The first aspect of the study was a descriptive analysis of patch occupancy and population size, in which we found a patch occupancy rate of 70% among 23 sites in the fragmented habitat compared to 100% among 48 sites with the same habitat characteristics in unfragmented habitat. Mark-recapture analyses yielded most-likely population size estimates of between 3 and 85 among the 16 occupied patches in the fragmented habitat. Hierarchical partitioning and model selection were used to identify geographic and habitat-related characteristics that influence patch occupancy and population size. Patch occupancy was primarily influenced by geographic isolation and habitat quality (vegetation basal area). The variance in population size among occupied sites was influenced primarily by forest type (dominant Eucalyptus species) and, to a lesser extent, by patch area and topographic context (gully sites had larger populations). A comparison of the sex ratios between the samples from the two habitat contexts revealed a significant deficiency of males in the fragmented habitat. We hypothesise that this is due to male-biased dispersal in an environment with increased dispersal-associated mortality. The population size and sex ratio data were incorporated into a simulation study to estimate the proportion of genetic diversity that would have been lost over the known timescale since fragmentation if the patch populations had been totally isolated. The observed difference in genetic diversity (gene diversity and allelic richness at microsatellite and mitochondrial markers) between 16 fragmented and 12 unfragmented sites was extremely low and inconsistent with the isolation of the patch populations. Our results show that although the remnant habitat patches comprise approximately 2% of the study area, they can support non-isolated populations. However, the distribution of agile antechinus populations in the fragmented system is dependent on habitat quality and patch connectivity. (C) 2004 Elsevier Ltd. All rights reserved.

Genetic resource impacts of habitat loss and degradation; reconciling empirical evidence and predicted theory for neotropical trees

The theoretical impacts of anthropogenic habitat degradation on genetic resources have been well articulated. Here we use a simulation approach to assess the magnitude of expected genetic change, and review 31 studies of 23 neotropical tree species to assess whether empirical case studies conform to theory. Major differences in the sensitivity of measures to detect the genetic health of degraded populations were obvious. Most studies employing genetic diversity (nine out of 13) found no significant consequences, yet most that assessed progeny inbreeding (six out of eight), reproductive output (seven out of 10) and fitness (all six) highlighted significant impacts. These observations are in line with theory, where inbreeding is observed immediately following impact, but genetic diversity is lost slowly over subsequent generations, which for trees may take decades. Studies also highlight the ecological, not just genetic, consequences of habitat degradation that can cause reduced seed set and progeny fitness. Unexpectedly, two studies examining pollen flow using paternity analysis highlight an extensive network of gene flow at smaller spatial scales (less than 10 km). Gene flow can thus mitigate against loss of genetic diversity and assist in long-term population viability, even in degraded landscapes. Unfortunately, the surveyed studies were too few and heterogeneous to examine concepts of population size thresholds and genetic resilience in relation to life history. Future suggested research priorities include undertaking integrated studies on a range of species in the same landscapes; better documentation of the extent and duration of impact; and most importantly, combining neutral marker, pollination dynamics, ecological consequences, and progeny fitness assessment within single studies.

Chloroplast and total genomic diversity in the endemic Costa Rican tree Lonchocarpus costaricensis (J.D. Smith) Pittier (Papilionaceae).

In Mesoamerica, tropical dry forest is a highly threatened habitat, and species endemic to this environment are under extreme pressure. The tree species, Lonchocarpus costaricensis is endemic to the dry northwest of Costa Rica and southwest Nicaragua. It is a locally important species but, as land has been cleared for agriculture, populations have experienced considerable reduction and fragmentation. To assess current levels and distribution of genetic diversity in the species, a combination of chloroplast-specific (cpDNA) and whole genome DNA markers (amplified fragment length polymorphism, AFLP) were used to fingerprint 121 individual trees in 6 populations. Two cpDNA haplotypes were identified, distributed among populations such that populations at the extremes of the distribution showed lowest diversity. A large number (487) of AFLP markers were obtained and indicated that diversity levels were highest in the two coastal populations (Cobano, Matapalo, H = 0.23, 0.28 respectively). Population differentiation was low overall, F-ST = 0.12, although Matapalo was strongly differentiated from all other populations (F-ST = 0.16-0.22), apart from Cobano (F., = 0.11). Spatial genetic structure was present in both datasets at different scales: cpDNA was structured at a range-wide distribution scale, whilst AFLP data revealed genetic neighbourhoods on a population scale. In general, the habitat degradation of recent times appears not to have yet impacted diversity levels in mature populations. However, although no data on seed or saplings were collected, it seems likely that reproductive mechanisms in the species will have been affected by land clearance. It is recommended that efforts should be made to conserve the extant genetic resource base and further research undertaken to investigate diversity levels in the progeny generation.

Contrasting quantitative traits and neutral genetic markers for genetic resource assessment of Mesoamerican Cedrela odorata

We compared within-population variability and degree of population differentiation for neutral genetic markers (RAPDS) and eight quantitative traits in Central American populations of the endangered tree, Cedrela odorata. Whilst population genetic diversity for neutral markers (Shannon index) and quantitative traits (heritability, coefficient of additive genetic variation) were uncorrelated, both marker types revealed strong differentiation between populations from the Atlantic coast of Costa Rica and the rest of the species' distribution. The degree of interpopulation differentiation was higher for RAPD markers (F-ST 0.67 for the sampled Mesoamerican range) than for quantitative traits (Q(ST) = 0.30). Hence, the divergence in quantitative traits was lower than could have been achieved by genetic drift alone, suggesting that balancing selection for similar phenotypes in different populations of this species. Nevertheless, a comparison of pair-wise estimates of population differentiation in neutral genetic markers and quantitative traits revealed a strong positive correlation (r = 0.66) suggesting that, for C. odorata, neutral marker divergence could be used as a surrogate for adaptive gene divergence for conservation planning. The utility of this finding and suggested further work are discussed.

Altitudinal variation in behavioural thermoregulation: Local adaptation vs. plasticity in California grasshoppers

We investigated the adaptive significance of behavioural thermoregulation in univoltine populations of the grasshopper Melanoplus sanguinipes along an altitudinal gradient in California using laboratory tests of animals raised under different temperatures. Trials consisted of continuous body temperature measurements with semi-implanted microprobes in a test arena, and observation and simultaneous recording of behavioural responses. These responses included mobility, basking and orientation of the body axes (aspect angle) towards a radiation source. Mobility and basking are determined by the altitudinal origin of the parental generation and not by the temperature treatments. With increasing altitude, individuals tend increasingly to raise body temperatures via mobility and increased basking. In contrast, body orientation towards the radiation source is influenced by the temperature treatments but not by the altitude of origin. Individuals experiencing higher temperatures during rearing show a lower tendency to lateral flanking. We conclude that body orientation responses are not adapted locally. In contrast other components of the behavioural syndrome that increase body temperature, such as mobility and basking, are adaptive in response to local selection pressure. The thermoregulatory syndrome of these grasshoppers is an important contribution to life-history adaptations that appropriately match season lengths.

Genetic uniformity of international isolates of Leifsonia xyli subsp xyli, causal agent of ratoon stunting disease of sugarcane

An international collection of the sugarcane ratoon stunting disease pathogen, Leifsonia xyli subsp. xyli, was analysed to assess genetic diversity. DNA fingerprinting using BOX primers was performed on 105 isolates, comprising 65 Australian isolates and an additional 40 isolates from Indonesia (n = 8), Japan (n = 1), USA (n = 3), Brazil (n = 2), Mali (n = 2), Zimbabwe (n = 13), South Africa (n = 9) and Reunion (n = 2). Sixty-two of these isolates were also screened using ERIC primers. No variation was found among any of the isolates. The intergenic spacer (IGS) region of the ribosomal RNA genes from 54 isolates was screened for sequence variation using single-stranded conformational polymorphism (SSCP), but none was observed. Direct sequencing of the IGS from a subset of nine isolates, representing all of the countries sampled in this study, confirmed the results of the SSCP analysis. Likewise, no sequence variation was found in the 16S ribosomal RNA genes of the same subset. Four Colombian isolates from sugarcane, morphologically similar to L. xyli subsp. xyli, were putatively shown to be an undescribed Agrococcus species of unknown pathogenicity. The lack of genetic variation among L. xyli subsp. xyli isolates, independent of time of sampling, cultivar of isolation, or country of origin, suggests the worldwide spread of a single pathogenic clone, and further suggests that sugarcane cultivars resistant to ratoon stunting disease in one area should retain this property in other regions.

Significant patterns of population genetic structure and limited gene flow in a threatened macropodid marsupial despite continuous habitat in southeast Queensland, Australia

Many endangered species worldwide are found in remnant populations, often within fragmented landscapes. However, when possible, an understanding of the natural extent of population structure and dispersal behaviour of threatened species would assist in their conservation and management. The brush-tailed rock-wallaby (Petrogale penicillata), a once abundant and widespread rock-wallaby species across southeastern Australia, has become nearly extinct across much of the southern part of its range. However, the northern part of the species' range still sustains many small colonies closely distributed across suitable habitat, providing a rare opportunity to investigate the natural population dynamics of a listed threatened species. We used 12 microsatellite markers to investigate genetic diversity, population structure and gene flow among brush-tailed rock-wallaby colonies within and among two valley regions with continuous habitat in southeast Queensland. We documented high and signifcant levels of population genetic structure between rock-wallaby colonies embedded in continuous escarpment habitat and forest. We found a strong and significant pattern of isolation-by-distance among colonies indicating restricted gene flow over a small geographic scale (< 10 km) and conclude that gene flow is more likely limited by intrinsic factors rather than environmental factors. In addition, we provide evidence that genetic diversity was significantly lower in colonies located in a more isolated valley region compared to colonies located in a valley region surrounded by continuous habitat. These findings shed light on the processes that have resulted in the endangered status of rock-wallaby species in Australia and they have strong implications for the conservation and management of both the remaining 'connected' brush-tailed rock-wallaby colonies in the northern parts of the species' range and the remnant endangered populations in the south.

Substantial changes in the genetic basis of tadpole morphology of Rana lessonae in the presence of predators

Predator-induced morphological plasticity is a model system for investigating phenotypic plasticity in an ecological context. We investigated the genetic basis of the predator-induced plasticity in Rana lessonae by determining the pattern of genetic covariation of three morphological traits that were found to be induced in a predatory environment. Body size decreased and tail dimensions increased when reared in the presence of preying dragonfly larvae. Genetic variance in body size increased by almost an order of magnitude in the predator environment, and the first genetic principal component was found to be highly significantly different between the two environments. The across environment genetic correlation for body size was significantly below 1 indicating that different genes contributed to this trait in the two environments. Body size may therefore be able to respond to selection independently in the two environments to some extent.

Molecular markers and human history: A tale of two haplogroup systems

To chronicle demographic movement across African Asian corridors, a variety of molecular (sequence analysis, restriction mapping and denaturing high performance liquid chromatography etc.) and statistical (correspondence analysis, AMOVA, calculation of diversity indices and phylogenetic inference, etc.) techniques were employed to assess the phylogeographic patterns of mtDNA control region and Y chromosomal variation among 14 sub-Saharan, North African and Middle Eastern populations. The patterns of genetic diversity revealed evidence of multiple migrations across several African Asian passageways as well within the African continent itself. The two-part analysis uncovered several interesting results which include the following: (1) a north (Egypt and Middle East Asia) to south (sub-Saharan Africa) partitioning of both mtDNA and Y chromosomal haplogroup diversity, (2) a genetic diversity gradient in sub-Saharan Africa from east to west, (3) evidence in favor of the Levantine Corridor over the Horn of Africa as the major genetic conduit since the Last Glacial Maximum, (4) a substantially higher mtDNA versus Y chromosomal sub-Saharan component in the Middle East collections, (5) a higher representation of East versus West African mtDNA haplotypes in the Arabian Peninsula populations versus no such bias in the Levant groups and lastly, (6) genetic remnants of the Bantu demographic expansion in sub-Saharan Africa. ^

High natural gene expression variation in the reef-building coral Acropora millepora: potential for acclimative and adaptive plasticity

Background: Ecosystems worldwide are suffering the consequences of anthropogenic impact. The diverse ecosystem of coral reefs, for example, are globally threatened by increases in sea surface temperatures due to global warming. Studies to date have focused on determining genetic diversity, the sequence variability of genes in a species, as a proxy to estimate and predict the potential adaptive response of coral populations to environmental changes linked to climate changes. However, the examination of natural gene expression variation has received less attention. This variation has been implicated as an important factor in evolutionary processes, upon which natural selection can act. Results: We acclimatized coral nubbins from six colonies of the reef-building coral Acropora millepora to a common garden in Heron Island (Great Barrier Reef, GBR) for a period of four weeks to remove any site-specific environmental effects on the physiology of the coral nubbins. By using a cDNA microarray platform, we detected a high level of gene expression variation, with 17% (488) of the unigenes differentially expressed across coral nubbins of the six colonies (jsFDR-corrected, p < 0.01). Among the main categories of biological processes found differentially expressed were transport, translation, response to stimulus, oxidation-reduction processes, and apoptosis. We found that the transcriptional profiles did not correspond to the genotype of the colony characterized using either an intron of the carbonic anhydrase gene or microsatellite loci markers. Conclusion: Our results provide evidence of the high inter-colony variation in A. millepora at the transcriptomic level grown under a common garden and without a correspondence with genotypic identity. This finding brings to our attention the importance of taking into account natural variation between reef corals when assessing experimental gene expression differences. The high transcriptional variation detected in this study is interpreted and discussed within the context of adaptive potential and phenotypic plasticity of reef corals. Whether this variation will allow coral reefs to survive to current challenges remains unknown.

Functional Evaluation of Causal Mutations Identified in Human Genetic Studies

Human genetics has been experiencing a wave of genetic discoveries thanks to the development of several technologies, such as genome-wide association studies (GWAS), whole-exome sequencing, and whole genome sequencing. Despite the massive genetic discoveries of new variants associated with human diseases, several key challenges emerge following the genetic discovery. GWAS is known to be good at identifying the locus associated with the patient phenotype. However, the actually causal variants responsible for the phenotype are often elusive. Another challenge in human genetics is that even the causal mutations are already known, the underlying biological effect might remain largely ambiguous. Functional evaluation plays a key role to solve these key challenges in human genetics both to identify causal variants responsible for the phenotype, and to further develop the biological insights from the disease-causing mutations.

We adopted various methods to characterize the effects of variants identified in human genetic studies, including patient genetic and phenotypic data, RNA chemistry, molecular biology, virology, and multi-electrode array and primary neuronal culture systems. Chapter 1 is a broader introduction for the motivation and challenges for functional evaluation in human genetic studies, and the background of several genetics discoveries, such as hepatitis C treatment response, in which we performed functional characterization.

Chapter 2 focuses on the characterization of causal variants following the GWAS study for hepatitis C treatment response. We characterized a non-coding SNP (rs4803217) of IL28B (IFNL3) in high linkage disequilibrium (LD) with the discovery SNP identified in the GWAS. In this chapter, we used inter-disciplinary approaches to characterize rs4803217 on RNA structure, disease association, and protein translation.

Chapter 3 describes another avenue of functional characterization following GWAS focusing on the novel transcripts and proteins identified near the IL28B (IFNL3) locus. It has been recently speculated that this novel protein, which was named IFNL4, may affect the HCV treatment response and clearance. In this chapter, we used molecular biology, virology, and patient genetic and phenotypic data to further characterize and understand the biology of IFNL4. The efforts in chapter 2 and 3 provided new insights to the candidate causal variant(s) responsible for the GWAS for HCV treatment response, however, more evidence is still required to make claims for the exact causal roles of these variants for the GWAS association.

Chapter 4 aims to characterize a mutation already known to cause a disease (seizure) in a mouse model. We demonstrate the potential use of multi-electrode array (MEA) system for the functional characterization and drug testing on mutations found in neurological diseases, such as seizure. Functional characterization in neurological diseases is relatively challenging and available systematic tools are relatively limited. This chapter shows an exploratory research and example to establish a system for the broader use for functional characterization and translational opportunities for mutations found in neurological diseases.

Overall, this dissertation spans a range of challenges of functional evaluations in human genetics. It is expected that the functional characterization to understand human mutations will become more central in human genetics, because there are still many biological questions remaining to be answered after the explosion of human genetic discoveries. The recent advance in several technologies, including genome editing and pluripotent stem cells, is also expected to make new tools available for functional studies in human diseases.

Dynamique de la variation génétique et épigénétique chez un vertébré kleptogène

La variation phénotypique est essentielle à la persistance des organismes dans le temps ainsi qu’à la colonisation de nouveaux habitats. Les principales sources de variation phénotypique sont la génétique et l'épigénétique. L'épigénétique a été proposé comme un atout important pour les organismes asexués pour compenser le manque de diversité génétique. L'objectif de cette étude est d'évaluer si l’absence de variation génétique est compensée par l'épigénétique en comparant les profils de méthylation d’individus gynogènes et kleptogènes des hybrides de salamandre à points bleus. Les individus échantillonnés s’organisent en cinq groupes génétiquement différenciés, provenant du même haplome paternel A. jeffersonianum. Deux des cinq groupes sont exclusivement gynogènes, pour des raisons écologiques ou génomiques. Les trois autres groupes sont formés d’individus parfois kleptogènes, car ils présentent une variation génétique plus élevée au sein d’un site qu’entre les sites, en plus de porter des allèles très divergents par rapport à la distribution globale des allèles hybrides, trouvés en haute fréquence dans les populations sympatriques de A. laterale. Les patrons épigénétiques sont variables et distincts entre les cinq groupes génétiques. Les groupes gynogènes sont les seuls à présenter un effet environnemental significatif sur leurs patrons épigénétiques, suggérant que ces individus clonaux doivent être en mesure de maximiser leur potentiel de variation épigénétique pour faire face à des variations environnementales.