Microsatellite markers derived from Quercus mongolica var. crispula (Fagaceae) inner bark expressed sequence tags

Asymmetric character displacement in mixed oak stands

Ecological character displacement (ECD) refers to a pattern of increased divergence at sites where species ranges overlap caused by competition for resources. Although ECD is believed to be common, there are few in-depth studies that clearly establish its existence, especially in plants. Thus, we have compared leaf traits in allopatric and sympatric populations of two East Asian deciduous oaks: Quercus dentata and Quercus aliena. In contrast to previous studies, we define sympatry and allopatry at a local scale, thereby comparing populations that can or cannot directly interact. Using genetic markers, we found greater genetic divergence between the two oak species growing in mixed stands and inferred that long-term gene flow has predominantly occurred asymmetrically from the cold-tolerant species (Q. dentata) to the warm-demanding later colonizing species (Q. aliena). Analysis of leaf traits revealed greater divergence in mixed than in pure oak stands. This was mostly due to the later colonizing species being characterized by more resource-conservative traits in the presence of the other species. Controlling for relevant environmental differences did not alter these conclusions. These results suggest that asymmetric trait divergence can take place where species coexist, possibly due to the imbalance in demographic history of species resulting in asymmetric inter-specific selection pressures.

Multiomics Molecular Research into the Recalcitrant and Orphan Quercus ilex Tree Species: Why, What for, and How

The holm oak (Quercus ilex L.) is the dominant tree species of the Mediterranean forest and the Spanish agrosilvopastoral ecosystem, “dehesa.” It has been, since the prehistoric period, an important part of the Iberian population from a social, cultural, and religious point of view, providing an ample variety of goods and services, and forming the basis of the economy in rural areas. Currently, there is renewed interest in its use for dietary diversification and sustainable food production. It is part of cultural richness, both economically (tangible) and environmentally (intangible), and must be preserved for future generations. However, a worrisome degradation of the species and associated ecosystems is occurring, observed in an increase in tree decline and mortality, which requires urgent action. Breeding programs based on the selection of elite genotypes by molecular markers is the only plausible biotechnological approach. To this end, the authors’ group started, in 2004, a research line aimed at characterizing the molecular biology of Q. ilex. It has been a challenging task due to its biological characteristics (long life cycle, allogamous, high phenotypic variability) and recalcitrant nature. The biology of this species has been characterized following the central dogma of molecular biology using the omics cascade. Molecular responses to biotic and abiotic stresses, as well as seed maturation and germination, are the two main objectives of our research. The contributions of the group to the knowledge of the species at the level of DNA-based markers, genomics, epigenomics, transcriptomics, proteomics, and metabolomics are discussed here. Moreover, data are compared with those reported for Quercus spp. All omics data generated, and the genome of Q. ilex available, will be integrated with morphological and physiological data in the systems biology direction. Thus, we will propose possible molecular markers related to resilient and productive genotypes to be used in reforestation programs. In addition, possible markers related to the nutritional value of acorn and derivate products, as well as bioactive compounds (peptides and phenolics) and allergens, will be suggested. Subsequently, the selected molecular markers will be validated by both genome-wide association and functional genomic analyses.

Population and Landscape Genetics Provide Insights Into Species Conservation of Two Evergreen Oaks in Qinghai-Tibet Plateau and Adjacent Regions

The combination of population and landscape genetics can facilitate the understanding of conservation strategy under the changing climate. Here, we focused on the two most diverse and ecologically important evergreen oaks: Quercus aquifolioides and Quercus spinosa in Qinghai-Tibetan Plateau (QTP), which is considered as world's biodiversity hotspot. We genotyped 1,657 individuals of 106 populations at 15 nuclear microsatellite loci throughout the species distribution range. Spatial patterns of genetic diversity were identified by mapping the allelic richness (AR) and locally common alleles (LCA) according to the circular neighborhood methodology. Migration routes from QTP were detected by historical gene flow estimation. The response pattern of genetic variation to environmental gradient was assessed by the genotype-environment association (GEA) analysis. The overall genetic structure showed a high level of intra-species genetic divergence of a strong west-east pattern. The West-to-East migration route indicated the complex demographic history of two oak species. We found evidence of isolation by the environment in Q. aqu-East and Q. spi-West lineage but not in Q. aqu-West and Q. spi-East lineage. Furthermore, priority for conservation should be given to populations that retain higher spatial genetic diversity or isolated at the edge of the distribution range. Our findings indicate that knowledge of spatial diversity and migration route can provide valuable information for the conservation of existing populations. This study provides an important guide for species conservation for two oak species by the integration of population and landscape genetic methods.

Different population size change and migration histories created genetic diversity of three oaks in Tokai region, central Japan

To understand genetic diversity in focal species, it is important to consider the possibility of speciation with gene flow, especially in species with porous genomes such as oaks. We studied genetic diversity and structure in three oaks, Quercus mongolica var. mongolicoides (QM), Q. mongolica var. crispula (QC) and Q. serrata (QS), growing in the Tokai region, central Japan. QM is semi-endemic to the region while the others are common taxa. We also conducted demographic modeling to infer their population size change and migration histories using an approximate Bayesian computation (ABC) approach. The three taxa showed distinct genetic structures but there was genetic admixture among the taxa, especially between QM and QC. ABC analysis of population size change revealed that the population size of QM was stable during and after the last glacial period, while QC and QS showed population expansion after the last glacial maximum. ABC analysis of population divergence and migration revealed that continuous gene flow between QM and QC after their divergence was supported, while between QM and QS, and between QC and QS, secondary contact after sufficient isolation was supported. These historical migration patterns among the three taxa indicate that QM and QC are currently in the early stage or gray zone of speciation, whereas speciation of the other two taxon pairs is considered to have almost been established. Observed gene flow patterns and strength between QM and QC, and between QM and QS, were explained by both flowering patterns and historical distributions, but those between QC and QS were not.

An admixture of Quercus dentata in the coastal ecotype of Q. mongolica var. crispula in northern Hokkaido and genetic and environmental effects on their traits

In northern Japan, coastal oak forests consist of Quercus dentata (Qd) on the coastal side and Q. mongolica var. crispula (Qc) on the inland side. In the forests of northern Hokkaido, Qd is rare, and a coastal ecotype of Qc with some Qd-like traits grows on the coastal side. To reveal the genetic background of this ecotype, nuclear microsatellite genotypes in closely related oak taxa were obtained from the Eurasian continent, Sakhalin, and Hokkaido. The clustering of these genotypes suggests an admixture of Qd in the coastal ecotype of Qc. Next, we evaluated the effects of admixture and coastal stress on the leaf and shoot traits of Qc and Qd along coastal-inland gradients in northern Hokkaido. The admixture of Qd in Qc was quantified by the Qd ancestry proportions. Coastal stress causes bud mortality in the upper parts of shoots and was quantified by the survival patterns of buds in shoots. The genetic and environmental effects on the traits at Qd-abundant and Qd-rare sites were estimated using linear mixed models. The genetic effect was detected in all traits. Both genetic and environmental effects were detected in most traits. Some traits differed between Qd-abundant and Qd-rare sites in addition to these effects, indicating more Qd-like traits at Qd-rare sites. The findings suggest that an admixture of Qd characterizes the genetic background of the coastal ecotype of Qc and that not only the coastal stress but also the genetic background is responsible for the leaf and shoot traits of Qc and Qd in northern Hokkaido.

Species Boundaries Between Three Sympatric Oak Species: Quercus aliena, Q. dentata, and Q. variabilis at the Northern Edge of Their Distribution in China

Oaks are important timber trees with wide distributions in China, but few genetic studies have been conducted on a fine scale. In this study, we seek to investigate the genetic diversity and differentiation of three sympatric oak species (Quercus aliena Blume, Quercus dentata Thunb. ex Murray, and Quercus variabilis Blume) in their northern distribution in China using 17 bi-parentally inherited nSSRs markers and five maternally inherited chloroplast DNA (cpDNA) fragments. Both the cpDNA and the nSSRs show a high level of genetic differentiation between different oak sections. The chloroplast haplotypes are clustered into two lineages. Clear species boundaries are detected between Q. variabilis and either Q. aliena or Q. dentata. The sharing of chloroplast haplotype H1 between Q. aliena and Q. dentata suggests very recent speciation and incomplete lineage sorting or introgression of H1 from one species to another. The nSSRs data indicate a complete fixation of variation within sites for all three oak species, and that extensive gene flow occurs within species whereas only limited gene flow is detected between Q. aliena and Q. dentata and nearly no gene flow can be detected between Q. aliena and Q. variabilis and between Q. dentata and Q. variabilis. Prezygotic isolation may have contributed to the species boundaries of these three sympatric oak species.

Positively selected amino acid replacements within the RuBisCO enzyme of oak trees are associated with ecological adaptations

Phylogenetic analysis by maximum likelihood (PAML) has become the standard approach to study positive selection at the molecular level, but other methods may provide complementary ways to identify amino acid replacements associated with particular conditions. Here, we compare results of the decision tree (DT) model method with ones of PAML using the key photosynthetic enzyme RuBisCO as a model system to study molecular adaptation to particular ecological conditions in oaks (Quercus). We sequenced the chloroplast rbcL gene encoding RuBisCO large subunit in 158 Quercus species, covering about a third of the global genus diversity. It has been hypothesized that RuBisCO has evolved differentially depending on the environmental conditions and leaf traits governing internal gas diffusion patterns. Here, we show, using PAML, that amino acid replacements at the residue positions 95, 145, 251, 262 and 328 of the RuBisCO large subunit have been the subject of positive selection along particular Quercus lineages associated with the leaf traits and climate characteristics. In parallel, the DT model identified amino acid replacements at sites 95, 219, 262 and 328 being associated with the leaf traits and climate characteristics, exhibiting partial overlap with the results obtained using PAML.

Introgression Threatens the Genetic Diversity of Quercus austrocochinchinensis (Fagaceae), an Endangered Oak: A Case Inferred by Molecular Markers

Natural introgression can cause negative effects where rare species experience genetic assimilation and invade by their abundant congeners. Quercus austrocochinchinensis and Q. kerrii (subgenus Cyclobalanopsis) are a pair of closely related species in the Indo-China area. Morphological intermediates of the two species have been reported in this region. In this study, we used AFLP, SSR and two key leaf morphological diagnostic traits to study the two Q. austrocochinchinensis populations, two pure Q. kerrii and two putative hybrid populations in China. Rates of individual admixture were examined using the Bayesian clustering programs STRUCTURE and NewHybrids, with no a priori species assignment. In total, we obtained 151 SSR alleles and 781 polymorphic loci of AFLP markers. Population differentiation inferred by SSR and AFLP was incoherent with recognized species boundaries. Bayesian admixture analyses and principal coordinate analysis identified more hybrids and backcrossed individuals than morphological intermediates in the populations. SSR inferred a wide genetic assimilation in Q. austrocochinchinensis, except for subpopulation D2 in the core area of Xi-Shuang-Ban-Na Nature Reserve (XSBN). However, AFLP recognized more Q. austrocochinchinensis purebreds than SSR. Analysis using NewHybrids on AFLP data indicated that these hybridized individuals were few F2 and predominantly backcrosses with both parental species. All these evidences indicate the formation of a hybrid swarm at XSBN where the two species co-exist. Both AFLP and SSR recognized that the core protected area of XSBN (D2) has a high percentage of Q. austrocochinchinensis purebreds and a unique germplasm. The Hainan population and the other subpopulations of XSBN of the species might have lost their genetic integrity. Our results revealed a clear genetic differentiation in the populations and subpopulations of Q. austrocochinchinensis and ongoing introgression between Q. austrocochinchinensis and Q. kerrii at the disturbed contact areas. Combining the results from genetic and morphological analyses, the conservation of subpopulation D2 should be prioritized. Conservation and restoration of the integrity of tropical ravine rainforest is an important long-term goal for the successful conservation of Q. austrocochinchinensis. The fine-scale landscape might play an essential role in shaping the spatial patterns of hybridization. Further studies are needed to evaluate these patterns and dynamics.

Footprints of divergent selection in natural populations of Castanopsis fargesii (Fagaceae)

Given predicted rapid climate change, an understanding of how environmental factors affect genetic diversity in natural populations is important. Future selection pressures are inherently unpredictable, so forest management policies should maintain both overall diversity and identify genetic markers associated with the environmental factors expected to change most rapidly, like temperature and rainfall. In this study, we genotyped 648 individuals in 28 populations of Castanopsis fargesii (Fagaceae) using 32 expressed sequence tag (EST)-derived microsatellite markers. After removing six loci that departed from Hardy-Weinberg equilibrium, we measured genetic variation, population structure and identified candidate loci putatively under selection by temperature and precipitation. We found that C. fargesii populations possessed high genetic diversity and moderate differentiation among them, indicating predominant outcrossing and few restrictions to gene flow. These patterns reduce the possible impact of stochastic effects or the influence of genetic isolation. Clear footprints of divergent selection at four loci were discovered. Frequencies of five alleles at these loci were strongly correlated with environmental factors, particularly extremes in precipitation. These alleles varied from being near fixation at one end of the gradient to being completely absent at the other. Our study species is an important forest tree in the subtropical regions of China and could have a major role in future management and reforestation plans. Our results demonstrate that the gene flow is widespread and abundant in natural populations, maintaining high diversity, while diversifying selection is acting on specific genomic regions.

Genetic differentiation and genetic diversity of Castanopsis (Fagaceae), the dominant tree species in Japanese broadleaved evergreen forests, revealed by analysis of EST-associated microsatellites

The broadleaved evergreen forests of the East Asian warm temperate zone are characterised by their high biodiversity and endemism, and there is therefore a need to extend our understanding of its genetic diversity and phylogeographic patterns. Castanopsis (Fagaceae) is one of the dominant tree species in the broadleaved evergreen forests of Japan. In this study we investigate the genetic diversity, genetic structure and leaf epidermal morphology of 63 natural populations of C. sieboldii and C. cuspidata, using 32 Expressed Sequence Tag associated microsatellites. The overall genetic differentiation between populations was low (G_ST = 0.069 in C. sieboldii and G_ST = 0.057 in C. cuspidata). Neighbor-joining tree and Bayesian clustering analyses revealed that the populations of C. sieboldii and C. cuspidata were genetically clearly differentiated, a result which is consistent with the morphology of their epidermal cell layers. This suggests that C. sieboldii and C. cuspidata should be treated as independent species, although intermediate morphologies are often observed, especially at sites where the two species coexist. The higher level of genetic diversity observed in the Kyushu region (for both species) and the Ryukyu Islands (for C. sieboldii) is consistent with the available fossil pollen data for Castanopsis-type broadleaved evergreen trees during the Last Glacial Maximum and suggests the existence of refugia for Castanopsis forests in southern Japan. Within the C. sieboldii populations, Bayesian clustering analyses detected three clusters, in the western and eastern parts of the main islands and in the Ryukyu Islands. The west-east genetic differentiation observed for this species in the main islands, a pattern which is also found in several plant and animal species inhabiting Castanopsis forests in Japan, suggests that they have been isolated from each other in the western and eastern populations for an extended period of time, and may imply the existence of eastern refugia.

Revisiting an important component of plant genomes: microsatellites

Microsatellites are some of the most highly variable repetitive DNA tracts in genomes. Few studies focus on whether the characteristic instability of microsatellites is linked to phenotypic effects in plants. We summarise recent data to investigate how microsatellite variations affect gene expression and hence phenotype. We discuss how the basic characteristics of microsatellites may contribute to phenotypic effects. In summary, microsatellites in plants are universal and highly mutable, they coexist and coevolve with transposable elements, and are under selective pressure. The number of motif nucleotides, the type of motif and transposon activity all contribute to the nonrandom generation and decay of microsatellites, and to conservation and distribution biases. Although microsatellites are generated by accident, they mature through responses to environmental change before final decay. This process is mediated by organism adjustment mechanisms, which maintain a balance between birth versus death and growth versus decay in microsatellites. Close relationships also exist between the physical structure, variation and functionality of microsatellites: in most plant species, sequences containing microsatellites are associated with catalytic activity and binding functions, are expressed in the membrane and organelles, and participate in the developmental and metabolic processes. Microsatellites contribute to genome structure and functional plasticity, and may be considered to promote species evolution in plants in response to environmental changes. In conclusion, the generation, loss, functionality and evolution of microsatellites can be related to plant gene expression and functional alterations. The effect of microsatellites on phenotypic variation may be as significant in plants as it is in animals.

Three EST-SSR markers associated with QTL for the growth of the clam Meretrix meretrix revealed by selective genotyping

The clam Meretrix meretrix is a member of widely cultured, commercially important clams. A marker-trait association analysis was performed using expressed sequence tag (EST) simple sequence repeat (SSR) markers for marker-assisted selection in M. meretrix. Three markers, MM1272, MM2034, and MM7721, were found to be significantly associated with quantitative trait loci (QTLs) controlling shell length (P < 0.0001) in clams of a fast-growing population (JSF) and a control population (JSC). The 144-bp allele of MM1272, the 154-bp allele of MM2034, and the 152- and 165-bp alleles of MM7721 showed a significantly higher frequency in the JSF population (17.65, 36.41, 28.67, and 29.33 %) than in the JSC population (4.65, 8.33, 3.47, and 5.56 %). The three markers showed lower values for the number of alleles and observed heterozygosity as well as a higher proportion of homozygotes in JSF than in JSC population. The three markers have been further confirmed in the high and low tails of another population (09G₃SPSB); similarly, lower values for the number of alleles and observed heterozygosity as well as a higher proportion of homozygotes were found in 09G₃SPSB(H). The putative functions of the three gene fragments containing MM1272, MM2034, and MM7721 also suggested that the three SSR-containing genes might be involved in growth of M. meretrix. All the results suggest that the three EST-SSR markers associated with growth QTLs would be useful for marker-assisted selection in M. meretrix breeding.

A second generation framework for the analysis of microsatellites in expressed sequence tags and the development of EST-SSR markers for a conifer, Cryptomeria japonica

Background

Microsatellites or simple sequence repeats (SSRs) in expressed sequence tags (ESTs) are useful resources for genome analysis because of their abundance, functionality and polymorphism. The advent of commercial second generation sequencing machines has lead to new strategies for developing EST-SSR markers, necessitating the development of bioinformatic framework that can keep pace with the increasing quality and quantity of sequence data produced. We describe an open scheme for analyzing ESTs and developing EST-SSR markers from reads collected by Sanger sequencing and pyrosequencing of sugi (Cryptomeria japonica).

Results

We collected 141,097 sequence reads by Sanger sequencing and 1,333,444 by pyrosequencing. After trimming contaminant and low quality sequences, 118,319 Sanger and 1,201,150 pyrosequencing reads were passed to the MIRA assembler, generating 81,284 contigs that were analysed for SSRs. 4,059 SSRs were found in 3,694 (4.54%) contigs, giving an SSR frequency lower than that in seven other plant species with gene indices (5.4–21.9%). The average GC content of the SSR-containing contigs was 41.55%, compared to 40.23% for all contigs. Tri-SSRs were the most common SSRs; the most common motif was AT, which was found in 655 (46.3%) di-SSRs, followed by the AAG motif, found in 342 (25.9%) tri-SSRs. Most (72.8%) tri-SSRs were in coding regions, but 55.6% of the di-SSRs were in non-coding regions; the AT motif was most abundant in 3′ untranslated regions. Gene ontology (GO) annotations showed that six GO terms were significantly overrepresented within SSR-containing contigs. Forty–four EST-SSR markers were developed from 192 primer pairs using two pipelines: read2Marker and the newly-developed CMiB, which combines several open tools. Markers resulting from both pipelines showed no differences in PCR success rate and polymorphisms, but PCR success and polymorphism were significantly affected by the expected PCR product size and number of SSR repeats, respectively. EST-SSR markers exhibited less polymorphism than genomic SSRs.

Conclusions

We have created a new open pipeline for developing EST-SSR markers and applied it in a comprehensive analysis of EST-SSRs and EST-SSR markers in C. japonica. The results will be useful in genomic analyses of conifers and other non-model species.

Characterization and comparison of gene-based simple sequence repeats across Brassica species

Simple sequence repeats (SSRs) are important components of eukaryotic genomes and may play an important role in regulating gene expression. However, the characteristics of genic SSRs and the effect of interspecific hybridization and polyploidization on genic SSRs seem not to have received desired attention in terms of scientific investigations. To determine the features of genic SSRs and elucidate their role in polyploidization process of the Brassica family, we identified SSRs in Plant Genome Database-assembled unique transcripts (PUTs) of Brassica species. A higher density of SSRs and a greater number of compound motif SSRs and mononucleotide motif types with large average number of repeats were detected in allotetraploid Brassica napus than in the diploid parental species (Brassica rapa and Brassica oleracea). In addition, a greater proportion of SSR-PUTs were found to be associated with the stress response and developmental processes in B. napus than in the parents. A negative correlation between the repeat number and the motif type and the total length, and a positive correlation between the repeat number and the total length of SSRs were observed. PUT-SSR might be generated from A/T-rich regions. The successful development of 123 pairs of SSR primers for Brassica PUTs showed that SSR-PUTs could be exploited as gene-based SSR functional markers for application in Brassica breeding. These results indicate that interspecific hybridization and polyploidization could trigger the amplification of SSRs, and long SSRs might become shorter to enable the plant to adapt to environmental and artificial selection.

Phylogeography of Phytophagous Weevils and Plant Species in Broadleaved Evergreen Forests: A Congruent Genetic Gap between Western and Eastern Parts of Japan

The Quaternary climate cycles played an important role in shaping the distribution of biodiversity among current populations, even in warm-temperate zones, where land was not covered by ice sheets. We focused on the Castanopsis-type broadleaved evergreen forest community in Japan, which characterizes the biodiversity and endemism of the warm-temperate zone. A comparison of the phylogeographic patterns of three types of phytophagous weevils associated with Castanopsis (a host-specific seed predator, a generalist seed predator, and a host-specific leaf miner) and several other plant species inhabiting the forests revealed largely congruent patterns of genetic differentiation between western and eastern parts of the main islands of Japan. A genetic gap was detected in the Kii Peninsula to Chugoku-Shikoku region, around the Seto Inland Sea. The patterns of western-eastern differentiation suggest past fragmentation of broadleaved evergreen forests into at least two separate refugia consisting of the southern parts of Kyushu to Shikoku and of Kii to Boso Peninsula. Moreover, the congruent phylogeographic patterns observed in Castanopsis and the phytophagous insect species imply that the plant-herbivore relationship has been largely maintained since the last glacial periods. These results reinforce the robustness of the deduced glacial and postglacial histories of Castanopsis-associated organisms.

Developing new SSR markers from ESTs of pea (Pisum sativum L.)

The development of expressed sequence tags (ESTs) from pea has provided a useful source for mining novel simple sequence repeat (SSR) markers. In the present research, in order to find EST-derived SSR markers, 18 552 pea ESTs from the National Center for Biotechnology Information (NCBI) database were downloaded and assembled into 10 086 unigenes. A total of 586 microsatellites in 530 unigenes were identified, indicating that merely 5.25% of sequences contained SSRs. The most abundant SSRs within pea were tri-nucleotide repeat motifs, and among all the tri-nucleotide repeats, the motif GAA was the most abundant type. In total, 49 SSRs were used for primer design. EST-SSR loci were subsequently screened on 10 widely adapted varieties in China. Of these, nine loci showed polymorphic profiles that revealed two to three alleles per locus. The polymorphism information content value ranged from 0.18 to 0.58 with an average of 0.41. Furthermore, transferable analysis revealed that some of these loci showed transferability to faba bean. Because of their polymorphism and transferability, these nine novel EST-SSRs will be valuable tools for marker-assisted breeding and comparative mapping of pea in the future.