Patterns and effects of gene flow on adaptation across spatial scales: implications for management

Gene flow can have rapid effects on adaptation and is an important evolutionary tool available when undertaking biological conservation and restoration. This tool is underused partly because of the perceived risk of outbreeding depression and loss of mean fitness when different populations are crossed. In this article, we briefly review some theory and empirical findings on how genetic variation is distributed across species ranges, describe known patterns of gene flow in nature with respect to environmental gradients, and highlight the effects of gene flow on adaptation in small or stressed populations in challenging environments (e.g., at species range limits). We then present a case study involving crosses at varying spatial scales among mountain populations of a trigger plant (Stylidium armeria: Stylidiaceae) in the Australian Alps to highlight how some issues around gene flow effects can be evaluated. We found evidence of outbreeding depression in seed production at greater geographic distances. Nevertheless, we found no evidence of maladaptive gene flow effects in likelihood of germination, plant performance (size), and performance variance, suggesting that gene flow at all spatial scales produces offspring with high adaptive potential. This case study demonstrates a path to evaluating how increasing sources of gene flow in managed wild and restored populations could identify some offspring with high fitness that could bolster the ability of populations to adapt to future environmental changes. We suggest further ways in which managers and researchers can act to understand and consider adaptive gene flow in natural and conservation contexts under rapidly changing conditions.

Genetic variation and structure of endemic and endangered wild celery (Kelussia odoratissima Mozaff.) quantified using novel microsatellite markers developed by next-generation sequencing

Kelussia odoratissima Mozaff. (Apiaceae) is a native plant that has been traditionally consumed in Iran's food and pharmaceutical industries. Overharvesting of the taxon, especially at the beginning of the growing season, due to its considerable medicinal and economic value, is believed to be the main reason for the extirpating of this plant. The consequences of the severe anthropogenic impacts on the genetic diversity of populations are poorly known. In order to investigate the level of genetic variation and patterns of the genetic structure of K. odoratissima, we developed novel microsatellite markers using the 454 Roche next-generation sequencing (NGS) platform for the first time. Out of 1,165 microsatellite markers bioinformatically confirmed, twenty-five were tested, of which 23 were used to screen genetic variation across 12 natural populations. Our results showed that the average number of alleles per locus and the polymorphic information content (PIC) were 10.87 (range 7 to 27), and 0.81 (range 0.67 to 0.94), respectively. The mean observed and expected heterozygosities (± SD) across all populations were 0.80 ± 0.31 and 0.72 ± 0.14, respectively. The average pairwise FST among the populations was 0.37 (range 0.04 to 0.81). Bayesian and distance-based clustering, and principal coordinate analyses revealed at least four major genetic clusters. Although high level of structure can be explained by landscape topography and geographic distance, presence of admixed populations can be associated to seed or pollen dispersal. Contrary to expectations, the high level of genetic variation and lack of inbreeding suggest that overexploitation has not yet significantly purged the allelic variability within the natural populations in protected areas.

Genome survey and genetic characterization of Acacia pachyceras O. Schwartz

Acacia pachyceras O. Schwartz (Leguminoseae), a woody tree growing in Kuwait is critically endangered. High throughput genomic research is immediately needed to formulate effective conservation strategies for its rehabilitation. We therefore, performed a genome survey analysis of the species. Whole genome sequencing generated ~97 Gb of raw reads (92x coverage) with a per base quality score above Q30. The k-mer analysis (17 mer) revealed its genome to be 720Mb in size with an average guanine-cytosine (GC) ratio of 35%. The assembled genome was analyzed for repeat regions (45.4%-interspersed repeats; 9%-retroelements; 2%-DNA transposons). BUSCO assessment of completeness of genome identified 93% of assembly to be complete. Gene alignments in BRAKER2 yielded 34,374 transcripts corresponding to 33,650 genes. Average length of coding sequences and protein sequences were recorded as 1,027nts and 342aa, respectively. GMATA software filtered a total of 901,755 simple sequence repeats (SSRs) regions against which 11,181 unique primers were designed. A subset of 110 SSR primers were PCR validated and demonstrated for its application in genetic diversity analysis of Acacia. The SSR primers successfully amplified A. gerrardii seedlings DNA depicting cross transferability among species. The principal coordinate analysis and the split decomposition tree (bootstrapping runs of 1000 replicates) distributed the Acacia genotypes into two clusters. The flow cytometry analysis revealed the A. pachyceras genome to be polyploid (6x). The DNA content was predicted as 2.46 pg, 1.23 pg, and 0.41 pg corresponding to 2C DNA, 1C DNA and 1Cx DNA, respectively. The results provide a base for further high throughput genomic studies and molecular breeding for its conservation.

A genome resource for Acacia, Australia's largest plant genus

Acacia (Leguminosae, Caesalpinioideae, mimosoid clade) is the largest and most widespread genus of plants in the Australian flora, occupying and dominating a diverse range of environments, with an equally diverse range of forms. For a genus of its size and importance, Acacia currently has surprisingly few genomic resources. Acacia pycnantha, the golden wattle, is a woody shrub or tree occurring in south-eastern Australia and is the country's floral emblem. To assemble a genome for A. pycnantha, we generated long-read sequences using Oxford Nanopore Technology, 10x Genomics Chromium linked reads, and short-read Illumina sequences, and produced an assembly spanning 814 Mb, with a scaffold N50 of 2.8 Mb, and 98.3% of complete Embryophyta BUSCOs. Genome annotation predicted 47,624 protein-coding genes, with 62.3% of the genome predicted to comprise transposable elements. Evolutionary analyses indicated a shared genome duplication event in the Caesalpinioideae, and conflict in the relationships between Cercis (subfamily Cercidoideae) and subfamilies Caesalpinioideae and Papilionoideae (pea-flowered legumes). Comparative genomics identified a suite of expanded and contracted gene families in A. pycnantha, and these were annotated with both GO terms and KEGG functional categories. One expanded gene family of particular interest is involved in flowering time and may be associated with the characteristic synchronous flowering of Acacia. This genome assembly and annotation will be a valuable resource for all studies involving Acacia, including the evolution, conservation, breeding, invasiveness, and physiology of the genus, and for comparative studies of legumes.

A range-wide analysis of population structure and genomic variation within the critically endangered spiny daisy (Acanthocladium dockeri)

Understanding population structure and genetic diversity is important for designing effective conservation strategies. As a critically endangered shrub, the six remaining extant populations of spiny daisy (Acanthocladium dockeri) are restricted to country roadsides in the mid-north of South Australia, where the species faces many ongoing abiotic and biotic threats to survival. Currently the spiny daisy is managed by selecting individuals from the extant populations and translocating them to establish insurance populations. However, there is little information available on the genetic differentiation between populations and diversity within source populations, which are essential components of planning translocations. To help fill this knowledge gap, we analysed population structure within and among all six of its known wild populations using 7,742 SNPs generated by a genotyping-by-sequencing approach. Results indicated that each population was strongly differentiated, had low levels of genetic diversity, and there was no evidence of inter-population gene flow. Individuals within each population were generally closely related, however, the Melrose population consisted entirely of clones. Our results suggest genetic rescue should be applied to wild spiny daisy populations to increase genetic diversity that will subsequently lead to greater intra-population fitness and adaptability. As a starting point, we suggest focussing on improving seed viability via inter-population crosses such as through hand pollination experiments to experimentally assess their sexual compatibility with the hope of increasing spiny daisy sexual reproduction and long-term reproductive fitness.

Comparison and Optimization of DNA Isolation Protocols for High Throughput Genomic Studies of Acacia pachyceras Schwartz

We describe the optimization and validation of six DNA isolation protocols from fresh leaves of the rare tree Acacia pachyceras. The first four protocols employed three commercial kits (Sigma, Nucleospin1, Nucleospin 2, Promega) whereas the remaining two were based on the traditional sodium dodecyl sulfate (SDS) and cetyltrimethylammonium bromide CTAB methods. Each protocol provided significantly different results concerning DNA concentration (p < 0.032), yield (p < 0.000), contaminant carry over, protocol duration, cost per sample, and comprehensive cost. We demonstrated the applicability of all the tested protocols in DNA barcoding. The protocol yielded maximum amounts (92.85 µg) of DNA in a rapid turnaround time (8 h). The quantity and purity surpassed all the other tested methods. DNA extracted by the CTAB method was the best for NGS (Phred score >Q30). These protocols will be useful tools for molecular research of Acacia pachyceras and other closely related tree species.

Draft genome sequence and SSR mining data of Acacia pachyceras Schwartz

Acacia tree population is declining in several countries of the world especially in the Arabian peninsula due to human-induced activities. The tree has potential medicinal and economic benefits as a source of fuel and timber. It can fix nitrogen, a significant property that assists in desert rehabilitation. However, the lack of genomic information of Acacia pachyceras hampers its genetic study and breeding process. We performed paired-end sequencing of A. pachyceras at a depth of 120X to obtain raw sequences of 108.9 GB with a per base quality >Q30. Filtered raw data was assembled into a fasta file of 4 GB. The assembled genomic sequences consisted of 901,755 single sequence repeats (SSRs). In total 11,596 primer pairs were designed against these SSR motifs. The data generated provides baseline genomic information about the species and formulates a base for further sequencing of A. pachyceras through PACBio and HiC technologies. The novel developed SSR markers will facilitate genetic diversity and conservation studies for Acacia species.

Genetic diversity and population structure of six South African Acacia mearnsii breeding populations based on SSR markers

Black wattle (Acacia mearnsii) has great economic value as a commercial source of tannins, timber and a source of firewood for local and international markets. It has been suggested that to maximize the genetic gain of A. mearnsii plantations in South Africa, the gene pool that exist within ICFR needs to be broadened via introduction of new genotypes with diverse traits. In this work, 282 A. mearnsii samples sourced from the ICFR breeding program were genotyped using 11 cross-species SSR markers. Our results showed low to moderate genetic differentiation (F_ST) among the six breeding subpopulations, with positive inbreeding (F_IS) values that could be attributed to an historical inbreeding event. Low levels of relatedness could however indicate some mechanism of inbreeding avoidance. The effects from a recent supplementation of genetic material from two native Australian populations were observed through genetic structuring analyses. Analysis of molecular variance (AMOVA) revealed that significant genetic variation was mainly distributed within populations (75%) and among individuals (23%). The results provide significant information on A. mearnsii population genetic diversity and structure, which can be used for conservation of the current subpopulations and future tree improvement programs.

Genomic, Habitat, and Leaf Shape Analyses Reveal a Possible Cryptic Species and Vulnerability to Climate Change in a Threatened Daisy

Olearia pannosa is a plant species listed as vulnerable in Australia. Two subspecies are currently recognised (O. pannosa subsp. pannosa (silver daisy) and O. pannosa subsp. cardiophylla (velvet daisy)), which have overlapping ranges but distinct leaf shape. Remnant populations face threats from habitat fragmentation and climate change. We analysed range-wide genomic data and leaf shape variation to assess population diversity and divergence and to inform conservation management strategies. We detected three distinct genetic groupings and a likely cryptic species. Samples identified as O. pannosa subsp. cardiophylla from the Flinders Ranges in South Australia were genetically distinct from all other samples and likely form a separate, range-restricted species. Remaining samples formed two genetic clusters, which aligned with leaf shape differences but not fully with current subspecies classifications. Levels of genetic diversity and inbreeding differed between the three genetic groups, suggesting each requires a separate management strategy. Additionally, we tested for associations between genetic and environmental variation and carried out habitat suitability modelling for O. pannosa subsp. pannosa populations. We found mean annual maximum temperature explained a significant proportion of genomic variance. Habitat suitability modelling identified mean summer maximum temperature, precipitation seasonality and mean annual rainfall as constraints on the distribution of O. pannosa subsp. pannosa, highlighting increasing aridity as a threat for populations located near suitability thresholds. Our results suggest maximum temperature is an important agent of selection on O. pannosa subsp. pannosa and should be considered in conservation strategies. We recommend taxonomic revision of O. pannosa and provide conservation management recommendations.