Preliminary Genomic Characterization of Ten Hardwood Tree Species from Multiplexed Low Coverage Whole Genome Sequencing

Applying molecular and genetic methods to trees and their fungal communities

Forests provide invaluable economic, ecological, and social services. At the same time, they are exposed to several threats, such as fragmentation, changing climatic conditions, or increasingly destructive pests and pathogens. Trees, the inherent species of forests, cannot be viewed as isolated organisms. Manifold (micro)organisms are associated with trees playing a pivotal role in forest ecosystems. Of these organisms, fungi may have the greatest impact on the life of trees. A multitude of molecular and genetic methods are now available to investigate tree species and their associated organisms. Due to their smaller genome sizes compared to tree species, whole genomes of different fungi are routinely compared. Such studies have only recently started in forest tree species. Here, we summarize the application of molecular and genetic methods in forest conservation genetics, tree breeding, and association genetics as well as for the investigation of fungal communities and their interrelated ecological functions. These techniques provide valuable insights into the molecular basis of adaptive traits, the impacts of forest management, and changing environmental conditions on tree species and fungal communities and can enhance tree-breeding cycles due to reduced time for field testing. It becomes clear that there are multifaceted interactions among microbial species as well as between these organisms and trees. We demonstrate the versatility of the different approaches based on case studies on trees and fungi. KEY POINTS: • Current knowledge of genetic methods applied to forest trees and associated fungi. • Genomic methods are essential in conservation, breeding, management, and research. • Important role of phytobiomes for trees and their ecosystems.

Development of Microsatellite Markers for Tanacetum cinerariifolium (Trevis.) Sch. Bip., a Plant with a Large and Highly Repetitive Genome

Dalmatian pyrethrum (Tanacetum cinerariifolium (Trevis.) Sch. Bip.) is an outcrossing plant species (2n = 18) endemic to the eastern Adriatic coast and source of the natural insecticide pyrethrin. Due to the high repeatability and large genome (1C-value = 9.58 pg) our previous attempts to develop microsatellite markers using the traditional method were unsuccessful. Now we have used Illumina paired-end whole genome sequencing and developed a specific procedure to obtain useful microsatellite markers. A total of 796,130,142 high-quality reads (approx. 12.5× coverage) were assembled into 6,909,675 contigs using two approaches (de novo assembly and joining of overlapped pair-end reads). A total of 31,380 contigs contained one or more microsatellite sequences, of which di-(59.7%) and trinucleotide (25.9%) repeats were the most abundant. Contigs containing microsatellites were filtered according to various criteria to achieve better yield of functional markers. After two rounds of testing, 17 microsatellite markers were developed and characterized in one natural population. Twelve loci were selected for preliminary genetic diversity analysis of three natural populations. Neighbor-joining tree, based on the proportion of shared alleles distances, grouped individuals into clusters according to population affiliation. The availability of codominant SSR markers will allow analysis of genetic diversity and structure of natural Dalmatian pyrethrum populations as well as identification of breeding lines and cultivars.

Strategies of tolerance reflected in two North American maple genomes

The first chromosome‐scale assemblies for North American members of the Acer genus, sugar maple (Acer saccharum) and boxelder (Acer negundo), as well as transcriptomic evaluation of the abiotic stress response in A. saccharum are reported. This integrated study describes in‐depth aspects contributing to each species' approach to tolerance and applies current knowledge in many areas of plant genome biology with Acer physiology to help convey the genomic complexities underlying tolerance in broadleaf tree species.

RNAi of AGAMOUS genes in sweetgum alters reproductive organ identity and decreases fruit persistence

Sweetgums (Liquidambar), members of the family Altingiaceae (Altingiales), have inflorescences and floral organs that are distinctive in structure compared with other angiosperms in which the roles of floral homeotic genes have been studied. To begin to understand the role of AGAMOUS (AG)-a floral homeotic gene that has a major role in stamen and carpel development-in development of the monosexual flowers of sweetgum, we used RNAi to reduce the expression of two members of the AG subfamily. Because AG suppression should induce floral sterility, RNAi might also provide a tool to mitigate the risks of invasiveness-and to reduce the production of its nuisance fruits or allergenic pollen-when sweetgum is used as an exotic shade or forest tree. We tested 33 independent transgenic events and non-transgenic controls during 10 years in the field. The RNAi-AG sweetgum trees maintained normal growth, phenology, and vivid fall coloration during the 10 years of study, but 8 insertion events had highly modified inflorescence and floral morphology. The modified flowers had anthers and carpels that were converted to flat leaf-like structures lacking pollen grains and ovules, respectively. The female inflorescences developed into dry papery structures that failed to produce seeds. These infructescences were smaller than control infructescences, and lost a greater percentage of biomass in a controlled decay assay. RNAi against AG genes was highly effective at impairing fertility and modifying reproductive development without significant vegetative effects in sweetgum and gave phenotypes distinct from, but similar to, that of AG loss of function in other angiosperms.

Repetitive genomic elements in Campomanesia xanthocarpa: prospection, characterization and cross amplification of molecular markers

Repetitive genomic elements were prospected in Campomanesia xanthocarpa, aiming to characterize these elements in a non-model plant species and to develop species-specific microsatellite markers. Approximately 4.12% of the partial genome of C. xanthocarpa is composed of repetitive elements, being retrotransposons the most widely represented. A total of nine polymorphic microsatellite markers were obtained: four nuclear-neutral, two nuclear EST, two plastidial and one mitochondrial. Levels of population genetic diversity of four natural populations of C. xanthocarpa were characterized using these markers. In addition, the cross-species amplification of the microsatellite markers was tested in seven species of tribe Myrteae (Myrtaceae). The characterized microsatellite markers revealed low to moderate levels of genetic diversity (expected heterozygosity range: 0.33-0.57; observed heterozygosity: 0.26-0.74 and number of alleles: 2.25-4.25). Cross-species amplification was successful for all loci, except Cxant76. These nine markers will contribute for studies on genetic diversity, gene flow, plant selection and breeding of this species, towards the conservation of natural populations, as well as its commercial use.

High-throughput sequencing analysis of Eugenia uniflora: insights into repetitive DNA, gene content and potential biotechnological applications

In this study, the Ion Torrent PGM platform was employed to conduct a whole-genome sequencing analysis of Eugenia uniflora L., aiming to generate new genomic information for this non-model species. The whole-genome set of assembled sequences (WGSAS) generated 2601 contigs, with a cumulative length of 3.15 Mbp. About 2.60% of the WGSAS was characterized by repetitive sequences, while 71.66% of the WGSAS was represented by predicted genes. A total of 147 metabolic pathways related to enzymes characterized in the E. uniflora WGSAS were detected. Proteins related to antifungal activity and to bacterium and fungi defense response as well as transcriptions factors related to abiotic stress responses were also identified. This sequencing provides important genomic resources for future breeding and selection regarding plants with higher production of such metabolites, as well as source of genes for genetic engineering of cultivated species, supporting a wide application of marker-assisted and genomic selection, towards the establishment of commercial orchards with improved cultivars of E. uniflora.

Discovery and characterization of SSR markers in Eugenia uniflora L. (Myrtaceae) using low coverage genome sequencing

Eugenia uniflora L. (Myrtaceae) is a tree species widely distributed in South America suffering the effects of the exploitation of natural populations. In this study, we employed low coverage sequencing of the E. uniflora genome for mining of SSR markers. The de novo assembly generated 2,601 contigs with an average length of 1139 bp and spans 3.15 Mb. A total of 76 dimer, 33 trimer and two compound SSR loci were identified. Twelve selected SSR loci were employed to genotype 30 individuals from two natural populations. A total of 73 alleles were detected (mean A= 6.1) were observed, the mean effective number of alleles was Ae = 3.91, mean Ho was 0.23 and mean HE was 0.70). The mean Wright's within population fixation index was FIS = 0.66 and significant deviation of HWE was observed in all loci, except one. The FST between populations equaled 0.27. The levels of genetic diversity and structure estimated with these 12 SSR markers are in accordance with data from genetics studies performed on other tree species of the Pampa biome, presenting moderate to high polymorphism and may be employed in studies of species conservation measures and breeding programs.

Microsatellite Borders and Micro-sequence Conservation in Juglans

Walnuts (Juglans spp.) are economically important nut and timber species with a worldwide distribution. Using the published Persian walnut genome as a reference for the assembly of short reads from six Juglans species and several interspecific hybrids, we identified simple sequence repeats in 12 Juglans nuclear and organellar genomes. The genome-wide distribution and polymorphisms of nuclear and organellar microsatellites (SSRs) for most Juglans genomes have not been previously studied. We compared the frequency of nuclear SSR motifs and their lengths across Juglans, and identified section-specific chloroplast SSR motifs. Primer pairs were designed for more than 60,000 SSR-containing sequences based on alignment against assembled scaffold sequences. Of the >60,000 loci, 39,000 were validated by e-PCR using unique primer pairs. We identified primers containing 100% sequence identity in multiple species. Across species, sequence identity in the SSR-flanking regions was generally low. Although SSRs are common and highly dispersed in the genome, their flanking sequences are conserved at about 90 to 95% identity within Juglans and within species. In a few rare cases, flanking sequences are identical across species of Juglans. This comprehensive report of nuclear and organellar SSRs in Juglans and the generation of validated SSR primers will be a useful resource for future genetic analyses, walnut breeding programs, high-level taxonomic evaluations, and genomic studies in Juglandaceae.

The first genetic linkage map for Fraxinus pennsylvanica and syntenic relationships with four related species

The genetic linkage map for green ash (Fraxinus pennsylvanica) contains 1201 DNA markers in 23 linkage groups spanning 2008.87cM. The green ash map shows stronger synteny with coffee than tomato. Green ash (Fraxinus pennsylvanica) is an outcrossing, diploid (2n = 46) hardwood tree species, native to North America. Native ash species in North America are being threatened by the rapid spread of the emerald ash borer (EAB, Agrilus planipennis), an invasive pest from Asia. Green ash, the most widely distributed ash species, is severely affected by EAB infestation, yet few genomic resources for genetic studies and improvement of green ash are available. In this study, a total of 5712 high quality single nucleotide polymorphisms (SNPs) were discovered using a minimum allele frequency of 1% across the entire genome through genotyping-by-sequencing. We also screened hundreds of genomic- and EST-based microsatellite markers (SSRs) from previous de novo assemblies (Staton et al., PLoS ONE 10:e0145031, 2015; Lane et al., BMC Genom 17:702, 2016). A first genetic linkage map of green ash was constructed from 90 individuals in a full-sib family, combining 2719 SNP and 84 SSR segregating markers among the parental maps. The consensus SNP and SSR map contains a total of 1201 markers in 23 linkage groups spanning 2008.87 cM, at an average inter-marker distance of 1.67 cM with a minimum logarithm of odds of 6 and maximum recombination fraction of 0.40. Comparisons of the organization the green ash map with the genomes of asterid species coffee and tomato, and genomes of the rosid species poplar and peach, showed areas of conserved gene order, with overall synteny strongest with coffee.

Characterization of Plastidial and Nuclear SSR Markers for Understanding Invasion Histories and Genetic Diversity of Schinus molle L

Invasive plant species are expected to display high dispersal capacity but low levels of genetic diversity due to the founder effect occurring at each invasion episode. Understanding the history of invasions and the levels of genetic diversity of such species is an important task for planning management and monitoring strategy for these events. Peruvian Peppertree (Schinus molle L.) is a pioneer tree species native from South America which was introduced in North America, Europe and Africa, becoming a threat to these non-native habitats. In this study, we report the discovery and characterization of 17 plastidial (ptSSR) and seven nuclear (nSSR) markers for S. molle based on low-coverage whole-genome sequencing data acquired through next-generation sequencing. The markers were tested in 56 individuals from two natural populations sampled in the Brazilian Caatinga and Pampa biomes. All loci are moderately to highly polymorphic and revealed to be suitable for genetic monitoring of new invasions, for understanding the history of old invasions, as well as for genetic studies of native populations in their natural occurrence range and of orchards established with commercial purposes.

A reassessment of the genome size-invasiveness relationship in reed canarygrass (Phalaris arundinacea)

Background and Aims

Genome size is hypothesized to affect invasiveness in plants. Key evidence comes from a previous study of invasive eastern North American populations of the grass Phalaris arundinacea: invasive genotypes with smaller genomes had higher growth rates, and genome sizes were smaller in the invasive vs. native range. This study aimed to re-investigate those patterns by examining a broader range of North American populations and by employing the modern best-practice protocol for plant genome size estimation in addition to the previously used protocol.

Methods

Genome sizes were measured using both internal and pseudo-internal standardization protocols for 20 invasive and nine native range accessions of P. arundinacea. After a round of vegetative propagation to reduce maternal environmental effects, growth (stem elongation) rates of these accessions were measured in the greenhouse.

Key Results

Using the best-practice protocol, there was no evidence of a correlation between genome size and growth rates (P = 0.704), and no evidence for differences in genome sizes of invasive and native range accessions (P > 0.353). However, using the older genome size estimation protocol, both relationships were significant (reproducing the results of the previous study).

Conclusions

Genome size reduction has not driven increased invasiveness in a broad sample of North American P. arundinacea. Further, inappropriate genome size estimation techniques can create spurious correlations between genome size and plant traits such as growth rate. Valid estimation is vital to progress in understanding the potentially widespread effects of genome size on biological processes and patterns.