Four chromosome scale genomes and a pan-genome annotation to accelerate pecan tree breeding

Systematic expression analysis of pecan GATA gene family during graft healing reveals that CiGATA8b and CiGATA12a are involved in stress responses

GATA transcription factors are type IV zinc-finger proteins that could bind to the GATA motif within the promoters of downstream genes, thus influencing plant development and stress responses. Presently, pecan GATA (CiGATA) gene family has yet to be systematically characterized. In this study, 33 CiGATA members were identified and grouped into four classes, with genes within the same class exhibiting structural commonality. Cis-Elements in the promoters of CiGATAs were predicted to be mainly associated with light, abscisic acid, methyl jasmonate, and anaerobic induction. Four members including CiGATA8b/12a/1b/3b were highly responsive to graft healing, among which CiGATA8b and CiGATA12a were likely related to the stress responses during graft healing, as revealed by the annotation of their co-expressed genes. CiGATA8b and CiGATA12a were both located in the nucleus and acted as transcriptional suppressor and activator, respectively. Yeast one-hybrid indicated that CiGATA8b and CiGATA12a could bind to the promoters of CiNLR and CiNAC30, respectively. Functional characterization via virus-induced gene silencing and overexpression revealed that CiGATA8b could increase disease resistance, and CiGATA12a was able to alleviate oxidative damage. Our results suggest that CiGATA8b and CiGATA12a are associated with stress responses, laying a foundation for understanding the molecular mechanisms of graft healing in pecan.

Chromosome-level genome assembly of Jaguar guapote (Parachromis manguensis) by massive parallel sequencing

Parachromis managuensis is a native cichlid fish from Central America and is the most commonly traded species within its genus. This study presents the first chromosome-scale genome assembly of P. managuensis using PacBio HiFi long reads and Hi-C sequencing data. The size of the P. managuensis genome is approximately 896.66 Mb, with a scaffold N50 of 38.19 Mb. The assembled genome demonstrates high quality in terms of completeness and accuracy, with a BUSCO score of 98.85% and a quality value (QV) of 50.95. A total of 888.60 Mb (99.10%) sequences were anchored to 24 pseudochromosomes. Additionally, 21,145 protein-coding genes and 325.58 Mb (~36.31%) repetitive sequences were identified. This chromosome-level genome assembly provides a crucial reference for studying the evolution and ecological adaptability of P. managuensis.

Weighted Gene Co-Expression Network Analysis Uncovers Core Drought Responsive Genes in Pecan (Carya illinoinensis)

Drought severely affects the growth and production of pecan (Carya illinoinensis), while genes conferred drought adaptation are yet to be fully elucidated. Here, an in-depth exploration of the two different RNA-seq projects regarding drought stress (designated as P1 and P2) was performed via weighted gene co-expression network analysis. For the two projects, there existed one pair of modules (P1 turquoise module and P2 blue module) that was probably associated with drought resistance, as the paired modules both exhibited an increased expression profile with increasing water shortage stress and were annotated to be involved in oxidative stress response and the signaling pathways of abscisic acid and jasmonic acid. There were 441 and 1258 hub genes in the P1 turquoise module and P2 blue module, respectively, among which, 140 were overlapped and thus were recognized as core drought responsive genes. An additional drought stress experiment was conducted for RT-qPCR validation, and the results showed that the 20 core genes selected for detection were highly responsive to water deficit. Together, our results will be helpful for understanding the molecular mechanism of drought response and improving drought resistance in pecan.

Comparative analysis using a chromosome-scale genome assembly for Funaria hygrometrica suggests greater collinearity in mosses than in seed plants

Mosses, the largest lineage of seed-free plants, have smaller and less variable genome sizes than flowering plants. Nevertheless, whether this difference results from divergent genome dynamics is poorly known. Here, we use newly generated chromosome-scale genome assemblies for Funaria hygrometrica and comparative analysis with other moss and seed plant genomes to investigate moss genome dynamics. Although some aspects of moss genome dynamics are seed plant-like, such as the mechanism of genome size change and de novo gain/loss of genes, moss genomes retain higher synteny, and collinearity over evolutionary time than seed plant genomes. Furthermore, transposable elements and genes are more evenly distributed along chromosomes in mosses than in seed plants, a feature shared with other sequenced seed-free plant genomes. Overall, our findings support the hypothesis that large-scale genome structure and dynamics of mosses and seed plants differ. In particular, our data suggest a lower rate of gene order reshuffling along chromosomes in mosses compared to seed plants. We speculate that such lower rate of structural genomic variation and unique chromosome structure in mosses may contribute to their relatively smaller and less variable genome sizes.

Comprehensive Genomic Analysis of the CDPK Gene Family in Pecan (Carya illinoinensis) and Their Potential Roles in Salt Stress Response

Calcium-dependent protein kinases (CDPKs) are crucial for plant development and stress responses. In this study, we performed a comprehensive genomic analysis of the CDPK gene family in pecan (Carya illinoinensis) and evaluated their potential roles in salt stress responses. A total of 31 CiCDPK genes were identified and classified into four subgroups through phylogenetic analysis. Structural and promoter analyses revealed conserved motifs and regulatory elements linked to stress responses. Gene duplication analysis showed that WGD and DSD events were primary drivers of CiCDPK expansion, shaped by purifying selection. GO and KEGG annotations highlighted roles in kinase activity, calcium binding, and signal transduction, while interaction networks suggested involvement in ROS regulation and ATP-dependent phosphorylation. Tissue-specific expression patterns indicated distinct roles of CiCDPKs, with CiCDPK20 and CiCDPK31 predominantly expressed in male flowers and seeds, respectively. Transcriptome data showed that CiCDPKs exhibited distinct responses to abiotic and biotic stress, highlighting their functional specialization under various conditions. qRT-PCR analysis further confirmed the involvement of 16 CiCDPKs in salt stress adaptation, supporting their critical roles in signal transduction pathways during salinity stress. This study provides insights into CiCDPK functions, offering potential applications in breeding pecan varieties with enhanced salt tolerance.

Pan-genome analyses of 11 Fraxinus species provide insights into salt adaptation in ash trees

Ash trees (Fraxinus) exhibit rich genetic diversity and wide adaptation to various ecological environments, and several species are highly salt tolerant. Dissecting the genomic basis of salt adaptation in Fraxinus is vital for its resistance breeding. Here, we present 11 high-quality chromosome-level genome assemblies for Fraxinus species, which reveal two unequal subgenome compositions and two recent whole-genome triplication events in their evolutionary history. A Fraxinus pan-genome was constructed on the basis of structural variations and revealed that presence-absence variations (PAVs) of transmembrane transport genes have likely contributed to salt adaptation in Fraxinus. Through whole-genome resequencing of an F1 population from an interspecies cross of F. velutina 'Lula 3' (salt tolerant) with F. pennsylvanica 'Lula 5' (salt sensitive), we mapped salt-tolerance PAV-based quantitative trait loci (QTLs) and pinpointed two PAV-QTLs and candidate genes associated with Fraxinus salt tolerance. Mechanistically, FvbHLH85 enhances salt tolerance by mediating reactive oxygen species and Na+/K+ homeostasis, whereas FvSWEET5 enhances salt tolerance by mediating osmotic homeostasis. Collectively, these findings provide valuable genomic resources for Fraxinus salt-resistance breeding and the research community.

Ancient structural variants control sex-specific flowering time morphs in walnuts and hickories

Balanced mating type polymorphisms offer a distinct window into the forces shaping sexual reproduction strategies. Multiple hermaphroditic genera in Juglandaceae, including walnuts (Juglans) and hickories (Carya), show a 1:1 genetic dimorphism for male versus female flowering order (heterodichogamy). We map two distinct Mendelian inheritance mechanisms to ancient (>37 million years old) genus-wide structural DNA polymorphisms. The dominant haplotype for female-first flowering in Juglans contains tandem repeats of the 3' untranslated region of a gene putatively involved in trehalose-6-phosphate metabolism and is associated with increased cis gene expression in developing male flowers, possibly mediated by small RNAs. The Carya locus contains ~20 syntenic genes and shows molecular signatures of sex chromosome-like evolution. Inheritance mechanisms for heterodichogamy are deeply conserved, yet may occasionally turn over, as in sex determination.

Metabolomic analysis combined with machine learning algorithms enables the evaluation of postharvest pecan color stability

Nut kernel color is a crucial quality indicator affecting the consumers first impression of the product. While growing evidence suggests that plant phenolics and their derivatives are linked to nut kernel color, the compounds (biomarkers) responsible for kernel color stability during storage remain elusive. Here, pathway-based metabolomics with machine learning algorithms were employed to identify key metabolites of postharvest pecan color stability. Metabolites in phenylpropanoid, flavonoid, and anthocyanin biosynthetic pathways were analyzed in the testa of nine pecan cultivars using liquid chromatography-mass spectrometry. With color measurements, different machine learning models were compared to find relevant biomarkers of pecan color phenotypes. Results revealed potential marker compounds that included flavonoid precursors and anthocyanidins as well as anthocyanins (e.g., peonidin, delphinidin-3-O-glucoside). Our findings provide a foundation for future research in the area, and will help select genes/proteins for the breeding of pecans with stable and desirable kernel color.

Unlocking diversity from wild relatives of perennial fruit crops in the pan-genomics era

Crop wild relatives of perennial fruit crops have a wealth of untapped genetic diversity that can be utilized for cultivar development. However, barriers such as linkage drag, long juvenility, and high heterozygosity have hindered their utilization. Advancements in genome sequencing technologies and assembly methods, combined with the integration of chromosome conformation capture have made it possible to construct high-quality reference genomes. These genome assemblies can be combined into pan-genomes, capturing inter- and intraspecific variations across coding and non-coding regions. Pan-genomes of perennial fruit crops are being developed to identify the genetic basis of traits. This will help overcome breeding challenges, enabling faster and more targeted development of new cultivars with novel traits through breeding and biotechnology.

Transcriptome analysis under pecan scab infection reveals the molecular mechanisms of the defense response in pecans

Pecan scab, caused by the fungal pathogen Venturia effusa, is the most devastating disease of pecan (Carya illinoinensis) in the southeastern United States. Resistance to this pathogen is determined by a complex interaction between host genetics and disease pathotype with even field-susceptible cultivars being resistant to most scab isolates. To understand the underlying molecular mechanisms of scab resistance in pecan, we performed a transcriptome analysis of the pecan cultivar, 'Desirable', in response to inoculation with a pathogenic and a non-pathogenic scab isolate at three different time points (24, 48, and 96 hrs. post-inoculation). Differential gene expression and gene ontology enrichment analyses showed contrasting gene expression patterns and pathway enrichment in response to the contrasting isolates with varying pathogenicity. The weighted gene co-expression network analysis of differentially expressed genes detected 11 gene modules. Among them, two modules had significant enrichment of genes involved with defense responses. These genes were particularly upregulated in the resistant reaction at the early stage of fungal infection (24 h) compared to the susceptible reaction. Hub genes in these modules were predominantly related to receptor-like protein kinase activity, signal reception, signal transduction, biosynthesis and transport of plant secondary metabolites, and oxidoreductase activity. Results of this study suggest that the early response of pathogen-related signal transduction and development of cellular barriers against the invading fungus are likely defense mechanisms employed by pecan cultivars against non-virulent scab isolates. The transcriptomic data generated here provide the foundation for identifying candidate resistance genes in pecan against V. effusa and for exploring the molecular mechanisms of disease resistance.

Molecular Mechanisms of Heterosis and Its Applications in Tree Breeding: Progress and Perspectives

Heterosis, or hybrid vigor, refers to the phenomenon where hybrid progenies outperform their parents in traits such as yield and resistance. This phenomenon has been widely applied in plant breeding. Recent advances in high-throughput genomics have significantly advanced our understanding of heterosis. This review systematically summarizes the genetic, molecular, and epigenetic mechanisms underlying heterosis. Furthermore, we discuss recent advances in predictive methods for heterosis and their applications in improving growth rate, resistance to abiotic stresses, and wood yield in tree species. We also explore the role of tree genomics in unraveling the mechanisms underlying heterosis, emphasizing the potential of integrating high-resolution genomics, single-cell sequencing, and spatial transcriptomics to achieve a comprehensive understanding of heterosis from the molecular to spatial levels. Building on this, CRISPR-based gene-editing technologies can be employed to precisely edit heterotic loci, enabling the study of allele function. Additionally, molecular marker-assisted selection (MAS) can be utilized to identify heterotic loci in parental lines, facilitating the selection of optimal hybrid combinations and significantly reducing the labor and time costs of hybrid breeding. Finally, we review the utilization of heterosis in tree breeding and provide a forward-looking perspective on future research directions, highlighting the potential of integrating multi-omics approaches and emerging gene-editing tools to revolutionize tree hybrid breeding.

Comprehensive genomic analysis of CiPawPYL-PP2C-SnRK family genes in pecan (Carya illinoinensis) and functional characterization of CiPawSnRK2.1 under salt stress responses

Abscisic acid (ABA) is a pivotal regulator of plant growth, development, and responses to environmental stresses. The ABA signaling pathway involves three key components: ABA receptors known as PYLs, PP2Cs, and SnRK2s, which are conserved across higher plants. This study comprehensively investigated the PYL-PP2C-SnRK gene family in pecan, identifying 14 PYL genes, 97 PP2C genes, and 44 SnRK genes, which were categorized into subgroups through phylogenetic and sequence structure analysis. Whole-genome duplication (WGD) and dispersed duplication (DSD) were identified as major drivers of family expansion, and purifying selection was the primary evolutionary force. Tissue-specific expression analysis suggested diverse functions in different pecan tissues. qRT-PCR validation confirmed the involvement of CiPawPYLs, CiPawPP2CAs, and CiPawSnRK2s in salt stress response. Subcellular localization analysis revealed CiPawPP2C1 in the nucleus and CiPawPYL1 and CiPawSnRK2.1 in both the nucleus and the plasma membrane. In addition, VIGS indicated that CiPawSnRK2.1-silenced pecan seedling leaves display significantly reduced salt tolerance. Y2H and LCI assays verified that CiPawPP2C3 can interact with CiPawPYL5, CiPawPYL8, and CiPawSnRK2.1. This study characterizes the role of CiPawSnRK2.1 in salt stress and lays the groundwork for exploring the CiPawPYL-PP2C-SnRK module, highlighting the need to investigate the roles of other components in the pecan ABA signaling pathway.

Metabolome and Transcriptome Analyses Reveal Metabolomic Variations and Key Transcription Factors Involved in Lipid Biosynthesis During Seed Development in Carya illinoinensis

Plant oils are a large group of neutral lipids that play a vital role in the food and oleochemical industries. The pecan (Carya illinoinensis) is a promising woody oil crop known for its high-quality sources of essential fatty acids and various bioactive compounds that may aid in preventing heart diseases. However, there is still a lack of understanding regarding the accumulation of lipids and the molecular mechanism of lipid biosynthesis during seed development. This study aims to analyze the metabolite variations and molecular mechanisms of lipid biosynthesis by integrating untargeted metabolomics and transcriptomics during pecan seed development. A total of 293 differentially accumulated metabolites were identified and further categorized into 13 groups, with lipids and lipid-like molecules constituting the largest group. The oil content and fatty acid compositions of pecan embryos were assessed at three stages of seed development. Oleic acid (c18:1) and linoleic acid (c18:2n6) were found to be the most abundant unsaturated fatty acid components in pecan embryos. Additionally, a comprehensive analysis revealed 15,990 differentially expressed genes, with a focus on the key genes related to lipid metabolism. Furthermore, the study identified 1201 transcription factors from differentially expressed genes. These transcription factors were divided into 65 families, with different members in the same family exhibiting different expression patterns during seed development. The expression patterns of ten transcription factor genes during seed development were verified by qRT-PCR. Two key genes, CiABI3 and CiFUS3 were further cloned and found to be localized in the nucleus. This study used metabolome and transcriptome analysis during key periods of pecan seed development to identify the key genes associated with seed development and fatty acid biosynthesis.

Phenotypic diversity and population structure of Pecan (Carya illinoinensis) collections reveals geographic patterns

Pecan (Carya illinoinensis) is an economically important nut crop known for its genetic diversity and adaptability to various climates. Understanding the growth variability, phenological traits, and population structure of pecan populations is crucial for breeding programs and conservation. In this study, plant growth and phenological traits were evaluated over three consecutive seasons (2015-2017) for 550 genotypes from 26 provenances. Significant variations in plant height, stem diameter, and budbreak were observed among provenances, with Southern provenances exhibiting faster growth and earlier budbreak compared to Northern provenances. Population structure analysis using SNP markers revealed eight distinct subpopulations, reflecting genetic differentiation among provenances. Notably, Southern Mexico collections formed two separate clusters, while Western collections, such as 'Allen 3', 'Allen 4', and 'Riverside', were distinguished from others. 'Burkett' and 'Apache' were grouped together due to their shared maternal parentage. Principal component analysis and phylogenetic tree analysis further supported subpopulation differentiation. Genetic differentiation among the 26 populations was evident, with six clusters highly in agreement with the subpopulations identified by STRUCTURE and fastSTRUCTURE. Principal components analysis (PCA) revealed distinct groups, corresponding to subpopulations identified by genetic analysis. Discriminant analysis of PCA (DAPC) based on provenance origin further supported the genetic structure, with clear separation of provenances into distinct clusters. These findings provide valuable insights into the genetic diversity and growth patterns of pecan populations. Understanding the genetic basis of phenological traits and population structure is essential for selecting superior cultivars adapted to diverse environments. The identified subpopulations can guide breeding efforts to develop resilient rootstocks and contribute to the sustainable management of pecan genetic resources. Overall, this study enhances our understanding of pecan genetic diversity and informs conservation and breeding strategies for the long-term viability of pecan cultivation.

Detection of colinear blocks and synteny and evolutionary analyses based on utilization of MCScanX

As different taxa evolve, gene order often changes slowly enough that chromosomal 'blocks' with conserved gene orders (synteny) are discernible. The MCScanX toolkit ( https://github.com/wyp1125/MCScanX ) was published in 2012 as freely available software for the detection of such 'colinear blocks' and subsequent synteny and evolutionary analyses based on genome-wide gene location and protein sequence information. Owing to its simplicity and high efficiency for colinear block detection, MCScanX provides a powerful tool for conducting diverse synteny and evolutionary analyses. Moreover, the detection of colinear blocks has been embraced as an integral step for pangenome graph construction. Here, new application trends of MCScanX are explored, striving to better connect this increasingly used tool to other tools and accelerate insight generation from exponentially growing sequence data. We provide a detailed protocol that covers how to install MCScanX on diverse platforms, tune parameters, prepare input files from data from the National Center for Biotechnology Information, run MCScanX and its visualization and evolutionary analysis tools, and connect MCScanX with external tools, including MCScanX-transposed, Circos and SynVisio. This protocol is easily implemented by users with minimal computational background and is adaptable to new data of interest to them. The data and utility programs for this protocol can be obtained from http://bdx-consulting.com/mcscanx-protocol .

Standard Area Diagrams for Pecan Leaf Scab: Effect of Rater Experience, Location, and Leaf Size on Reliability and Accuracy of Visual Estimates

Assessments of the severity of scab (Venturia effusa), an economically significant disease of pecan, are critical for determining pecan cultivar susceptibility, disease epidemiology, and integrated disease management approaches. We developed a standard area diagram (SAD) set to aid in assessments of pecan leaflet scab. Leaflets with scab lesions were harvested and scanned using a flatbed scanner at 600 dpi, and Fiji (ImageJ) was used to determine the actual percent disease severity. The SADs had 10 leaflets ranging in severity from 0.2 to 48.9%. Forty "small" (1.34 to 7.43 cm2) and 40 "large" (7.67 to 25.9 cm2) leaflet images were randomized for rater assessments. The images were assessed twice by 36 raters, first without and then with the SADs as a guide. Data were subjected to analysis using Lin's concordance correlation coefficient (LCC, pc) to determine the accuracy of ratings and by intraclass correlation coefficient (ICC) analysis to determine interrater reliability. The effects of rater experience, rater location, and leaflet size were also determined. The SADs significantly improved the agreement between raters and the actual values (LCC, pc = 0.70 and 0.84 without and with the SADs, respectively). The reliability of estimates was improved (ICC = 0.54 and 0.82 without and with the SADs, respectively). The effect of rater location on overall concordance was significant without and with the SADs based on an analysis of variance using a generalized linear model and lsmeans separation (P < 0.05). A generalized linear mixed model analysis revealed that there was a significant interaction between rater location, experience, and the use of the SADs, with some raters having greater improvement in generalized bias and concordance. Raters had a significantly better accuracy when rating "small" leaves (LCC, pc = 0.86) compared with "large" leaves (LCC, pc = 0.82) when using the SADs, highlighting the impact of psychophysics on field evaluations of plant disease severity. The proposed SADs will serve as an improved tool for performing pecan leaflet scab assessments by the pecan research community.

Transcriptional dynamics reveals the asymmetrical events underlying graft union formation in pecan (Carya illinoinensis)

Grafting is a widely used technique for pecan propagation; however, the background molecular events underlying grafting are still poorly understood. In our study, the graft partners during pecan [Carya illinoinensis (Wangenh.) K. Koch] graft union formation were separately sampled for RNA-seq, and the transcriptional dynamics were described via weighted gene co-expression network analysis. To reveal the main events underlying grafting, the correlations between modules and grafting traits were analyzed. Functional annotation showed that during the entire graft process, signal transduction was activated in the scion, while messenger RNA splicing was induced in the rootstock. At 2 days after grafting, the main processes occurring in the scion were associated with protein synthesis and processing, while the primary processes occurring in the rootstock were energy release-related. During the period of 7-14 days after grafting, defense response was a critical process taking place in the scion; however, the main process functioning in the rootstock was photosynthesis. From 22 to 32 days after grafting, the principal processes taking place in the scion were jasmonic acid biosynthesis and defense response, whereas the highly activated processes associated with the rootstock were auxin biosynthesis and plant-type secondary cell wall biogenesis. To further prove that the graft partners responded asymmetrically to stress, hydrogen peroxide contents as well as peroxidase and β-1,3-glucanase activities were detected, and the results showed that their levels were increased in the scion not the rootstock at certain time points after grafting. Our study reveals that the scion and rootstock might respond asymmetrically to grafting in pecan, and the scion was likely associated with stress response, while the rootstock was probably involved in energy supply and xylem bridge differentiation during graft union formation.

Characteristic analysis of BZR genes family and their responses to hormone treatments and abiotic stresses in Carya illinoinensis

As the core of Brassinosteroids (BR) signaling pathway, BR-resistant (BZR) transcription factor regulates thousands of targeted genes mediating photomophogenesis, pollen sterility, cell expansion and stress response. Pecan (Carya illinoinensis) is a famous trees species of Carya, and its nut has high nutritional and economic values. However, there has no report on BZR genes family in pecan yet. Herein, totals of seven CiBZR members were identified in pecan genome, which were predicted to be hydrophilic unstable proteins and located in the nucleus. CiBZR genes had close evolutionary relationships with CcBZRs and JrBZRs in both Carya cathayensis and Juglans regia. These seven CiBZR genes were located independently on 7 chromosomes without doubling or tandem duplication. Based on the analysis of conserved motifs and gene structures, CiBZR genes were divided into three categories. More than 40 cis-acting elements were found in the 2 kb promoter regions of CiBZRs, which were mainly involved in hormone, light, and stress response, and plant growth and development. Notably, some of these CiBZR proteins were mainly located in the nucleus, had the self-activation ability and interaction relationship with BIN2 kinase, and negatively regulated the expression of CiCPD and CiDWF4. Gene expressions analysis further showed that CiBZR genes could express in many tissues and shared similar expression trends during embryo development. Moreover, most CiBZR genes responded to BR, Gibberellin (GA), Strigolactone (SL), salt, acid and osmotic stress. This study provides theoretical basis for the subsequent study on the role of CiBZR family genes in plant growth, development and stress responses.

Chromosome-level genome assembly of Platycarya strobilacea

Platycarya strobilacea belongs to the walnut family (Juglandaceae), is commonly known as species endemic to East Asia, and is an ecologically important, wind pollinated, woody deciduous tree. To facilitate this ancient tree for the ecological value and conservation of this ancient tree, we report a new high-quality genome assembly of P. strobilacea. The genome size was 677.30 Mb, with a scaffold N50 size of 45,791,698 bp, and 98.43% of the assembly was anchored to 15 chromosomes. We annotated 32,246 protein-coding genes in the genome, of which 96.30% were functionally annotated in six databases. This new high-quality assembly of P. strobilacea provide valuable resource for the phylogenetic and evolutionary analysis of the walnut family and angiosperm.

Transcriptome profile of pecan scab resistant and susceptible trees from a pecan provenance collection

Pecan scab is a devastating disease that causes damage to pecan (Carya illinoinensis (Wangenh.) K. Koch) fruit and leaves. The disease is caused by the fungus Venturia effusa (G. Winter) and the main management practice for controlling the disease is by application of fungicides at 2-to-3-week intervals throughout the growing season. Besides disease-related yield loss, application of fungicides can result in considerable cost and increases the likelihood of fungicide resistance developing in the pathogen. Resistant cultivars are available for pecan growers; although, in several cases resistance has been overcome as the pathogen adapts to infect resistant hosts. Despite the importance of host resistance in scab management, there is little information regarding the molecular basis of genetic resistance to pecan scab.The purpose of this study was to elucidate mechanisms of natural pecan scab resistance by analyzing transcripts that are differentially expressed in pecan leaf samples from scab resistant and susceptible trees. The leaf samples were collected from trees in a provenance collection orchard that represents the natural range of pecan in the US and Mexico. Trees in the orchard have been exposed to natural scab infections since planting in 1989, and scab ratings were collected over three seasons. Based on this data, ten susceptible trees and ten resistant trees were selected for analysis. RNA-seq data was collected and analyzed for diseased and non-diseased parts of susceptible trees as well as for resistant trees. A total of 313 genes were found to be differentially expressed when comparing resistant and susceptible trees without disease. For susceptible samples showing scab symptoms, 1,454 genes were identified as differentially expressed compared to non-diseased susceptible samples. Many genes involved in pathogen recognition, defense responses, and signal transduction were up-regulated in diseased samples of susceptible trees, whereas differentially expressed genes in pecan scab resistant samples were generally down-regulated compared to non-diseased susceptible samples.Our results provide the first account of candidate genes involved in resistance/susceptibility to pecan scab under natural conditions in a pecan orchard. This information can be used to aid pecan breeding programs and development of biotechnology-based approaches for generating pecan cultivars with more durable scab resistance.

Distinct ancient structural polymorphisms control heterodichogamy in walnuts and hickories

The maintenance of stable mating type polymorphisms is a classic example of balancing selection, underlying the nearly ubiquitous 50/50 sex ratio in species with separate sexes. One lesser known but intriguing example of a balanced mating polymorphism in angiosperms is heterodichogamy - polymorphism for opposing directions of dichogamy (temporal separation of male and female function in hermaphrodites) within a flowering season. This mating system is common throughout Juglandaceae, the family that includes globally important and iconic nut and timber crops - walnuts (Juglans), as well as pecan and other hickories (Carya). In both genera, heterodichogamy is controlled by a single dominant allele. We fine-map the locus in each genus, and find two ancient (>50 Mya) structural variants involving different genes that both segregate as genus-wide trans-species polymorphisms. The Juglans locus maps to a ca. 20 kb structural variant adjacent to a probable trehalose phosphate phosphatase (TPPD-1), homologs of which regulate floral development in model systems. TPPD-1 is differentially expressed between morphs in developing male flowers, with increased allele-specific expression of the dominant haplotype copy. Across species, the dominant haplotype contains a tandem array of duplicated sequence motifs, part of which is an inverted copy of the TPPD-1 3' UTR. These repeats generate various distinct small RNAs matching sequences within the 3' UTR and further downstream. In contrast to the single-gene Juglans locus, the Carya heterodichogamy locus maps to a ca. 200-450 kb cluster of tightly linked polymorphisms across 20 genes, some of which have known roles in flowering and are differentially expressed between morphs in developing flowers. The dominant haplotype in pecan, which is nearly always heterozygous and appears to rarely recombine, shows markedly reduced genetic diversity and is over twice as long as its recessive counterpart due to accumulation of various types of transposable elements. We did not detect either genetic system in other heterodichogamous genera within Juglandaceae, suggesting that additional genetic systems for heterodichogamy may yet remain undiscovered.

Pan-genome analysis of 13 Malus accessions reveals structural and sequence variations associated with fruit traits

Structural variations (SVs) and copy number variations (CNVs) contribute to trait variations in fleshy-fruited species. Here, we assemble 10 genomes of genetically diverse Malus accessions, including the ever-green cultivar 'Granny Smith' and the widely cultivated cultivar 'Red Fuji'. Combining with three previously reported genomes, we assemble the pan-genome of Malus species and identify 20,220 CNVs and 317,393 SVs. We also observe CNVs that are positively correlated with expression levels of the genes they are associated with. Furthermore, we show that the noncoding RNA generated from a 209 bp insertion in the intron of mitogen-activated protein kinase homology encoding gene, MMK2, regulates the gene expression and affects fruit coloration. Moreover, we identify overlapping SVs associated with fruit quality and biotic resistance. This pan-genome uncovers possible contributions of CNVs to gene expression and highlights the role of SVs in apple domestication and economically important traits.

GET_PANGENES: calling pangenes from plant genome alignments confirms presence-absence variation

Crop pangenomes made from individual cultivar assemblies promise easy access to conserved genes, but genome content variability and inconsistent identifiers hamper their exploration. To address this, we define pangenes, which summarize a species coding potential and link back to original annotations. The protocol get_pangenes performs whole genome alignments (WGA) to call syntenic gene models based on coordinate overlaps. A benchmark with small and large plant genomes shows that pangenes recapitulate phylogeny-based orthologies and produce complete soft-core gene sets. Moreover, WGAs support lift-over and help confirm gene presence-absence variation. Source code and documentation: https://github.com/Ensembl/plant-scripts .

Genomic Analysis of Plastid-Nuclear Interactions and Differential Evolution Rates in Coevolved Genes across Juglandaceae Species

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions.

The Progression in Developing Genomic Resources for Crop Improvement

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.

High-quality genome assembly and multi-omics analysis of pigment synthesis pathway in Auricularia cornea

Owing to its great market potential for food and health care, white Auricularia cornea, a rare edible fungus, has received increased attention in recent years. This study presents a high-quality genome assembly of A. cornea and multi-omics analysis of its pigment synthesis pathway. Continuous Long Reads libraries, combined with Hi-C-assisted assembly were used to assemble of white A. cornea. Based on this data, we analyzed the transcriptome and metabolome of purple and white strains during the mycelium, primordium, and fruiting body stages. Finally, we obtained the genome of A.cornea assembled from 13 clusters. Comparative and evolutionary analysis suggests that A.cornea is more closely related to Auricularia subglabra than to Auricularia heimuer. The divergence of white/purple A.cornea occurred approximately 40,000 years ago, and there were numerous inversions and translocations between homologous regions of the two genomes. Purple strain synthesized pigment via the shikimate pathway. The pigment in the fruiting body of A. cornea was γ-glutaminyl-3,4-dihydroxy-benzoate. During pigment synthesis, α-D-glucose-1P, citrate, 2-Oxoglutarate, and glutamate were four important intermediate metabolites, whereas polyphenol oxidase and other 20 enzyme genes were the key enzymes. This study sheds light on the genetic blueprint and evolutionary history of the white A.cornea genome, revealing the mechanism of pigment synthesis in A.cornea. It has important theoretical and practical implications for understanding the evolution of basidiomycetes, molecular breeding of white A.cornea, and deciphering the genetic regulations of edible fungi. Additionally, it provides valuable insights for the study of phenotypic traits in other edible fungi.

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae)

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes and its only congeneric species, Platycarya strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole-genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Transcriptome Analysis of Resistant and Susceptible Pecan (Carya illinoinensis) Reveals the Mechanism of Resistance to Black Spot Disease (Colletotrichum fioriniae)

Pecan, Carya illinoinensis (Wangenh.) K. Koch, is an important dried fruit and woody oil tree species grown worldwide. With continuous expansion of pecan cultivation, the frequency and scope of diseases, especially black spot disease, are increasing, damaging trees and reducing yields. In this study, the key factors in resistance to black spot disease (Colletotrichum fioriniae) were investigated between the high-resistance pecan variety "Kanza" and the low-resistance variety "Mahan". Leaf anatomy and antioxidase activities confirmed much stronger resistance to black spot disease in "Kanza" than in "Mahan". Transcriptome analysis indicated that the increased expression of genes associated with defense response, oxidation-reduction, and catalytic activity was involved in disease resistance. A connection network identified a highly expressed hub gene CiFSD2 (CIL1242S0042), which might participate in redox reactions to affect disease resistance. Overexpression of CiFSD2 in tobacco inhibited enlargement of necrotic spots and increased disease resistance. Overall, the expression of differentially expressed genes differed in pecan varieties with different levels of resistance to C. fioriniae infection. In addition, the hub genes associated with black spot resistance were identified and the functions clarified. The in-depth understanding of resistance to black spot disease provides new insights for early screening of resistant varieties and molecular-assisted breeding in pecan.

Pan-genome and transcriptome analyses provide insights into genomic variation and differential gene expression profiles related to disease resistance and fatty acid biosynthesis in eastern black walnut (Juglans nigra)

Walnut (Juglans) species are used as nut crops worldwide. Eastern black walnut (EBW, Juglans nigra), a diploid, horticultural important woody species is native to much of eastern North America. Although it is highly valued for its wood and nut, there are few resources for understanding EBW genetics. Here, we present a high-quality genome assembly of J. nigra based on Illumina, Pacbio, and Hi-C technologies. The genome size was 540.8 Mb, with a scaffold N50 size of 35.1 Mb, and 99.0% of the assembly was anchored to 16 chromosomes. Using this genome as a reference, the resequencing of 74 accessions revealed the effective population size of J. nigra declined during the glacial maximum. A single whole-genome duplication event was identified in the J. nigra genome. Large syntenic blocks among J. nigra, Juglans regia, and Juglans microcarpa predominated, but inversions of more than 600 kb were identified. By comparing the EBW genome with those of J. regia and J. microcarpa, we detected InDel sizes of 34.9 Mb in J. regia and 18.3 Mb in J. microcarpa, respectively. Transcriptomic analysis of differentially expressed genes identified five presumed NBS-LRR (NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT) genes were upregulated during the development of walnut husks and shells compared to developing embryos. We also identified candidate genes with essential roles in seed oil synthesis, including FAD (FATTY ACID DESATURASE) and OLE (OLEOSIN). Our work advances the understanding of fatty acid bioaccumulation and disease resistance in nut crops, and also provides an essential resource for conducting genomics-enabled breeding in walnut.

Proteomic Analysis of Pecan (Carya illinoinensis) Nut Development

Pecan (Carya illinoinensis) nuts are an economically valuable crop native to the United States and Mexico. A proteomic summary from two pecan cultivars at multiple time points was used to compare protein accumulation during pecan kernel development. Patterns of soluble protein accumulation were elucidated using qualitative gel-free and label-free mass-spectrometric proteomic analyses and quantitative (label-free) 2-D gel electrophoresis. Two-dimensional (2-D) gel electrophoresis distinguished a total of 1267 protein spots and shotgun proteomics identified 556 proteins. Rapid overall protein accumulation occurred in mid-September during the transition to the dough stage as the cotyledons enlarge within the kernel. Pecan allergens Car i 1 and Car i 2 were first observed to accumulate during the dough stage in late September. While overall protein accumulation increased, the presence of histones diminished during development. Twelve protein spots accumulated differentially based on 2-D gel analysis in the weeklong interval between the dough stage and the transition into a mature kernel, while eleven protein spots were differentially accumulated between the two cultivars. These results provide a foundation for more focused proteomic analyses of pecans that may be used in the future to identify proteins that are important for desirable traits, such as reduced allergen content, improved polyphenol or lipid content, increased tolerance to salinity, biotic stress, seed hardiness, and seed viability.

Comparative transcriptome analyses reveal insights into catkin bloom patterns in pecan protogynous and protandrous cultivars

In perennial plants such as pecan, once reproductive maturity is attained, there are genetic switches that are regulated and required for flower development year after year. Pecan trees are heterodichogamous with both pistillate and staminate flowers produced on the same tree. Therefore, defining genes exclusively responsible for pistillate inflorescence and staminate inflorescence (catkin) initiation is challenging at best. To understand these genetic switches and their timing, this study analyzed catkin bloom and gene expression of lateral buds collected from a protogynous (Wichita) and a protandrous (Western) pecan cultivar in summer, autumn and spring. Our data showed that pistillate flowers in the current season on the same shoot negatively impacted catkin production on the protogynous 'Wichita' cultivar. Whereas fruit production the previous year on 'Wichita' had a positive effect on catkin production on the same shoot the following year. However, fruiting the previous year nor current year pistillate flower production had no significant effect on catkin production on 'Western' (protandrous cultivar) cultivar. The RNA-Seq results present more significant differences between the fruiting and non-fruiting shoots of the 'Wichita' cultivar compared to the 'Western' cultivar, revealing the genetic signals likely responsible for catkin production. Our data presented here, indicates the genes showing expression for the initiation of both types of flowers the season before bloom.

Genome structure-based Juglandaceae phylogenies contradict alignment-based phylogenies and substitution rates vary with DNA repair genes

In lineages of allopolyploid origin, sets of homoeologous chromosomes may coexist that differ in gene content and syntenic structure. Presence or absence of genes and microsynteny along chromosomal blocks can serve to differentiate subgenomes and to infer phylogenies. We here apply genome-structural data to infer relationships in an ancient allopolyploid lineage, the walnut family (Juglandaceae), by using seven chromosome-level genomes, two of them newly assembled. Microsynteny and gene-content analyses yield identical topologies that place Platycarya with Engelhardia as did a 1980s morphological-cladistic study. DNA-alignment-based topologies here and in numerous earlier studies instead group Platycarya with Carya and Juglans, perhaps misled by past hybridization. All available data support a hybrid origin of Juglandaceae from extinct or unsampled progenitors nested within, or sister to, Myricaceae. Rhoiptelea chiliantha, sister to all other Juglandaceae, contains proportionally more DNA repair genes and appears to evolve at a rate 2.6- to 3.5-times slower than the remaining species.

Influence of Geographical Orchard Location on the Microbiome from the Progeny of a Pecan Controlled Cross

Carya&amp;nbsp;illinoinensis (Wangenh.) K.Koch production has expanded beyond the native distribution as the genetic diversity of the species, in part, has allowed the trees to grow under broad geographic and climatic ranges. Research in other plant species has demonstrated that the phytobiome enhances their ability to survive and thrive in specific environments and, conversely, is influenced by the prevailing environment and plant genetics, among other factors. We sought to analyze the microbiota of pecan seedlings from the controlled cross 'Lakota' × 'Oaxaca' that were made in Georgia and Texas, respectively, to determine if the maternal geographical origin influences the microbiome of the resulting progeny. No significant differences in bacterial communities were observed between the seeds obtained from the two different states (p = 0.081). However, seed origin did induce significant differences in leaf fungal composition (p = 0.012). Results suggest that, in addition to some environmental, epigenetics, or host genetic components, ecological processes, such as dispersal mechanisms of the host, differentially impact the pecan microbiome, which may have ramifications for the health of trees grown in different environments. Future studies on the role of the microbiome in plant health and productivity will aid in the development of sustainable agriculture for improved food security.

Identification and Expression Analysis of MPK and MKK Gene Families in Pecan (Carya illinoinensis)

Mitogen-activated protein kinases consist of three kinase modules composed of MPKs, MKKs, and MPKKKs. As members of the protein kinase (PK) superfamily, they are involved in various processes, such as developmental programs, cell division, hormonal progression, and signaling responses to biotic and abiotic stresses. In this study, a total of 18 MPKs and 10 MKKs were annotated on the pecan genome, all of which could be classified into four subgroups, respectively. The gene structures and conserved sequences of family members in the same branch were relatively similar. All MPK proteins had a conserved motif TxY, and D(L/I/V)K and VGTxxYMSPER existed in all MKK proteins. Duplication events contributed largely to the expansion of the pecan MPK and MKK gene families. Phylogenetic analysis of protein sequences from six plants indicated that species evolution occurred in pecan. Organ-specific expression profiles of MPK and MKK showed functional diversity. Ka/Ks values indicated that all genes with duplicated events underwent strong negative selection. Seven CiPawMPK and four CiPawMKK genes with high expression levels were screened by transcriptomic data from different organs, and these candidates were validated by qRT-PCR analysis of hormone-treated and stressed samples.

How do emerging long-read sequencing technologies function in transforming the plant pathology research landscape?

Long-read sequencing technologies are revolutionizing the sequencing and analysis of plant and pathogen genomes and transcriptomes, as well as contributing to emerging areas of interest in plant-pathogen interactions, disease management techniques, and the introduction of new plant varieties or cultivars. Long-read sequencing (LRS) technologies are progressively being implemented to study plants and pathogens of agricultural importance, which have substantial economic effects. The variability and complexity of the genome and transcriptome affect plant growth, development and pathogen responses. Overcoming the limitations of second-generation sequencing, LRS technology has significantly increased the length of a single contiguous read from a few hundred to millions of base pairs. Because of the longer read lengths, new analysis methods and tools have been developed for plant and pathogen genomics and transcriptomics. LRS technologies enable faster, more efficient, and high-throughput ultralong reads, allowing direct sequencing of genomes that would be impossible or difficult to investigate using short-read sequencing approaches. These benefits include genome assembly in repetitive areas, creating more comprehensive and exact genome determinations, assembling full-length transcripts, and detecting DNA and RNA alterations. Furthermore, these technologies allow for the identification of transcriptome diversity, significant structural variation analysis, and direct epigenetic mark detection in plant and pathogen genomic regions. LRS in plant pathology is found efficient for identifying and characterization of effectors in plants as well as known and unknown plant pathogens. In this review, we investigate how these technologies are transforming the landscape of determination and characterization of plant and pathogen genomes and transcriptomes efficiently and accurately. Moreover, we highlight potential areas of interest offered by LRS technologies for future study into plant-pathogen interactions, disease control strategies, and the development of new plant varieties or cultivars.

Chromosome-level genome assembly provides new insights into Japanese chestnut (Castanea crenata) genomes

Japanese chestnut (Castanea crenata Sieb. et Zucc) is an economically and ecologically important chestnut species in East Asia. Here, we presented a high-quality chromosome-level reference genome of the Japanese chestnut cultivar 'Tsukuba' by combining Nanopore long reads and Hi-C sequencing. The final assembly has a size of 718.30 Mb and consists of 12 pseudochromosomes ranging from 41.03 to 92.03 Mb, with a BUSCO complete gene percentage of 97.6%. A total of 421.37 Mb repetitive sequences and 46,744 gene models encoding 46,463 proteins were predicted in the genome. Genome evolution analysis showed that Japanese chestnut is closely related to Chinese chestnut and these species shared a common ancestor ~6.5 million years ago. This high-quality Japanese chestnut genome represents an important resource for the chestnut genomics community and will improve our understanding of chestnut biology and evolution.

Population Genetic Characteristics and Mating Type Frequency of Venturia effusa from Pecan in South America

Scab, caused by the plant-pathogenic fungus Venturia effusa, is a major disease of pecan in South America, resulting in loss of quantity and quality of nut yield. Characteristics of the populations of V. effusa in South America are unknown. We used microsatellites to describe the genetic diversity and population structure of V. effusa in South America, and determined the mating type status of the pathogen. The four hierarchically sampled orchard populations from Argentina (AR), Brazil (BRC and BRS), and Uruguay (UR) had moderate to high genotypic and gene diversity. There was evidence of population differentiation (F_st = 0.196) but the correlation between geographic distance and genetic distance was not statistically significant. Genetic differentiation was minimal between the UR, BRC, and BRS populations, and these populations were more clearly differentiated from the AR population. The MAT1-1 and MAT1-2 mating types occurred in all four orchards and their frequencies did not deviate from the 1:1 ratio expected under random mating; however, multilocus linkage equilibrium was rejected in three of the four populations. The population genetics of South American populations of V. effusa has many similarities to the population genetics of V. effusa previously described in the United States. Characterizing the populations genetics and reproductive systems of V. effusa are important to establish the evolutionary potential of the pathogen and, thus, its adaptability-and can provide a basis for informed approaches to utilizing available host resistance and determining phytosanitary needs.

GENESPACE tracks regions of interest and gene copy number variation across multiple genomes

The development of multiple chromosome-scale reference genome sequences in many taxonomic groups has yielded a high-resolution view of the patterns and processes of molecular evolution. Nonetheless, leveraging information across multiple genomes remains a significant challenge in nearly all eukaryotic systems. These challenges range from studying the evolution of chromosome structure, to finding candidate genes for quantitative trait loci, to testing hypotheses about speciation and adaptation. Here, we present GENESPACE, which addresses these challenges by integrating conserved gene order and orthology to define the expected physical position of all genes across multiple genomes. We demonstrate this utility by dissecting presence–absence, copy-number, and structural variation at three levels of biological organization: spanning 300 million years of vertebrate sex chromosome evolution, across the diversity of the Poaceae (grass) plant family, and among 26 maize cultivars. The methods to build and visualize syntenic orthology in the GENESPACE R package offer a significant addition to existing gene family and synteny programs, especially in polyploid, outbred, and other complex genomes.

The genome is the complete DNA sequence of an individual. It is a crucial foundation for many studies in medicine, agriculture, and conservation biology. Advances in genetics have made it possible to rapidly sequence, or read out, the genome of many organisms. For closely related species, scientists can then do detailed comparisons, revealing similar genes with a shared past or a common role, but comparing more distantly related organisms remains difficult.

One major challenge is that genes are often lost or duplicated over evolutionary time. One way to be more confident is to look at ‘synteny’, or how genes are organized or ordered within the genome. In some groups of species, synteny persists across millions of years of evolution. Combining sequence similarity with gene order could make comparisons between distantly related species more robust.

To do this, Lovell et al. developed GENESPACE, a software that links similarities between DNA sequences to the order of genes in a genome. This allows researchers to visualize and explore related DNA sequences and determine whether genes have been lost or duplicated. To demonstrate the value of GENESPACE, Lovell et al. explored evolution in vertebrates and flowering plants. The software was able to highlight the shared sequences between unique sex chromosomes in birds and mammals, and it was able to track the positions of genes important in the evolution of grass crops including maize, wheat, and rice.

Exploring the genetic code in this way could lead to a better understanding of the evolution of important sections of the genome. It might also allow scientists to find target genes for applications like crop improvement. Lovell et al. have designed the GENESPACE software to be easy for other scientists to use, allowing them to make graphics and perform analyses with few programming skills.

Representing sex chromosomes in genome assemblies

Sex chromosomes have evolved hundreds of independent times across eukaryotes. As genome sequencing, assembly, and scaffolding techniques rapidly improve, it is now feasible to build fully phased sex chromosome assemblies. Despite technological advances enabling phased assembly of whole chromosomes, there are currently no standards for representing sex chromosomes when publicly releasing a genome. Furthermore, most computational analysis tools are unable to efficiently investigate their unique biology relative to autosomes. We discuss a diversity of sex chromosome systems and consider the challenges of representing sex chromosome pairs in genome assemblies. By addressing these issues now as technologies for full phasing of chromosomal assemblies are maturing, we can collectively ensure that future genome analysis toolkits can be broadly applied to all eukaryotes with diverse types of sex chromosome systems. Here we provide best practice guidelines for presenting a genome assembly that contains sex chromosomes. These guidelines can also be applied to other non-recombining genomic regions, such as S-loci in plants and mating-type loci in fungi and algae.

Insight into the CBL and CIPK gene families in pecan (Carya illinoinensis): identification, evolution and expression patterns in drought response

BACKGROUND

Calcium (Ca²⁺) serves as a ubiquitous second messenger and plays a pivotal role in signal transduction. Calcineurin B-like proteins (CBLs) are plant-specific Ca²⁺ sensors that interact with CBL-interacting protein kinases (CIPKs) to transmit Ca²⁺ signals. CBL-CIPK complexes have been reported to play pivotal roles in plant development and response to drought stress; however, limited information is available about the CBL and CIPK genes in pecan, an important nut crop.

RESULTS

In the present study, a total of 9 CBL and 30 CIPK genes were identified from the pecan genome and divided into four and five clades based on phylogeny, respectively. Gene structure and distribution of conserved sequence motif analysis suggested that family members in the same clade commonly exhibited similar exon-intron structures and motif compositions. The segmental duplication events contributed largely to the expansion of pecan CBL and CIPK gene families, and Ka/Ks values revealed that all of them experienced strong negative selection. Phylogenetic analysis of CIPK proteins from 14 plant species revealed that CIPKs in the intron-poor clade originated in seed plants. Tissue-specific expression profiles of CiCBLs and CiCIPKs were analysed, presenting functional diversity. Expression profiles derived from RNA-Seq revealed distinct expression patterns of CiCBLs and CiCIPKs under drought treatment in pecan. Moreover, coexpression network analysis helped to elucidate the relationships between these genes and identify potential candidates for the regulation of drought response, which were verified by qRT-PCR analysis.

CONCLUSIONS

The characterization and analysis of CBL and CIPK genes in pecan genome could provide a basis for further functional analysis of CiCBLs and CiCIPKs in the drought stress response of pecan.

Are we there yet? Driving the road to evolutionary graph-pangenomics

With increase in the number of sequenced genomes, it is now recognized that graph-based pangenomes can provide a comprehensive platform to study diversity in a population or species, from point mutations to large chromosomal rearrangements. By incorporating concepts from graph theory, a graph-pangenome can be studied directly to identify genomic regions and genes that underlie important evolutionary processes and traits. Here, I discuss how basic concepts in graph theory can be implemented to address questions in evolutionary genomics and guide future breeding efforts. Despite its compelling versatility, a graph-pangenome assembly is still challenging especially in species with large complex genomes. As technology is rapidly improving, the graph-pangenome is expected to become a central platform in genomics studies and applications. Thus, development of tools and methods that exploit the graph structure are urged to pave the route to evolutionary graph-pangenomics.

Chromosome-level assembly of the common vetch (Vicia sativa) reference genome

Vicia sativa L. (common vetch, n = 6) is an annual, herbaceous, climbing legume, originating in the Fertile Crescent of the Middle East and now widespread in the Mediterranean basin, West, Central and Eastern Asia, North and South America. V. sativa is of economic importance as a forage legume in countries such as Australia, China, and the USA, and contributes valuable nitrogen to agricultural rotation cropping systems. To accelerate precision genome breeding and genomics-based selection of this legume, we present a chromosome-level reference genome sequence for V. sativa, constructed using a combination of long-read Oxford Nanopore sequencing, short-read Illumina sequencing, and high-throughput chromosome conformation data (CHiCAGO and Hi-C) analysis. The chromosome-level assembly of six pseudo-chromosomes has a total genome length of 1.65 Gbp, with a median contig length of 684 Kbp. BUSCO analysis of the assembly demonstrated very high completeness of 98% of the dicotyledonous orthologs. RNA-seq analysis and gene modelling enabled the annotation of 53,218 protein-coding genes. This V. sativa assembly will provide insights into vetch genome evolution and be a valuable resource for genomic breeding, genetic diversity and for understanding adaption to diverse arid environments.

Genetic Diversity Revealed by Microsatellites in Genus Carya

The genus Carya consists of 17 species divided into 3 sections: Carya or the true hickories, Apocarya or the pecan hickories, and Sinocarya or the Asian hickories. Interspecific hybrids exist and have been used in pecan cultivar development. Nuclear and plastid microsatellite or SSR markers have been useful in distinguishing species, sections, and populations. They provide evidence for hybridity between species and can confirm heredity within crosses. As more sophisticated methods of genomic evaluation are cooperatively developed for use in pecan breeding and selection, the use of these methods will be supplemented and informed by the lessons provided by microsatellite markers, as interpreted across broad germplasm collections. In this study, over 400 Carya accessions from diverse diploid and tetraploid taxa and their interspecific hybrids, maintained at the USDA National Collection of Genetic Resources for Carya (NCGR-Carya), were analyzed using 14 nuclear and 3 plastid microsatellite markers. Principal coordinate analysis showed clear taxonomic classifications at multiple taxonomic levels along with patterns of interspecific hybridity. Evidence was also found for genetic differences associated with geographic distribution. The results indicate that this group of markers is useful in examining and characterizing populations and hybrids in the genus Carya and may help delineate the composition of a core collection to help characterize the NCGR-Carya repository collection for use in its pecan breeding program. The SSR fingerprints of the inventories of the USDA NCGR-Carya repository can also be used as a reference for identifying unknown pecan trees for growers.

Green plant genomes: What we know in an era of rapidly expanding opportunities

Green plants play a fundamental role in ecosystems, human health, and agriculture. As de novo genomes are being generated for all known eukaryotic species as advocated by the Earth BioGenome Project, increasing genomic information on green land plants is essential. However, setting standards for the generation and storage of the complex set of genomes that characterize the green lineage of life is a major challenge for plant scientists. Such standards will need to accommodate the immense variation in green plant genome size, transposable element content, and structural complexity while enabling research into the molecular and evolutionary processes that have resulted in this enormous genomic variation. Here we provide an overview and assessment of the current state of knowledge of green plant genomes. To date fewer than 300 complete chromosome-scale genome assemblies representing fewer than 900 species have been generated across the estimated 450,000 to 500,000 species in the green plant clade. These genomes range in size from 12 Mb to 27.6 Gb and are biased toward agricultural crops with large branches of the green tree of life untouched by genomic-scale sequencing. Locating suitable tissue samples of most species of plants, especially those taxa from extreme environments, remains one of the biggest hurdles to increasing our genomic inventory. Furthermore, the annotation of plant genomes is at present undergoing intensive improvement. It is our hope that this fresh overview will help in the development of genomic quality standards for a cohesive and meaningful synthesis of green plant genomes as we scale up for the future.

From Genome Sequencing to CRISPR-Based Genome Editing for Climate-Resilient Forest Trees

Due to the economic and ecological importance of forest trees, modern breeding and genetic manipulation of forest trees have become increasingly prevalent. The CRISPR-based technology provides a versatile, powerful, and widely accepted tool for analyzing gene function and precise genetic modification in virtually any species but remains largely unexplored in forest species. Rapidly accumulating genetic and genomic resources for forest trees enabled the identification of numerous genes and biological processes that are associated with important traits such as wood quality, drought, or pest resistance, facilitating the selection of suitable gene editing targets. Here, we introduce and discuss the latest progress, opportunities, and challenges of genome sequencing and editing for improving forest sustainability.

Comparative plastomes of Carya species provide new insights into the plastomes evolution and maternal phylogeny of the genus

Carya, in the Juglandiodeae subfamily, is to a typical temperate-subtropical forest-tree genus for studying the phylogenetic evolution and intercontinental disjunction between eastern Asia (EA) and North America (NA). Species of the genus have high economic values worldwide for their high-quality wood and the rich healthy factors of their nuts. Although previous efforts based on multiple molecular markers or genome-wide SNPs supported the monophyly of Carya and its two EA and NA major subclades, the maternal phylogeny of Carya still need to be comprehensively evaluated. The variation of Carya plastome has never been thoroughly characterized. Here, we novelly present 19 newly generated plastomes of congeneric Carya species, including the recently rediscovered critically endangered C. poilanei. The overall assessment of plastomes revealed highly conservative in the general structures. Our results indicated that remarkable differences in several plastome features are highly consistent with the EA-NA disjunction and showed the relatively diverse matrilineal sources among EA Carya compared to NA Carya. The maternal phylogenies were conducted with different plastome regions and full-length plastome datasets from 30 plastomes, representing 26 species in six genera of Juglandoideae and Myrica rubra (as root). Six out of seven phylogenetic topologies strongly supported the previously reported relationships among genera of Juglandoideae and the two subclades of EA and NA Carya, but displayed significant incongruencies between species within the EA and NA subclades. The phylogenetic tree generated from full-length plastomes demonstrated the optimal topology and revealed significant geographical maternal relationships among Carya species, especially for EA Carya within overlapping distribution areas. The full-length plastome-based phylogenetic topology also strongly supported the taxonomic status of five controversial species as separate species of Carya. Historical and recent introgressive hybridization and plastid captures might contribute to plastome geographic patterns and inconsistencies between topologies built from different datasets, while incomplete lineage sorting could account for the discordance between maternal topology and the previous nuclear genome data-based phylogeny. Our findings highlight full-length plastomes as an ideal tool for exploring maternal relationships among the subclades of Carya, and potentially in other outcrossing perennial woody plants, for resolving plastome phylogenetic relationships.

Characterization and expression analysis of the SPL gene family during floral development and abiotic stress in pecan (Carya illinoinensis)

SQUAMOSA promoter binding protein-like (SPL) genes are a type of plant-specific transcription factors that play crucial roles in the regulation of phase transition, floral transformation, fruit development, and various stresses. Although SPLs have been characterized in several model species, no systematic analysis has been studied in pecans, an important woody oil tree species. In this study, a total of 32 SPL genes (CiSPLs) were identified in the pecan genome. After conducting phylogenetic analysis of the conserved SBP proteins from Arabidopsis, rice, and poplar, the CiSPLs were separated into eight subgroups. The CiSPL genes within the same subgroup contained very similar exon-intron structures and conserved motifs. Nine segmentally duplicated gene pairs in the pecan genome and 16 collinear gene pairs between the CiSPL and AtSPL genes were identified. Cis-element analysis showed that CiSPL genes may regulate plant meristem differentiation and seed development, participate in various biological processes, and respond to plant hormones and environmental stresses. Therefore, we focused our study on the expression profiles of CiSPL genes during flower and fruit development. Most of the CiSPL genes were predominantly expressed in buds and/or female flowers. Additionally, quantitative real time PCR (qRT-PCR) analyses confirmed that CiSPL genes showed distinct spatiotemporal expression patterns in response to drought and salt treatments. The study provides foundation for the further exploration of the function and evolution of SPL genes in pecan.

Chromosome-scale assembly reveals asymmetric paleo-subgenome evolution and targets for the acceleration of fungal resistance breeding in the nut crop, pecan

Pecan (Carya illinoinensis) is a tree nut crop of worldwide economic importance that is rich in health-promoting factors. However, pecan production and nut quality are greatly challenged by environmental stresses such as the outbreak of severe fungal diseases. Here, we report a high-quality, chromosome-scale genome assembly of the controlled-cross pecan cultivar 'Pawnee' constructed by integrating Nanopore sequencing and Hi-C technologies. Phylogenetic and evolutionary analyses reveal two whole-genome duplication (WGD) events and two paleo-subgenomes in pecan and walnut. Time estimates suggest that the recent WGD event and considerable genome rearrangements in pecan and walnut account for expansions in genome size and chromosome number after the divergence from bayberry. The two paleo-subgenomes differ in size and protein-coding gene sets. They exhibit uneven ancient gene loss, asymmetrical distribution of transposable elements (especially LTR/Copia and LTR/Gypsy), and expansions in transcription factor families (such as the extreme pecan-specific expansion in the far-red impaired response 1 family), which are likely to reflect the long evolutionary history of species in the Juglandaceae. A whole-genome scan of resequencing data from 86 pecan scab-associated core accessions identified 47 chromosome regions containing 185 putative candidate genes. Significant changes were detected in the expression of candidate genes associated with the chitin response pathway under chitin treatment in the scab-resistant and scab-susceptible cultivars 'Excell' and 'Pawnee'. These findings enable us to identify key genes that may be important susceptibility factors for fungal diseases in pecan. The high-quality sequences are valuable resources for pecan breeders and will provide a foundation for the production and quality improvement of tree nut crops.

Mitochondrial Phylogenomics of Fagales Provides Insights Into Plant Mitogenome Mosaic Evolution

Fagales are an order of woody plants and comprise more than 1,100 species, most of which produce economically important timbers, nuts, and fruits. Their nuclear and plastid genomes are well-sequenced and provided valuable resources to study their phylogeny, breeding, resistance, etc. However, little is known about the mitochondrial genomes (mitogenomes), which hinder a full understanding of their genome evolution. In this study, we assembled complete mitogenomes of 23 species, covering five of the seven families of Fagales. These mitogenomes had similar gene sets but varied 2.4 times in size. The mitochondrial genes were highly conserved, and their capacity in phylogeny was challenging. The mitogenomic structure was extremely dynamic, and synteny among species was poor. Further analyses of the Fagales mitogenomes revealed extremely mosaic characteristics, with horizontal transfer (HGT)-like sequences from almost all seed plant taxa and even mitoviruses. The largest mitogenome, Carpinus cordata, did not have large amounts of specific sequences but instead contained a high proportion of sequences homologous to other Fagales. Independent and unequal transfers of third-party DNA, including nuclear genome and other resources, may partially account for the HGT-like fragments and unbalanced size expansions observed in Fagales mitogenomes. Supporting this, a mitochondrial plasmid-like of nuclear origin was found in Carpinus. Overall, we deciphered the last genetic materials of Fagales, and our large-scale analyses provide new insights into plant mitogenome evolution and size variation.