Genome-Wide Comparative Analysis of Flowering-Time Genes; Insights on the Gene Family Expansion and Evolutionary Perspective

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Identification and analysis of MAPK cascade gene families of Camellia oleifera and their roles in response to cold stress

BACKGROUND

Low-temperature severely limits the growth and development of Camellia oleifera (C. oleifera). The mitogen-activated protein kinase (MAPK) cascade plays a key role in the response to cold stress.

METHODS AND RESULTS

Our study aims to identify MAPK cascade genes in C. oleifera and reveal their roles in response to cold stress. In our study, we systematically identified and analyzed the MAPK cascade gene families of C. oleifera, including their physical and chemical properties, conserved motifs, and multiple sequence alignments. In addition, we characterized the interacting networks of MAPKK kinase (MAPKKK)-MAPK kinase (MAPKK)-MAPK in C. oleifera. The molecular mechanism of cold stress resistance of MAPK cascade genes in wild C. oleifera was analyzed by differential gene expression and real-time quantitative reverse transcription-PCR (qRT-PCR).

CONCLUSION

In this study, 21 MAPKs, 4 MAPKKs and 55 MAPKKKs genes were identified in the leaf transcriptome of C. oleifera. According to the phylogenetic results, MAPKs were divided into 4 groups (A, B, C and D), MAPKKs were divided into 3 groups (A, B and D), and MAPKKKs were divided into 2 groups (MEKK and Raf). Motif analysis showed that the motifs in each subfamily were conserved, and most of the motifs in the same subfamily were basically the same. The protein interaction network based on Arabidopsis thaliana (A. thaliana) homologs revealed that MAPK, MAPKK, and MAPKKK genes were widely involved in C. oleifera growth and development and in responses to biotic and abiotic stresses. Gene expression analysis revealed that the CoMAPKKK5/CoMAPKKK43/CoMAPKKK49-CoMAPKK4-CoMAPK8 module may play a key role in the cold stress resistance of wild C. oleifera at a high-elevation site in Lu Mountain (LSG). This study can facilitate the mining and utilization of genetic resources of C. oleifera with low-temperature tolerance.

Multistage temporal transcriptomic atlas unveils major contributor to reproductive phase in upland cotton

Flowering is a major developmental transition in plants, but asynchronous flowering hinders the utilization of wild cotton relatives in breeding programs. We performed comparative transcriptomic profiling of early- and late-flowering Gossypium hirsutum genotypes to elucidate genetic factors influencing reproductive timing. Shoot apices were sampled from the photoperiod-sensitive landrace G. hirsutum purpurascens (GhP) and early-maturing variety ZhongMianSuo (ZMS) at five time points following the emergence of sympodial nodes. RNA-sequencing revealed extensive transcriptional differences during floral transition. Numerous flowering-associated genes exhibited genotype-specific expression, including FLOWERING LOCUS T (FT) homologs upregulated in ZMS. FT-interacting factors like SOC1 and CO-like also showed higher expression in ZMS, implicating florigen pathways in early flowering. Additionally, circadian clock and light signalling components were misregulated between varieties, suggesting altered photoperiod responses in GhP. Weighted co-expression network analysis specifically linked a module enriched for circadian-related genes to GhP's late flowering. Through an integrated transcriptome analysis, we defined a regulatory landscape of reproductive phase change in cotton. Differentially expressed genes related to photoperiod, circadian clock, and light signalling likely contribute to delayed flowering in wild cottons. Characterization of upstream flowering regulators will enable modifying photoperiod sensitivity and expand germplasm use for cotton improvement. This study provides candidate targets for elucidating interactive mechanisms that control cotton flowering time across diverse genotypes.

Chromosome-level genome assembly provides insights into the genome evolution and functional importance of the phenylpropanoid-flavonoid pathway in Thymus mongolicus

BACKGROUND

Thymus mongolicus (family Lamiaceae) is a Thyme subshrub with strong aroma and remarkable environmental adaptability. Limited genomic information limits the use of this plant.

RESULTS

Chromosome-level 605.2 Mb genome of T. mongolicus was generated, with 96.28% anchored to 12 pseudochromosomes. The repetitive sequences were dominant, accounting for 70.98%, and 32,593 protein-coding genes were predicted. Synteny analysis revealed that Lamiaceae species generally underwent two rounds of whole genome duplication; moreover, species-specific genome duplication was identified. A recent LTR retrotransposon burst and tandem duplication might play important roles in the formation of the Thymus genome. Using comparative genomic analysis, phylogenetic tree of seven Lamiaceae species was constructed, which revealed that Thyme plants evolved recently in the family. Under the phylogenetic framework, we performed functional enrichment analysis of the genes on nodes that contained the most gene duplication events (> 50% support) and of relevant significant expanded gene families. These genes were highly associated with environmental adaptation and biosynthesis of secondary metabolites. Combined transcriptome and metabolome analyses revealed that Peroxidases, Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferases, and 4-coumarate-CoA ligases genes were the essential regulators of the phenylpropanoid-flavonoid pathway. Their catalytic products (e.g., apigenin, naringenin chalcone, and several apigenin-related compounds) might be responsible for the environmental tolerance and aromatic properties of T. mongolicus.

CONCLUSION

This study enhanced the understanding of the genomic evolution of T. mongolicus, enabling further exploration of its unique traits and applications, and contributed to the understanding of Lamiaceae genomics and evolutionary biology.

Genome-wide identification and characterization of NBLRR genes in finger millet (Eleusine coracana L.) and their expression in response to Magnaporthe grisea infection

BACKGROUND

The nucleotide binding site leucine rich repeat (NBLRR) genes significantly regulate defences against phytopathogens in plants. The genome-wide identification and analysis of NBLRR genes have been performed in several species. However, the detailed evolution, structure, expression of NBLRRs and functional response to Magnaporthe grisea are unknown in finger millet (Eleusine coracana (L.) Gaertn.).

RESULTS

The genome-wide scanning of the finger millet genome resulted in 116 NBLRR (EcNBLRRs1-116) encompassing 64 CC-NB-LRR, 47 NB-LRR and 5 CCR-NB-LRR types. The evolutionary studies among the NBLRRs of five Gramineae species, viz., purple false brome (Brachypodium distachyon (L.) P.Beauv.), finger millet (E. coracana), rice (Oryza sativa L.), sorghum (Sorghum bicolor L. (Moench)) and foxtail millet (Setaria italica (L.) P.Beauv.) showed the evolution of NBLRRs in the ancestral lineage of the target species and subsequent divergence through gene-loss events. The purifying selection (Ka/Ks < 1) shaped the expansions of NBLRRs paralogs in finger millet and orthologs among the target Gramineae species. The promoter sequence analysis showed various stress- and phytohormone-responsive cis-acting elements besides growth and development, indicating their potential role in disease defence and regulatory mechanisms. The expression analysis of 22 EcNBLRRs in the genotypes showing contrasting responses to Magnaporthe grisea infection revealed four and five EcNBLRRs in early and late infection stages, respectively. The six of these nine candidate EcNBLRRs proteins, viz., EcNBLRR21, EcNBLRR26, EcNBLRR30, EcNBLRR45, EcNBLRR55 and EcNBLRR76 showed CC, NB and LRR domains, whereas the EcNBLRR23, EcNBLRR32 and EcNBLRR83 showed NB and LRR somains.

CONCLUSION

The identification and expression analysis of EcNBLRRs showed the role of EcNBLRR genes in assigning blast resistance in finger millet. These results pave the foundation for in-depth and targeted functional analysis of EcNBLRRs through genome editing and transgenic approaches.

Identification of MYB Transcription Factors Involving in Fruit Quality Regulation of Fragaria × ananassa Duch

The cultivated strawberry (Fragaria × ananassa Duch.) is an important horticultural crop. The economic values of strawberry cultivars are decided by their fruit qualities including taste, color and aroma. The important role of MYB transcription factors in fruit quality regulation is recognized increasingly with the identification of MYB genes involved in metabolism. A total of 407 MYB genes of F. × ananassa (FaMYBs) were identified in the genome-wide scale and named according to subgenome locations. The 407 FaMYBs were clustered into 36 groups based on phylogenetic analysis. According to synteny analysis, whole genome duplication and segmental duplication contributed over 90% of the expansion of the FaMYBs family. A total of 101 FaMYB loci with 1-6 alleles were identified by the homologous gene groups on homologous chromosomes. The differentially expressed FaMYB profiles of three cultivars with different fruit quality and fruit ripe processes provided the 8 candidate loci involved in fruit quality regulation. In this experiment, 7, 5, and 4 FaMYBs were screeded as candidate genes involved in the regulation of metabolism/transportation of anthocyanins, sugars or organic acids and 4-hydroxy-2, 5-dimethyl-3(2H)-furanone, respectively. These results pointed out the key FaMYBs for further functional analysis of gene regulation of strawberry fruit quality and would be helpful in the clarification on ofe roles of MYBs in the metabolism of fruit crops.

FLOWERING LOCUS T indel variants confer vernalization-independent and photoperiod-insensitive flowering of yellow lupin (Lupinus luteus L.)

Ongoing climate change has considerably reduced the seasonal window for crop vernalization, concurrently expanding cultivation area into northern latitudes with long-day photoperiod. To address these changes, cool season legume breeders need to understand molecular control of vernalization and photoperiod. A key floral transition gene integrating signals from these pathways is the Flowering locus T (FT). Here, a recently domesticated grain legume, yellow lupin (Lupinus luteus L.), was explored for potential involvement of FT homologues in abolition of vernalization and photoperiod requirements. Two FTa (LlutFTa1a and LlutFTa1b) and FTc (LlutFTc1 and LlutFTc2) homologues were identified and sequenced for two contrasting parents of a reference recombinant inbred line (RIL) population, an early-flowering cultivar Wodjil and a late-flowering wild-type P28213. Large deletions were detected in the 5' promoter regions of three FT homologues. Quantitative trait loci were identified for flowering time and vernalization response in the RIL population and in a diverse panel of wild and domesticated accessions. A 2227 bp deletion found in the LlutFTc1 promoter was linked with early phenology and vernalization independence, whereas LlutFTa1a and LlutFTc2 indels with photoperiod responsiveness. Comparative mapping highlighted convergence of FTc1 indel evolution in two Old World lupin species, addressing both artificial selection during domestication and natural adaptation to short season environmental conditions. We concluded that rapid flowering in yellow lupin is associated with the de-repression of the LlutFTc1 homologue from the juvenile phase, putatively due to the elimination of all binding sites in the promoter region for the AGAMOUS-like 15 transcription factor.

The Antarctic Moss Pohlia nutans Genome Provides Insights Into the Evolution of Bryophytes and the Adaptation to Extreme Terrestrial Habitats

The Antarctic continent has extreme natural environment and fragile ecosystem. Mosses are one of the dominant floras in the Antarctic continent. However, their genomic features and adaptation processes to extreme environments remain poorly understood. Here, we assembled the high-quality genome sequence of the Antarctic moss (Pohlia nutans) with 698.20 Mb and 22 chromosomes. We found that the high proportion of repeat sequences and a recent whole-genome duplication (WGD) contribute to the large size genome of P. nutans when compared to other bryophytes. The genome of P. nutans harbors the signatures of massive segmental gene duplications and large expansions of gene families, likely facilitating neofunctionalization. Genomic characteristics that may support the Antarctic lifestyle of this moss comprise expanded gene families involved in phenylpropanoid biosynthesis, unsaturated fatty acid biosynthesis, and plant hormone signal transduction. Additional contributions include the significant expansion and upregulation of several genes encoding DNA photolyase, antioxidant enzymes, flavonoid biosynthesis enzymes, possibly reflecting diverse adaptive strategies. Notably, integrated multi-omic analyses elucidate flavonoid biosynthesis may function as the reactive oxygen species detoxification under UV-B radiation. Our studies provide insight into the unique features of the Antarctic moss genome and their molecular responses to extreme terrestrial environments.