The NLRs or NBS-LRRs (nucleotide-binding, leucine-rich-repeat) form the largest resistance gene family in plants, with lineage-specific contingents of TNL, CNL and RNL subfamilies and a central role in resilience to stress. The origin, evolution and distribution of NLR sequences has been unclear owing in part to the variable size and diversity of the RNL subfamily and a lack of data in Gymnosperms. We developed, searched and annotated transcriptomes assemblies of seven conifers and identified a resource of 3816 expressed NLR sequences. Our analyses encompassed sequences data spanning the major groups of land plants and determinations of NLR transcripts levels in response to drought in white spruce. We showed that conifers have among the most diverse and numerous RNLs in tested land plants. We report an evolutionary swap in the formation of RNLs, which emerged from the fusion of an RPW8 domain to a NB-ARC domain of CNL. We uncovered a quantitative relationship between RNLs and TNLs across all land plants investigated, with an average ratio of 1:10. The conifer RNL repertoire harbours four distinct groups, with two that differ from Angiosperms, one of which contained several upregulated sequences in response to drought while the majority of responsive NLRs are downregulated.

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

The gene encoding the monolignol biosynthetic enzyme cinnamyl alcohol dehydrogenase (CAD, E.C. 1.1.1.195) can be expressed in response to different developmental and environmental cues. Control of Cad gene expression could involve either differential regulation of more than one Cad gene or, alternatively combinatorial regulation of a single Cad gene. In loblolly pine (Pinus taeda L.), we found several electrophoretic variants (allozymes) of CAD and a high level of heterozygosity (he = 0.46). Analysis of inheritance patterns of pine CAD allozymes gave segregation ratios that were consistent with Mendelian expectations for a single functional gene. The identity of the full-length Cad cDNA sequence was confirmed by alignment with peptide sequences obtained from purified active enzyme and by extensive similarity to Cad sequences from other species. Southern blot analysis of genomic DNA using the Cad cDNA as a hybridization probe gave simple patterns, consistent with our interpretation that pine Cad is a single-copy gene. Phylogenetic analysis and evolution rate estimates showed that Cad sequences are diverging less rapidly in the gymnosperms than in the angiosperms. The Cad mRNA was present in both lignifying tissues and a non lignifying tissue (the megagametophyte) of pine. The presence of a single gene suggests that different regulatory mechanisms for a single Cad gene, rather than differential regulation of several genes, can account for its expression in response to different cues.