Genome Improvement and Genetic Map Construction for Aethionema arabicum, the First Divergent Branch in the Brassicaceae Family

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.

The dimorphic diaspore model Aethionema arabicum (Brassicaceae): Distinct molecular and morphological control of responses to parental and germination temperatures

Plants in habitats with unpredictable conditions often have diversified bet-hedging strategies that ensure fitness over a wider range of variable environmental factors. A striking example is the diaspore (seed and fruit) heteromorphism that evolved to maximize species survival in Aethionema arabicum (Brassicaceae) in which external and endogenous triggers allow the production of two distinct diaspores on the same plant. Using this dimorphic diaspore model, we identified contrasting molecular, biophysical, and ecophysiological mechanisms in the germination responses to different temperatures of the mucilaginous seeds (M+ seed morphs), the dispersed indehiscent fruits (IND fruit morphs), and the bare non-mucilaginous M- seeds obtained by pericarp (fruit coat) removal from IND fruits. Large-scale comparative transcriptome and hormone analyses of M+ seeds, IND fruits, and M- seeds provided comprehensive datasets for their distinct thermal responses. Morph-specific differences in co-expressed gene modules in seeds, as well as in seed and pericarp hormone contents, identified a role of the IND pericarp in imposing coat dormancy by generating hypoxia affecting abscisic acid (ABA) sensitivity. This involved expression of morph-specific transcription factors, hypoxia response, and cell wall remodeling genes, as well as altered ABA metabolism, transport, and signaling. Parental temperature affected ABA contents and ABA-related gene expression and altered IND pericarp biomechanical properties. Elucidating the molecular framework underlying the diaspore heteromorphism can provide insight into developmental responses to globally changing temperatures.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

Single-molecule real-time sequencing of the full-length transcriptome of purple garlic (Allium sativum L. cv. Leduzipi) and identification of serine O-acetyltransferase family proteins involved in cysteine biosynthesis

BACKGROUND

Garlic (Allium sativum L.), whose bioactive components are mainly organosulfur compounds (OSCs), is a herbaceous perennial widely consumed as a green vegetable and a condiment. Yet, the metabolic enzymes involved in the biosynthesis of OSCs are not identified in garlic.

RESULTS

Here, a full-length transcriptome of purple garlic was generated via PacBio and Illumina sequencing, to characterize the garlic transcriptome and identify key proteins mediating the biosynthesis of OSCs. Overall, 22.56 Gb of clean data were generated, resulting in 454 698 circular consensus sequence (CCS) reads, of which 83.4% (379 206) were identified as being full-length non-chimeric reads - their further transcript clustering facilitated identification of 36 571 high-quality consensus reads. Once corrected, their genome-wide mapping revealed that 6140 reads were novel isoforms of known genes, and 2186 reads were novel isoforms from novel genes. We detected 1677 alternative splicing events, finding 2902 genes possessing either two or more poly(A) sites. Given the importance of serine O-acetyltransferase (SERAT) in cysteine biosynthesis, we investigated the five SERAT homologs in garlic. Phylogenetic analysis revealed a three-tier classification of SERAT proteins, each featuring a serine acetyltransferase domain (N-terminal) and one or two hexapeptide transferase motifs. Template-based modeling showed that garlic SERATs shared a common homo-trimeric structure with homologs from bacteria and other plants. The residues responsible for substrate recognition and catalysis were highly conserved, implying a similar reaction mechanism. In profiling the five SERAT genes' transcript levels, their expression pattern varied significantly among different tissues.

CONCLUSION

This study's findings deepen our knowledge of SERAT proteins, and provide timely genetic resources that could advance future exploration into garlic's genetic improvement and breeding. © 2021 Society of Chemical Industry.

A tale of two morphs: developmental patterns and mechanisms of seed coat differentiation in the dimorphic diaspore model Aethionema arabicum (Brassicaceae)

The developmental transition from a fertilized ovule to a dispersed diaspore (seed or fruit) involves complex differentiation processes of the ovule's integuments leading to the diversity in mature seed coat structures in angiosperms. In this study, comparative imaging and transcriptome analysis were combined to investigate the morph-specific developmental differences during outer seed coat differentiation and mucilage production in Aethionema arabicum, the Brassicaceae model for diaspore dimorphism. One of the intriguing adaptations of this species is the production and dispersal of morphologically distinct, mucilaginous and non-mucilaginous diaspores from the same plant (dimorphism). The dehiscent fruit morph programme producing multiple mucilaginous seed diaspores was used as the default trait combination, similar to Arabidopsis thaliana, and was compared with the indehiscent fruit morph programme leading to non-mucilaginous diaspores. Synchrotron-based radiation X-ray tomographic microscopy revealed a co-ordinated framework of morph-specific early changes in internal anatomy of developing A. arabicum gynoecia including seed abortion in the indehiscent programme and mucilage production by the mucilaginous seed coat. The associated comparative analysis of the gene expression patterns revealed that the unique seed coat dimorphism of Ae. arabicum provides an excellent model system for comparative study of the control of epidermal cell differentiation and mucilage biosynthesis by the mucilage transcription factor cascade and their downstream cell wall and mucilage remodelling genes. Elucidating the underlying molecular framework of the dimorphic diaspore syndrome is key to understanding differential regulation of bet-hedging survival strategies in challenging environments, timely in the face of global climatic change.

Aethionema arabicum genome annotation using PacBio full-length transcripts provides a valuable resource for seed dormancy and Brassicaceae evolution research

Aethionema arabicum is an important model plant for Brassicaceae trait evolution, particularly of seed (development, regulation, germination, dormancy) and fruit (development, dehiscence mechanisms) characters. Its genome assembly was recently improved but the gene annotation was not updated. Here, we improved the Ae. arabicum gene annotation using 294 RNA-seq libraries and 136 307 full-length PacBio Iso-seq transcripts, increasing BUSCO completeness by 11.6% and featuring 5606 additional genes. Analysis of orthologs showed a lower number of genes in Ae. arabicum than in other Brassicaceae, which could be partially explained by loss of homeologs derived from the At-α polyploidization event and by a lower occurrence of tandem duplications after divergence of Aethionema from the other Brassicaceae. Benchmarking of MADS-box genes identified orthologs of FUL and AGL79 not found in previous versions. Analysis of full-length transcripts related to ABA-mediated seed dormancy discovered a conserved isoform of PIF6-β and antisense transcripts in ABI3, ABI4 and DOG1, among other cases found of different alternative splicing between Turkey and Cyprus ecotypes. The presented data allow alternative splicing mining and proposition of numerous hypotheses to research evolution and functional genomics. Annotation data and sequences are available at the Ae. arabicum DB (https://plantcode.online.uni-marburg.de/aetar_db).

New model species for arctic-alpine plant molecular ecology

Arctic and alpine, high latitude and high elevation environments are one of the most stressful environments for species to inhabit. This harshness manifests itself in lower species richness in comparison to more southern vegetation zones (Francis & Currie, 2003). Furthermore, the climatic oscillations-past and predicted-have the most dramatic effect on these ecosystems. For example, in regions of continental ice sheets-the northernmost part of Western Europe and North America-the Arctic species assemblages are no older than a few thousands of years, which is a relatively short period from an evolutionary perspective. Although similar environments may have existed further south during the Ice Age, allowing some preadaptation for the Arctic species, the current habitat is a unique combination of environmental factors such as the climate, soil, bedrock, and photoperiod. Hence, understanding the evolutionary forces shaping Arctic-alpine species will be important for predicting these vulnerable environments' population viability and adaptive potential in the future. In this issue of Molecular Ecology Resources, Nowak et al. (Molecular Ecology Resources) present extensive genome-wide resources for an Arctic-alpine plant Draba nivalis. This adds a valuable new member into the cabbage family models for evolutionary genetics and adaptation studies, to accompany e.g., Arabidopsis (Nature Genetics, 43, 476; Nature, 408, 796), Arabis (Nature Plants, 1, 14023) and Capsella (Nature Genetics, 45, 831). A whole new avenue will open up for molecular ecological studies not only for D. nivalis, but the whole large Draba genus with its diverse ecological and evolutionary characteristics.

Single Nucleotide Polymorphism Charting of P. patens Reveals Accumulation of Somatic Mutations During in vitro Culture on the Scale of Natural Variation by Selfing

Introduction: Physcomitrium patens (Hedw.) Mitten (previously known as Physcomitrella patens) was collected by H.L.K. Whitehouse in Gransden Wood (Huntingdonshire, United Kingdom) in 1962 and distributed across the globe starting in 1974. Hence, the Gransden accession has been cultured in vitro in laboratories for half a century. Today, there are more than 13 different pedigrees derived from the original accession. Additionally, accessions from other sites worldwide were collected during the last decades. Methods and Results: In this study, 250 high throughput RNA sequencing (RNA-seq) samples and 25 gDNA samples were used to detect single nucleotide polymorphisms (SNPs). Analyses were performed using five different P. patens accessions and 13 different Gransden pedigrees. SNPs were overlaid with metadata and known phenotypic variations. Unique SNPs defining Gransden pedigrees and accessions were identified and experimentally confirmed. They can be successfully employed for PCR-based identification. Conclusion: We show independent mutations in different Gransden laboratory pedigrees, demonstrating that somatic mutations occur and accumulate during in vitro culture. The frequency of such mutations is similar to those observed in naturally occurring populations. We present evidence that vegetative propagation leads to accumulation of deleterious mutations, and that sexual reproduction purges those. Unique SNP sets for five different P. patens accessions were isolated and can be used to determine individual accessions as well as Gransden pedigrees. Based on that, laboratory methods to easily determine P. patens accessions and Gransden pedigrees are presented.

High-Density Genetic Linkage Mapping of Lepidium Based on Genotyping-by-Sequencing SNPs and Segregating Contig Tag Haplotypes

Lepidium campestre has been targeted for domestication as future oilseed and catch crop. Three hundred eighty plants comprising genotypes of L. campestre, Lepidium heterophyllum, and their interspecific F2 mapping population were genotyped using genotyping by sequencing (GBS), and the generated polymorphic markers were used for the construction of high-density genetic linkage map. TASSEL-GBS, a reference genome-based pipeline, was used for this analysis using a draft L. campestre whole genome sequence. The analysis resulted in 120,438 biallelic single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) above 0.01. The construction of genetic linkage map was conducted using MSTMap based on phased SNPs segregating in 1:2:1 ratio for the F2 individuals, followed by genetic mapping of segregating contig tag haplotypes as dominant markers against the linkage map. The final linkage map consisted of eight linkage groups (LGs) containing 2,330 SNP markers and spanned 881 Kosambi cM. Contigs (10,302) were genetically mapped to the eight LGs, which were assembled into pseudomolecules that covered a total of ∼120.6 Mbp. The final size of the pseudomolecules ranged from 9.4 Mbp (LG-4) to 20.4 Mpb (LG-7). The following major correspondence between the eight Lepidium LGs (LG-1 to LG-8) and the five Arabidopsis thaliana (At) chromosomes (Atx-1-Atx-5) was revealed through comparative genomics analysis: LG-1&2_Atx-1, LG-3_Atx-2&3, LG-4_Atx-2, LG-5_Atx-2&Atx-3, LG-6_Atx-4&5, LG-7_Atx-4, and LG-8_Atx-5. This analysis revealed that at least 66% of the sequences of the LGs showed high collinearity with At chromosomes. The sequence identity between the corresponding regions of the LGs and At chromosomes ranged from 80.6% (LG-6) to 86.4% (LG-8) with overall mean of 82.9%. The map positions on Lepidium LGs of the homologs of 24 genes that regulate various traits in A. thaliana were also identified. The eight LGs revealed in this study confirm the previously reported (1) haploid chromosome number of eight in L. campestre and L. heterophyllum and (2) chromosomal fusion, translocation, and inversion events during the evolution of n = 8 karyotype in ancestral species shared by Lepidium and Arabidopsis to n = 5 karyotype in A. thaliana. This study generated highly useful genomic tools and resources for Lepidium that can be used to accelerate its domestication.

Genomic Blocks in Aethionema arabicum Support Arabideae as Next Diverging Clade in Brassicaceae

The tribe Aethionemeae is sister to all other crucifers, making it a crucial group for unraveling genome evolution and phylogenetic relationships within the crown group Brassicaceae. In this study, we extend the analysis of Brassicaceae genomic blocks (GBs) to Aethionema whereby we identified unique block boundaries shared only with the tribe Arabideae. This was achieved using bioinformatic methods to analyze synteny between the recently updated genome sequence of Aethionema arabicum and other high-quality Brassicaceae genome sequences. We show that compared to the largely conserved genomic structure of most non-polyploid Brassicaceae lineages, GBs are highly rearranged in Aethionema. Furthermore, we detected similarities between the genomes of Aethionema and Arabis alpina, in which also a high number of genomic rearrangements compared to those of other Brassicaceae was found. These similarities suggest that tribe Arabideae, a clade showing conflicting phylogenetic position between studies, may have diverged before diversification of the other major lineages, and highlight the potential of synteny information for phylogenetic inference.