The RUBY reporter enabled the evaluation of different transgene expression constructs in lettuce, revealing that the lettuce ubiquitin promoter and terminator had strong expression that was stable over multiple generations. Nearly four decades after the first transgenic lettuce was reported, constructs for stable transgene expression remain limited. Notably, the 35S promoter from the Cauliflower Mosaic Virus (35S), which drives strong expression of transgenes in several plant species, has often shown silencing and instability in lettuce. Other promoter/terminator combinations that are commonly used in plant expression vectors have not been extensively studied in lettuce. In this study, we evaluated three different expression constructs in two different horticultural types of lettuce using the non-invasive RUBY reporter, which allowed for the monitoring of transgene expression throughout the process of regeneration during tissue culture, throughout development of the primary transgenics, and in two subsequent sexual generations. The LsUBI promoter/terminator combination resulted in strong, uniform expression throughout regeneration, during growth of the primary transgenics, and in both subsequent generations. The AtUBI promoter/tRBCS combination showed slightly lower levels of expression and intermediate levels of silencing, while the 35S promoter/tHSP combination showed both initial strong expression and frequent silencing. Therefore, our data show that the LsUBI promoter/terminator combination provides strong, uniform expression that is unlikely to result in silencing and that the AtUBI promoter/tRBCS combination is an additional option for stable expression of transgenes in lettuce, especially if an intermediate expression level is desired.

Nearly four decades after the first transgenic lettuce was reported, constructs for stable transgene expression remain limited. Notably, the 35S promoter from the Cauliflower Mosaic Virus (35S), which drives strong expression of transgenes in several plant species, has often shown silencing and instability in lettuce. Other promoter/terminator combinations that are commonly used in plant expression vectors have not been extensively studied in lettuce. In this study, we evaluated three different expression constructs in two different horticultural types of lettuce using the non-invasive RUBY reporter, which allowed for the monitoring of transgene expression throughout the process of regeneration during tissue culture, throughout development of the primary transgenics, and in two subsequent sexual generations. The LsUBI promoter/terminator combination resulted in strong, uniform expression throughout regeneration, during growth of the primary transgenics, and in both subsequent generations. The AtUBI promoter/tRBCS combination showed slightly lower levels of expression and intermediate levels of silencing, while the 35S promoter/tHSP combination showed both initial strong expression and frequent silencing. Therefore, our data show that the LsUBI promoter/terminator combination provides strong, uniform expression that is unlikely to result in silencing and that the AtUBI promoter/tRBCS combination is an additional option for stable expression of transgenes in lettuce, especially if an intermediate expression level is desired. ### Competing Interest Statement The authors have declared no competing interest.

Same Cas9 protein, two different jobs: the CRISPR-Combo genome engineering strategy enables simultaneous gene activation and genome editing for different targets through changes to the guide RNA structure.

The change in allele frequencies within a population over time represents a fundamental process of evolution. By monitoring allele frequencies, we can analyze the effects of natural selection and genetic drift on populations. To efficiently track time-resolved genetic change, large experimental or wild populations can be sequenced as pools of individuals sampled over time using high-throughput genome sequencing (called the Evolve & Resequence approach, E&R). Here, we present a set of experiments using hundreds of natural genotypes of the model plant Arabidopsis thaliana to showcase the power of this approach to study rapid evolution at large scale. First, we validate that sequencing DNA directly extracted from pools of flowers from multiple plants -- organs that are relatively consistent in size and easy to sample -- produces comparable results to other, more expensive state-of-the-art approaches such as sampling and sequencing of individual leaves. Sequencing pools of flowers from 25-50 individuals at ∼40X coverage recovers genome-wide frequencies in diverse populations with accuracy r > 0.95. Secondly, to enable analyses of evolutionary adaptation using E&R approaches of plants in highly replicated environments, we provide open source tools that streamline sequencing data curation and calculate various population genetic statistics two orders of magnitude faster than current software. To directly demonstrate the usefulness of our method, we conducted a two-year outdoor evolution experiment with A. thaliana to show signals of rapid evolution in multiple genomic regions. We demonstrate how these laboratory and computational Pool-seq-based methods can be scaled to study hundreds of populations across many climates.

During meiosis, DNA double‐strand breaks undergo interhomolog repair to yield crossovers between homologous chromosomes. To investigate how interhomolog sequence polymorphism affects crossovers, we sequenced multiple recombinant populations of the model plant Arabidopsis thaliana. Crossovers were elevated in the diverse pericentromeric regions, showing a local preference for polymorphic regions. We provide evidence that crossover association with elevated diversity is mediated via the Class I crossover formation pathway, although very high levels of diversity suppress crossovers. Interhomolog polymorphism causes mismatches in recombining molecules, which can be detected by MutS homolog (MSH) mismatch repair protein heterodimers. Therefore, we mapped crossovers in a msh2 mutant, defective in mismatch recognition, using multiple hybrid backgrounds. Although total crossover numbers were unchanged in msh2 mutants, recombination was remodelled from the diverse pericentromeres towards the less‐polymorphic sub‐telomeric regions. Juxtaposition of megabase heterozygous and homozygous regions causes crossover remodelling towards the heterozygous regions in wild type Arabidopsis, but not in msh2 mutants. Immunostaining showed that MSH2 protein accumulates on meiotic chromosomes during prophase I, consistent with MSH2 regulating meiotic recombination. Our results reveal a pro‐crossover role for MSH2 in regions of higher sequence diversity in A. thaliana. Meiotic DNA double‐strand breaks undergo interhomolog repair to yield crossovers between homologous chromosomes. Sequencing of polymorphisms in wild‐type and mutant plants reveals an unexpected pro‐crossover role of the mismatch repair protein MSH2 in regions of higher relative sequence diversity. Sequencing of polymorphisms between homologous chromosomes in wild‐type and mutant plants reveals an unexpected pro‐crossover role of the mismatch repair protein MSH2 in regions of higher relative sequence diversity.

An important component of cellular biochemistry is the concentration of proteins and nucleic acids in non-membranous compartments(1,2). These biomolecular condensates are formed from processes that include liquid-liquid phase separation. The multivalent interactions necessary for liquid-liquid phase separation have been extensively studied in vitro(1,3). However, the regulation of this process in vivo is poorly understood. Here we identify an in vivo regulator of liquid-liquid phase separation through a genetic screen targeting factors required for Arabidopsis RNA-binding protein FCA function. FCA contains prion-like domains that phase-separate in vitro, and exhibits behaviour in vivo that is consistent with phase separation. The mutant screen identified a functional requirement for FLL2, a coiled-coil protein, in the formation of FCA nuclear bodies. FCA reduces transcriptional read-through by promoting proximal polyadenylation at many sites in the Arabidopsis genome(3,4). FLL2 was required to promote this proximal polyadenylation, but not the binding of FCA to target RNA. Ectopic expression of FLL2 increased the size and number of FCA nuclear bodies. Crosslinking with formaldehyde captured in vivo interactions between FLL2, FCA and the polymerase and nuclease modules of the RNA 3'-end processing machinery. These 3' RNA-processing components colocalized with FCA in the nuclear bodies in vivo, which indicates that FCA nuclear bodies compartmentalize 3'-end processing factors to enhance polyadenylation at specific sites. Our findings show that coiled-coil proteins can promote liquid-liquid phase separation, which expands our understanding of the principles that govern the in vivo dynamics of liquid-like bodies.

Rowan et al. generated a dataset of over 17,000 meiotic crossovers (COs) from over 2000 F2 individuals from a single Arabidopsis thaliana cross. The unprecedented density of COs and the high-quality reference genomes of the two... Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana. COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

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.

Meiotic crossovers (COs) ensure proper chromosome segregation and redistribute the genetic variation that is transmitted to the next generation. Large populations and the demand for genome-wide, fine-scale resolution challenge existing methods for CO identification. Taking advantage of linked-read sequencing, we develop a highly efficient method for genome-wide identification of COs at kilobase resolution in pooled recombinants. We first test this method using a pool of Arabidopsis F2 recombinants, and recapitulate results obtained from the same plants using individual whole-genome sequencing. By applying this method to a pool of pollen DNA from an F1 plant, we establish a highly accurate CO landscape without generating or sequencing a single recombinant plant. The simplicity of this approach enables the simultaneous generation and analysis of multiple CO landscapes, accelerating the pace at which mechanisms for the regulation of recombination can be elucidated through efficient comparisons of genotypic and environmental effects on recombination.

Hybrids occasionally exhibit genetic interactions resulting in reduced fitness in comparison to their parents. Studies of Arabidopsis thaliana have highlighted the role of immune conflicts, but less is known about the role of other factors in hybrid incompatibility in plants. Here, we present a new hybrid incompatibility phenomenon in this species. We have characterized a new case of F1 hybrid incompatibility from a cross between the A. thaliana accessions Krotzenburg-0 (Kro-0) and BG-5, by conducting transcript, metabolite and hormone analyses, and identified the causal loci through genetic mapping. The F1 hybrids showed arrested growth of the main stem, altered shoot architecture, and altered concentrations of hormones in comparison to parents. The F1 phenotype could be rescued in a developmental-stage-dependent manner by shifting to a higher growth temperature. These F1 phenotypes were linked to two loci, one on chromosome 2 and one on chromosome 3. The F2 generation segregated plants with more severe phenotypes which were linked to the same loci as those in the F1 . This study provides novel insights into how previously unknown mechanisms controlling shoot branching and stem growth can result in hybrid incompatibility.

The growing collection of sequenced or genotyped Arabidopsis thaliana accessions includes mostly individuals from the native Eurasian and N. African range and introduced North American populations. Here, we describe the genetic and phenotypic diversity, along with habitats and life history, of A. thaliana plants collected at the southernmost end of its worldwide distribution. Seed samples were harvested from plants growing in four sites within a similar to 3500-km(2)-area in Patagonia, Argentina, and represent the first germplasm to be collected in South America for this species. Whole-genome resequencing revealed that plants from the four sites and a Patagonia herbarium specimen collected in 1967 formed a single haplogroup (Pat), indicating that the phenotypic variation observed in the field reflected plastic responses to the environment. ADMIXTURE and principal components analyses suggest that the ancestor of the Pat haplogroup either came from Italy or the Balkan/Caucasus regions of Eurasia. In the laboratory, plants from the Pat haplogroup were hyposensitive to continuous red (Rc) and shade light, with corresponding changes in the expression of phytochrome signalling genes. Pat had higher PIF3 and PIF5 and lower HY5 expression under Rc light; and lower expression of PIL1, ATHB2 and HFR1 under shade compared to Col-0. In addition, Pat plants had a strong vernalization requirement associated with high levels of FLC expression. We conclude that including Pat in studies of natural variation and in comparison with other introduced populations will provide additional information for association studies and allow for a more detailed assessment of the demographic events following colonization.

The gold standard for studying natural selection and adaptation in the wild is to quantify lifetime fitness of individuals from natural populations that have been grown together in a common garden, or that have been reciprocally transplanted. By combining fitness values with species traits and genome sequences, one can infer selection coefficients at the genetic level. Here we present a rainfall-manipulation experiment with 517 whole-genome sequenced natural accessions of the plant Arabidopsis thaliana spanning the global distribution of the species. The experiments were conducted in two field stations in contrasting climates, in the Mediterranean and in Central Europe, where we built rainout shelters and simulated high and low rainfall. Using custom image analysis we quantified fitness- and phenology-related traits for 23,154 pots, which contained about 14,500 plants growing independently, and over 310,000 plants growing in small populations (max. 30 plants). This large field experiment dataset, which associates fitness and ecologically-relevant traits with genomes, will provide an important resource to test eco-evolutionary genetic theories and to understand the potential evolutionary impacts of future climates on an important plant model species.

The redox imbalanced 6 mutant (rimb6) of Arabidopsis thaliana was isolated in a genetic screening approach for mutants with defects in chloroplast-to-nucleus redox signaling. It has an atypically low activation status of the 2-Cys peroxiredoxin-A promoter in the seedling stage. rimb6 shows wildtype-like germination, seedling development and greening, but slower growth and reduced biomass in the rosette stage. Mapping of the casual mutation revealed that rimb6 carries a single nucleotide polymorphism in the gene encoding CONSTITUTIVE EXPRESSER OF PATHOGENESIS RELATED (PR) GENES 1, CPR1 (At4g12560), leading to a premature stop codon. CPR1 is known as a repressor of pathogen signaling and regulator of microtubule organization. Allelism of rimb6 and cpr1 revealed a function of CPR1 in chloroplast stress protection. Expression studies in pathogen signaling mutants demonstrated that CPR1-mediated activation of genes for photosynthesis and chloroplast antioxidant protection is, in contrast to activation of pathogen responses, regulated independently from PAD4-controlled salicylic acid (SA) accumulation. We conclude that the support of plastid function is a basic, SA-independent function of CPR1.

A major goal for biologists is to understand the connection between genes and phenotypic traits, and genetic mapping in experimental populations remains a powerful approach for discovering the causal genes underlying phenotypes. For genetic mapping, the process of genotyping was previously a major rate-limiting step. Modern sequencing technology has greatly improved the resolution and speed of genetic mapping by reducing the time, labor, and cost per genotyping marker. In addition, the ability to perform genotyping-by-sequencing (GBS) has facilitated large-scale population genetic analyses by providing a simpler way to survey segregating genetic variation in natural populations. Here we present two protocols for GBS, using the Illumina platform, that can be applied to a wide range of genotyping projects in different species. The first protocol is for genotyping a subset of marker positions genome-wide using restriction digestion, and the second is for preparing inexpensive paired-end whole-genome libraries. We discuss the suitability of each approach for different genotyping applications and provide notes for adapting these protocols for use with a liquid-handling robot.

Arabidopsis thaliana serves as a model organism for the study of fundamental physiological, cellular, and molecular processes. It has also greatly advanced our understanding of intraspecific genome variation. We present a detailed map of variation in 1,135 high-quality re-sequenced natural inbred lines representing the native Eurasian and North African range and recently colonized North America. We identify relict populations that continue to inhabit ancestral habitats, primarily in the Iberian Peninsula. They have mixed with a lineage that has spread to northern latitudes from an unknown glacial refugium and is now found in a much broader spectrum of habitats. Insights into the history of the species and the fine-scale distribution of genetic diversity provide the basis for full exploitation of A. thaliana natural variation through integration of genomes and epigenomes with molecular and non-molecular phenotypes.

The reshuffling of existing genetic variation during meiosis is important both during evolution and in breeding. The reassortment of genetic variants relies on the formation of crossovers (COs) between homologous chromosomes. The pattern of genome-wide CO distributions can be rapidly and precisely established by the short-read sequencing of individuals from F-2 populations, which in turn are useful for quantitative trait locus (QTL) mapping. Although sequencing costs have decreased precipitously in recent years, the costs of library preparation for hundreds of individuals have remained high. To enable rapid and inexpensive CO detection and QTL mapping using low-coverage whole-genome sequencing of large mapping populations, we have developed a new method for library preparation along with Trained Individual GenomE Reconstruction, a probabilistic method for genotype and CO predictions for recombinant individuals. In an example case with hundreds of F-2 individuals from two Arabidopsis thaliana accessions, we resolved most CO breakpoints to within 2 kb and reduced a major flowering time QTL to a 9-kb interval. In addition, an extended region of unusually low recombination revealed a 1.8-Mb inversion polymorphism on the long arm of chromosome 4. We observed no significant differences in the frequency and distribution of COs between F-2 individuals with and without a functional copy of the DNA helicase gene RECQ4A. In summary, we present a new, cost-efficient method for large-scale, high-precision genotyping-by-sequencing.

Motivation: Genetic heterogeneity, the fact that several sequence variants give rise to the same phenotype, is a phenomenon that is of the utmost interest in the analysis of complex phenotypes. Current approaches for finding regions in the genome that exhibit genetic heterogeneity suffer from at least one of two shortcomings: (i) they require the definition of an exact interval in the genome that is to be tested for genetic heterogeneity, potentially missing intervals of high relevance, or (ii) they suffer from an enormous multiple hypothesis testing problem due to the large number of potential candidate intervals being tested, which results in either many false positives or a lack of power to detect true intervals. Results: Here, we present an approach that overcomes both problems: it allows one to automatically find all contiguous sequences of single nucleotide polymorphisms in the genome that are jointly associated with the phenotype. It also solves both the inherent computational efficiency problem and the statistical problem of multiple hypothesis testing, which are both caused by the huge number of candidate intervals. We demonstrate on Arabidopsis thaliana genome-wide association study data that our approach can discover regions that exhibit genetic heterogeneity and would be missed by single-locus mapping. Conclusions: Our novel approach can contribute to the genome-wide discovery of intervals that are involved in the genetic heterogeneity underlying complex phenotypes. Availability and implementation: The code can be obtained at: http://www.bsse.ethz.ch/mlcb/research/bioinformatics-and-computational-biology/sis.html. Contact: felipe.llinares@bsse.ethz.ch Supplementary information : Supplementary data are available at Bioinformatics online.

In their commentary, [Golczyk et al. (2014)][1] state that the level of plastid DNA (ptDNA) remains approximately constant during mature, aging, and senescent stages of leaf development, in apparent contrast with our conclusion of ptDNA decline from previous studies. However, the ptDNA decline that