The Arabidopsis lyrata genome sequence and the basis of rapid genome size change

Natural Variation in Sexual Traits and Gene Expression between Selfing and Outcrossing Arabidopsis lyrata Suggests Sexual Selection at Work

Flowering plants show significant diversity in sexual strategies, profoundly impacting the evolution of sexual traits and associated genes. Sexual selection is one of the primary evolutionary forces driving sexual trait variation, particularly evident during pollen-pistil interactions, where pollen grains compete for fertilization and females select mating partners. Multiple mating may intensify competition among pollen donors for siring, while in contrast, self-fertilization reduces sire-sire competition, relaxing the sexual selection pressure. Traits involved in male-male competition and female choice are well described, and molecular mechanisms underlying pollen development and pollen-pistil interactions have been extensively studied in the model species Arabidopsis thaliana. However, whether these molecular mechanisms are involved in sexual selection in nature remains unclear. To address this gap, we measured intrinsic pollen performance and its interaction with female choice and investigated the associated gene expression patterns in a selfing and an outcrossing population of Arabidopsis lyrata. We found that pollen germination and pollen tube growth were significantly higher in outcrossers than selfers, and this difference was accompanied by changes in the expression of genes involved in vesicle transport and cytoskeleton. Outcrosser mother plants showed a negative impact on pollen tube growth compared to selfer mother plants, together with a difference of expression for genes involved in auxin and stress response, suggesting a potential mechanism for female choice through molecular cross talk at the post-pollination stage. Our study provides insight into the impact of sexual selection on the evolution of sexual gene expression in plants.

Repeat competition and ecological shifts drive the evolution of the mobilome in Rhynchospora Vahl (Cyperaceae), the holocentric beaksedges

BACKGROUND AND AIMS

Genomic changes triggered by polyploidy, chromosomal rearrangements and/ or environmental stress are among factors that affect the activity of mobile elements, particularly long terminal repeat retrotransposons (LTR-RTs) and DNA transposons. Because these elements can proliferate and move throughout host genomes, altering the genetic, epigenetic and nucleotypic landscape, they have been recognized as a relevant evolutionary force. Beaksedges (Rhynchospora) stand out for their wide cosmopolitan distribution, high diversity (~400 species) and holocentric chromosomes related to high karyotypic diversity and a centromere-specific satDNA, Tyba. This makes the genus an interesting model to investigate the interactions between repetitive elements, phylogenetic relationships and ecological variables.

METHODS

Here we used comparative phylogenetic methods to investigate the forces driving the evolution of the entire set of mobile elements (mobilome) in the holocentric genus Rhynchospora. We statistically tested the impact of phylogenetic relationships, abundance of holocentromeric satDNA Tyba, diversity of repeatome composition, ecological variables and chromosome number in mobile element diversification.

KEY RESULTS

Tyba abundance was found to be inversely correlated with LTR-RT content. Decrease of LTR abundance and diversity was also related to increase in chromosome number (likely due to fission events) and colonization of dry environments in the northern hemisphere. In contrast, we found constant LTR insertions throughout time in species with lower chromosome numbers in rainier environments in South America. A multivariate model showed that different traits drive LTR abundance, especially repeat diversity and Tyba abundance. Other mobile elements, such as non-LTR RTs and DNA transposons, had insufficient abundance to be included in our models.

CONCLUSIONS

Our findings suggest that LTR evolution is strongly impacted by the holocentric characteristics of Rhynchospora chromosomes, correlating with species diversification and biome shifts, and supporting a holokinetic drive model of evolution and a competitive scenario with Tyba. Altogether, our results present evidence of multi-trait influence on LTR-RT dynamics and provide a broader understanding of transposable element evolution in a macroevolutionary context.

The Demographic History of Populations and Genomic Imprinting have Shaped the Transposon Patterns in Arabidopsis lyrata

Purifying selection is expected to prevent the accumulation of transposable elements (TEs) within their host, especially when located in and around genes and if affected by epigenetic silencing. However, positive selection may favor the spread of TEs, causing genomic imprinting under parental conflict, as genomic imprinting allows parent-specific influence over resource accumulation to the progeny. Concomitantly, the number and frequency of TE insertions in natural populations are conditioned by demographic events. In this study, we aimed to test how demography and selective forces interact to affect the accumulation of TEs around genes, depending on their epigenetic silencing, with a particular focus on imprinted genes. To this aim, we compared the frequency and distribution of TEs in Arabidopsis lyrata from Europe and North America. Generally, we found that TE insertions showed a lower frequency when they were inserted in or near genes, especially TEs targeted by epigenetic silencing, suggesting purifying selection at work. We also found that many TEs were lost or got fixed in North American populations during the colonization and the postglacial range expansion from refugia of the species in North America, as well as during the transition to selfing, suggesting a potential "TE load." Finally, we found that silenced TEs increased in frequency and even tended to reach fixation when they were linked to imprinted genes. We conclude that in A. lyrata, genomic imprinting has spread in natural populations through demographic events and positive selection acting on silenced TEs, potentially under a parental conflict scenario.

Ectopic Gene Conversion Causing Quantitative Trait Variation

Why is there so much non-neutral genetic variation segregating in natural populations? We dissect function and evolution of a near-cryptic quantitative trait locus (QTL) for defense metabolites in Arabidopsis using the CRISPR/Cas9 system and nucleotide polymorphism patterns. The QTL is explained by genetic variation in a family of 4 tightly linked indole-glucosinolate O-methyltransferase genes. Some of this variation appears to be maintained by balancing selection, some appears to be generated by non-reciprocal transfer of sequence, also known as ectopic gene conversion (EGC), between functionally diverged gene copies. Here, we elucidate how EGC, as an inevitable consequence of gene duplication, could be a general mechanism for generating genetic variation for fitness traits.

A chromosome-level genome assembly of the varied leaved jewelflower, Streptanthus diversifolius, reveals a recent whole genome duplication

The Streptanthoid complex, a clade of primarily Streptanthus and Caulanthus species in the Thelypodieae (Brassicaceae) is an emerging model system for ecological and evolutionary studies. This complex spans the full range of the California Floristic Province including desert, foothill, and mountain environments. The ability of these related species to radiate into dramatically different environments makes them a desirable study subject for exploring how plant species expand their ranges and adapt to new environments over time. Ecological and evolutionary studies for this complex have revealed fascinating variation in serpentine soil adaptation, defense compounds, germination, flowering, and life history strategies. Until now a lack of publicly available genome assemblies has hindered the ability to relate these phenotypic observations to their underlying genetic and molecular mechanisms. To help remedy this situation, we present here a chromosome-level genome assembly and annotation of Streptanthus diversifolius, a member of the Streptanthoid Complex, developed using Illumina, Hi-C, and HiFi sequencing technologies. Construction of this assembly also provides further evidence to support the previously reported recent whole genome duplication unique to the Thelypodieae. This whole genome duplication may have provided individuals in the Streptanthoid Complex the genetic arsenal to rapidly radiate throughout the California Floristic Province and to occupy commonly inhospitable environments including serpentine soils.

SNP-RFLP Markers for the Study of Arabidopsis lyrata

Arabidopsis lyrata has become a useful system for the study of comparative genomics, hybridization, polyploidization, and evolutionary transitions from outcrossing to selfing. Previous studies of its mating system have used microsatellite loci, but low allelic diversity, particularly in self-compatible populations characterized by low levels of outcrossing, reduces the utility of these markers for more detailed studies. Here, we aimed to develop population-level SNP markers for A. lyrata ssp. lyrata sampled from a self-compatible population at Rondeau Provincial Park, Ontario, Canada. We performed de novo SNP discovery and identified 6808 putative SNPs from genome-wide sequences of 22 individuals originating from a highly selfing population. Further filtering and marker validation enabled the development of 17 SNP marker loci that can be visualized using standard PCR-RFLP protocols. These markers had average minor-allele frequencies of 0.40 in the target population, and four of seven markers were variable in a small sample from nine other A. lyrata populations. These PCR-RFLP markers have the potential to be useful for the analysis of mating patterns within and beyond the inbred self-compatible populations of A. lyrata studied here and enable the continued development of A. lyrata as a model for studying evolutionary transitions from outcrossing to selfing.

Genome-wide identification and comparative evolution of 14-3-3 gene family members in five Brassicaceae species

BACKGROUND

The 14-3-3 proteins are highly conserved regulatory eukaryotic proteins, which are crucial in growth, development, and stress responses. However, systematic characterization of the 14-3-3 gene family in Brassicaceae species and their evolutionary relationships have not been comprehensively reported.

RESULTS

This study conducted genome-wide identification, structural characteristics, and comparative evolutionary analysis of 14-3-3 gene family members in Arabidopsis thaliana, A. lyrata, A. pumila, Camelina sativa, and Brassica oleracea using comparative genomics. Overall, a total of 108 14-3-3 genes, which were phylogenetically classified into ε and non-ε groups were identified in the five species, with the non-ε members exhibiting more similar exon-intron structures and conserved motif patterns. Collinearity analysis revealed that the Brassicaceae 14-3-3 gene family members underwent varying degrees of expansion following whole-genome duplication (WGD) events. Notably, the number of 14-3-3 gene family members between A. lyrata and A. thaliana remained similar despite the former having approximately 1.66-fold larger genome size. In contrast, the number of 14-3-3 gene family members in A. pumila and C. sativa increased in proportionately to their genome size, while gene members in the more distantly related species to A. thaliana, B. oleracea, showed irregular expansion patterns. Selection pressure analysis revealed that 14-3-3 homologs in all the five species underwent purifying selection, with the group ε members experiencing relatively weaker purifying selection. Cloning of ApGRF6-2 gene from A. pumila indicated that the ApGRF6-2 protein was localized in the cell membrane and cytoplasm, while ectopic overexpression of ApGRF6-2 in A. thaliana could promote early flowering by upregulating the expression of floral meristem identity genes.

CONCLUSION

This study provides a comprehensive and systematic identification of the 14-3-3 gene family members in five Brassicaceae species using updated genome sequences, and the results could form a basis for further validation of functional and molecular mechanisms of 14-3-3 genes in plant growth, development, abiotic stress responses, as well as flowering regulation.

Incorporating genetic load contributes to predicting Arabidopsis thaliana's response to climate change

Understanding how species respond to climate change can facilitate species conservation and crop breeding. Current prediction frameworks about population vulnerability focused on predicting range shifts or local adaptation but ignored genetic load, which is also crucial for adaptation. By analyzing 1115 globally distributed Arabidopsis thaliana natural accessions, we find that effective population size (Ne) is the major contributor of genetic load variation, both along genome and among populations, and can explain 74-94% genetic load variation in natural populations. Intriguingly, Ne affects genetic load by changing both effectiveness of purifying selection and GC biased gene conversion strength. In particular, by incorporating genetic load, genetic offset and species distribution models (SDM), we predict that, the populations at species' range edge are generally at higher risk. The populations at the eastern range perform poorer in all aspects, southern range have higher genetic offset and lower SDM suitability, while northern range have higher genetic load. Among the diverse natural populations, the Yangtze River basin population is the most vulnerable population under future climate change. Overall, here we deciphered the driving forces of genetic load in A. thaliana, and incorporated SDM, local adaptation and genetic load to predict the fate of populations under future climate change.

Genome-wide identification and analysis of phosphate utilization related genes (PURs) reveal their roles involved in low phosphate responses in Brassica napus L

BACKGROUND

Phosphorus (P) is an essential macronutrient for Brassica napus L. growth and development, and is mainly acquired from the soil as phosphate (Pi). However, there is no research on the system analysis of Pi utilization related genes (PURs) in B. napus yet.

RESULTS

In this study, 285 PURs were identified in B. napus genome, including 4 transcription factor (TF) gene families (83 genes) and 17 structural gene families (202 genes). Subcellular localization analysis showed that the proteins encoded by B. napus PURs were mainly located in the nucleus (~ 46.0%) and cell membrane (~ 36.5%). Chromosome localization analysis suggested that B. napus PURs were distributed on An (131) and Cn (149) subgenomes without bias. Analysis of 35 representative species confirmed that PURs were widely present in plants ranging from Chlorophyta to angiosperms with a rapid expansion trend. Collinearity analysis revealed that allopolyploidization and small-scale duplication events resulted in the large expansion of B. napus PURs. For each gene pair of B. napus PURs, the sequence identity of promoter was significantly lower than that of CDS, proving the significant difference in promoter region that might be related to the divergence of PURs expression and function. Transcription factor (TF) binding site prediction, cis-element analysis, and microRNA prediction suggested that the expressions of B. napus PURs are regulated by multiple factors including 32 TF gene families (362), 108 types of CRE (29,770) and 25 types of miRNAs (66). Spatiotemporal expression analysis demonstrated that B. napus PURs were widely expressed during the whole developmental stages, and most synteny-gene pairs (76.42%) shared conserved expression patterns. RNA-seq analyses revealed that most B. napus PURs were induced by low Pi stress, and the hub genes were generally the Pi transporter (PHT) family members. qRT-PCR analysis proved that the expression levels of four B. napus PURs were positively correlated with the root system architecture of three B. napus varieties under low Pi supply at the seedling stage.

CONCLUSION

The 285 PURs were identified from B. napus with strong LP inducible expression profile. Our findings regarding the evolution, transcriptional regulation, and expression of B. napus PURs provide valuable information for further functional research.

Green revolution gene drives adaptation of Arabidopsis to the extremely high altitude

To elucidate the process of adaptation, particularly the traits subject to natural selection and the molecular mechanisms underlying their natural variation, is one of the primary objectives of evolutionary biology. The uplifted landscape offers an excellent framework for understanding how organisms adapt to dramatic climatic gradients. To investigate the genetic basis of plant adaptation to the extremely high altitude, we first compared the genomic and phenotypic variations of two closely related Arabidopsis thaliana accessions from high altitude (Xizang, also known as "Tibet") and low altitude (Yunnan), respectively. The Xizang population represents a relict group characterized by a small effective population size. Notably, the Xizang genome has more transposable elements (TEs) and more gene loss-of-function (LoF) mutations. Differentially expressed genes were enriched in biological processes of cellular response to oxygen-containing compound, regulation of defense response, and response to light intensity. Intriguingly, the phenotypic selection analysis revealed that silique density was under natural selection. Furthermore, we genetically mapped and validated that the LoF mutation of GA20ox1, the homologous gene of green revolution in rice, resulted in a higher silique density in Xizang Arabidopsis. Given that GA20ox1 is linked to Arabidopsis adaptation to the Alps Mountains, its parallel evolution plays an important role in the adaptation to Alpine habitats. Overall, our results highlight that identifying adaptive traits and elucidating the molecular mechanisms underlying natural variation of these traits is crucial for unraveling the mystery of adaptive evolution and has significant implications for crop breeding.

Molecular mechanism of interaction between SHORT VEGETATIVE PHASE and APETALA1 in Arabidopsis thaliana

Point mutations were introduced into specific leucine (L) amino acids within the K domain of SHORT VEGETATIVE PHASE (SVP), and their effects on the SVP-AP1 interaction were assessed. Yeast two-hybrid experiments and β-galactosidase activity assays demonstrated that SVP maintained its capacity to interact with APETALA1 (AP1) despite point mutations at the 108th, 116th, 119th, and 127th leucine residues, where leucine was substituted with alanine (A). However, the mutation of the leucine residue at position 124 to alanine abolished the interaction between SVP and AP1 regardless of whether the mutation was singular or combined with others. Pull-down experiments confirmed that the leucine residue at position 124 is particularly critical for the SVP-AP1 interaction. Arabidopsis plants overexpressing 35S::AtSVP-L124A exhibited a delayed flowering phenotype compared to wild-type Col-0 Arabidopsis plants, but showed early-flowering phenotype compared to SVP overexpressing plants. SVP binds to the promoters of AP1, APETALA3 (AP3), PISTILLATA (PI), and SEPALLATA3 (SEP3), as well as to the intron of AGAMOUS (AG). Through the formation of heterodimers with AP1, SVP regulates the expression of B-class and C-class floral homeotic genes, thereby modulating floral organ development. The leucine residue at position 124 of SVP is essential for its interaction with AP1, and 35S::AtSVP-L124A transgenic plants exhibited an extended period of vegetative growth.

Genomic implications of the repeated shift to self-fertilization across a species' geographic distribution

The ability to self-fertilize often varies among closely related hermaphroditic plant species, though, variation can also exist within species. In the North American Arabidopsis lyrata, the shift from self-incompatibility (SI) to selfing established in multiple regions independently, mostly since recent postglacial range expansion. This has made the species an ideal model for the investigation of the genomic basis of the breakdown of SI and its population genetic consequences. By comparing nearby selfing and outcrossing populations across the entire species' geographic distribution, we investigated variation at the self-incompatibility (S-)locus and across the genome. Furthermore, a diallel crossing experiment on one mixed-mating population was performed to gain insight into the inheritance of mating system variation. We confirmed that the breakdown of SI had evolved in several S-locus backgrounds. The diallel suggested the involvement of biparental contributions with dominance relations. Though, the population-level genome-wide association study did not single out clear-cut candidate genes but several regions with one near the S-locus. On the implication side, selfing as compared to outcrossing populations had less than half of the genomic diversity, while the number and length of runs of homozygosity (ROHs) scaled with the degree of inbreeding. Selfing populations with a history of long expansion had the longest ROHs. The results highlight that mating system shift to selfing, its genetic underpinning and the likely negative genomic consequences for evolutionary potential can be strongly interlinked with past range dynamics.

Centrophilic retrotransposon integration via CENH3 chromatin in Arabidopsis

In organisms ranging from vertebrates to plants, major components of centromeres are rapidly evolving repeat sequences, such as tandem repeats (TRs) and transposable elements (TEs), which harbour centromere-specific histone H3 (CENH3)1,2. Complete centromere structures recently determined in human and Arabidopsis suggest frequent integration and purging of retrotransposons within the TR regions of centromeres3-5. Despite the high impact of 'centrophilic' retrotransposons on the paradox of rapid centromere evolution, the mechanisms involved in centromere targeting remain poorly understood in any organism. Here we show that both Ty3 and Ty1 long terminal repeat retrotransposons rapidly turnover within the centromeric TRs of Arabidopsis species. We demonstrate that the Ty1/Copia element Tal1 (Transposon of Arabidopsis lyrata 1) integrates de novo into regions occupied by CENH3 in Arabidopsis thaliana, and that ectopic expansion of the CENH3 region results in spread of Tal1 integration regions. The integration spectra of chimeric TEs reveal the key structural variations responsible for contrasting chromatin-targeting specificities to centromeres versus gene-rich regions, which have recurrently converted during the evolution of these TEs. Our findings show the impact of centromeric chromatin on TE-mediated rapid centromere evolution, with relevance across eukaryotic genomes.

Fixation of Expression Divergences by Natural Selection in Arabidopsis Coding Genes

Functional divergences of coding genes can be caused by divergences in their coding sequences and expression. However, whether and how expression divergences and coding sequence divergences coevolve is not clear. Gene expression divergences in differentiated cells and tissues recapitulate developmental models within a species, while gene expression divergences between analogous cells and tissues resemble traditional phylogenies in different species, suggesting that gene expression divergences are molecular traits that can be used for evolutionary studies. Using transcriptomes and evolutionary proxies to study gene expression divergences among differentiated cells and tissues in Arabidopsis, expression divergences of coding genes are shown to be strongly anti-correlated with phylostrata (gene ages), indicators of selective constraint Ka/Ks (nonsynonymous replacement rate/synonymous substitution rate) and indicators of positive selection (frequency of loci with Ka/Ks > 1), but only weakly or not correlated with indicators of neutral selection (Ks). Our results thus suggest that expression divergences largely coevolve with coding sequence divergences, suggesting that expression divergences of coding genes are selectively fixed by natural selection but not neutral selection, which provides a molecular framework for trait diversification, functional adaptation and speciation. Our findings therefore support that positive selection rather than negative or neutral selection is a major driver for the origin and evolution of Arabidopsis genes, supporting the Darwinian theory at molecular levels.

Heat shock responsive genes in Brassicaceae: genome-wide identification, phylogeny, and evolutionary associations within and between genera

Heat stress affects the growth and development of Brassicaceae crops. Plant breeders aim to mitigate the effects of heat stress by selecting for heat stress tolerance, but the genes responsible for heat stress in Brassicaceae remain largely unknown. During heat stress, heat shock proteins (HSPs) function as molecular chaperones to aid in protein folding, and heat shock transcription factors (HSFs) serve as transcriptional regulators of HSP expression. We identified 5002 heat shock related genes, including HSPs and HSFs, across 32 genomes in Brassicaceae. Among these, 3347 genes were duplicated, with segmented duplication primarily contributing to their expansion. We identified 466 physical gene clusters, including 240 homogenous clusters and 226 heterogeneous clusters, shedding light on the organization of heat shock related genes. Notably, 37 genes were co-located with published thermotolerance quantitative trait loci, which supports their functional role in conferring heat stress tolerance. This study provides a comprehensive resource for the identification of functional Brassicaceae heat shock related genes, elucidates their clustering and duplication patterns and establishes the genomic foundation for future heat tolerance research. We hypothesise that genetic variants in HSP and HSF genes in certain species have potential for improving heat stress tolerance in Brassicaceae crops.

Conserved carotenoid pigmentation in reproductive organs of Charophyceae

Sexual reproduction in Charophyceae abounds in complex traits. Their gametangia develop as intricate structures, with oogonia spirally surrounded by envelope cells and richly pigmented antheridia. The red-probably protectant-pigmentation of antheridia is conserved across Charophyceae. Chara tomentosa is, however, unique in exhibiting this pigmentation and also in vegetative tissue. Here, we investigated the two sympatric species, C. tomentosa and Chara baltica, and compared their molecular chassis for pigmentation. Using reversed phase C30 high performance liquid chromatography (RP-C30-HPLC), we uncover that the major pigments are β-carotene, δ-carotene and γ-carotene; using headspace solid-phase microextraction coupled to gas chromatography equipped with a mass spectrometer (HS-SPME-GC-MS), we pinpoint that the unusually large carotenoid pool in C. tomentosa gives rise to diverse volatile apocarotenoids, including abundant 6-methyl-5-hepten-2-one. Based on transcriptome analyses, we uncover signatures of the unique biology of Charophycaee and genes for pigment production, including monocyclized carotenoids. The rich carotenoid pool probably serves as a substrate for diverse carotenoid-derived metabolites, signified not only by (i) the volatile apocarotenoids we detected but (ii) the high expression of a gene coding for a cytochrome P450 enzyme related to land plant proteins involved in the biosynthesis of carotenoid-derived hormones. Overall, our data shed light on a key protection strategy of sexual reproduction in the widespread group of macroalgae. The genetic underpinnings of this are shared across hundreds of millions of years of plant and algal evolution. This article is part of the theme issue 'The evolution of plant metabolism'.

Exploring the complexity of genome size reduction in angiosperms

The genome sizes of angiosperms decreased significantly more than the genome sizes of their ancestors (pteridophytes and gymnosperms). Decreases in genome size involve a highly complex process, with remnants of the genome size reduction scattered across the genome and not directly linked to specific genomic structures. This is because the associated mechanisms operate on a much smaller scale than the mechanisms mediating increases in genome size. This review thoroughly summarizes the available literature regarding the molecular mechanisms underlying genome size reductions and introduces Utricularia gibba and Arabidopsis thaliana as model species for the examination of the effects of these molecular mechanisms. Additionally, we propose that phosphorus deficiency and drought stress are the major external factors contributing to decreases in genome size. Considering these factors affect almost all land plants, angiosperms likely gained the mechanisms for genome size reductions. These environmental factors may affect the retention rates of deletions, while also influencing the mutation rates of deletions via the functional diversification of the proteins facilitating double-strand break repair. The biased retention and mutation rates of deletions may have synergistic effects that enhance deletions in intergenic regions, introns, transposable elements, duplicates, and repeats, leading to a rapid decrease in genome size. We suggest that these selection pressures and associated molecular mechanisms may drive key changes in angiosperms during recurrent cycles of genome size decreases and increases.

Butterfly eggs prime anti-herbivore defense in an annual but not perennial Arabidopsis species

MAIN CONCLUSION

Unlike Arabidopsis thaliana, defenses of Arabidopsis lyrata against Pieris brassicae larval feeding are not primable by P. brassicae eggs. Thus, egg primability of plant anti-herbivore defenses is not phylogenetically conserved in the genus Arabidopsis. While plant anti-herbivore defenses of the annual species Arabidopsis thaliana were shown to be primable by Pieris brassicae eggs, the primability of the phylogenetically closely related perennial Arabidopsis lyrata has not yet been investigated. Previous studies revealed that closely related wild Brassicaceae plant species, the annual Brassica nigra and the perennial B. oleracea, exhibit an egg-primable defense trait, even though they have different life spans. Here, we tested whether P. brassicae eggs prime anti-herbivore defenses of the perennial A. lyrata. We exposed A. lyrata to P. brassicae eggs and larval feeding and assessed their primability by (i) determining the biomass of P. brassicae larvae after feeding on plants with and without prior P. brassicae egg deposition and (ii) investigating the plant transcriptomic response after egg deposition and/or larval feeding. For comparison, these studies were also conducted with A. thaliana. Consistent with previous findings, A. thaliana's response to prior P. brassicae egg deposition negatively affected conspecific larvae feeding upon A. thaliana. However, this was not observed in A. lyrata. Arabidopsis thaliana responded to P. brassicae eggs with strong transcriptional reprogramming, whereas A. lyrata responses to eggs were negligible. In response to larval feeding, A. lyrata exhibited a greater transcriptome change compared to A. thaliana. Among the strongly feeding-induced A. lyrata genes were those that are egg-primed in feeding-induced A. thaliana, i.e., CAX3, PR1, PR5, and PDF1.4. These results suggest that A. lyrata has evolved a robust feeding response that is independent from prior egg exposure.

Differences in activity and stability drive transposable element variation in tropical and temperate maize

Much of the profound interspecific variation in genome content has been attributed to transposable elements (TEs). To explore the extent of TE variation within species, we developed an optimized open-source algorithm, panEDTA, to de novo annotate TEs in a pangenome context. We then generated a unified TE annotation for a maize pangenome derived from 26 reference-quality genomes, which reveals an excess of 35.1 Mb of TE sequences per genome in tropical maize relative to temperate maize. A small number (n = 216) of TE families, mainly LTR retrotransposons, drive these differences. Evidence from the methylome, transcriptome, LTR age distribution, and LTR insertional polymorphisms reveals that 64.7% of the variability is contributed by LTR families that are young, less methylated, and more expressed in tropical maize, whereas 18.5% is driven by LTR families with removal or loss in temperate maize. Additionally, we find enrichment for Young LTR families adjacent to nucleotide-binding and leucine-rich repeat (NLR) clusters of varying copy number across lines, suggesting TE activity may be associated with disease resistance in maize.

The genetic architecture of the load linked to dominant and recessive self-incompatibility alleles in Arabidopsis halleri and Arabidopsis lyrata

The long-term balancing selection acting on mating types or sex-determining genes is expected to lead to the accumulation of deleterious mutations in the tightly linked chromosomal segments that are locally 'sheltered' from purifying selection. However, the factors determining the extent of this accumulation are poorly understood. Here, we took advantage of variations in the intensity of balancing selection along a dominance hierarchy formed by alleles at the sporophytic self-incompatibility system of the Brassicaceae to compare the pace at which linked deleterious mutations accumulate among them. We first experimentally measured the phenotypic manifestation of the linked load at three different levels of the dominance hierarchy. We then sequenced and phased polymorphisms in the chromosomal regions linked to 126 distinct copies of S-alleles in two populations of Arabidopsis halleri and three populations of Arabidopsis lyrata. We find that linkage to the S-locus locally distorts phylogenies over about 10-30 kb along the chromosome. The more intense balancing selection on dominant S-alleles results in greater fixation of linked deleterious mutations, while recessive S-alleles accumulate more linked deleterious mutations that are segregating. Hence, the structure rather than the overall magnitude of the linked genetic load differs between dominant and recessive S-alleles. Our results have consequences for the long-term evolution of new S-alleles, the evolution of dominance modifiers between them, and raise the question of why the non-recombining regions of some sex and mating type chromosomes expand over evolutionary times while others, such as the S-locus of the Brassicaceae, remain restricted to small chromosomal regions.

Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons

BACKGROUND

Centromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution.

RESULTS

Based on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche.

CONCLUSION

Our data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

Genome-wide identification and characterization of SNW/SKIP domain-containing proteins in plants

Sessile organisms, such as plants, developed various ways to sense and respond to external and internal stimuli to maximize their fitness through evolutionary time. Transcripts and protein regulation are, among many, the main mechanisms that plants use to respond to environmental changes. SKIP protein is one such, presenting an SNKW interacting domain, which is highly conserved among eukaryotes, where SKI interacting protein acts in regulating key processes. In the present work, many bioinformatics tools, such as phylogenetic relationships, gene structure, physical-chemical properties, conserved motifs, prediction of regulatory cis-elements, chromosomal localization, and protein-protein interaction network, were used to better understand the genome-wide SNW/SKIP domain-containing proteins. In total, 28 proteins containing the SNW/SKIP domain were identified in different plant species, including plants of agronomic interest. Two main protein clusters were formed in phylogenetic analysis, and gene structure analysis revealed that, in general, the coding region had no introns. Also, expression of these genes is possibly induced by abiotic stress stimuli. Primary structure analysis of the proteins revealed the existence of an evolutionarily conserved functional unit. But physicochemical properties show that proteins containing the SNW/SKIP domain are commonly unstable under in vivo conditions. In addition, the protein network, demonstrated that SKIP homologues could act by modulating plant fitness through gene expression regulation at the transcriptional and post-transcriptional levels. This could be corroborated by the expression number of gene copies of SKIP proteins in many species, highlighting it's crucial role in plant development and tolerance through the course of evolution.

Two imprinted genes primed by DEMETER in the central cell and activated by WRKY10 in the endosperm

The early development of the endosperm is crucial for balancing the allocation of maternal nutrients to offspring. This process is believed to be evolutionarily associated with genomic imprinting, resulting in parentally biased allelic gene expression. Beyond FertilizationIndependentSeed (FIS) genes, the number of imprinted genes involved in early endosperm development and seed size determination remains limited. This study introduces early endosperm-expressed HAIKU (IKU) downstream Candidate F-box 1 (ICF1) and ICF2 as maternally expressed imprinted genes (MEGs) in Arabidopsis thaliana. Although these genes are also demethylated by DEMETER (DME) in the central cell, their activation differs from the direct DME-mediated activation seen in classical MEGs such as the FIS genes. Instead, ICF maternal alleles carry pre-established hypomethylation in their promoters, priming them for activation by the WRKY10 transcription factor in the endosperm. On the contrary, paternal alleles are predominantly suppressed by CG methylation. Furthermore, we find that ICF genes partially contribute to the small seed size observed in iku mutants. Our discovery reveals a two-step regulatory mechanism that highlights the important role of conventional transcription factors in the activation of imprinted genes, which was previously not fully recognized. Therefore, the mechanism provides a new dimension to understand the transcriptional regulation of imprinting in plant reproduction and development.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

Plastome evolution of Engelhardia facilitates phylogeny of Juglandaceae

BACKGROUND

Engelhardia (Juglandaceae) is a genus of significant ecological and economic importance, prevalent in the tropics and subtropics of East Asia. Although previous efforts based on multiple molecular markers providing profound insights into species delimitation and phylogeography of Engelhardia, the maternal genome evolution and phylogeny of Engelhardia in Juglandaceae still need to be comprehensively evaluated. In this study, we sequenced plastomes from 14 samples of eight Engelhardia species and the outgroup Rhoiptelea chiliantha, and incorporated published data from 36 Juglandaceae and six outgroup species to test phylogenetic resolution. Moreover, comparative analyses of the plastomes were conducted to investigate the plastomes evolution of Engelhardia and the whole Juglandaceae family.

RESULTS

The 13 Engelhardia plastomes were highly similar in genome size, gene content, and order. They exhibited a typical quadripartite structure, with lengths from 161,069 bp to 162,336 bp. Three mutation hotspot regions (TrnK-rps16, ndhF-rpl32, and ycf1) could be used as effective molecular markers for further phylogenetic analyses and species identification. Insertion and deletion (InDels) may be an important driving factor for the evolution of plastomes in Juglandoideae and Engelhardioideae. A total of ten codons were identified as the optimal codons in Juglandaceae. The mutation pressure mostly contributed to shaping codon usage. Seventy-eight protein-coding genes in Juglandaceae experienced relaxed purifying selection, only rpl22 and psaI genes showed positive selection (Ka/Ks > 1). Phylogenetic results fully supported Engelhardia as a monophyletic group including two sects and the division of Juglandaceae into three subfamilies. The Engelhardia originated in the Late Cretaceous and diversified in the Late Eocene, and Juglandaceae originated in the Early Cretaceous and differentiated in Middle Cretaceous. The phylogeny and divergence times didn't support rapid radiation occurred in the evolution history of Engelhardia.

CONCLUSION

Our study fully supported the taxonomic treatment of at the section for Engelhardia species and three subfamilies for Juglandaceae and confirmed the power of phylogenetic resolution using plastome sequences. Moreover, our results also laid the foundation for further studying the course, tempo and mode of plastome evolution of Engelhardia and the whole Juglandaceae family.

Evolutionary assembly of the plant terrestrialization toolkit from protein domains

Land plants (embryophytes) came about in a momentous evolutionary singularity: plant terrestrialization. This event marks not only the conquest of land by plants but also the massive radiation of embryophytes into a diverse array of novel forms and functions. The unique suite of traits present in the earliest land plants is thought to have been ushered in by a burst in genomic novelty. Here, we asked the question of how these bursts were possible. For this, we explored: (i) the initial emergence and (ii) the reshuffling of domains to give rise to hallmark environmental response genes of land plants. We pinpoint that a quarter of the embryophytic genes for stress physiology are specific to the lineage, yet a significant portion of this novelty arises not de novo but from reshuffling and recombining of pre-existing domains. Our data suggest that novel combinations of old genomic substrate shaped the plant terrestrialization toolkit, including hallmark processes in signalling, biotic interactions and specialized metabolism.

Exploring the emerging role of long non-coding RNAs (lncRNAs) in plant biology: Functions, mechanisms of action, and future directions

Long non-coding RNAs (lncRNAs) are a class of RNA transcripts that surpass 200 nucleotides in length and lack discernible coding potential. LncRNAs that have been functionally characterized have pivotal functions in several plant processes, including the regulation of flowering, and development of lateral roots. It also plays a crucial role in the plant's response to abiotic stressors and exhibits vital activities in environmental adaptation. The progress in NGS (next-generation sequencing) and functional genomics technology has facilitated the discovery of lncRNA in plant species. This review is a brief explanation of lncRNA genomics, its molecular role, and the mechanism of action in plants. The review also addresses the challenges encountered in this field and highlights promising molecular and computational methodologies that can aid in the comparative and functional analysis of lncRNAs.

Genome-wide identification and evolution of the tubulin gene family in Camelina sativa

BACKGROUND

Tubulins play crucial roles in numerous fundamental processes of plant development. In flowering plants, tubulins are grouped into α-, β- and γ-subfamilies, while α- and β-tubulins possess a large isotype diversity and gene number variations among different species. This circumstance leads to insufficient recognition of orthologous isotypes and significantly complicates extrapolation of obtained experimental results, and brings difficulties for the identification of particular tubulin isotype function. The aim of this research is to identify and characterize tubulins of an emerging biofuel crop Camelina sativa.

RESULTS

We report comprehensive identification and characterization of tubulin gene family in C. sativa, including analyses of exon-intron organization, duplicated genes comparison, proper isotype designation, phylogenetic analysis, and expression patterns in different tissues. 17 α-, 34 β- and 6 γ-tubulin genes were identified and assigned to a particular isotype. Recognition of orthologous tubulin isotypes was cross-referred, involving data of phylogeny, synteny analyses and genes allocation on reconstructed genomic blocks of Ancestral Crucifer Karyotype. An investigation of expression patterns of tubulin homeologs revealed the predominant role of N6 (A) and N7 (B) subgenomes in tubulin expression at various developmental stages, contrarily to general the dominance of transcripts of H7 (C) subgenome.

CONCLUSIONS

For the first time a complete set of tubulin gene family members was identified and characterized for allohexaploid C. sativa species. The study demonstrates the comprehensive approach of precise inferring gene orthology. The applied technique allowed not only identifying C. sativa tubulin orthologs in model Arabidopsis species and tracking tubulin gene evolution, but also uncovered that A. thaliana is missing orthologs for several particular isotypes of α- and β-tubulins.

The Carthamus tinctorius L. genome sequence provides insights into synthesis of unsaturated fatty acids

Domesticated safflower (Carthamus tinctorius L.) is a widely cultivated edible oil crop. However, despite its economic importance, the genetic basis underlying key traits such as oil content, resistance to biotic and abiotic stresses, and flowering time remains poorly understood. Here, we present the genome assembly for C. tinctorius variety Jihong01, which was obtained by integrating Oxford Nanopore Technologies (ONT) and BGI-SEQ500 sequencing results. The assembled genome was 1,061.1 Mb, and consisted of 32,379 protein-coding genes, 97.71% of which were functionally annotated. Safflower had a recent whole genome duplication (WGD) event in evolution history and diverged from sunflower approximately 37.3 million years ago. Through comparative genomic analysis at five seed development stages, we unveiled the pivotal roles of fatty acid desaturase 2 (FAD2) and fatty acid desaturase 6 (FAD6) in linoleic acid (LA) biosynthesis. Similarly, the differential gene expression analysis further reinforced the significance of these genes in regulating LA accumulation. Moreover, our investigation of seed fatty acid composition at different seed developmental stages unveiled the crucial roles of FAD2 and FAD6 in LA biosynthesis. These findings offer important insights into enhancing breeding programs for the improvement of quality traits and provide reference resource for further research on the natural properties of safflower.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

Lineage-specific gene duplication and expansion of DUF1216 gene family in Brassicaceae

Proteins containing domain of unknown function (DUF) are prevalent in eukaryotic genome. The DUF1216 proteins possess a conserved DUF1216 domain resembling to the mediator protein of Arabidopsis RNA polymerase II transcriptional subunit-like protein. The DUF1216 family are specifically existed in Brassicaceae, however, no comprehensive evolutionary analysis of DUF1216 genes have been performed. We performed a first comprehensive genome-wide analysis of DUF1216 proteins in Brassicaceae. Totally 284 DUF1216 genes were identified in 27 Brassicaceae species and classified into four subfamilies on the basis of phylogenetic analysis. The analysis of gene structure and conserved motifs revealed that DUF1216 genes within the same subfamily exhibited similar intron/exon patterns and motif composition. The majority members of DUF1216 genes contain a signal peptide in the N-terminal, and the ninth position of the signal peptide in most DUF1216 is cysteine. Synteny analysis revealed that segmental duplication is a major mechanism for expanding of DUF1216 genes in Brassica oleracea, Brassica juncea, Brassica napus, Lepidium meyneii, and Brassica carinata, while in Arabidopsis thaliana and Capsella rubella, tandem duplication plays a major role in the expansion of the DUF1216 gene family. The analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios for DUF1216 paralogous indicated that most of gene pairs underwent purifying selection. DUF1216 genes displayed a specifically high expression in reproductive tissues in most Brassicaceae species, while its expression in Brassica juncea was specifically high in root. Our studies offered new insights into the phylogenetic relationships, gene structures and expressional patterns of DUF1216 members in Brassicaceae, which provides a foundation for future functional analysis.

Differential expression and evolutionary diversification of RNA helicases in Boechera sexual and apomictic reproduction

In higher plants, sexual reproduction is characterized by meiosis of the first cells of the germlines, and double fertilization of the egg and central cell after gametogenesis. In contrast, in apomicts of the genus Boechera, meiosis is omitted or altered and only the central cell requires fertilization, while the embryo forms parthenogenetically from the egg cell. To deepen the understanding of the transcriptional basis underlying these differences, we applied RNA-seq to compare expression in reproductive tissues of different Boechera accessions. This confirmed previous evidence of an enrichment of RNA helicases in plant germlines. Furthermore, few RNA helicases were differentially expressed in female reproductive ovule tissues harboring mature gametophytes from apomictic and sexual accessions. For some of these genes, we further found evidence for a complex recent evolutionary history. This included a homolog of Arabidopsis thaliana FASCIATED STEM4 (FAS4). In contrast to AtFAS4, which is a single-copy gene, FAS4 is represented by three homologs in Boechera, suggesting a potential for subfunctionalization to modulate reproductive development. To gain first insights into functional roles of FAS4, we studied Arabidopsis lines carrying mutant alleles. This identified the crucial importance of AtFAS4 for reproduction, as we observed developmental defects and arrest during male and female gametogenesis.

Forces driving transposable element load variation during Arabidopsis range expansion

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Machine learning on alignment features for parent-of-origin classification of simulated hybrid RNA-seq

BACKGROUND

Parent-of-origin allele-specific gene expression (ASE) can be detected in interspecies hybrids by virtue of RNA sequence variants between the parental haplotypes. ASE is detectable by differential expression analysis (DEA) applied to the counts of RNA-seq read pairs aligned to parental references, but aligners do not always choose the correct parental reference.

RESULTS

We used public data for species that are known to hybridize. We measured our ability to assign RNA-seq read pairs to their proper transcriptome or genome references. We tested software packages that assign each read pair to a reference position and found that they often favored the incorrect species reference. To address this problem, we introduce a post process that extracts alignment features and trains a random forest classifier to choose the better alignment. On each simulated hybrid dataset tested, our machine-learning post-processor achieved higher accuracy than the aligner by itself at choosing the correct parent-of-origin per RNA-seq read pair.

CONCLUSIONS

For the parent-of-origin classification of RNA-seq, machine learning can improve the accuracy of alignment-based methods. This approach could be useful for enhancing ASE detection in interspecies hybrids, though RNA-seq from real hybrids may present challenges not captured by our simulations. We believe this is the first application of machine learning to this problem domain.

Deleterious phenotypes in wild Arabidopsis arenosa populations are common and linked to runs of homozygosity

In this study, we aimed to systematically assess the frequency at which potentially deleterious phenotypes appear in natural populations of the outcrossing model plant Arabidopsis arenosa, and to establish their underlying genetics. For this purpose, we collected seeds from wild A. arenosa populations and screened over 2,500 plants for unusual phenotypes in the greenhouse. We repeatedly found plants with obvious phenotypic defects, such as small stature and necrotic or chlorotic leaves, among first-generation progeny of wild A. arenosa plants. Such abnormal plants were present in about 10% of maternal sibships, with multiple plants with similar phenotypes in each of these sibships, pointing to a genetic basis of the observed defects. A combination of transcriptome profiling, linkage mapping and genome-wide runs of homozygosity patterns using a newly assembled reference genome indicated a range of underlying genetic architectures associated with phenotypic abnormalities. This included evidence for homozygosity of certain genomic regions, consistent with alleles that are identical by descent being responsible for these defects. Our observations suggest that deleterious alleles with different genetic architectures are segregating at appreciable frequencies in wild A. arenosa populations.

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica

BACKGROUND

Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds.

RESULTS

In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay.

CONCLUSIONS

This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

Sequence similarity networks bear out hierarchical relationships of green cytochrome P450

Land plants have diversified enzyme families. One of the most prominent is the cytochrome P450 (CYP or CYP450) family. With over 443,000 CYP proteins sequenced across the tree of life, CYPs are ubiquitous in archaea, bacteria, and eukaryotes. Here, we focused on land plants and algae to study the role of CYP diversification. CYPs, acting as monooxygenases, catalyze hydroxylation reactions crucial for specialized plant metabolic pathways, including detoxification and phytohormone production; the CYPome consists of one enormous superfamily that is divided into clans and families. Their evolutionary history speaks of high substrate promiscuity; radiation and functional diversification have yielded numerous CYP families. To understand the evolutionary relationships within the CYPs, we employed sequence similarity network analyses. We recovered distinct clusters representing different CYP families, reflecting their diversified sequences that we link to the prediction of functionalities. Hierarchical clustering and phylogenetic analysis further elucidated relationships between CYP clans, uncovering their shared deep evolutionary history. We explored the distribution and diversification of CYP subfamilies across plant and algal lineages, uncovering novel candidates and providing insights into the evolution of these enzyme families. This identified unexpected relationships between CYP families, such as the link between CYP82 and CYP74, shedding light on their roles in plant defense signaling pathways. Our approach provides a methodology that brings insights into the emergence of new functions within the CYP450 family, contributing to the evolutionary history of plants and algae. These insights can be further validated and implemented via experimental setups under various external conditions.

tRNA-Cys gene clusters exhibit high variability in Arabidopsis thaliana

Although most of the genes encoding tRNAs in plants are dispersed throughout the genome, a fraction of them form tRNA gene clusters. In Arabidopsis thaliana, the smallest of tRNA clusters on chromosome 5 consists of four tRNA-Cys-GCA genes placed within repeating units of 0.4 kbp. A systematic analysis of the genomic sequences of syntenic regions from various ecotypes of A. thaliana showed that the general structure of the cluster, consisting of a tRNA-Cys pseudogene followed by repeating units containing tRNA-Cys genes, is well conserved. However, there is significant heterogeneity in the number of repeating units between different ecotypes. A unique feature of this cluster is the presence of putative transposable elements (Helitron). In addition, two further tRNA-Cys gene mini-clusters (gene pairs) in A. thaliana were identified. RNA-seq-based evaluation of expression of tRNA-Cys-GCA genes showed a positive signal for 11 out of 13 unique transcripts. An analysis of the conservation of the tRNA-Cys clusters from A. thaliana with the corresponding regions from four other Arabidopsis species suggests a sequence of events that led to the divergence of these regions.

Structural evidence for MADS-box type I family expansion seen in new assemblies of Arabidopsis arenosa and A. lyrata

Arabidopsis thaliana diverged from A. arenosa and A. lyrata at least 6 million years ago. The three species differ by genome-wide polymorphisms and morphological traits. The species are to a high degree reproductively isolated, but hybridization barriers are incomplete. A special type of hybridization barrier is based on the triploid endosperm of the seed, where embryo lethality is caused by endosperm failure to support the developing embryo. The MADS-box type I family of transcription factors is specifically expressed in the endosperm and has been proposed to play a role in endosperm-based hybridization barriers. The gene family is well known for its high evolutionary duplication rate, as well as being regulated by genomic imprinting. Here we address MADS-box type I gene family evolution and the role of type I genes in the context of hybridization. Using two de-novo assembled and annotated chromosome-level genomes of A. arenosa and A. lyrata ssp. petraea we analyzed the MADS-box type I gene family in Arabidopsis to predict orthologs, copy number, and structural genomic variation related to the type I loci. Our findings were compared to gene expression profiles sampled before and after the transition to endosperm cellularization in order to investigate the involvement of MADS-box type I loci in endosperm-based hybridization barriers. We observed substantial differences in type-I expression in the endosperm of A. arenosa and A. lyrata ssp. petraea, suggesting a genetic cause for the endosperm-based hybridization barrier between A. arenosa and A. lyrata ssp. petraea.

Role of phenotypic and transcriptomic plasticity in alpine adaptation of Arabidopsis arenosa

Plasticity is an important component of the response of organism to environmental changes, but whether plasticity facilitates adaptation is still largely debated. Using transcriptomic and phenotypic data, we explored the evolution of ancestral plasticity during alpine colonization in Arabidopsis arenosa. We leveraged naturally replicated adaptation in four distinct mountain regions in Central Europe. We sampled seeds from ancestral foothill and independently formed alpine populations in each region and raised them in growth chambers under conditions approximating their natural environments. We gathered RNA-seq and genetic data of 48 and 63 plants and scored vegetative and flowering traits in 203 and 272 plants respectively. Then, we compared gene expression and trait values over two treatments differing in temperature and irradiance and elevations of origin and quantified the extent of ancestral and derived plasticity. At the transcriptomic level, initial plastic changes tended to be more reinforced than reversed in adapted alpine populations. Genes showing reinforcement were involved in the stress response, developmental processes and morphogenesis and those undergoing reversion were related to the stress response (light and biotic stress). At the phenotypic level, initial plastic changes in all but one trait were also reinforced supporting a facilitating role of phenotypic plasticity during colonization of an alpine environment. Our results contrasted with previous studies that showed generally higher reversion than reinforcement and supported the idea that ancestral plasticity tends to be reinforced in the context of alpine adaptation. However, plasticity may also be the source of potential maladaptation, especially at the transcriptomic level.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

High-quality Cymbidium mannii genome and multifaceted regulation of crassulacean acid metabolism in epiphytes

Epiphytes with crassulacean acid metabolism (CAM) photosynthesis are widespread among vascular plants, and repeated evolution of CAM photosynthesis is a key innovation for micro-ecosystem adaptation. However, we lack a complete understanding of the molecular regulation of CAM photosynthesis in epiphytes. Here, we report a high-quality chromosome-level genome assembly of a CAM epiphyte, Cymbidium mannii (Orchidaceae). The 2.88-Gb orchid genome with a contig N50 of 22.7 Mb and 27 192 annotated genes was organized into 20 pseudochromosomes, 82.8% of which consisted of repetitive elements. Recent expansions of long terminal repeat retrotransposon families have made a major contribution to the evolution of genome size in Cymbidium orchids. We reveal a holistic scenario of molecular regulation of metabolic physiology using high-resolution transcriptomics, proteomics, and metabolomics data collected across a CAM diel cycle. Patterns of rhythmically oscillating metabolites, especially CAM-related products, reveal circadian rhythmicity in metabolite accumulation in epiphytes. Genome-wide analysis of transcript and protein level regulation revealed phase shifts during the multifaceted regulation of circadian metabolism. Notably, we observed diurnal expression of several core CAM genes (especially βCA and PPC) that may be involved in temporal fixation of carbon sources. Our study provides a valuable resource for investigating post-transcription and translation scenarios in C. mannii, an Orchidaceae model for understanding the evolution of innovative traits in epiphytes.

Transition to Self-compatibility Associated With Dominant S-allele in a Diploid Siberian Progenitor of Allotetraploid Arabidopsis kamchatica Revealed by Arabidopsis lyrata Genomes

A transition to selfing can be beneficial when mating partners are scarce, for example, due to ploidy changes or at species range edges. Here, we explain how self-compatibility evolved in diploid Siberian Arabidopsis lyrata, and how it contributed to the establishment of allotetraploid Arabidopsis kamchatica. First, we provide chromosome-level genome assemblies for two self-fertilizing diploid A. lyrata accessions, one from North America and one from Siberia, including a fully assembled S-locus for the latter. We then propose a sequence of events leading to the loss of self-incompatibility in Siberian A. lyrata, date this independent transition to ∼90 Kya, and infer evolutionary relationships between Siberian and North American A. lyrata, showing an independent transition to selfing in Siberia. Finally, we provide evidence that this selfing Siberian A. lyrata lineage contributed to the formation of the allotetraploid A. kamchatica and propose that the selfing of the latter is mediated by the loss-of-function mutation in a dominant S-allele inherited from A. lyrata.

Short Interrupted Repeat Cassette (SIRC)-Novel Type of Repetitive DNA Element Found in Arabidopsis thaliana

Short interrupted repeat cassette (SIRC)-a novel DNA element found throughout the A. thaliana nuclear genome. SIRCs are represented by short direct repeats interrupted by diverse DNA sequences. The maxima of SIRC's distribution are located within pericentromeric regions. We suggest that originally SIRC was a special case of the complex internal structure of the miniature inverted repeat transposable element (MITE), and further MITE amplification, transposition, and loss of terminal inverted repeats gave rise to SIRC as an independent DNA element. SIRC sites were significantly enriched with several histone modifications associated with constitutive heterochromatin and mobile genetic elements. The majority of DNA-binding proteins, strongly associated with SIRC, are related to histone modifications for transcription repression. A part of SIRC was found to overlap highly inducible protein-coding genes, suggesting a possible regulatory role for these elements, yet their definitive functions need further investigation.

Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences

Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].

Genome-wide analysis revealed the stepwise origin and functional diversification of HSDs from lower to higher plant species

Hydroxysteroid dehydrogenase (HSDs) is an oil-body sterol protein (steroleosin) with an NADP(H) binding domain that belongs to the short-chain dehydrogenase/reductase (SDR) superfamily. There are numerous studies on the characterization of HSDs in plants. However, thus far, the evolutionary differentiation and divergence analysis of these genes remain to be explored. The current study used an integrated method to elucidate the sequential evolution of HSDs in 64 sequenced plant genomes. Analyses were conducted on their origins, distribution, duplication, evolutionary paths, domain functions, motif composition, properties, and cis-elements. Results indicate that except for algae, HSD1 was widely distributed in plant species ranging from lower to higher plants, while HSD5 was restricted to terrestrial plants, and HSD2 was identified in fewer monocots and several dicot plants. Phylogenetic analysis of HSD proteins revealed that monocotyledonous HSD1 in moss and ferns appeared closest to the outgroup, V. carteri HSD-like, M. musculus HSD1, and H. sapiens HSD1. These data support the hypothesis that HSD1 originated in bryophytes and then in non-vascular and vascular plants, followed by HSD5 only in land plants. Gene structure analysis suggests that HSDs in plant species came up with a fixed number of six exons, and the intron phase was primarily 0, 1, 0, 0, and 0. Similarly, duplication analysis revealed that segmental duplications were the main reason for HSDs in plant species. Physicochemical properties suggest that dicotyledonous HSD1s and HSD5s were mainly acidic. The monocotyledonous HSD1s and HSD2s and the dicotyledonous HSD2s, HSD3s, HSD4s, and HSD6s were mainly basic, implying that HSDs in plants may have a variety of functions. Cis-regulatory elements and expression analysis revealed that HSDs in plants might have roles in several abiotic stresses. Due to the high expression of HSD1s and HSD5s in seeds, these HSDs in plants may have roles in fatty acid accumulation and degradation.

Cycles of satellite and transposon evolution in Arabidopsis centromeres

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation^1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox^3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.