Species in the genus Armillaria (fungi, basidiomycota) are well-known as saprophytes and pathogens on plants. Many of them cause white-rot root disease in diverse woody plants worldwide. Mitochondrial genomes (mitogenomes) are widely used in evolutionary and population studies, but despite the importance and wide distribution of Armillaria, the complete mitogenomes have not previously been reported for this genus. Meanwhile, the well-supported phylogeny of Armillaria species provides an excellent framework in which to study variation in mitogenomes and how they have evolved over time. Here we completely sequenced, assembled, and annotated the circular mitogenomes of four species: A. borealis, A. gallica, A. sinapina, and A. solidipes (116,443, 98,896, 103,563, and 122,167 bp, respectively). The variation in mitogenome size can be explained by variable numbers of mobile genetic elements, introns, and plasmid-related sequences. Most Armillaria introns contained open reading frames (ORFs) that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Insertions of mobile elements were also evident as fragments of plasmid-related sequences in Armillaria mitogenomes. We also found several truncated gene duplications in all four mitogenomes. Our study showed that fungal mitogenomes have a high degree of variation in size, gene content, and genomic organization even among closely related species of Armillara. We suggest that mobile genetic elements invading introns and intergenic sequences in the Armillaria mitogenomes have played a significant role in shaping their genome structure. The mitogenome changes we describe here are consistent with widely accepted phylogenetic relationships among the four species.

The genetic diversity of SSR markers was studied in six populations of Siberian larch, two popu-lations of Gmelin larch, and four populations of Cajander larch. Seven pairs of nuclear microsatellite loci were used for this analysis. 103 allelic variants were detected in 365 individuals of three species of larch. According to the AMOVA results, the variability proportion that characterizes the differences between three Larix species was twice as higher (15%) than the proportion that accounts for among-population differences within the species (6%). The differentiation of the populations of three species of larch based on SSR markers exceeded 12% (F _ST = 0.121). A significant correlation of the genetic distances with the geographic distances between populations was found (r = 0.835, P = 0.01).

We evaluated the population structure of the bog and dry land populations of the Siberian pine Pinus sibirica ( P. sibrica) in Western Siberia using nuclear genome markers. Six pairs of nuclear...

Forest trees growing in high altitude conditions offer a convenient model for studying adaptation processes. They are subject to a whole range of adverse factors that are likely to cause local adaptation and related genetic changes. Siberian larch (Larix sibirica Ledeb.), whose distribution covers different altitudes, makes it possible to directly compare lowland with highland populations. This paper presents for the first time the results of studying the genetic differentiation of Siberian larch populations, presumably associated with adaptation to the altitudinal gradient of climatic conditions, based on a joint analysis of altitude and six other bioclimatic variables, together with a large number of genetic markers, single nucleotide polymorphisms (SNPs), obtained from double digest restriction-site-associated DNA sequencing (ddRADseq). In total, 25,143 SNPs were genotyped in 231 trees. In addition, a dataset of 761 supposedly selectively neutral SNPs was assembled by selecting SNPs located outside coding regions in the Siberian larch genome and mapped to different contigs. The analysis using four different methods (PCAdapt, LFMM, BayeScEnv and RDA) revealed 550 outlier SNPs, including 207 SNPs whose variation was significantly correlated with the variation of some of environmental factors and presumably associated with local adaptation, including 67 SNPs that correlated with altitude based on either LFMM or BayeScEnv and 23 SNPs based on both of them. Twenty SNPs were found in the coding regions of genes, and 16 of them represented non-synonymous nucleotide substitutions. They are located in genes involved in the processes of macromolecular cell metabolism and organic biosynthesis associated with reproduction and development, as well as organismal response to stress. Among these 20 SNPs, nine were possibly associated with altitude, but only one of them was identified as associated with altitude by all four methods used in the study, a nonsynonymous SNP in scaffold_31130 in position 28092, a gene encoding a cell membrane protein with uncertain function. Among the studied populations, at least two main groups (clusters), the Altai populations and all others, were significantly genetically different according to the admixture analysis based on any of the three SNP datasets as follows: 761 supposedly selectively neutral SNPs, all 25,143 SNPs and 550 adaptive SNPs. In general, according to the AMOVA results, genetic differentiation between transects or regions or between population samples was relatively low, although statistically significant, based on 761 neutral SNPs (F_ST = 0.036) and all 25,143 SNPs (F_ST = 0.017). Meanwhile, the differentiation based on 550 adaptive SNPs was much higher (F_ST = 0.218). The data showed a relatively weak but highly significant linear correlation between genetic and geographic distances (r = 0.206, p = 0.001).

The recent release of the nuclear, chloroplast and mitochondrial genome assemblies of Siberian larch (Larix sibirica Ledeb.), one of the most cold-resistant tree species in the only deciduous genus of Pinaceae, with seasonal senescence and a rot-resistant valuable timber widely used in construction, greatly contributed to the development of genomic resources for the larch genus. Here, we present an extensive repeatome analysis and the first annotation of the draft nuclear Siberian larch genome assembly. About 66% of the larch genome consists of highly repetitive elements (REs), with the likely wave of retrotransposons insertions into the larch genome estimated to occur 4–5 MYA. In total, 39,370 gene models were predicted, with 87% of them having homology to the Arabidopsis-annotated proteins and 78% having at least one GO term assignment. The current state of the genome annotations allows for the exploration of the gymnosperm and angiosperm species for relative gene abundance in different functional categories. Comparative analysis of functional gene categories across different angiosperm and gymnosperm species finds that the Siberian larch genome has an overabundance of genes associated with programmed cell death (PCD), autophagy, stress hormone biosynthesis and regulatory pathways; genes that may play important roles in seasonal senescence and stress response to extreme cold in larch. Despite being incomplete, the draft assemblies and annotations of the conifer genomes are at a point of development where they now represent a valuable source for further genomic, genetic and population studies.

Dendrogenomics is a new interdisciplinary approach that allows joint analysis of dendrological and genomic data and opens up new ways to study the temporal dynamics of forest treelines, delineate spatial and temporal population structures, decipher individual tree responses to abiotic and biotic stresses, and evaluate the adaptive genetic potential of forest tree populations. These data are needed for the prediction of climate change effects and mitigation of the negative effects. We present here an association analysis of the variation of 27 individual tree traits, including adaptive dendrophenotypes reflecting the individual responses of trees to drought stress, such as the resistance (Rt), recovery (Rc), resilience (Rs), and relative resilience (RRs) indexes measured in 136 Siberian larch trees in 5 populations in the foothills of the Batenevsky Ridge (Kuznetsk Alatau, Republic of Khakassia, Russia), with variation of 9742 SNPs genotyped using ddRADseq in the same trees. The population structure of five closely located Siberian larch populations was relatively weak (F_ST = 0.018). We found that the level of individual heterozygosity positively correlated with the Rc and RR indices for the five studied drought periods and partly with the Rs indices for three drought periods. It seems that higher individual heterozygosity improves the adaptive capabilities of the tree. We also discovered a significant negative relationship between individual heterozygosity and the Rt index in four out of five periods, which means that growth slows down during droughts more in trees with higher individual heterozygosity and is likely associated with energy and internal resource reallocation toward more efficient water and energy usage and optimization of larch growth during drought years. We found 371 SNPs with potentially adaptive variations significantly associated with the variation of adaptive dendrophenotypes based on all three different methods of association analysis. Among them, 26 SNPs were located in genomic regions carrying functional genes: 21 in intergenic regions and 5 in gene-coding regions. Based on the obtained results, it can be assumed that these populations of Siberian larch have relatively high standing adaptive genetic variation and adaptive potential underlying the adaptations of larch to various climatic conditions.