Search

Lettuce (Lactuca sativaL.) is a leafy vegetable crop with ongoing breeding efforts related to quality, resilience, and innovative production systems. Genetic variation of important traits in close relatives is necessary to meet lettuce breeding goals.Lactuca virosa(2x=2n=18), a wild relative assigned to the tertiary lettuce gene pool, has a much larger genome (3.7 Gbp) thanLactuca sativa(2.5 Gbp). It has been used in interspecific crosses and is a donor to modern crisphead lettuce cultivars. Here, we present ade novoreference assembly ofL. virosawith high continuity and complete gene space. This assembly facilitated comparisons to the genome ofL. sativaand to that of the wild speciesL. saligna, a representative of the secondary lettuce gene pool. To assess the diversity in gene content, we classified the genes of the threeLactucaspecies as core, accessory and unique. In addition, we identified three interspecific chromosomal inversions compared toL. sativa, which each may cause recombination suppression and thus hamper future introgression breeding. Using three-way comparisons in both reference-based and reference-free manners, we show that the proliferation of long-terminal repeat elements has driven the genome expansion ofL. virosa. Further, we performed a genome-wide comparison of immune genes, nucleotide-binding leucine-rich repeat, and receptor-like kinases amongLactucaspp. and indicate the evolutionary patterns and mechanisms behind their expansions. These genome analyses greatly facilitate the understanding of genetic variation inL. virosa, which is beneficial for the breeding of improved lettuce varieties.

Lactuca saligna L. is a wild relative of cultivated lettuce (Lactuca sativa L.), with which it is partially interfertile. Hybrid progeny suffer from hybrid incompatibilities (HI), resulting in reduced fertility and distorted transmission ratios. Lactuca saligna displays broad spectrum resistance against lettuce downy mildew caused by Bremia lactucae Regel and is considered a non-host species. This phenomenon of resistance in L. saligna is called non-host resistance (NHR). One possible mechanism behind this NHR is through the plant–pathogen interaction triggered by pathogen-recognition receptors, including nucleotide-binding leucin-rich repeats (NLRs) and receptor-like kinases (RLKs). We report a chromosome-level genome assembly of L. saligna (accession CGN05327), leading to the identification of two large paracentric inversions (>50 Mb) between L. saligna and L. sativa. Genome-wide searches delineated the major resistance clusters as regions enriched in NLRs and RLKs. Three of the enriched regions co-locate with previously identified NHR intervals. RNA-seq analysis of Bremia infected lettuce identified several differentially expressed RLKs in NHR regions. Three tandem wall-associated kinase-encoding genes (WAKs) in the NHR8 interval display particularly high expression changes at an early stage of infection. We propose RLKs as strong candidate(s) for determinants for the NHR phenotype of L. saligna.

SummaryLactuca salignaL. is a wild relative of cultivated lettuce (Lactuca sativaL.), with which it is partially interfertile. Hybrid progeny suffer from hybrid incompatibilities (HI), resulting in reduced fertility and distorted transmission ratios.Lactuca salignadisplays broad spectrum resistance against lettuce downy mildew caused byBremia lactucaeRegel and is considered a non-host species. This phenomenon of resistance inL. salignais called non-host resistance (NHR). One possible mechanism behind this NHR is through the plant–pathogen interaction triggered by pathogen-recognition receptors, including nucleotide-binding leucin-rich repeats (NLRs) and receptor-like kinases (RLKs). We report a chromosome-level genome assembly ofL. saligna(accession CGN05327), leading to the identification of two large paracentric inversions (>50 Mb) betweenL. salignaandL. sativa. Genome-wide searches delineated the major resistance clusters as regions enriched inNLRs andRLKs. Three of the enriched regions co-locate with previously identified NHR intervals. RNA-seq analysis ofBremiainfected lettuce identified several differentially expressedRLKs in NHR regions. Three tandem wall-associated kinase-encoding genes (WAKs) in the NHR8 interval display particularly high expression changes at an early stage of infection. We proposeRLKs as strong candidate(s) for determinants for the NHR phenotype ofL. saligna.

1AbstractThe fungus Fusarium oxysporum is infamous for its devastating effects on economically important crops worldwide. F. oxysporum isolates are grouped into formae speciales based on their ability to cause disease on different hosts. Assigning F. oxysporum strains to formae speciales using non-experimental procedures has proven to be challenging due to their genetic heterogeneity and polyphyletic nature. However, genetically diverse isolates of the same forma specialis encode similar repertoires of effectors, proteins that are secreted by the fungus and contribute to the establishment of compatibility with the host. Based on this observation, we previously designed the F. oxysporum Effector Clustering (FoEC) pipeline which is able to classify F. oxysporum strains by forma specialis based on hierarchical clustering of the presence of predicted putative effector sequences, solely using genome assemblies as input. Here we present the updated FoEC2 pipeline which is more user friendly, customizable and, due to multithreading, has improved scalability. It is designed as a Snakemake pipeline and incorporates a new interactive visualization app. We showcase FoEC2 by clustering 537 publicly available F. oxysporum genomes and further analysis of putative effector families as multiple sequence alignments. We confirm classification of isolates into formae speciales and are able to further identify their subtypes. The pipeline is available on github: https://github.com/pvdam3/FoEC2.

KNOX genes encode transcription factors (TFs), several of which act non-cell-autonomously. KNOX genes evolved in algae, and two classes, class I KNOX and class II KNOX genes, were already present in charophytes. In tracheophytes, class I KNOX genes are expressed in shoot apical meristems (SAMs) and thought to inhibit cell differentiation, whereas class II KNOX genes are expressed in mature organs regulating differentiation. In this review, we summarize the data available on gene families and expression patterns of class I and class II KNOX genes in embryophytes. The expression patterns of class I KNOX genes should be seen in the context of SAM structure and of leaf primordium development where the inhibition of cell differentiation needs to be lifted. Although the SAMs of angiosperms and gnetophytes almost always belong to the duplex type, several other types are distributed in gymnosperms, ferns, lycopods and bryophytes. KNOX gene families remained small (maximally five genes) in the representatives of bryophytes, lycopods and ferns examined thus far; however, they expanded to some extent in gymnosperms and, independently and much more strongly, in angiosperms. The growing sophistication of mechanisms to repress and re-induce class KNOX I expression played a major role in the evolution of leaf shape.

Background Bacterial plant pathogens of the Pectobacterium genus are responsible for a wide spectrum of diseases in plants, including important crops such as potato, tomato, lettuce, and banana. Investigation of the genetic diversity underlying virulence and host specificity can be performed at genome level by using a comprehensive comparative approach called pangenomics. A pangenomic approach, using newly developed functionalities in PanTools, was applied to analyze the complex phylogeny of the Pectobacterium genus. We specifically used the pangenome to investigate genetic differences between virulent and avirulent strains of P. brasiliense, a potato blackleg causing species dominantly present in Western Europe. Results Here we generated a multilevel pangenome for Pectobacterium, comprising 197 strains across 19 species, including type strains, with a focus on P. brasiliense. The extensive phylogenetic analysis of the Pectobacterium genus showed robust distinct clades, with most detail provided by 452,388 parsimony-informative single-nucleotide polymorphisms identified in single-copy orthologs. The average Pectobacterium genome consists of 47% core genes, 1% unique genes, and 52% accessory genes. Using the pangenome, we zoomed in on differences between virulent and avirulent P. brasiliense strains and identified 86 genes associated to virulent strains. We found that the organization of genes is highly structured and linked with gene conservation, function, and transcriptional orientation. Conclusion The pangenome analysis demonstrates that evolution in Pectobacteria is a highly dynamic process, including gene acquisitions partly in clusters, genome rearrangements, and loss of genes. Pectobacterium species are typically not characterized by a set of species-specific genes, but instead present themselves using new gene combinations from the shared gene pool. A multilevel pangenomic approach, fusing DNA, protein, biological function, taxonomic group, and phenotypes, facilitates studies in a flexible taxonomic context.

Plant mitochondrial genomes are usually assembled and displayed as circular maps based on the widely-held view across the broad community of life scientists that circular genome-sized molecules are the primary form of plant mitochondrial DNA, despite the understanding by plant mitochondrial researchers that this is an inaccurate and outdated concept. Many plant mitochondrial genomes have one or more pairs of large repeats that can act as sites for inter- or intramolecular recombination, leading to multiple alternative arrangements (isoforms). Most mitochondrial genomes have been assembled using methods unable to capture the complete spectrum of isoforms within a species, leading to an incomplete inference of their structure and recombinational activity. To document and investigate underlying reasons for structural diversity in plant mitochondrial DNA, we used long-read (PacBio) and short-read (Illumina) sequencing data to assemble and compare mitochondrial genomes of domesticated (Lactuca sativa) and wild (L. saligna and L. serriola) lettuce species. We characterized a comprehensive, complex set of isoforms within each species and compared genome structures between species. Physical analysis of L. sativa mtDNA molecules by fluorescence microscopy revealed a variety of linear, branched, and circular structures. The mitochondrial genomes for L. sativa and L. serriola were identical in sequence and arrangement and differed substantially from L. saligna, indicating that the mitochondrial genome structure did not change during domestication. From the isoforms in our data, we infer that recombination occurs at repeats of all sizes at variable frequencies. The differences in genome structure between L. saligna and the two other Lactuca species can be largely explained by rare recombination events that rearranged the structure. Our data demonstrate that representations of plant mitochondrial genomes as simple, circular molecules are not accurate descriptions of their true nature and that in reality plant mitochondrial DNA is a complex, dynamic mixture of forms., Author summary Plant mitochondrial genomes are commonly depicted in research articles and textbooks as circular molecules that are the size of the genome. Although research on mitochondrial DNA (mtDNA) over the past few decades has revealed that genome-sized circles are exceedingly rare and that alternative forms of mtDNA are more common, many biologists still perceive circular maps as representing one or more physical chromosomes. This misconception can lead to biases in how mitochondrial genomes are assembled and misinterpretation of their evolutionary relationships, synteny, and histories. In this study, we present an assembly methodology that uses short- and long-read sequencing data to determine the mitochondrial genome structures of three lettuce species. We show that these mitochondrial genomes are fluid and dynamic, with multiple sequence arrangements of the genome coexisting within individuals of the same species. Differences in sequence arrangements between species can be explained by rare recombination events. Inspection of physical molecules of mtDNA reveals primarily non-circular forms. We demonstrate that plant mitochondrial genomes are a complex mixture of physical forms and sequence arrangements. Our data suggest that plant mitochondrial genomes should be presented as multiple sequence units showing their variable and dynamic connections, rather than as circles.

Lycopodiophyta—consisting of three orders, Lycopodiales, Isoetales and Selaginellales, with different types of shoot apical meristems (SAMs)—form the earliest branch among the extant vascular plants. They represent a sister group to all other vascular plants, from which they differ in that their leaves are microphylls—that is, leaves with a single, unbranched vein, emerging from the protostele without a leaf gap—not megaphylls. All leaves represent determinate organs originating on the flanks of indeterminate SAMs. Thus, leaf formation requires the suppression of indeterminacy, that is, of KNOX transcription factors. In seed plants, this is mediated by different groups of transcription factors including ARP and YABBY. We generated a shoot tip transcriptome of Huperzia selago (Lycopodiales) to examine the genes involved in leaf formation. Our H. selago transcriptome does not contain any ARP homolog, although transcriptomes of Selaginella spp. do. Surprisingly, we discovered a YABBY homolog, although these transcription factors were assumed to have evolved only in seed plants. The existence of a YABBY homolog in H. selago suggests that YABBY evolved already in the common ancestor of the vascular plants, and subsequently was lost in some lineages like Selaginellales, whereas ARP may have been lost in Lycopodiales. The presence of YABBY in the common ancestor of vascular plants would also support the hypothesis that this common ancestor had a simplex SAM. Furthermore, a comparison of the expression patterns of ARP in shoot tips of Selaginella kraussiana (Harrison CJ, etal. 2005. Independent recruitment of a conserved developmental mechanism during leaf evolution. Nature 434(7032):509–514.) and YABBY in shoot tips of H. selago implies that the development of microphylls, unlike megaphylls, does not seem to depend on the combined activities of ARP and YABBY. Altogether, our data show that Lycopodiophyta are a diverse group; so, in order to understand the role of Lycopodiophyta in evolution, representatives of Lycopodiales, Selaginellales, as well as of Isoetales, have to be examined.

How plants relate their requirements for energy with the reducing power necessary to fuel growth is not understood. The activated glucose forms and NADPH are key precursors in pathways yielding, respectively, energy and reducing power for anabolic metabolism. Moreover, they are substrates or allosteric regulators of trehalose-phosphate synthase (TPS1) in fungi and probably also in plants. TPS1 synthesizes the signalling metabolite trehalose-6-phosphate (T6P) and, therefore, has the potential to relate reducing power with energy metabolism to fuel growth. A working model is discussed where trehalose-6-phosphate (T6P) inhibition of SnRK1 is part of a growth-regulating loop in young and metabolically active heterotrophic plant tissues. SnRK1 is the Snf1 Related Kinase 1 and the plant homologue of the AMP-dependent protein kinase of animals, a central energy gauge. T6P accumulation in response to high sucrose levels in a cell inhibits SnRK1 activity, thus promoting anabolic processes and growth. When T6P levels drop due to low glucose-6-phosphate, uridine-diphosphoglucose, and altered NADPH or due to restricted TPS1 activity, active SnRK1 promotes catabolic processes required to respond to energy and carbon deprivation. The model explains why too little or too much T6P has been found to be growth inhibitory: Arabidopsis thaliana embryos and seedlings without TPS1 are growth arrested and Arabidopsis seedlings accumulating T6P on a trehalose medium are growth arrested. Finally, the insight gained with respect to the possible role of T6P metabolism, where it is known to alter developmental and environmental responses of plants, is discussed.

The first large dataset of hop EST sequences has recently become available to the hop research community. In this work we assembled all single pass hop EST sequences into tentative consensus sequences and used this information for development of molecular markers. A total of 14,944 EST sequences were acquired in a total length of 7,786,981 bp. After vector clipping, polyA clipping and removal of low-quality and short sequences, the remaining 14,731 sequences were clustered using the TGICL tool. Clustering yielded 1,318 contigs and only 76 singletons in a total length of 1,067,770 bp. A total of 13 clusters contained more than 100 EST sequences each. The 20 contigs with the highest number of ESTs clustered together were compared with blastn or blastx algorithms for simple annotation. The most abundant were ribosomal DNA and chloroplast transcripts. All sequences were subsequently screened for the presence of microsatellite motifs (2-10 bp motifs) and 235 contained such motifs. The most common type of repeat was dinucleotide (131), followed by trinucleotide (118). PCR primers were developed for chalcone synthase (CHS) and valerophenone synthase (VPS) genes and the sequences were checked for the presence of single nucleotide polymorphisms (SNPs). Markers were found in both genes and will be further used to map the two genes in the hop linkage map.

Organ formation in animals and plants relies on precise control of cell state transitions to turn stem cell daughters into fully differentiated cells. In plants, cells cannot rearrange due to shared cell walls. Thus, differentiation progression and the accompanying cell expansion must be tightly coordinated across tissues. PLETHORA (PLT) transcription factor gradients are unique in their ability to guide the progression of cell differentiation at different positions in the growing Arabidopsis thaliana root, which contrasts with well-described transcription factor gradients in animals specifying distinct cell fates within an essentially static context. To understand the output of the PLT gradient, we studied the gene set transcriptionally controlled by PLTs. Our work reveals how the PLT gradient can regulate cell state by region-specific induction of cell proliferation genes and repression of differentiation. Moreover, PLT targets include major patterning genes and autoregulatory feedback components, enforcing their role as master regulators of organ development.

Background: Following gene duplication, retained paralogs undergo functional divergence, which is reflected in changes in DNA sequence and expression patterns. The extent of divergence is influenced by several factors, including protein function. We examine whether an epigenetic modification, trimethylation of histone H3 at lysine 27 (H3K27me3), could be a factor in the evolution of expression patterns after gene duplication. Whereas in animals this repressive mark for transcription is deposited on long regions of DNA, in plants its localization is gene-specific. Because of this and a well-annotated recent whole-genome duplication, Arabidopsis thaliana is uniquely suited for studying the potential association of H3K27me3 with the evolutionary fate of genes. Results: Paralogous pairs with H3K27me3 show the highest coding sequence divergence, which can be explained by their low expression levels. Interestingly, they also show the highest similarity in expression patterns and upstream regulatory regions, while paralogous pairs where only one gene is an H3K27me3 target show the highest divergence in expression patterns and upstream regulatory sequence. These trends in divergence of expression and upstream regions are especially pronounced for transcription factors. Conclusions: After duplication, a histone modification can be associated with a particular fate of paralogs: H3K27me3 is linked to lower expression divergence yet higher coding sequence divergence. Our results show that H3K27me3 constrains expression divergence after duplication. Moreover, its association with higher conservation of upstream regions provides a potential mechanism for the conserved H3K27me3 targeting of the paralogs.