DNA sequence variation and selection of tag single-nucleotide polymorphisms at candidate genes for drought-stress response in Pinus taeda L

Colonization of root endophytic fungus Serendipita indica improves drought tolerance of Pinus taeda seedlings by regulating metabolome and proteome

Pinus taeda is an important forest tree species for plantations because of its rapid growth and high yield of oleoresins. Although P. taeda plantations distribute in warm and wet southern China, drought, sometime serious and long time, often occurs in the region. To explore drought tolerance of P. taeda and usage of beneficial microorganisms, P. taeda seedlings were planted in pots and were inoculated with root endophytic fungus Serendipita indica and finally were treated with drought stress for 53 d. Metabolome and proteome of their needles were analyzed. The results showed that S. indica inoculation of P. taeda seedlings under drought stress caused great changes in levels of some metabolites in their needles, especially some flavonoids and organic acids. Among them, the levels of eriocitrin, trans-aconitic acid, vitamin C, uric acid, alpha-ketoglutaric acid, vitamin A, stachydrine, coumalic acid, itaconic acid, calceolarioside B, 2-oxoglutaric acid, and citric acid were upregulated more than three times in inoculated seedlings under drought stress, compared to those of non-inoculated seedlings under drought stress. KEGG analysis showed that some pathways were enriched in inoculated seedlings under drought stress, such as flavonoid biosynthesis, ascorbate and aldarate metabolism, C5-branched dibasic acid metabolism. Proteome analysis revealed some specific differential proteins. Two proteins, namely, H9X056 and H9VDW5, only appeared in the needles of inoculated seedlings under drought stress. The protein H9VNE7 was upregulated more than 11.0 times as that of non-inoculated seedlings under drought stress. In addition, S. indica inoculation increased enrichment of water deficient-inducible proteins (such as LP3-1, LP3-2, LP3-3, and dehydrins) and those involved in ribosomal structures (such as A0A385JF23). Meanwhile, under drought stress, the inoculation caused great changes in biosynthesis and metabolism pathways, mainly including phenylpropanoid biosynthesis, cutin, suberine and wax biosynthesis, and 2-oxocarboxylic acid metabolism. In addition, there were positive relationships between accumulation of some metabolites and enrichment of proteins in P. taeda under drought stress. Altogether, our results showed great changes in metabolome and proteome in inoculated seedlings under drought stress and provided a guideline to further study functions of metabolites and proteins, especially those related to drought stress.

NLR diversity and candidate fusiform rust resistance genes in loblolly pine

Resistance to fusiform rust disease in loblolly pine (Pinus taeda) is a classic gene-for-gene system. Early resistance gene mapping in the P. taeda family 10-5 identified RAPD markers for a major fusiform rust resistance gene, Fr1. More recently, single nucleotide polymorphism (SNP) markers associated with resistance were mapped to a full-length gene model in the loblolly pine genome encoding for a nucleotide-binding site leucine-rich repeat (NLR) protein. NLR genes are one of the most abundant gene families in plant genomes and are involved in effector-triggered immunity. Inter- and intraspecies studies of NLR gene diversity and expression have resulted in improved disease resistance. To characterize NLR gene diversity and discover potential resistance genes, we assembled de novo transcriptomes from 92 loblolly genotypes from across the natural range of the species. In these transcriptomes, we identified novel NLR transcripts that are not present in the loblolly pine reference genome and found significant geographic diversity of NLR genes providing evidence of gene family evolution. We designed capture probes for these NLRs to identify and map SNPs that stably cosegregate with resistance to the SC20-21 isolate of Cronartium quercuum f.sp. fusiforme (Cqf) in half-sib progeny of the 10-5 family. We identified 10 SNPs and 2 quantitative trait loci associated with resistance to SC20-21 Cqf. The geographic diversity of NLR genes provides evidence of NLR gene family evolution in loblolly pine. The SNPs associated with rust resistance provide a resource to enhance breeding and deployment of resistant pine seedlings.

Resolving the origins of invertebrate colonists in the Yangtze River Estuary with molecular markers: Implications for ecological connectivity

Understanding connectivity over different spatial and temporal scales is fundamental for managing of ecological systems. However, controversy exists for wintertime ecological connectivity between the Yangtze River Estuary (YRE) and inner southwestern Yellow Sea. Here, we investigated ecological connectivity between the YRE and inner southwestern Yellow Sea in wintertime by precisely pinpointing the source of the newly colonized populations of a winter-spawning rocky intertidal invertebrate, Littorina brevicula (Philippi, 1844), on artificial structures along the coast of the Yangtze River Delta (YRD) using mitochondrial ND6 sequences and microsatellite data. Clear phylogeographic and genetic differentiation were detected between natural rocky populations south and north of the YRE, which resulted from the lack of hard substrate for rocky invertebrates in the large YRD coast. For the newly colonized populations on the coast of YRD, most individuals (98%) to the south of ~33.5°N were from natural rocky populations to the south of the YRE and most of those (94%) to the north of ~33.5°N were from the northern natural rocky populations, which demonstrated strong ecological connectivity between the inner southwestern Yellow Sea and the YRE in winter time. We presented the first genetic evidence that demonstrated a northward wintertime coastal current in the inner southwestern Yellow Sea, and precisely illustrated the boundary of the coastal current recently proposed by numerical experiment. These results indicated that the YRE serves as an important source of materials and energy for the inner southwestern Yellow Sea in winter, which can be crucial for the function of the Yellow Sea ecosystem.

Gene Frequency Shift in Relict Abies pinsapo Forests Associated with Drought-Induced Mortality: Preliminary Evidence of Local-Scale Divergent Selection

Current climate change constitutes a challenge for the survival of several drought-sensitive forests. The study of the genetic basis of adaptation offers a suitable way to understand how tree species may respond to future climatic conditions, as well as to design suitable conservation and management strategies. Here, we focus on selected genetic signatures of the drought-sensitive relict fir, Abies pinsapo Boiss. Field sampling of 156 individuals was performed in two elevation ecotones, characterized by widespread A. pinsapo decline and mortality. The DNA from dead trees was investigated and compared to living individuals, accounting for different ages and elevations. We studied the genes gated outwardly-rectifying K+ (GORK) channel and Plasma membrane Intrinsic Protein (PIP1) aquaporin, previously related to drought response in plant model species, to test whether drought was the main abiotic factor driving the decline of A. pinsapo forests. A combination of linear regression and factor models were used to test these selection signatures, as well as a fixation index (Fst), used here to analyze the genetic structure. The results were consistent among these approaches, supporting a statistically significant association of the GORK gene with survival in one of the A. pinsapo populations. These results provide preliminary evidence for the potential role of the GORK gene in the resilience to drought of A. pinsapo.

Molecular approaches reveal weak sibship aggregation and a high dispersal propensity in a non-native fish parasite

Inferring parameters related to the aggregation pattern of parasites and to their dispersal propensity are important for predicting their ecological consequences and evolutionary potential. Nonetheless, it is notoriously difficult to infer these parameters from wildlife parasites given the difficulty in tracking these organisms. Molecular-based inferences constitute a promising approach that has yet rarely been applied in the wild. Here, we combined several population genetic analyses including sibship reconstruction to document the genetic structure, patterns of sibship aggregation, and the dispersal dynamics of a non-native parasite of fish, the freshwater copepod ectoparasite Tracheliastes polycolpus. We collected parasites according to a hierarchical sampling design, with the sampling of all parasites from all host individuals captured in eight sites spread along an upstream-downstream river gradient. Individual multilocus genotypes were obtained from 14 microsatellite markers, and used to assign parasites to full-sib families and to investigate the genetic structure of T. polycolpus among both hosts and sampling sites. The distribution of full-sibs obtained among the sampling sites was used to estimate individual dispersal distances within families. Our results showed that T. polycolpus sibs tend to be aggregated within sites but not within host individuals. We detected important upstream-to-downstream dispersal events of T. polycolpus between sites (modal distance: 25.4 km; 95% CI [22.9, 27.7]), becoming scarcer as the geographic distance from their family core location increases. Such a dispersal pattern likely contributes to the strong isolation-by-distance observed at the river scale. We also detected some downstream-to-upstream dispersal events (modal distance: 2.6 km; 95% CI [2.2-23.3]) that likely result from movements of infected hosts. Within each site, the dispersal of free-living infective larvae among hosts likely contributes to increasing genetic diversity on hosts, possibly fostering the evolutionary potential of T. polycolpus.

The GenTree Platform: growth traits and tree-level environmental data in 12 European forest tree species

BACKGROUND

Progress in the field of evolutionary forest ecology has been hampered by the huge challenge of phenotyping trees across their ranges in their natural environments, and the limitation in high-resolution environmental information.

FINDINGS

The GenTree Platform contains phenotypic and environmental data from 4,959 trees from 12 ecologically and economically important European forest tree species: Abies alba Mill. (silver fir), Betula pendula Roth. (silver birch), Fagus sylvatica L. (European beech), Picea abies (L.) H. Karst (Norway spruce), Pinus cembra L. (Swiss stone pine), Pinus halepensis Mill. (Aleppo pine), Pinus nigra Arnold (European black pine), Pinus pinaster Aiton (maritime pine), Pinus sylvestris L. (Scots pine), Populus nigra L. (European black poplar), Taxus baccata L. (English yew), and Quercus petraea (Matt.) Liebl. (sessile oak). Phenotypic (height, diameter at breast height, crown size, bark thickness, biomass, straightness, forking, branch angle, fructification), regeneration, environmental in situ measurements (soil depth, vegetation cover, competition indices), and environmental modeling data extracted by using bilinear interpolation accounting for surrounding conditions of each tree (precipitation, temperature, insolation, drought indices) were obtained from trees in 194 sites covering the species' geographic ranges and reflecting local environmental gradients.

CONCLUSION

The GenTree Platform is a new resource for investigating ecological and evolutionary processes in forest trees. The coherent phenotyping and environmental characterization across 12 species in their European ranges allow for a wide range of analyses from forest ecologists, conservationists, and macro-ecologists. Also, the data here presented can be linked to the GenTree Dendroecological collection, the GenTree Leaf Trait collection, and the GenTree Genomic collection presented elsewhere, which together build the largest evolutionary forest ecology data collection available.

Imprints of selection in peripheral and ecologically marginal central-eastern European Scots pine populations

Knowledge of the molecular mechanisms underlying the stress response in plants is essential to understand evolutionary processes that result in long-term persistence of populations. Populations inhabiting marginal ecological conditions at the distribution range periphery may have preserved imprints of natural selection that have shaped functional genetic variation of the species. Our aim was to evaluate the extent of selection processes in the extremely fragmented, peripheral and isolated populations of Scots pine in central-eastern Europe. Autochthonous populations of the Carpathian Mts. and the Pannonian Basin were sampled and drought stress-related candidate genes were re-sequenced. Neutrality tests and outlier detection approaches were applied to infer the effect and direction of selection. Populations retained high genetic diversity by preserving a high number of alleles and haplotypes, many of them being population specific. Neutrality tests and outlier detection highlighted nucleotide positions that are under divergent selection and may be involved in local adaptation. The detected genetic pattern confirms that natural selection has played an important role in shaping modern-day genetic variation in marginal Scots pine populations, allowing for the long-term persistence of populations. Selection detected at functional regions possibly acts to maintain diversity and counteract the effect of genetic erosion.

Spatial genetic structure and landscape connectivity in black bears: Investigating the significance of using different land cover datasets and classifications in landscape genetics analyses

Landscape genetic analyses allow detection of fine-scale spatial genetic structure (SGS) and quantification of effects of landscape features on gene flow and connectivity. Typically, analyses require generation of resistance surfaces. These surfaces characteristically take the form of a grid with cells that are coded to represent the degree to which landscape or environmental features promote or inhibit animal movement. How accurately resistance surfaces predict association between the landscape and movement is determined in large part by (a) the landscape features used, (b) the resistance values assigned to features, and (c) how accurately resistance surfaces represent landscape permeability. Our objective was to evaluate the performance of resistance surfaces generated using two publicly available land cover datasets that varied in how accurately they represent the actual landscape. We genotyped 365 individuals from a large black bear population (Ursus americanus) in the Northern Lower Peninsula (NLP) of Michigan, USA at 12 microsatellite loci, and evaluated the relationship between gene flow and landscape features using two different land cover datasets. We investigated the relative importance of land cover classification and accuracy on landscape resistance model performance. We detected local spatial genetic structure in Michigan's NLP black bears and found roads and land cover were significantly correlated with genetic distance. We observed similarities in model performance when different land cover datasets were used despite 21% dissimilarity in classification between the two land cover datasets. However, we did find the performance of land cover models to predict genetic distance was dependent on the way the land cover was defined. Models in which land cover was finely defined (i.e., eight land cover classes) outperformed models where land cover was defined more coarsely (i.e., habitat/non-habitat or forest/non-forest). Our results show that landscape genetic researchers should carefully consider how land cover classification changes inference in landscape genetic studies.

Heat stress responses and population genetics of the kelp Laminaria digitata (Phaeophyceae) across latitudes reveal differentiation among North Atlantic populations

To understand the thermal plasticity of a coastal foundation species across its latitudinal distribution, we assess physiological responses to high temperature stress in the kelp Laminaria digitata in combination with population genetic characteristics and relate heat resilience to genetic features and phylogeography. We hypothesize that populations from Arctic and cold-temperate locations are less heat resilient than populations from warm distributional edges. Using meristems of natural L. digitata populations from six locations ranging between Kongsfjorden, Spitsbergen (79°N), and Quiberon, France (47°N), we performed a common-garden heat stress experiment applying 15°C to 23°C over eight days. We assessed growth, photosynthetic quantum yield, carbon and nitrogen storage, and xanthophyll pigment contents as response traits. Population connectivity and genetic diversity were analyzed with microsatellite markers. Results from the heat stress experiment suggest that the upper temperature limit of L. digitata is nearly identical across its distribution range, but subtle differences in growth and stress responses were revealed for three populations from the species' ecological range margins. Two populations at the species' warm distribution limit showed higher temperature tolerance compared to other populations in growth at 19°C and recovery from 21°C (Quiberon, France), and photosynthetic quantum yield and xanthophyll pigment responses at 23°C (Helgoland, Germany). In L. digitata from the northernmost population (Spitsbergen, Norway), quantum yield indicated the highest heat sensitivity. Microsatellite genotyping revealed all sampled populations to be genetically distinct, with a strong hierarchical structure between southern and northern clades. Genetic diversity was lowest in the isolated population of the North Sea island of Helgoland and highest in Roscoff in the English Channel. All together, these results support the hypothesis of moderate local differentiation across L. digitata's European distribution, whereas effects are likely too weak to ameliorate the species' capacity to withstand ocean warming and marine heatwaves at the southern range edge.

Lack of conspicuous sex-biased dispersal patterns at different spatial scales in an Asian endemic goose species breeding in unpredictable steppe wetlands

Dispersal affects the spatial distribution and population structure of species. Dispersal is often male-biased in mammals while female-biased in birds, with the notable exception of the Anatidae. In this study, we tested genetic evidence for sex-biased dispersal (SBD) in the Swan Goose Anser cygnoides, an Asian endemic and IUCN vulnerable species, which has been increasingly restricted to breeding on Mongolian steppe wetlands. We analyzed the genotypes of 278 Swan Geese samples from 14 locations at 14 microsatellite loci. Results from assignment indices, analysis of molecular variance, and five other population descriptors all failed to support significant SBD signals for the Swan Goose at the landscape level. Although overall results showed significantly high relatedness within colonies (suggesting high levels of philopatry in both sexes), local male genetic structure at the 1,050 km distance indicated greater dispersal distance for females from the eastern sector of the breeding range. Hence, local dispersal is likely scale-dependent and female-biased within the eastern breeding range. These findings are intriguing considering the prevailing expectation for there to be female fidelity in most goose species. We suggest that while behavior-related traits may have facilitated the local genetic structure for the Swan Goose, several extrinsic factors, including the decreasing availability of the nesting sites and the severe fragmentation of breeding habitats, could have contributed to the absence of SBD at the landscape level. The long-distance molt migration that is typical of goose species such as the Swan Goose may also have hampered our ability to detect SBD. Hence, we urge further genetic sampling from other areas in summer to extend our results, complemented by field observations to confirm our DNA analysis conclusions about sex-specific dispersal patterns at different spatial scales in this species.

Environmental Genome-Wide Association Reveals Climate Adaptation Is Shaped by Subtle to Moderate Allele Frequency Shifts in Loblolly Pine

Understanding the genomic basis of local adaptation is crucial to determine the potential of long-lived woody species to withstand changes in their natural environment. In the past, efforts to dissect the genomic architecture in gymnosperms species have been limited due to the absence of reference genomes. Recently, the genomes of some commercially important conifers, such as loblolly pine, have become available, allowing whole-genome studies of these species. In this study, we test for associations between 87k SNPs, obtained from whole-genome resequencing of loblolly pine individuals, and 270 environmental variables and combinations of them. We determine the geographic location of significant loci and identify their genomic location using our newly constructed ultradense 26k SNP linkage map. We found that water availability is the main climatic variable shaping local adaptation of the species, and found 821 SNPs showing significant associations with climatic variables or combinations of them based on the consistent results of three different genotype–environment association methods. Our results suggest that adaptation to climate in the species might have occurred by many changes in the frequency of alleles with moderate to small effect sizes, and by the smaller contribution of large effect alleles in genes related to moisture deficit, temperature and precipitation. Genomic regions of low recombination and high population differentiation harbored SNPs associated with groups of environmental variables, suggesting climate adaptation might have evolved as a result of different selection pressures acting on groups of genes associated with an aspect of climate rather than on individual environmental variables.

Development of 32 novel microsatellite loci in Juglans sigillata using genomic data

PREMISE

A novel set of microsatellite markers was developed for Juglans sigillata (Juglandaceae), an endemic walnut species in southwestern China, to facilitate cultivar identification and future investigations into the genetic structure and domestication history of this species and its close relatives.

METHODS AND RESULTS

We developed 32 microsatellite loci for J. sigillata using genomic data and used them to examine 60 individuals from three natural populations. A high level of polymorphism was detected by these primers, with up to eight alleles observed per locus, and an average of four alleles across populations. The levels of observed and expected heterozygosity ranged from 0.000-1.000 and 0.000-0.785, respectively. All but two of the loci were also successfully amplified in three closely related Eurasian Juglans species (J. regia, J. cathayensis, and J. mandshurica).

CONCLUSIONS

The microsatellite loci identified here provide a powerful resource for examining the genetic structure and domestication history of Juglans, as well as identification of its cultivars.

Genome-Wide Association Study of Wood Anatomical and Morphological Traits in Populus trichocarpa

To understand the genetic mechanisms underlying wood anatomical and morphological traits in Populus trichocarpa, we used 869 unrelated genotypes from a common garden in Clatskanie, Oregon that were previously collected from across the distribution range in western North America. Using GEMMA mixed model analysis, we tested for the association of 25 phenotypic traits and nine multitrait combinations with 6.741 million SNPs covering the entire genome. Broad-sense trait heritabilities ranged from 0.117 to 0.477. Most traits were significantly correlated with geoclimatic variables suggesting a role of climate and geography in shaping the variation of this species. Fifty-seven SNPs from single trait GWAS and 11 SNPs from multitrait GWAS passed an FDR threshold of 0.05, leading to the identification of eight and seven nearby candidate genes, respectively. The percentage of phenotypic variance explained (PVE) by the significant SNPs for both single and multitrait GWAS ranged from 0.01% to 6.18%. To further evaluate the potential roles of candidate genes, we used a multi-omic network containing five additional data sets, including leaf and wood metabolite GWAS layers and coexpression and comethylation networks. We also performed a functional enrichment analysis on coexpression nearest neighbors for each gene model identified by the wood anatomical and morphological trait GWAS analyses. Genes affecting cell wall composition and transport related genes were enriched in wood anatomy and stomatal density trait networks. Signaling and metabolism related genes were also common in networks for stomatal density. For leaf morphology traits (leaf dry and wet weight) the networks were significantly enriched for GO terms related to photosynthetic processes as well as cellular homeostasis. The identified genes provide further insights into the genetic control of these traits, which are important determinants of the suitability and sustainability of improved genotypes for lignocellulosic biofuel production.

The distribution and spread of naturally occurring Medea selfish genetic elements in the United States

Selfish genetic elements (SGEs) are DNA sequences that are transmitted to viable offspring in greater than Mendelian frequencies. Medea SGEs occur naturally in some populations of red flour beetle (Tribolium castaneum) and are expected to increase in frequency within populations and spread among populations. The large-scale U.S. distributions of Medea-4 (M4) had been mapped based on samples from 1993 to 1995. We sampled beetles in 2011-2014 and show that the distribution of M4 in the United States is dynamic and has shifted southward. By using a genetic marker of Medea-1 (M1), we found five unique geographic clusters with high and low M1 frequencies in a pattern not predicted by microsatellite-based analysis of population structure. Our results indicate the absence of rigid barriers to Medea spread in the United States, so assessment of what factors have limited its current distribution requires further investigation. There is great interest in using synthetic SGEs, including synthetic Medea, to alter or suppress pest populations, but there is concern about unpredicted spread of these SGEs and potential for populations to become resistant to them. The finding of patchy distributions of Medea elements suggests that released synthetic SGEs cannot always be expected to spread uniformly, especially in target species with limited dispersal.

Population genetics of fruit bat reservoir informs the dynamics, distribution and diversity of Nipah virus

The structure and connectivity of wildlife host populations may influence zoonotic disease dynamics, evolution and therefore spillover risk to people. Fruit bats in the genus Pteropus, or flying foxes, are the primary natural reservoir for henipaviruses-a group of emerging paramyxoviruses that threaten livestock and public health. In Bangladesh, Pteropus medius is the reservoir for Nipah virus-and viral spillover has led to human fatalities nearly every year since 2001. Here, we use mitochondrial DNA and nuclear microsatellite markers to measure the population structure, demographic history and phylogeography of P. medius in Bangladesh. We combine this with a phylogeographic analysis of all known Nipah virus sequences and strains currently available to better inform the dynamics, distribution and evolutionary history of Nipah virus. We show that P. medius is primarily panmictic, but combined analysis of microsatellite and morphological data shows evidence for differentiation of two populations in eastern Bangladesh, corresponding to a divergent strain of Nipah virus also found in bats from eastern Bangladesh. Our demographic analyses indicate that a large, expanding population of flying foxes has existed in Bangladesh since the Late Pleistocene, coinciding with human population expansion in South Asia, suggesting repeated historical spillover of Nipah virus likely occurred. We present the first evidence of mitochondrial introgression, or hybridization, between P. medius and flying fox species found in South-East Asia (P. vampyrus and P. hypomelanus), which may help to explain the distribution of Nipah virus strains across the region.

To remain or leave: Dispersal variation and its genetic consequences in benthic freshwater invertebrates

Variation in dispersal capacity may influence population genetic variation and relatedness of freshwater animals thus demonstrating how life-history traits influence patterns and processes that in turn influence biodiversity. The majority of studies have focused on the consequences of dispersal variation in taxa inhabiting riverine systems whose dendritic nature and upstream/downstream gradients facilitate characterizing populations along networks. We undertook extensive, large-scale investigations of the impacts of hydrological connectivity on population genetic variation in two freshwater bryozoan species whose dispersive propagules (statoblasts) are either attached to surfaces (Fredericella sultana) or are released as buoyant stages (Cristatella mucedo) and that live primarily in either lotic (F. sultana) or lentic environments (C. mucedo). Describing population genetic structure in multiple sites characterized by varying degrees of hydrological connectivity within each of three (or four) UK regions enabled us to test the following hypotheses: (1) genetic diversity and gene flow will be more influenced by hydrological connectivity in populations of C. mucedo (because F. sultana dispersal stages are retained); (2) populations of F. sultana will be characterized by greater genetic divergence than those of C. mucedo (reflecting their relative dispersal capacities); and (3) genetic variation will be greatest in F. sultana (reflecting a propensity for genetic divergence as a result of its low dispersal potential). We found that hydrological connectivity enhanced genetic diversity and gene flow among C. mucedo populations but not in F. sultana while higher overall measures of clonal diversity and greater genetic divergence characterized populations of F. sultana. We suggest that genetic divergence over time within F. sultana populations reflects a general constraint of releasing propagules that might eventually be swept to sea when taxa inhabit running waters. In contrast, taxa that primarily inhabit lakes and ponds may colonize across hydrologically connected regions, establishing genetically related populations. Our study contributes more nuanced views about drivers of population genetic structures in passively dispersing freshwater invertebrates as outlined by the Monopolization Hypothesis (Acta Oecologica, 23, 2002, 121) by highlighting how a range of demographic and evolutionary processes reflect life-history attributes of benthic colonial invertebrates (bryozoans) and cyclically parthenogenetic zooplankton. In addition, growing evidence that genetic divergence may commonly characterize populations of diverse groups of riverine taxa suggests that organisms inhabiting lotic systems may be particularly challenged by environmental change. Such change may predispose riverine populations to extinction as a result of genetic divergence combined with limited dispersal and gene flow.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.1tm8705.

Seeing-good-gene-based mate choice: From genes to behavioural preferences

Although vertebrates have been reported to gain higher reproductive outputs by choosing mates, few studies have been conducted on threatened species. However, species recovery should benefit if natural mate choice could improve reproductive output (i.e. pair performance related to offspring number, such as increased clutch size, numbers of fertilized egg and fledglings). We assessed the evidence for major histocompatibility complex (MHC)‐based mate preference in the endangered crested ibis (Nipponia nippon) and quantified the impacts of such choice on reproductive output.

We tested the hypothesis that crested ibis advertise “good genes” through external traits, by testing whether nuptial plumage characteristics and body morphology mediate mate choice for underlying genetic MHC variation.

We found differences between males and females in preferred MHC genotypes, external traits used in mate choice and contributions to reproductive outputs. Females preferred MHC‐heterozygous males, which had darker [i.e. lower total reflectance and ultraviolet (UV) reflectance] nuptial plumage. Males preferred females lacking the DAB*d allele at the MHC class II DAB locus, which had higher average body mass. DAB*d‐free females yielded heavier eggs and more fledglings, while MHC‐heterozygous males contributed to more fertilized eggs and fledglings. Fledging rate was highest when both parents had the preferred MHC genotypes (i.e. MHC‐heterozygous father and DAB*d‐free mother). Comparisons showed that free‐mating wild and semi‐natural pairs yielded more fertilized eggs and more fledglings, with a higher fledging rate, than captive pairs matched artificially based on pedigree.

Conservation programmes seldom apply modern research results to population management, which could hinder recovery of threatened species. Our results show that mate choice can play an important role in improving reproductive output, with an example in which an endangered bird selects mates using UV visual capability. Despite the undoubted importance of pedigree‐based matching of mates in conservation programmes, we show that free mating can be a better alternative strategy.

Estimating seed dispersal distance: A comparison of methods using animal movement and plant genetic data on two primate-dispersed Neotropical plant species

Seed dispersal distance (SDD) critically influences the survival of seedlings, spatial patterns of genetic diversity within plant populations, and gene flow among plant populations. In animal-dispersed species, foraging behavior and movement patterns determine SDD. Direct observations of seed dispersal events by animals in natural plant populations are mostly constrained by the high mobility and low visibility of seed dispersers. Therefore, diverse alternative methods are used to estimate seed dispersal distance, but direct comparisons of these approaches within the same seed dispersal system are mostly missing.We investigated two plant species with different life history traits, Leonia cymosa and Parkia panurensis, exclusively dispersed by two tamarin species, Saguinus mystax and Leontocebus nigrifrons. We compared SDD estimates obtained from direct observations, genetic identification of mother plants from seed coats, parentage analysis of seedlings/saplings, and phenomenological and mechanistic modeling approaches.SDD derived from the different methods ranged between 158 and 201 m for P. panurensis and between 178 and 318 m for L. cymosa. In P. panurensis, the modeling approaches resulted in moderately higher estimates than observations and genotyping of seed coats. In L. cymosa, parentage analysis resulted in a lower estimate than all other methods. Overall, SDD estimates for P. panurensis (179 ± 16 m; mean ± SD) were significantly lower than for L. cymosa (266 ± 59 m; mean ± SD).Differences among methods were related to processes of the seed dispersal loop integrated by the respective methods (e.g., seed deposition or seedling distribution). We discuss the merits and limitations of each method and highlight the aspects to be considered when comparing SDD derived from different methodologies. Differences among plant species were related to differences in reproductive traits influencing gut passage time and feeding behavior, highlighting the importance of plant traits on animal-mediated seed dispersal distance.

Phylogeography and paleodistribution models of a widespread birch (Betula platyphylla Suk.) across East Asia: Multiple refugia, multidirectional expansion, and heterogeneous genetic pattern

Widespread tree species cover large geographical areas and play important roles in various vegetation types. Understanding how these species responded to historical climatic changes is important for understanding community assembly mechanisms with evolutionary and conservation implications. However, the location of refugial areas and postglacial history of widespread trees in East Asia remain poorly known. We combined microsatellite data (63 populations, 1756 individuals) and ecological niche modeling to examine the range-wide population diversity, genetic structure, and historical demography of a pioneer tree species, Asian white birch (Betula platyphylla Suk.) across East Asia. We found a north-to-south trend of declining genetic diversity and five clusters, corresponding to geographical regions. Different clusters were inferred to have diverged through Pleistocene climatic oscillations and have different expansion routes, leading to genetic admixture in some populations. Ecological niche models indicated that the distribution of B. platyphylla during the last glacial maximum still had a large latitude span with slight shifts toward southeast, and northern populations had more variable distribution ranges than those in the south during later climatic oscillations. Our results reflect the relatively stable distribution through the last glacial-interglacial cycles and recent multidirectional expansion of B. platyphylla, providing new hypotheses for the response pattern of widespread tree species to climate change. The gradual genetic pattern from northeast to southwest and alternative distribution dynamics possibly resulted from environmental differences caused by latitude and topographic heterogeneity.

Assessment of flux through oleoresin biosynthesis in epithelial cells of loblolly pine resin ducts

The shoot system of pines contains abundant resin ducts, which harbor oleoresins that play important roles in constitutive and inducible defenses. In a pilot study, we assessed the chemical diversity of oleoresins obtained from mature tissues of loblolly pine trees (Pinus taeda L.). Building on these data sets, we designed experiments to assess oleoresin biosynthesis in needles of 2-year-old saplings. Comparative transcriptome analyses of single cell types indicated that genes involved in the biosynthesis of oleoresins are significantly enriched in isolated epithelial cells of resin ducts, compared with those expressed in mesophyll cells. Simulations using newly developed genome-scale models of epithelial and mesophyll cells, which incorporate our data on oleoresin yield and composition as well as gene expression patterns, predicted that heterotrophic metabolism in epithelial cells involves enhanced levels of oxidative phosphorylation and fermentation (providing redox and energy equivalents). Furthermore, flux was predicted to be more evenly distributed across the metabolic network of mesophyll cells, which, in contrast to epithelial cells, do not synthesize high levels of specialized metabolites. Our findings provide novel insights into the remarkable specialization of metabolism in epithelial cells.

Are Survivors Different? Genetic-Based Selection of Trees by Mountain Pine Beetle During a Climate Change-Driven Outbreak in a High-Elevation Pine Forest

Increased mortality of forest trees, driven directly or indirectly by climate change, is occurring around the world. In western North America, whitebark pine, a high elevation keystone species, and lodgepole pine, a widespread ecologically and economically important tree, have experienced extensive mortality in recent climate-driven outbreaks of the mountain pine beetle. However, even in stands experiencing high levels of mortality, some mature trees have survived. We hypothesized that the outbreak acted as a natural selection event, removing trees most susceptible to the beetle and least adapted to warmer drier conditions. If this was the case, genetic change would be expected at loci underlying beetle resistance. Given we did not know the basis for resistance, we used inter-simple sequence repeats to compare the genetic profiles of two sets of trees, survivors (mature, living trees) and general population (trees just under the diameter preferred by the beetles and expected to approximate the genetic structure of each tree species at the site without beetle selection). This method detects high levels of polymorphism and has often been able to detect patterns associated with phenotypic traits. For both whitebark and lodgepole pine, survivors and general population trees mostly segregated independently indicating a genetic basis for survivorship. Exceptions were a few general population trees that segregated with survivors in proportions roughly reflecting the proportion of survivors versus beetle-killed trees. Our results indicate that during outbreaks, beetle choice may result in strong selection for trees with greater resistance to attack. Our findings suggest that survivorship is genetically based and, thus, heritable. Therefore, retaining survivors after outbreaks to act as primary seed sources could act to promote adaptation. Further research will be needed to characterize the actual mechanism(s) of resistance.

The genetics of drought tolerance in conifers

Contents 1034 I. 1034 II. 1035 III. 1037 IV. 1038 V. 1042 VI. 1043 VII. 1045 References 1045 SUMMARY: As temperatures warm and precipitation patterns shift as a result of climate change, interest in the identification of tree genotypes that will thrive under more arid conditions has grown. In this review, we discuss the multiple definitions of 'drought tolerance' and the biological processes involved in drought responses. We describe the three major approaches taken in the study of genetic variation in drought responses, the advantages and shortcomings of each, and what each of these approaches has revealed about the genetic basis of adaptation to drought in conifers. Finally, we discuss how a greater knowledge of the genetics of drought tolerance may aid forest management, and provide recommendations for how future studies may overcome the limitations of past approaches. In particular, we urge a more direct focus on survival, growth and the traits that directly predict them (rather than on proxies, such as water use efficiency), combining research approaches with complementary strengths and weaknesses, and the inclusion of a wider range of taxa and life stages.

Differential physiological and genetic responses of five European Scots pine provenances to induced water stress

Pinus sylvestris L. (Scots pine) is the conifer with widest natural distribution area. Portugal constitutes its westernmost limit of distribution. Most of the Portuguese populations were planted but two autochthonous populations were recently ascribed to 'Serra do Gerês' (NW Portugal), and seem to be well adapted to the temperate climate. However, the ongoing climate changes may compromise their survival. With this study we intend to evaluate the anatomic-physiological and genetic responses of Scots pine from five European provenances ('Gerês', 'Puebla de Lillo', 'Montes Universales', Germany and Sweden) to three water availability regimes, in order to determine which one(s) present higher resistance to drought. Individuals from 'Gerês' presented the highest stability in photosynthetic reactions as well as the better photochemical and metabolic behaviours under drought (T3). Hence, the relative expression ratio of three water stress-responsive genes during drought was lower and gradual in 'Gerês', compared to all other provenances, followed by Germany. The results achieved in 'Gerês' and Germany provenances are very interesting since they reflected that the native populations of 'Gerês' along with the Portuguese Scots pine planted populations with a probable German provenance, have ability and high adaptive potential to respond to situations of water deficit. Moreover, the present genetic and physiological data demonstrated the urgent demand for the conservation of Portuguese Scots pine genetic resources as well as its use in plantation/afforestation of areas where the warming and drought has been affecting the survival of this species.

Different Alleles of a Gene Encoding Leucoanthocyanidin Reductase (PaLAR3) Influence Resistance against the Fungus Heterobasidion parviporum in Picea abies

Despite the fact that fungal diseases are a growing menace for conifers in modern silviculture, only a very limited number of molecular markers for pathogen resistance have been validated in conifer species. A previous genetic study indicated that the resistance of Norway spruce (Picea abies) to Heterobasidion annosum s.l., a pathogenic basidiomycete species complex, is linked to a quantitative trait loci that associates with differences in fungal growth in sapwood (FGS) that includes a gene, PaLAR3, which encodes a leucoanthocyanidin reductase. In this study, gene sequences showed the presence of two PaLAR3 allelic lineages in P. abies. Higher resistance was associated with the novel allele, which was found in low frequency in the four P. abies populations that we studied. Norway spruce plants carrying at least one copy of the novel allele showed a significant reduction in FGS after inoculation with Heterobasidion parviporum compared to their half-siblings carrying no copies, indicating dominance of this allele. The amount of (+) catechin, the enzymatic product of PaLAR3, was significantly higher in bark of trees homozygous for the novel allele. Although we observed that the in vitro activities of the enzymes encoded by the two alleles were similar, we could show that allele-specific transcript levels were significantly higher for the novel allele, indicating that regulation of gene expression is responsible for the observed effects in resistance, possibly caused by differences in cis-acting elements that we observe in the promoter region of the two alleles.

Genetic architecture of wood properties based on association analysis and co-expression networks in white spruce

Association studies are widely utilized to analyze complex traits but their ability to disclose genetic architectures is often limited by statistical constraints, and functional insights are usually minimal in nonmodel organisms like forest trees. We developed an approach to integrate association mapping results with co-expression networks. We tested single nucleotide polymorphisms (SNPs) in 2652 candidate genes for statistical associations with wood density, stiffness, microfibril angle and ring width in a population of 1694 white spruce trees (Picea glauca). Associations mapping identified 229-292 genes per wood trait using a statistical significance level of P < 0.05 to maximize discovery. Over-representation of genes associated for nearly all traits was found in a xylem preferential co-expression group developed in independent experiments. A xylem co-expression network was reconstructed with 180 wood associated genes and several known MYB and NAC regulators were identified as network hubs. The network revealed a link between the gene PgNAC8, wood stiffness and microfibril angle, as well as considerable within-season variation for both genetic control of wood traits and gene expression. Trait associations were distributed throughout the network suggesting complex interactions and pleiotropic effects. Our findings indicate that integration of association mapping and co-expression networks enhances our understanding of complex wood traits.

Conflicting genomic signals affect phylogenetic inference in four species of North American pines

Adaptive evolutionary processes in plants may be accompanied by episodes of introgression, parallel evolution and incomplete lineage sorting that pose challenges in untangling species evolutionary history. Genus Pinus (pines) is one of the most abundant and most studied groups among gymnosperms, and a good example of a lineage where these phenomena have been observed. Pines are among the most ecologically and economically important plant species. Some, such as the pines of the southeastern USA (southern pines in subsection Australes), are subjects of intensive breeding programmes. Despite numerous published studies, the evolutionary history of Australes remains ambiguous and often controversial. We studied the phylogeny of four major southern pine species: shortleaf (Pinus echinata), slash (P. elliottii), longleaf (P. palustris) and loblolly (P. taeda), using sequences from 11 nuclear loci and maximum likelihood and Bayesian methods. Our analysis encountered resolution difficulties similar to earlier published studies. Although incomplete lineage sorting and introgression are two phenomena presumptively underlying our results, the phylogenetic inferences seem to be also influenced by the genes examined, with certain topologies supported by sets of genes sharing common putative functionalities. For example, genes involved in wood formation supported the clade echinata-taeda, genes linked to plant defence supported the clade echinata-elliottii and genes linked to water management properties supported the clade echinata-palustris The support for these clades was very high and consistent across methods. We discuss the potential factors that could underlie these observations, including incomplete lineage sorting, hybridization and parallel or adaptive evolution. Our results likely reflect the relatively short evolutionary history of the subsection that is thought to have begun during the middle Miocene and has been influenced by climate fluctuations.

Landscape genomic analysis of candidate genes for climate adaptation in a California endemic oak, Quercus lobata

PREMISE OF THE STUDY

The ability of California tree populations to survive anthropogenic climate change will be shaped by the geographic structure of adaptive genetic variation. Our goal is to test whether climate-associated candidate genes show evidence of spatially divergent selection in natural populations of valley oak, Quercus lobata, as preliminary indication of local adaptation.

METHODS

Using DNA from 45 individuals from 13 localities across the species' range, we sequenced portions of 40 candidate genes related to budburst/flowering, growth, osmotic stress, and temperature stress. Using 195 single nucleotide polymorphisms (SNPs), we estimated genetic differentiation across populations and correlated allele frequencies with climate gradients using single-locus and multivariate models.

RESULTS

The top 5% of FST estimates ranged from 0.25 to 0.68, yielding loci potentially under spatially divergent selection. Environmental analyses of SNP frequencies with climate gradients revealed three significantly correlated SNPs within budburst/flowering genes and two SNPs within temperature stress genes with mean annual precipitation, after controlling for multiple testing. A redundancy model showed a significant association between SNPs and climate variables and revealed a similar set of SNPs with high loadings on the first axis. In the RDA, climate accounted for 67% of the explained variation, when holding climate constant, in contrast to a putatively neutral SSR data set where climate accounted for only 33%.

CONCLUSIONS

Population differentiation and geographic gradients of allele frequencies in climate-associated functional genes in Q. lobata provide initial evidence of adaptive genetic variation and background for predicting population response to climate change.

The genetic regulatory network centered on Pto-Wuschela and its targets involved in wood formation revealed by association studies

Transcription factors (TFs) regulate gene expression and can strongly affect phenotypes. However, few studies have examined TF variants and TF interactions with their targets in plants. Here, we used genetic association in 435 unrelated individuals of Populus tomentosa to explore the variants in Pto-Wuschela and its targets to decipher the genetic regulatory network of Pto-Wuschela. Our bioinformatics and co-expression analysis identified 53 genes with the motif TCACGTGA as putative targets of Pto-Wuschela. Single-marker association analysis showed that Pto-Wuschela was associated with wood properties, which is in agreement with the observation that it has higher expression in stem vascular tissues in Populus. Also, SNPs in the 53 targets were associated with growth or wood properties under additive or dominance effects, suggesting these genes and Pto-Wuschela may act in the same genetic pathways that affect variation in these quantitative traits. Epistasis analysis indicated that 75.5% of these genes directly or indirectly interacted Pto-Wuschela, revealing the coordinated genetic regulatory network formed by Pto-Wuschela and its targets. Thus, our study provides an alternative method for dissection of the interactions between a TF and its targets, which will strength our understanding of the regulatory roles of TFs in complex traits in plants.

Association Analysis of Simple Sequence Repeat (SSR) Markers with Agronomic Traits in Tall Fescue (Festuca arundinacea Schreb.)

Tall fescue is widely used in temperate regions throughout the world as a dominant forage grass as well as a turfgrass, in pastoral and turf industry. However, the utilization of tall fescue was limited because of its leaf roughness, poor regeneration ability and poor stress resistance. New cultivars were desirable in modern pastoral industries exceed the potential of existing cultivars. Therefore, well understanding the agronomic traits and describing germplasms would help to overcome these constraints, and morphological evaluation of tall fescue germplasm is the key component in selecting rational parents for hybridization breeding. However, describing the morphological traits of tall fescue germplasm is costly and time-consuming. Fortunately, biotechnology approaches can supplement conventional breeding efforts for tall fescue improvement. Association mapping, as a powerful approach to identify association between agronomic traits and molecular markers has been widely used for enhancing the utilization, conservation and management of the tall fescue germplasms. Therefore, in the present research, 115 tall fescue accessions from different origins (25 accessions are cultivars; 31 accessions from America; 32 accessions from European; 7 accessions from Africa; 20 accessions from Asia), were evaluated for agronomic traits and genetic diversity with 90 simple sequence repeat (SSR) markers. The panel displayed significant variation in spike count per plant (SCP) and spike weight (SW). However, BCS performed the lowest CV among all the observed agronomic traits. Three subpopulations were identified within the collections but no obvious relative kinship (K) was found. The GLM model was used to describe the association between SSR and agronomic traits. Fifty-one SSR markers associated with agronomic traits were observed. Twelve single-associated markers were associated with PH; six single-associated markers were associated with BCS; eight single-associated markers were associated with SW; five single-associated markers were associated with SC; seven single-associated markers were associated with SCP; three single-associated markers were associated with SL. Especially, we observed that the genetic variation of SW was explained 11.6 % by M37 marker. It is interesting to observe that nine markers (M1, M2, M35, M54 marker was associated with both BCS and SC; M3, M4 markers were associated with BCS, SW, and SC; M19 marker was associated with both pH and PD, M40 marker was associated with both SCP and SW; and M193 marker was associated with both PH and SL) were associated with more than two agronomic traits. Notably, Branch count per spike (BCS) was explained by four markers (M1, M2, M3, and M4) exceeding 10 %. These identified marker alleles associated with agronomic traits could provide important information and markers for molecular-assisted breeding that facilitate the breeding process in tall fescue.

Transcript abundance patterns of Populus C-repeat binding factor2 orthologs and genetic association of PsCBF2 allelic variation with physiological and biochemical traits in response to abiotic stress

We conducted a candidate gene-based approach to search for genetic associations between 10 SNPs in PsCBF2 and 10 abiotic stress-related traits. The increasing incidence of abiotic stresses and the limitations of available treatments, particularly in trees, highlight the need to improve our understanding of the mechanisms of stress responses. In Arabidopsis, C-repeat binding factor 2 (CBF2) plays an important role in freezing tolerance and cold acclimation. Here, we isolated orthologs of CBF2 from five Populus species. Expression profiling revealed that the Populus CBF2s were preferentially induced in response to cold, with CBF2 transcript abundances ranging from 5.4- to 62-fold higher than in unstressed controls of the corresponding species. In addition, we used a candidate gene-based approach in Populus simonii Carr. to identify single nucleotide polymorphisms (SNPs) in PsCBF2 associated with physiological and biochemical traits. PsCBF2 showed high nucleotide diversity (π T = 0.00549, θ w = 0.01406) and low average linkage disequilibrium (r (2) = 0.061). Association studies in 528 individuals of an association population showed that nine SNPs (false discovery rate Q < 0.10) and one haplotype (Q < 0.10) were significantly associated with differences in four physiological and biochemical traits (P < 0.005), with each marker explaining 1.31-5.87 % of the total variance in the corresponding trait. PsCBF2 transcript levels differed significantly in abundance among genotypic classes for most of the significant SNPs. Identification of these significant associations will help reveal the molecular basis of physiological differences and provide a starting point for marker-assisted selection for traits involved in stress tolerance in P. simonii.

Exploring potential of pearl millet germplasm association panel for association mapping of drought tolerance traits

A pearl millet inbred germplasm association panel (PMiGAP) comprising 250 inbred lines, representative of cultivated germplasm from Africa and Asia, elite improved open-pollinated cultivars, hybrid parental inbreds and inbred mapping population parents, was recently established. This study presents the first report of genetic diversity in PMiGAP and its exploitation for association mapping of drought tolerance traits. For diversity and genetic structure analysis, PMiGAP was genotyped with 37 SSR and CISP markers representing all seven linkage groups. For association analysis, it was phenotyped for yield and yield components and morpho-physiological traits under both well-watered and drought conditions, and genotyped with SNPs and InDels from seventeen genes underlying a major validated drought tolerance (DT) QTL. The average gene diversity in PMiGAP was 0.54. The STRUCTURE analysis revealed six subpopulations within PMiGAP. Significant associations were obtained for 22 SNPs and 3 InDels from 13 genes under different treatments. Seven SNPs associations from 5 genes were common under irrigated and one of the drought stress treatments. Most significantly, an important SNP in putative acetyl CoA carboxylase gene showed constitutive association with grain yield, grain harvest index and panicle yield under all treatments. An InDel in putative chlorophyll a/b binding protein gene was significantly associated with both stay-green and grain yield traits under drought stress. This can be used as a functional marker for selecting high yielding genotypes with 'stay green' phenotype under drought stress. The present study identified useful marker-trait associations of important agronomics traits under irrigated and drought stress conditions with genes underlying a major validated DT-QTL in pearl millet. Results suggest that PMiGAP is a useful panel for association mapping. Expression patterns of genes also shed light on some physiological mechanisms underlying pearl millet drought tolerance.

Differential Gene Expression Reveals Candidate Genes for Drought Stress Response in Abies alba (Pinaceae)

Increasing drought periods as a result of global climate change pose a threat to many tree species by possibly outpacing their adaptive capabilities. Revealing the genetic basis of drought stress response is therefore implemental for future conservation strategies and risk assessment. Access to informative genomic regions is however challenging, especially for conifers, partially due to their large genomes, which puts constraints on the feasibility of whole genome scans. Candidate genes offer a valuable tool to reduce the complexity of the analysis and the amount of sequencing work and costs. For this study we combined an improved drought stress phenotyping of needles via a novel terahertz water monitoring technique with Massive Analysis of cDNA Ends to identify candidate genes for drought stress response in European silver fir (Abies alba Mill.). A pooled cDNA library was constructed from the cotyledons of six drought stressed and six well-watered silver fir seedlings, respectively. Differential expression analyses of these libraries revealed 296 candidate genes for drought stress response in silver fir (247 up- and 49 down-regulated) of which a subset was validated by RT-qPCR of the twelve individual cotyledons. A majority of these genes code for currently uncharacterized proteins and hint on new genomic resources to be explored in conifers. Furthermore, we could show that some traditional reference genes from model plant species (GAPDH and eIF4A2) are not suitable for differential analysis and we propose a new reference gene, TPC1, for drought stress expression profiling in needles of conifer seedlings.

Identification of additive, dominant, and epistatic variation conferred by key genes in cellulose biosynthesis pathway in Populus tomentosa†

Economically important traits in many species generally show polygenic, quantitative inheritance. The components of genetic variation (additive, dominant and epistatic effects) of these traits conferred by multiple genes in shared biological pathways remain to be defined. Here, we investigated 11 full-length genes in cellulose biosynthesis, on 10 growth and wood-property traits, within a population of 460 unrelated Populus tomentosa individuals, via multi-gene association. To validate positive associations, we conducted single-marker analysis in a linkage population of 1,200 individuals. We identified 118, 121, and 43 associations (P< 0.01) corresponding to additive, dominant, and epistatic effects, respectively, with low to moderate proportions of phenotypic variance (R(2)). Epistatic interaction models uncovered a combination of three non-synonymous sites from three unique genes, representing a significant epistasis for diameter at breast height and stem volume. Single-marker analysis validated 61 associations (false discovery rate, Q ≤ 0.10), representing 38 SNPs from nine genes, and its average effect (R(2) = 3.8%) nearly 2-fold higher than that identified with multi-gene association, suggesting that multi-gene association can capture smaller individual variants. Moreover, a structural gene-gene network based on tissue-specific transcript abundances provides a better understanding of the multi-gene pathway affecting tree growth and lignocellulose biosynthesis. Our study highlights the importance of pathway-based multiple gene associations to uncover the nature of genetic variance for quantitative traits and may drive novel progress in molecular breeding.

Association genetics and expression patterns of a CBF4 homolog in Populus under abiotic stress

New strategies for prevention and treatment of abiotic stress require an improved understanding of stress responses. Here, we examined response differences of a C-repeat binding factor gene (PsCBF4) between five species in the genus Populus. We also used a candidate gene-based approach to identify single nucleotide polymorphisms (SNPs) within PsCBF4 that were associated with physiological and biochemical traits in a natural population (528 unrelated individuals) of Populus simonii. We first isolated a 1,044-bp PsCBF4 cDNA encoding a polypeptide of 256 amino acids. Expression profiling revealed that CBF4 is differentially expressed under cold, heat, drought, and salt conditions among five Populus species. Cold stress is the most significant interspecific difference, and PsCBF4 transcript levels ranged from 6.5 to 379.5 times higher than in unstressed controls. A natural population of P. simonii showed high nucleotide diversity (π T = 0.00880, θ w = 0.01192) and low linkage disequilibrium (r (2) ≥ 0.1, within 700 bp) across PsCBF4. Association analysis showed that nine SNPs (false discovery rate Q < 0.10) and two haplotypes (Q < 0.10) were significantly associated with six physiological and biochemical traits, with each marker explaining 3.36-6.12 % of the phenotypic variance in the corresponding trait. Transcript analysis further detected significant differences among genotypic classes for all significant SNPs. Identification of these significant associations will help reveal the molecular bases of physiological and biochemical differences and provide a starting point for marker-assisted selection for traits involved in stress tolerance in P. simonii.

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra- and interspecific genetic variation at 25 climate-related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two-thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. FST outlier analysis identified six SNPs in five climate-related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation-with-migration model indicated that gene flow between the species has been more restricted for climate-related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species-specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.

Spatial and temporal genetic structure of a river-resident Atlantic salmon (Salmo salar) after millennia of isolation

The river-resident Salmo salar ("småblank") has been isolated from other Atlantic salmon populations for 9,500 years in upper River Namsen, Norway. This is the only European Atlantic salmon population accomplishing its entire life cycle in a river. Hydropower development during the last six decades has introduced movement barriers and changed more than 50% of the river habitat to lentic conditions. Based on microsatellites and SNPs, genetic variation within småblank was only about 50% of that in the anadromous Atlantic salmon within the same river. The genetic differentiation (F ST) between småblank and the anadromous population was 0.24. This is similar to the differentiation between anadromous Atlantic salmon in Europe and North America. Microsatellite analyses identified three genetic subpopulations within småblank, each with an effective population size Ne of a few hundred individuals. There was no evidence of reduced heterozygosity and allelic richness in contemporary samples (2005-2008) compared with historical samples (1955-56 and 1978-79). However, there was a reduction in genetic differentiation between sampling localities over time. SNP data supported the differentiation of småblank into subpopulations and revealed downstream asymmetric gene flow between subpopulations. In spite of this, genetic variation was not higher in the lower than in the upper areas. The meta-population structure of småblank probably maintains genetic variation better than one panmictic population would do, as long as gene flow among subpopulations is maintained. Småblank is a unique endemic island population of Atlantic salmon. It is in a precarious situation due to a variety of anthropogenic impacts on its restricted habitat area. Thus, maintaining population size and avoiding further habitat fragmentation are important.

Genetic control of functional traits related to photosynthesis and water use efficiency in Pinus pinaster Ait. drought response: integration of genome annotation, allele association and QTL detection for candidate gene identification

BACKGROUND

Understanding molecular mechanisms that control photosynthesis and water use efficiency in response to drought is crucial for plant species from dry areas. This study aimed to identify QTL for these traits in a Mediterranean conifer and tested their stability under drought.

RESULTS

High density linkage maps for Pinus pinaster were used in the detection of QTL for photosynthesis and water use efficiency at three water irrigation regimes. A total of 28 significant and 27 suggestive QTL were found. QTL detected for photochemical traits accounted for the higher percentage of phenotypic variance. Functional annotation of genes within the QTL suggested 58 candidate genes for the analyzed traits. Allele association analysis in selected candidate genes showed three SNPs located in a MYB transcription factor that were significantly associated with efficiency of energy capture by open PSII reaction centers and specific leaf area.

CONCLUSIONS

The integration of QTL mapping of functional traits, genome annotation and allele association yielded several candidate genes involved with molecular control of photosynthesis and water use efficiency in response to drought in a conifer species. The results obtained highlight the importance of maintaining the integrity of the photochemical machinery in P. pinaster drought response.

Signatures of volcanism and aridity in the evolution of an insular pine (Pinus canariensis Chr. Sm. Ex DC in Buch)

Oceanic islands of volcanic origin provide useful templates for the study of evolution because they are subjected to recurrent perturbations that generate steep environmental gradients that may promote adaptation. Here we combine population genetic data from nuclear genes with the analysis of environmental variation and phenotypic data from common gardens to disentangle the confounding effects of demography and selection to identify the factors of importance for the evolution of the insular pine P. canariensis. Eight nuclear genes were partially sequenced in a survey covering the entire species range, and phenotypic traits were measured in four common gardens from contrasting environments. The explanatory power of population substrate age and environmental indices were assessed against molecular and phenotypic diversity estimates. In addition, neutral genetic variability (FST) and the genetic differentiation of phenotypic variation (QST) were compared in order to identify the evolutionary forces acting on these traits. Two key factors in the evolution of the species were identified: (1) recurrent volcanic activity has left an imprint in the genetic diversity of the nuclear genes; (2) aridity in southern slopes promotes local adaptation in the driest localities of P. canariensis, despite high levels of gene flow among populations.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.

Genetic structure of pike (Esox lucius) reveals a complex and previously unrecognized colonization history of Ireland

AIM

We investigated genetic variation of Irish pike populations and their relationship with European outgroups, in order to elucidate the origin of this species to the island, which is largely assumed to have occurred as a human-mediated introduction over the past few hundred years. We aimed thereby to provide new insights into population structure to improve fisheries and biodiversity management in Irish freshwaters.

LOCATION

Ireland, Britain and continental Europe.

METHODS

A total of 752 pike (Esox lucius) were sampled from 15 locations around Ireland, and 9 continental European sites, and genotyped at six polymorphic microsatellite loci. Patterns and mechanisms of population genetic structure were assessed through a diverse array of methods, including Bayesian clustering, hierarchical analysis of molecular variance, and approximate Bayesian computation.

RESULTS

Varying levels of genetic diversity and a high degree of population genetic differentiation were detected. Clear substructure within Ireland was identified, with two main groups being evident. One of the Irish populations showed high similarity with British populations. The other, more widespread, Irish strain did not group with any European population examined. Approximate Bayesian computation suggested that this widespread Irish strain is older, and may have colonized Ireland independently of humans.

MAIN CONCLUSIONS

Population genetic substructure in Irish pike is high and comparable to the levels observed elsewhere in Europe. A comparison of evolutionary scenarios upholds the possibility that pike may have colonized Ireland in two 'waves', the first of which, being independent of human colonization, would represent the first evidence for natural colonization of a non-anadromous freshwater fish to the island of Ireland. Although further investigations using comprehensive genomic techniques will be necessary to confirm this, the present results warrant a reappraisal of current management strategies for this species.