High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar.

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Genetics of flight in spongy moths (Lymantria dispar ssp.): functionally integrated profiling of a complex invasive trait.

BACKGROUND: Flight can drastically enhance dispersal capacity and is a key trait defining the potential of exotic insect species to spread and invade new habitats. The phytophagous European spongy moths (ESM, Lymantria dispar dispar) and Asian spongy moths (ASM; a multi-species group represented here by L. d. asiatica and L. d. japonica), are globally invasive species that vary in adult female flight capability-female ASM are typically flight capable, whereas female ESM are typically flightless. Genetic markers of flight capability would supply a powerful tool for flight profiling of these species at any intercepted life stage. To assess the functional complexity of spongy moth flight and to identify potential markers of flight capability, we used multiple genetic approaches aimed at capturing complementary signals of putative flight-relevant genetic divergence between ESM and ASM: reduced representation genome-wide association studies, whole genome sequence comparisons, and developmental transcriptomics. We then judged the candidacy of flight-associated genes through functional analyses aimed at addressing the proximate demands of flight and salient features of the ecological context of spongy moth flight evolution. RESULTS: Candidate gene sets were typically non-overlapping across different genetic approaches, with only nine gene annotations shared between any pair of approaches. We detected an array of flight-relevant functional themes across gene sets that collectively suggest divergence in flight capability between European and Asian spongy moth lineages has coincided with evolutionary differentiation in multiple aspects of flight development, execution, and surrounding life history. Overall, our results indicate that spongy moth flight evolution has shaped or been influenced by a large and functionally broad network of traits. CONCLUSIONS: Our study identified a suite of flight-associated genes in spongy moths suited to exploration of the genetic architecture and evolution of flight, or validation for flight profiling purposes. This work illustrates how complementary genetic approaches combined with phenotypically targeted functional analyses can help to characterize genetically complex traits.

Does urbanisation lead to parallel demographic shifts across the world in a cosmopolitan plant?

Urbanisation is occurring globally, leading to dramatic environmental changes that are altering the ecology and evolution of species. In particular, the expansion of human infrastructure and the loss and fragmentation of natural habitats in cities is predicted to increase genetic drift and reduce gene flow by reducing the size and connectivity of populations. Alternatively, the 'urban facilitation model' suggests that some species will have greater gene flow into and within cities leading to higher diversity and lower differentiation in urban populations. These alternative hypotheses have not been contrasted across multiple cities. Here, we used the genomic data from the GLobal Urban Evolution project (GLUE), to study the effects of urbanisation on non-adaptive evolutionary processes of white clover (Trifolium repens) at a global scale. We found that white clover populations presented high genetic diversity and no evidence of reduced Ne linked to urbanisation. On the contrary, we found that urban populations were less likely to experience a recent decrease in effective population size than rural ones. In addition, we found little genetic structure among populations both globally and between urban and rural populations, which showed extensive gene flow between habitats. Interestingly, white clover displayed overall higher gene flow within urban areas than within rural habitats. Our study provides the largest comprehensive test of the demographic effects of urbanisation. Our results contrast with the common perception that heavily altered and fragmented urban environments will reduce the effective population size and genetic diversity of populations and contribute to their isolation.

Efficient genomics-based 'end-to-end' selective tree breeding framework.

Since their initiation in the 1950s, worldwide selective tree breeding programs followed the recurrent selection scheme of repeated cycles of selection, breeding (mating), and testing phases and essentially remained unchanged to accelerate this process or address environmental contingencies and concerns. Here, we introduce an "end-to-end" selective tree breeding framework that: (1) leverages strategically preselected GWAS-based sequence data capturing trait architecture information, (2) generates unprecedented resolution of genealogical relationships among tested individuals, and (3) leads to the elimination of the breeding phase through the utilization of readily available wind-pollinated (OP) families. Individuals' breeding values generated from multi-trait multi-site analysis were also used in an optimum contribution selection protocol to effectively manage genetic gain/co-ancestry trade-offs and traits' correlated response to selection. The proof-of-concept study involved a 40-year-old spruce OP testing population growing on three sites in British Columbia, Canada, clearly demonstrating our method's superiority in capturing most of the available genetic gains in a substantially reduced timeline relative to the traditional approach. The proposed framework is expected to increase the efficiency of existing selective breeding programs, accelerate the start of new programs for ecologically and environmentally important tree species, and address climate-change caused biotic and abiotic stress concerns more effectively.

Chromosome-scale genome assembly and insights into the metabolome and gene regulation of leaf color transition in an important oak species, Quercus dentata.

Quercus dentata Thunb., a dominant forest tree species in northern China, has significant ecological and ornamental value due to its adaptability and beautiful autumn coloration, with color changes from green to yellow into red resulting from the autumnal shifts in leaf pigmentation. However, the key genes and molecular regulatory mechanisms for leaf color transition remain to be investigated. First, we presented a high-quality chromosome-scale assembly for Q. dentata. This 893.54 Mb sized genome (contig N50 = 4.21 Mb, scaffold N50 = 75.55 Mb; 2n = 24) harbors 31 584 protein-coding genes. Second, our metabolome analyses uncovered pelargonidin-3-O-glucoside, cyanidin-3-O-arabinoside, and cyanidin-3-O-glucoside as the main pigments involved in leaf color transition. Third, gene co-expression further identified the MYB-bHLH-WD40 (MBW) transcription activation complex as central to anthocyanin biosynthesis regulation. Notably, transcription factor (TF) QdNAC (QD08G038820) was highly co-expressed with this MBW complex and may regulate anthocyanin accumulation and chlorophyll degradation during leaf senescence through direct interaction with another TF, QdMYB (QD01G020890), as revealed by our further protein-protein and DNA-protein interaction assays. Our high-quality genome assembly, metabolome, and transcriptome resources further enrich Quercus genomics and will facilitate upcoming exploration of ornamental values and environmental adaptability in this important genus.

SubPhaser: a robust allopolyploid subgenome phasing method based on subgenome-specific k-mers.

With advanced sequencing technology, dozens of complex polyploid plant genomes have been characterized. However, for many polyploid species, their diploid ancestors are unknown or extinct, making it impossible to unravel the subgenomes and genome evolution directly. We developed a novel subgenome-phasing algorithm, SubPhaser, specifically designed for a neoallopolyploid or a homoploid hybrid. SubPhaser first searches for the subgenome-specific sequence (k-mer), then assigns homoeologous chromosomes into subgenomes, and further provides tools to annotate and investigate specific sequences. SubPhaser works well on neoallopolyploids and homoploid hybrids containing subgenome-specific sequences like wheat, but fails on autopolyploids lacking subgenome-specific sequences like alfalfa, indicating that SubPhaser can phase neoallopolyploid/homoploid hybrids with high accuracy, sensitivity and performance. This highly accurate, highly sensitive, ancestral data free chromosome phasing algorithm, SubPhaser, offers significant application value for subgenome phasing in neoallopolyploids and homoploid hybrids, and for the subsequent exploration of genome evolution and related genetic/epigenetic mechanisms.

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits.

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism's potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species' evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module "leaf size gradient" was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under "global regulation" in vegetative bud establishment and development. Lamina basal shape and intercalary veins' number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits' ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

In-depth transcriptome characterization uncovers distinct gene family expansions for Cupressus gigantea important to this long-lived species' adaptability to environmental cues.

BACKGROUND: Cupressus gigantea, a rare and endangered tree species with remarkable medicinal value, is endemic to the Tibetan Plateau. Yet, little is known about the underlying genetics of the unique ecological adaptability of this extremely long-lived conifer with a large genome size. Here, we present its first de novo and multi-tissue transcriptome in-depth characterization. RESULTS: We performed Illumina paired-end sequencing and RNA libraries assembly derived from terminal buds, male and female strobili, biennial leaves, and cambium tissues taken from adult C. gigantea. In total, large-scale high-quality reads were assembled into 101,092 unigenes, with an average sequence length of 1029 bp, and 6848 unigenes (6.77%) were mapped against the KEGG databases to identify 292 pathways. A core set of 41,373 genes belonging to 2412 orthologous gene families shared between C. gigantea and nine other plants was revealed. In addition, we identified 2515 small to larger-size gene families containing in total 9223 genes specific to C. gigantea, and enriched for gene ontologies relating to biotic interactions. We identified an important terpene synthases gene family expansion with its 121 putative members. CONCLUSIONS: This study presents the first comprehensive transcriptome characterization of C. gigantea. Our results will facilitate functional genomic studies to support genetic improvement and conservation programs for this endangered conifer.

Gene copy number variations involved in balsam poplar (Populus balsamifera L.) adaptive variations.

Gene copy number variations (CNVs) involved in phenotypic variations have already been shown in plants, but genomewide testing of CNVs for adaptive variation was not doable until recent technological developments. Thus, reports of the genomic architecture of adaptation involving CNVs remain scarce to date. Here, we investigated F1 progenies of an intraprovenance cross (north-north cross, 58th parallel) and an interprovenances cross (north-south cross, 58th/49th parallels) for CNVs using comparative genomic hybridization on arrays of probes targeting gene sequences in balsam poplar (Populus balsamifera L.), a widespread North American forest tree. A total of 1,721 genes were found in varying copy numbers over the set of 19,823 tested genes. These gene CNVs presented an estimated average size of 8.3 kb and were distributed over poplar's 19 chromosomes including 22 hotspot regions. Gene CNVs number was higher for the interprovenance progeny in accordance with an expected higher genetic diversity related to the composite origin of this family. Regression analyses between gene CNVs and seven adaptive trait variations resulted in 23 significant links; among these adaptive gene CNVs, 30% were located in hotspots. One-to-five gene CNVs were found related to each of the measured adaptive traits and annotated for both biotic and abiotic stress responses. These annotations can be related to the occurrence of a higher pathogenic pressure in the southern parts of balsam poplar's distribution, and higher photosynthetic assimilation rates and water-use efficiency at high latitudes. Overall, our findings suggest that gene CNVs typically having higher mutation rates than SNPs may in fact represent efficient adaptive variations against fast-evolving pathogens.

Genomic selection of juvenile height across a single-generational gap in Douglas-fir.

Here, we perform cross-generational GS analysis on coastal Douglas-fir (Pseudotsuga menziesii), reflecting trans-generational selective breeding application. A total of 1321 trees, representing 37 full-sib F1 families from 3 environments in British Columbia, Canada, were used as the training population for (1) EBVs (estimated breeding values) of juvenile height (HTJ) in the F1 generation predicting genomic EBVs of HTJ of 136 individuals in the F2 generation, (2) deregressed EBVs of F1 HTJ predicting deregressed genomic EBVs of F2 HTJ, (3) F1 mature height (HT35) predicting HTJ EBVs in F2, and (4) deregressed F1 HT35 predicting genomic deregressed HTJ EBVs in F2. Ridge regression best linear unbiased predictor (RR-BLUP), generalized ridge regression (GRR), and Bayes-B GS methods were used and compared to pedigree-based (ABLUP) predictions. GS accuracies for scenarios 1 (0.92, 0.91, and 0.91) and 3 (0.57, 0.56, and 0.58) were similar to their ABLUP counterparts (0.92 and 0.60, respectively) (using RR-BLUP, GRR, and Bayes-B). Results using deregressed values fell dramatically for both scenarios 2 and 4 which approached zero in many cases. Cross-generational GS validation of juvenile height in Douglas-fir produced predictive accuracies almost as high as that of ABLUP. Without capturing LD, GS cannot surpass the prediction of ABLUP. Here we tracked pedigree relatedness between training and validation sets. More markers or improved distribution of markers are required to capture LD in Douglas-fir. This is essential for accurate forward selection among siblings as markers that track pedigree are of little use for forward selection of individuals within controlled pollinated families.

Partial correlation analysis of transcriptomes helps detangle the growth and defense network in spruce.

In plants, there can be a trade-off between resource allocations to growth vs defense. Here, we use partial correlation analysis of gene expression to make inferences about the nature of this interaction. We studied segregating progenies of Interior spruce subject to weevil attack. In a controlled experiment, we measured pre-attack plant growth and post-attack damage with several morphological measures, and profiled transcriptomes of 188 progeny. We used partial correlations of individual transcripts (expressed sequence tags, ESTs) with pairs of growth/defense traits to identify important nodes and edges in the inferred underlying gene network, for example, those pairs of growth/defense traits with high mutual correlation with a single EST transcript. We give a method to identify such ESTs. A terpenoid ABC transporter gene showed strongest correlations (P = 0.019); its transcript represented a hub within the compact 166-member gene-gene interaction network (P = 0.004) of the negative genetic correlations between growth and subsequent pest attack. A small 21-member interaction network (P = 0.004) represented the uncovered positive correlations. Our study demonstrates partial correlation analysis identifies important gene networks underlying growth and susceptibility to the weevil in spruce. In particular, we found transcripts that strongly modify the trade-off between growth and defense, and allow identification of networks more central to the trade-off.

Genomic prediction accuracies in space and time for height and wood density of Douglas-fir using exome capture as the genotyping platform.

BACKGROUND: Genomic selection (GS) can offer unprecedented gains, in terms of cost efficiency and generation turnover, to forest tree selective breeding; especially for late expressing and low heritability traits. Here, we used: 1) exome capture as a genotyping platform for 1372 Douglas-fir trees representing 37 full-sib families growing on three sites in British Columbia, Canada and 2) height growth and wood density (EBVs), and deregressed estimated breeding values (DEBVs) as phenotypes. Representing models with (EBVs) and without (DEBVs) pedigree structure. Ridge regression best linear unbiased predictor (RR-BLUP) and generalized ridge regression (GRR) were used to assess their predictive accuracies over space (within site, cross-sites, multi-site, and multi-site to single site) and time (age-age/ trait-trait). RESULTS: The RR-BLUP and GRR models produced similar predictive accuracies across the studied traits. Within-site GS prediction accuracies with models trained on EBVs were high (RR-BLUP: 0.79-0.91 and GRR: 0.80-0.91), and were generally similar to the multi-site (RR-BLUP: 0.83-0.91, GRR: 0.83-0.91) and multi-site to single-site predictive accuracies (RR-BLUP: 0.79-0.92, GRR: 0.79-0.92). Cross-site predictions were surprisingly high, with predictive accuracies within a similar range (RR-BLUP: 0.79-0.92, GRR: 0.78-0.91). Height at 12 years was deemed the earliest acceptable age at which accurate predictions can be made concerning future height (age-age) and wood density (trait-trait). Using DEBVs reduced the accuracies of all cross-validation procedures dramatically, indicating that the models were tracking pedigree (family means), rather than marker-QTL LD. CONCLUSIONS: While GS models' prediction accuracies were high, the main driving force was the pedigree tracking rather than LD. It is likely that many more markers are needed to increase the chance of capturing the LD between causal genes and markers.

A transcriptomics investigation into pine reproductive organ development.

The development of reproductive structures in gymnosperms is still poorly studied because of a lack of genomic information and useful genetic tools. The hermaphroditic reproductive structure derived from unisexual gymnosperms is an even less studied aspect of seed plant evolution. To extend our understanding of the molecular mechanism of hermaphroditism and the determination of sexual identity of conifer reproductive structures in general, unisexual and bisexual cones from Pinus tabuliformis were profiled for gene expression using 60K microarrays. Expression patterns of genes during progression of sexual cone development were analysed using RNA-seq. The results showed that, overall, the transcriptomes of male structures in bisexual cones were more similar to those of female cones. However, the expression of several MADS-box genes in the bisexual cones was similar to that of male cones at the more juvenile developmental stage, while despite these expression shifts, male structures of bisexual cones and normal male cones were histologically indistinguishable and cone development was continuous. This study represents a starting point for in-depth analysis of the molecular regulation of cone development and also the origin of hermaphroditism in pine.

Prediction accuracies for growth and wood attributes of interior spruce in space using genotyping-by-sequencing.

BACKGROUND: Genomic selection (GS) in forestry can substantially reduce the length of breeding cycle and increase gain per unit time through early selection and greater selection intensity, particularly for traits of low heritability and late expression. Affordable next-generation sequencing technologies made it possible to genotype large numbers of trees at a reasonable cost. RESULTS: Genotyping-by-sequencing was used to genotype 1,126 Interior spruce trees representing 25 open-pollinated families planted over three sites in British Columbia, Canada. Four imputation algorithms were compared (mean value (MI), singular value decomposition (SVD), expectation maximization (EM), and a newly derived, family-based k-nearest neighbor (kNN-Fam)). Trees were phenotyped for several yield and wood attributes. Single- and multi-site GS prediction models were developed using the Ridge Regression Best Linear Unbiased Predictor (RR-BLUP) and the Generalized Ridge Regression (GRR) to test different assumption about trait architecture. Finally, using PCA, multi-trait GS prediction models were developed. The EM and kNN-Fam imputation methods were superior for 30 and 60% missing data, respectively. The RR-BLUP GS prediction model produced better accuracies than the GRR indicating that the genetic architecture for these traits is complex. GS prediction accuracies for multi-site were high and better than those of single-sites while multi-site predictability produced the lowest accuracies reflecting type-b genetic correlations and deemed unreliable. The incorporation of genomic information in quantitative genetics analyses produced more realistic heritability estimates as half-sib pedigree tended to inflate the additive genetic variance and subsequently both heritability and gain estimates. Principle component scores as representatives of multi-trait GS prediction models produced surprising results where negatively correlated traits could be concurrently selected for using PCA2 and PCA3. CONCLUSIONS: The application of GS to open-pollinated family testing, the simplest form of tree improvement evaluation methods, was proven to be effective. Prediction accuracies obtained for all traits greatly support the integration of GS in tree breeding. While the within-site GS prediction accuracies were high, the results clearly indicate that single-site GS models ability to predict other sites are unreliable supporting the utilization of multi-site approach. Principle component scores provided an opportunity for the concurrent selection of traits with different phenotypic optima.

High-resolution genetic mapping of allelic variants associated with cell wall chemistry in Populus.

BACKGROUND: QTL cloning for the discovery of genes underlying polygenic traits has historically been cumbersome in long-lived perennial plants like Populus. Linkage disequilibrium-based association mapping has been proposed as a cloning tool, and recent advances in high-throughput genotyping and whole-genome resequencing enable marker saturation to levels sufficient for association mapping with no a priori candidate gene selection. Here, multiyear and multienvironment evaluation of cell wall phenotypes was conducted in an interspecific P. trichocarpa x P. deltoides pseudo-backcross mapping pedigree and two partially overlapping populations of unrelated P. trichocarpa genotypes using pyrolysis molecular beam mass spectrometry, saccharification, and/ or traditional wet chemistry. QTL mapping was conducted using a high-density genetic map with 3,568 SNP markers. As a fine-mapping approach, chromosome-wide association mapping targeting a QTL hot-spot on linkage group XIV was performed in the two P. trichocarpa populations. Both populations were genotyped using the 34 K Populus Infinium SNP array and whole-genome resequencing of one of the populations facilitated marker-saturation of candidate intervals for gene identification. RESULTS: Five QTLs ranging in size from 0.6 to 1.8 Mb were mapped on linkage group XIV for lignin content, syringyl to guaiacyl (S/G) ratio, 5- and 6-carbon sugars using the mapping pedigree. Six candidate loci exhibiting significant associations with phenotypes were identified within QTL intervals. These associations were reproducible across multiple environments, two independent genotyping platforms, and different plant growth stages. cDNA sequencing for allelic variants of three of the six loci identified polymorphisms leading to variable length poly glutamine (PolyQ) stretch in a transcription factor annotated as an ANGUSTIFOLIA C-terminus Binding Protein (CtBP) and premature stop codons in a KANADI transcription factor as well as a protein kinase. Results from protoplast transient expression assays suggested that each of the polymorphisms conferred allelic differences in the activation of cellulose, hemicelluloses, and lignin pathway marker genes. CONCLUSION: This study illustrates the utility of complementary QTL and association mapping as tools for gene discovery with no a priori candidate gene selection. This proof of concept in a perennial organism opens up opportunities for discovery of novel genetic determinants of economically important but complex traits in plants.

Extensive functional pleiotropy of REVOLUTA substantiated through forward genetics.

In plants, genes may sustain extensive pleiotropic functional properties by individually affecting multiple, distinct traits. We discuss results from three genome-wide association studies of approximately 400 natural poplar (Populus trichocarpa) accessions phenotyped for 60 ecological/biomass, wood quality, and rust fungus resistance traits. Single-nucleotide polymorphisms (SNPs) in the poplar ortholog of the class III homeodomain-leucine zipper transcription factor gene REVOLUTA (PtREV) were significantly associated with three specific traits. Based on SNP associations with fungal resistance, leaf drop, and cellulose content, the PtREV gene contains three potential regulatory sites within noncoding regions at the gene's 3' end, where alternative splicing and messenger RNA processing actively occur. The polymorphisms in this region associated with leaf abscission and cellulose content are suggested to represent more recent variants, whereas the SNP associated with leaf rust resistance may be more ancient, consistent with REV's primary role in auxin signaling and its functional evolution in supporting fundamental processes of vascular plant development.

Genome-wide association implicates numerous genes underlying ecological trait variation in natural populations of Populus trichocarpa.

In order to uncover the genetic basis of phenotypic trait variation, we used 448 unrelated wild accessions of black cottonwood (Populus trichocarpa) from much of its range in western North America. Extensive data from large-scale trait phenotyping (with spatial and temporal replications within a common garden) and genotyping (with a 34 K Populus single nucleotide polymorphism (SNP) array) of all accessions were used for gene discovery in a genome-wide association study (GWAS). We performed GWAS with 40 biomass, ecophysiology and phenology traits and 29,355 filtered SNPs representing 3518 genes. The association analyses were carried out using a Unified Mixed Model accounting for population structure effects among accessions. We uncovered 410 significant SNPs using a Bonferroni-corrected threshold (P<1.7×10(-6)). Markers were found across 19 chromosomes, explained 1-13% of trait variation, and implicated 275 unique genes in trait associations. Phenology had the largest number of associated genes (240 genes), followed by biomass (53 genes) and ecophysiology traits (25 genes). The GWAS results propose numerous loci for further investigation. Many traits had significant associations with multiple genes, underscoring their genetic complexity. Genes were also identified with multiple trait associations within and/or across trait categories. In some cases, traits were genetically correlated while in others they were not.

A cryptic syngameon within Betula shrubs revealed: Implications for conservation in changing subarctic environments.

Arctic and subarctic ecosystems are rapidly transforming due to global warming, emphasizing the need to understand the genetic diversity and adaptive strategies of northern plant species for effective conservation. This study focuses on Betula glandulosa, a native North American tundra shrub known as dwarf birch, which demonstrates an apparent capacity to adapt to changing climate conditions. To address the taxonomic challenges associated with shrub birches and logistical complexities of sampling in the northernmost areas where species' ranges overlap, we adopted a multicriteria approach. Incorporating molecular data, ploidy level assessment and leaf morphology, we aimed to distinguish B. glandulosa individuals from other shrub birch species sampled. Our results revealed three distinct species and their hybrids within the 537 collected samples, suggesting the existence of a shrub birch syngameon, a reproductive network of interconnected species. Additionally, we identified two discrete genetic clusters within the core species, B. glandulosa, that likely correspond to two different glacial lineages. A comparison between the nuclear and chloroplast SNP data emphasizes a long history of gene exchange between different birch species and genetic clusters. Furthermore, our results highlight the significance of incorporating interfertile congeneric species in conservation strategies and underscores the need for a holistic approach to conservation in the context of climate change, considering the complex dynamics of species interactions. While further research will be needed to describe this shrub birches syngameon and its constituents, this study is a first step in recognizing its existence and disseminating awareness among ecologists and conservation practitioners. This biological phenomenon, which offers evolutionary flexibility and resilience beyond what its constituent species can achieve individually, may have significant ecological implications.

Unraveling the transcriptional features and gene expression networks of pathogenic and saprotrophic Ophiostoma species during the infection of Ulmus americana.

American elm (Ulmus americana), highly prized for its ornamental value, has suffered two successive outbreaks of Dutch elm disease (DED) caused by ascomycete fungi belonging to the genus Ophiostoma. To identify the genes linked to the pathogenicity of different species and lineages of Ophiostoma, we inoculated 2-year-old U. americana saplings with six strains representing three species of DED fungi, and one strain of the saprotroph Ophiostoma quercus. Differential expression analyses were performed following RNA sequencing of fungal transcripts recovered at 3- and 10-days post-infection. Based on a total of 8,640 Ophiostoma genes, we observed a difference in fungal gene expression depending on the strain inoculated and the time of incubation in host tissue. Some genes overexpressed in the more virulent strains of Ophiostoma encode hydrolases that possibly act synergistically. A mutant of Ophiostoma novo-ulmi in which the gene encoding the ogf1 transcription factor had been deleted did not produce transcripts for the gene encoding the hydrophobin cerato-ulmin and was less virulent. Weighted gene correlation network analyses identified several candidate pathogenicity genes distributed among 13 modules of interconnected genes.IMPORTANCEOphiostoma is a genus of cosmopolitan fungi that belongs to the family Ophiostomataceae and includes the pathogens responsible for two devastating pandemics of Dutch elm disease (DED). As the mechanisms of action of DED agents remain unclear, we carried out the first comparative transcriptomic study including representative strains of the three Ophiostoma species causing DED, along with the phylogenetically close saprotrophic species Ophiostoma quercus. Statistical analyses of the fungal transcriptomes recovered at 3 and 10 days following infection of Ulmus americana saplings highlighted several candidate genes associated with virulence and host-pathogen interactions wherein each strain showed a distinct transcriptome. The results of this research underscore the importance of investigating the transcriptional behavior of different fungal taxa to understand their pathogenicity and virulence in relation to the timeline of infection.

Is there hybridization between 2 species of the same genus in sympatry?-The genetic relationships between Anoplophora glabripennis, Anoplophora chinensis, and putative hybrids.

Anoplophora glabripennis (Asian longhorn beetle, ALB) and Anoplophora chinensis (Citrus longhorn beetle, CLB) are native forest pests in China; they have become important international quarantine pests. They are found using the same Salix aureo-pendula host tree of Cixi, Zhejiang province, China. On this host tree, we collected additional beetles that appeared to be morphologically intermediate between ALB and CLB. By using a stereoscope, we observed that there were several bumps on the base of the elytra, which was inconsistent with ALB, which typically has a smooth elytral base, but was more like CLB, which has numerous short tubercles on the elytral base. Given their sympatry and intermediate morphology, we hypothesized that these may represent ALB × CLB hybrids. We studied the genomic profiles for 46 samples (ALB, CLB, and putative hybrids) using genotyping-by-sequencing (GBS) providing a reduced representation of the entire genome. Employing principal component analyses on the 163 GBS-derived single nucleotide polymorphism data, we found putative hybrids tightly clustered with ALB, but genetically distinct from the CLB individuals. Therefore, our initial hybrid hypothesis was not supported by genomic data. Further, while mating experiments between adult ALB and CLB were successful in 4 separate years (2017, 2018, 2020, and 2021), and oviposition behavior was observed, no progeny was produced. Having employed population genomic analysis and biological hybridization experiments, we conclude that the putative hybrids represent newly discovered morphological variants within ALB. Our approach further confirmed the advantage of genome-wide information for Anoplophora species assignment in certain ambiguous classification cases.