A simulation study comparing advanced marker-assisted selection with genomic selection in tree breeding programs.

Advances in DNA sequencing technologies allow the sequencing of whole genomes of thousands of individuals and provide several million single nucleotide polymorphisms (SNPs) per individual. These data combined with precise and high-throughput phenotyping enable genome-wide association studies (GWAS) and the identification of SNPs underlying traits with complex genetic architectures. The identified causal SNPs and estimated allelic effects could then be used for advanced marker-assisted selection (MAS) in breeding programs. But could such MAS compete with the broadly used genomic selection (GS)? This question is of particular interest for the lengthy tree breeding strategies. Here, with our new software "SNPscan breeder," we simulated a simple tree breeding program and compared the impact of different selection criteria on genetic gain and inbreeding. Further, we assessed different genetic architectures and different levels of kinship among individuals of the breeding population. Interestingly, apart from progeny testing, GS using gBLUP performed best under almost all simulated scenarios. MAS based on GWAS results outperformed GS only if the allelic effects were estimated in large populations (ca. 10,000 individuals) of unrelated individuals. Notably, GWAS using 3,000 extreme phenotypes performed as good as the use of 10,000 phenotypes. GS increased inbreeding and thus reduced genetic diversity more strongly compared to progeny testing and GWAS-based selection. We discuss the practical implications for tree breeding programs. In conclusion, our analyses further support the potential of GS for forest tree breeding and improvement, although MAS may gain relevance with decreasing sequencing costs in the future.

ARR17 controls dioecy in Populus by repressing B-class MADS-box gene expression.

The number of dioecious species for which the genetic basis of sex determination has been resolved is rapidly increasing. Nevertheless, the molecular mechanisms downstream of the sex determinants remain largely elusive. Here, by RNA-sequencing early-flowering isogenic aspen (Populus tremula) lines differing exclusively for the sex switch gene ARR17, we show that a narrowly defined genetic network controls differential development of female and male flowers. Although ARR17 encodes a type-A response regulator supposedly involved in cytokinin (CK) hormone signalling, clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9-mediated arr17 knockout only affected the expression of a strikingly small number of genes, indicating a specific role in the regulation of floral development rather than a generic function in hormone signalling. Notably, the UNUSUAL FLORAL ORGANS (UFO) gene, encoding an F-box protein acting as a transcriptional cofactor with LEAFY (LFY) to activate B-class MADS-box gene expression, and the B-class gene PISTILLATA (PI), necessary for male floral organ development, were strongly de-repressed in the arr17 CRISPR mutants. Our data highlight a CK-independent role of the poplar response regulator ARR17 and further emphasize the minimal differences between female and male individuals. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Sex determination and sex chromosome evolution in land plants.

Linnaeus's very first opus, written when he was 22 years old, dealt with the analogy that exists between plants and animals in how they 'propagate their species', and a revised version with a plate depicting the union of male and female Mercurialis annua plants became a foundational text on the sexuality of plants. The question how systems with separate males and females have evolved in sedentary organisms that appear ancestrally bisexual has fascinated biologists ever since. The phenomenon, termed dioecy, has important consequences for plant reproductive success and is of commercial interest since it affects seed quality and fruit production. This theme issue presents a series of articles that synthesize and challenge the current understanding of how plants achieve dioecy. The articles deal with a broad set of taxa, including Coccinia, Ginkgo, Mercurialis, Populus, Rumex and Silene, as well as overarching topics, such as the field's terminology, analogies with animal sex determination systems, evolutionary pathways to dioecy, dosage compensation, and the longevity of the two sexes. In this introduction, we focus on four topics, each addressed by several articles from different angles and with different conclusions. Our highlighting of unclear or controversial issues may help future studies to build on the current understanding and to ask new questions that will expand our knowledge of plant sexual systems. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

The Hda1 histone deacetylase limits divergent non-coding transcription and restricts transcription initiation frequency.

Nucleosome-depleted regions (NDRs) at gene promoters support initiation of RNA polymerase II transcription. Interestingly, transcription often initiates in both directions, resulting in an mRNA and a divergent non-coding (DNC) transcript of unclear purpose. Here, we characterized the genetic architecture and molecular mechanism of DNC transcription in budding yeast. Using high-throughput reverse genetic screens based on quantitative single-cell fluorescence measurements, we identified the Hda1 histone deacetylase complex (Hda1C) as a repressor of DNC transcription. Nascent transcription profiling showed a genome-wide role of Hda1C in repression of DNC transcription. Live-cell imaging of transcription revealed that mutations in the Hda3 subunit increased the frequency of DNC transcription. Hda1C contributed to decreased acetylation of histone H3 in DNC transcription regions, supporting DNC transcription repression by histone deacetylation. Our data support the interpretation that DNC transcription results as a consequence of the NDR-based architecture of eukaryotic promoters, but that it is governed by locus-specific repression to maintain genome fidelity.

Whole-genome sequencing and genome regions of special interest: Lessons from major histocompatibility complex, sex determination, and plant self-incompatibility.

Whole-genome sequencing of non-model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data-handling. Three such regions are the major histocompatibility complex (MHC), sex-determining regions (SDRs) and the plant self-incompatibility locus (S-locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short-read sequences against multiple alternative references improves genotyping comprehensiveness at the S-locus thereby contributing to more accurate assessments of allelic frequencies. Long-read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S-locus haplotypes. These advances are adding value to short-read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches.

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration.

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.

Default Sex and Single Gene Sex Determination in Dioecious Plants.

A well-established hypothesis for the evolution of dioecy involves two genes linked at a sex-determining region (SDR). Recently there has been increased interest in possible single gene sex determination. Work in Populus has finally provided direct experimental evidence for single gene sex determination in plants using CRISPR-Cas9 to knock out a single gene and convert individuals from female to male. In poplar, the feminizing factor popARR17 acts as a "master regulator", analogous to the mammalian masculinizing factor SRY. The production of fully functional males from females by a simple single gene knockout is experimental evidence that an antagonistic male-determining factor does not exist in Populus. Mammals have a "default sex" (female), as do poplar trees (Populus), although the default sex in poplars is male. The occurrence of single gene sex determination with a default sex may be much commoner in plants than hitherto expected, especially when dioecy evolved via monoecy. The master regulator does not even need to be at the SDR (although it may be). In most poplars the feminizing factor popARR17 is not at the SDR, but instead a negative regulator of it. So far there is little information on how high-level regulators are connected to floral phenotype. A model is presented of how sex-determining genes could lead to different floral morphologies via MADS-box floral developmental genes.

A single gene underlies the dynamic evolution of poplar sex determination.

Although hundreds of plant lineages have independently evolved dioecy (that is, separation of the sexes), the underlying genetic basis remains largely elusive1. Here we show that diverse poplar species carry partial duplicates of the ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17) orthologue in the male-specific region of the Y chromosome. These duplicates give rise to small RNAs apparently causing male-specific DNA methylation and silencing of the ARR17 gene. CRISPR-Cas9-induced mutations demonstrate that ARR17 functions as a sex switch, triggering female development when on and male development when off. Despite repeated turnover events, including a transition from the XY system to a ZW system, the sex-specific regulation of ARR17 is conserved across the poplar genus and probably beyond. Our data reveal how a single-gene-based mechanism of dioecy can enable highly dynamic sex-linked regions and contribute to maintaining recombination and integrity of sex chromosomes.

ENO regulates tomato fruit size through the floral meristem development network.

A dramatic evolution of fruit size has accompanied the domestication and improvement of fruit-bearing crop species. In tomato (Solanum lycopersicum), naturally occurring cis-regulatory mutations in the genes of the CLAVATA-WUSCHEL signaling pathway have led to a significant increase in fruit size generating enlarged meristems that lead to flowers with extra organs and bigger fruits. In this work, by combining mapping-by-sequencing and CRISPR/Cas9 genome editing methods, we isolated EXCESSIVE NUMBER OF FLORAL ORGANS (ENO), an AP2/ERF transcription factor which regulates floral meristem activity. Thus, the ENO gene mutation gives rise to plants that yield larger multilocular fruits due to an increased size of the floral meristem. Genetic analyses indicate that eno exhibits synergistic effects with mutations at the LOCULE NUMBER (encoding SlWUS) and FASCIATED (encoding SlCLV3) loci, two central players in the evolution of fruit size in the domestication of cultivated tomatoes. Our findings reveal that an eno mutation causes a substantial expansion of SlWUS expression domains in a flower-specific manner. In vitro binding results show that ENO is able to interact with the GGC-box cis-regulatory element within the SlWUS promoter region, suggesting that ENO directly regulates SlWUS expression domains to maintain floral stem-cell homeostasis. Furthermore, the study of natural allelic variation of the ENO locus proved that a cis-regulatory mutation in the promoter of ENO had been targeted by positive selection during the domestication process, setting up the background for significant increases in fruit locule number and fruit size in modern tomatoes.

Correction to: Long-term study of a subdioecious Populus ×canescens family reveals sex lability of females and reproduction behaviour of cosexual plants.

Table 4 in the original publication reports incomplete genotype names in the column "Cross" and wrong codes in the column "Generation".

Long-term study of a subdioecious Populus ×canescens family reveals sex lability of females and reproduction behaviour of cosexual plants.

KEY MESSAGE: Cosexual Populus ×canescens plants are inconstant females with life course plasticity of sex phenotype and can reproduce by selfing. Populus species are dioecious, but deviations from dioecy are reported in some cases. The objectives of this study were to investigate the phenotypic expression and the inheritance of subdioecy in a Populus ×canescens pedigree. The F1 progeny was monitored for sex during 14 years. Thirty per cent of individuals expressed deviations from dioecy and long-term plasticity of sex. Some plants started flowering as male, then became cosexual, and finally turned female. Two cosexual individuals were self-pollinated and generated a selfed progeny markedly impaired by inbreeding depression, but able to reproduce by outcrossing. Sex segregation of the F1 progeny statistically fitted the expected ratio 1:2:1 (female:male:cosexual). By analysis of DNA markers, the cosexual individuals were genetically clustered with the females. The segregation ratio and the genetic profile indicated that cosexual plants were female with altered sex phenotype. Linkage analysis identified a putative sex-determining region with suppressed recombination on chromosome 19 of the male Populus tremula parent. The male sex trait was linked to the pericentromeric region of the P. tremula chromosome 19, whereas the cosexual trait was linked to chromosome 19 of the female Populus alba parent. A genetic model is proposed to explain inheritance and phenotypic expression of sex.

RNA-seq of eight different poplar clones reveals conserved up-regulation of gene expression in response to insect herbivory.

BACKGROUND: Herbivorous insects can have a profound impact on plant growth performance. In some years, canopy damage in poplar plantations exceeds 50% of the total leaf surface, thereby possibly compromising carbon fixation and biomass yield. To assess the transcriptional response of elite poplar clones to insect feeding and to test whether this response varies between different genotypes, we performed an RNA-sequencing experiment. We deeply sequenced the transcriptomes of eight elite clones belonging to three poplar species (Populus trichocarpa, P. nigra and P. maximowiczii), under Phratora vitellinae feeding and control conditions. This allowed us to precisely quantify transcript levels of about 24,000 expressed genes. RESULTS: Our data reveal a striking overall up-regulation of gene expression under insect attack in all eight poplar clones studied. The up-regulated genes were markedly enriched for the biological process 'regulation of transcription' indicating a highly concerted restructuring of the transcriptome. A search for potential cis-regulatory elements (CREs) that may be involved in this process identified the G-box (CACGTG) as the most significant motif in the promoters of the induced genes. In line with the role of the G-box in jasmonate (JA)-mediated activation of gene expression by MYC2, several genes involved in JA biosynthesis and signaling were up-regulated in our dataset. A co-expression network analysis additionally highlighted WRKY transcription factors. Within the most prominent expression module, WRKYs were strongly overrepresented and occupied several network hubs. Finally, the insect-induced genes comprised several protein families known to be involved in plant defenses, e.g. cytochrome P450s, chitinases and protease inhibitors. CONCLUSIONS: Our data represent a comprehensive characterization of the transcriptional response of selected elite poplar clones to insect herbivory. Our results suggest that the concerted up-regulation of gene expression is controlled by JA signaling and WRKY transcription factors, and activates several defense mechanisms. Our data highlight potential targets of selection and may thus contribute to breeding insect-resistant poplar clones.

Mutations in EID1 and LNK2 caused light-conditional clock deceleration during tomato domestication.

Circadian period and phase of cultivated tomato (Solanum lycopersicum) were changed during domestication, likely adapting the species to its new agricultural environments. Whereas the delayed circadian phase is mainly caused by allelic variation of EID1, the genetic basis of the long circadian period has remained elusive. Here we show that a partial deletion of the clock gene LNK2 is responsible for the period lengthening in cultivated tomatoes. We use resequencing data to phylogenetically classify hundreds of tomato accessions and investigate the evolution of the eid1 and lnk2 mutations along successive domestication steps. We reveal signatures of selection across the genomic region of LNK2 and different patterns of fixation of the mutant alleles. Strikingly, LNK2 and EID1 are both involved in light input to the circadian clock, indicating that domestication specifically targeted this input pathway. In line with this, we show that the clock deceleration in the cultivated tomato is light-dependent and requires the phytochrome B1 photoreceptor. Such conditional variation in circadian rhythms may be key for latitudinal adaptation in a variety of species, including crop plants and livestock.

Variation in the flowering gene SELF PRUNING 5G promotes day-neutrality and early yield in tomato.

Plants evolved so that their flowering is triggered by seasonal changes in day length. However, day-length sensitivity in crops limits their geographical range of cultivation, and thus modification of the photoperiod response was critical for their domestication. Here we show that loss of day-length-sensitive flowering in tomato was driven by the florigen paralog and flowering repressor SELF-PRUNING 5G (SP5G). SP5G expression is induced to high levels during long days in wild species, but not in cultivated tomato because of cis-regulatory variation. CRISPR/Cas9-engineered mutations in SP5G cause rapid flowering and enhance the compact determinate growth habit of field tomatoes, resulting in a quick burst of flower production that translates to an early yield. Our findings suggest that pre-existing variation in SP5G facilitated the expansion of cultivated tomato beyond its origin near the equator in South America, and they provide a compelling demonstration of the power of gene editing to rapidly improve yield traits in crop breeding.

Analysis of Circadian Leaf Movements.

The circadian clock is a molecular timekeeper that controls a wide variety of biological processes. In plants, clock outputs range from the molecular level, with rhythmic gene expression and metabolite content, to physiological processes such as stomatal conductance or leaf movements. Any of these outputs can be used as markers to monitor the state of the circadian clock. In the model plant Arabidopsis thaliana, much of the current knowledge about the clock has been gained from time course experiments profiling expression of endogenous genes or reporter constructs regulated by the circadian clock. Since these methods require labor-intensive sample preparation or transformation, monitoring leaf movements is an interesting alternative, especially in non-model species and for natural variation studies. Technological improvements both in digital photography and image analysis allow cheap and easy monitoring of circadian leaf movements. In this chapter we present a protocol that uses an autonomous point and shoot camera and free software to monitor circadian leaf movements in tomato.

Domestication selected for deceleration of the circadian clock in cultivated tomato.

The circadian clock is a critical regulator of plant physiology and development, controlling key agricultural traits in crop plants. In addition, natural variation in circadian rhythms is important for local adaptation. However, quantitative modulation of circadian rhythms due to artificial selection has not yet been reported. Here we show that the circadian clock of cultivated tomato (Solanum lycopersicum) has slowed during domestication. Allelic variation of the tomato homolog of the Arabidopsis gene EID1 is responsible for a phase delay. Notably, the genomic region harboring EID1 shows signatures of a selective sweep. We find that the EID1 allele in cultivated tomatoes enhances plant performance specifically under long day photoperiods, suggesting that humans selected slower circadian rhythms to adapt the cultivated species to the long summer days it encountered as it was moved away from the equator.