Abstract:Genome-wide identification and comparative gene family analyses have been commonly performed to investigate species-specific evolution linked to various traits and molecular pathways. However, most previous studies were limited to gene screening in a single reference genome, failing to account for the gene presence/absence variations (gPAVs) in a species. Here, we propose an innovative pangenome-based approach of gene family analyses based on orthologous gene groups (OGGs). Using the basic helix-loop-helix (bHLH) transcription factor family in barley as an example, we identified 161 ∼ 176 bHLHs in 20 barley genomes, which could be classified into 201 OGGs. These 201 OGGs were further classified into 140 core, 12 soft-core, 29 shell, and 20 line-specific/cloud bHLHs, revealing a complete profile of bHLH in barley. Using a genome-scan approach, we overcome the genome annotation bias and identified on average 1.5 un-annotated core bHLHs per barley genome. We found that all core bHLHs belong to whole genome/segmental duplicates whilst dispensable bHLHs were more likely to result from small scale duplication events. Interestingly, we noticed that the dispensable bHLHs tended to enrich in specific subfamilies SF13, SF27, and SF28, implying the potential biased expansion of specific bHLHs in barley. We found that 50% of the bHLHs contain at least one intact transposon element within the 2kb upstream-to-downstream region. bHLHs with CNV have 1.48 TEs on average, significantly higher than 1.36 for core bHLH without CNV, supporting TEs' potential role in bHLH expansion. Selection pressure analyses showed that dispensable bHLHs had experienced clear relaxed selection compared to core bHLHs, consistent with their conservation patterns. We further integrate pangenome with recently available barley pantranscriptome data in 5 tissues and discovered apparent transcriptional divergence within and across bHLH subfamilies. We conclude that pangenome-based gene family analyses can better describe the genuine evolution status of bHLHs untapped before and provided novel insights into bHLH evolution in barley. We expect this study will inspire similar analyses in many other gene families and species.

Gene-set association and epistatic analyses reveal complex gene interaction networks affecting flowering time in a worldwide barley collection

Multi-Locus Genome-Wide Association Studies Reveal Novel Alleles for Flowering Time Under Vernalisation and Extended Photoperiod in a Barley MAGIC Population

Genome-wide Association Mapping Reveals Novel Genes Associated with Coleoptile Length in a Worldwide Collection of Barley

Identification of Epistasis Loci Underlying Rice Flowering Time by Controlling Population Stratification and Polygenic Effect

Genetic and Epigenetic Variations in Barley (hordeum Vulgare L.) Using SSR and MSAP Approaches

Qtl Underlying Some Agronomic Traits in Barley Detected by Snp Markers

Barley Genome, Pan Genome and Its Application to Barley Breeding

Dissecting the Genetic Basis of Grain Size and Weight in Barley (hordeum Vulgare L.) by QTL and Comparative Genetic Analyses

Genetic Dissection of the Interactions Between Semi-Dwarfing Genes Sdw1 and Ari-E and Their Effects on Agronomic Traits in a Barley MAGIC Population

Pangenome and Pantranscriptome As the New Reference for Gene Family Characterisation – a Case Study of Basic Helix-Loop-helix (bhlh) Genes in Barley

Hybridisation‐based target enrichment of phenology genes to dissect the genetic basis of yield and adaptation in barley

Genetic mapping reveals new loci and alleles for flowering time and plant height using the HvDRR population of barley

Genetic architecture underlying light and temperature mediated flowering in Arabidopsis, rice, and temperate cereals

Erratum To: Genome-wide Association Mapping of Agronomic and Morphologic Traits in Highly Structured Populations of Barley Cultivars

A multilocus association analysis method integrating phenotype and expression data reveals multiple novel associations to flowering time variation in wild‐collected Arabidopsis thaliana

Dual-trait genomic analysis in highly stratified Arabidopsis thaliana populations using genome-wide association summary statistics

Adaptive signals of flowering time pathways in wild barley from Israel over 28 generations

Integration of GWAS, linkage analysis and transcriptome analysis to reveal the genetic basis of flowering time-related traits in maize

Multi-Locus Genome-Wide Association Studies for 14 Main Agronomic Traits in Barley

Compressed variance component mixed model reveals epistasis associated with flowering in Arabidopsis