Variation in microbial feature perception in the Rutaceae family with immune receptor conservation in citrus.

Although much is known about the responses of model plants to microbial features, we still lack an understanding of the extent of variation in immune perception across members of a plant family. In this work, we analyzed immune responses in Citrus and wild relatives, surveying 86 Rutaceae genotypes with differing leaf morphologies and disease resistances. We found that responses to microbial features vary both within and between members. Species in 2 subtribes, the Balsamocitrinae and Clauseninae, can recognize flagellin (flg22), cold shock protein (csp22), and chitin, including 1 feature from Candidatus Liberibacter species (csp22CLas), the bacterium associated with Huanglongbing. We investigated differences at the receptor level for the flagellin receptor FLAGELLIN SENSING 2 (FLS2) and the chitin receptor LYSIN MOTIF RECEPTOR KINASE 5 (LYK5) in citrus genotypes. We characterized 2 genetically linked FLS2 homologs from "Frost Lisbon" lemon (Citrus ×limon, responsive) and "Washington navel" orange (Citrus ×aurantium, nonresponsive). Surprisingly, FLS2 homologs from responsive and nonresponsive genotypes were expressed in Citrus and functional when transferred to a heterologous system. "Washington navel" orange weakly responded to chitin, whereas "Tango" mandarin (C. ×aurantium) exhibited a robust response. LYK5 alleles were identical or nearly identical between the 2 genotypes and complemented the Arabidopsis (Arabidopsis thaliana) lyk4/lyk5-2 mutant with respect to chitin perception. Collectively, our data indicate that differences in chitin and flg22 perception in these citrus genotypes are not the results of sequence polymorphisms at the receptor level. These findings shed light on the diversity of perception of microbial features and highlight genotypes capable of recognizing polymorphic pathogen features.

Establishment and inheritance of minichromosomes from Arabidopsis haploid induction.

Minichromosomes are small, sometimes circular, rearranged chromosomes consisting of one centromere and short chromosomal arms formed by treatments that break DNA, including plant transformation. Minichromosomes have the potential to serve as vectors to quickly move valuable genes across a wide range of germplasm, including into adapted crop varieties. To realize this potential, minichromosomes must be reliably generated, easily manipulated, and stably inherited. Here we show a reliable method for minichromosome formation in haploids resulting from CENH3-mediated genome elimination, a process that generates genome instability and karyotypic novelty specifically on one parental genome. First, we identified 2 out of 260 haploids, each containing a single-copy minichromosome originating from centromeric regions of chromosomes 1 and 3, respectively. The chromosome 1 minichromosome we characterized did not pair at meiosis but displayed consistent transmission over nine selfing generations. Next, we demonstrated that CENH3-based haploid induction can produce minichromosomes in a targeted manner. Haploid inducers carrying a selectable pericentromeric marker were used to isolate additional chromosome-specific minichromosomes, which occurred in 3 out of 163 haploids. Our findings document the formation of heritable, rearranged chromosomes, and we provide a method for convenient minichromosome production.

Genetic modulation of Valencia sweet orange field performance by 50 rootstocks under huanglongbing-endemic conditions.

Although the citrus scion cultivar primarily determines the characteristics of the fruit, the rootstock cultivar of the graft combination has a major role in determining the horticultural performance of the tree. The disease huanglongbing (HLB) is particularly devastating to citrus, and the rootstock has been demonstrated to modulate tree tolerance. However, no existing rootstock is entirely suitable in the HLB-endemic environment, and citrus rootstocks are particularly challenging to breed because of a long life cycle and several biological characteristics that interfere with breeding and commercial use. This study with Valencia sweet orange scion documents the multi-season performance of 50 new hybrid rootstocks and commercial standards in one trial that forms the first wave of a new breeding strategy, with the aim of identifying superior rootstocks for commercial use now, and mapping important traits to be used in selection for the next generation of outstanding rootstocks. A large assortment of traits were quantified for all trees in the study, including traits associated with tree size, health, cropping, and fruit quality. Among the quantitative traits compared between rootstock clones, all except one were observed to have significant rootstock influence. Multiple progeny from eight different parental combinations were included in the trial study, and significant differences between parental combinations of the rootstocks were observed for 27 of the 32 traits compared. Pedigree information was integrated with quantitative trait measurements to dissect the genetic components of rootstock-mediated tree performance. Results suggest there is a significant genetic component underlying rootstock-mediated tolerance to HLB and other critical traits, and that integration of pedigree-based genetic information with quantitative phenotypic data from trials should enable marker-based breeding approaches for the rapid selection of next-generation rootstocks with superior combinations of traits that are needed for commercial success. The current generation of new rootstocks included in this trial is a step toward this goal. Based on results from this trial, the new hybrids US-1649, US-1688, US-1709, and US-2338 were considered the four most promising new rootstocks. Release of these rootstocks for commercial use is being considered, pending the evaluation of continuing performance in this trial and the results from other trials.

PICLN modulates alternative splicing and light/temperature responses in plants.

Plants undergo transcriptome reprograming to adapt to daily and seasonal fluctuations in light and temperature conditions. While most efforts have focused on the role of master transcription factors, the importance of splicing factors modulating these processes is now emerging. Efficient pre-mRNA splicing depends on proper spliceosome assembly, which in plants and animals requires the methylosome complex. Ion Chloride nucleotide-sensitive protein (PICLN) is part of the methylosome complex in both humans and Arabidopsis (Arabidopsis thaliana), and we show here that the human PICLN ortholog rescues phenotypes of Arabidopsis picln mutants. Altered photomorphogenic and photoperiodic responses in Arabidopsis picln mutants are associated with changes in pre-mRNA splicing that partially overlap with those in PROTEIN ARGININE METHYL TRANSFERASE5 (prmt5) mutants. Mammalian PICLN also acts in concert with the Survival Motor Neuron (SMN) complex component GEMIN2 to modulate the late steps of UsnRNP assembly, and many alternative splicing events regulated by PICLN but not PRMT5, the main protein of the methylosome, are controlled by Arabidopsis GEMIN2. As with GEMIN2 and SM PROTEIN E1/PORCUPINE (SME1/PCP), low temperature, which increases PICLN expression, aggravates morphological and molecular defects of picln mutants. Taken together, these results establish a key role for PICLN in the regulation of pre-mRNA splicing and in mediating plant adaptation to daily and seasonal fluctuations in environmental conditions.

Gene Body Methylation in Plants: Mechanisms, Functions, and Important Implications for Understanding Evolutionary Processes.

Gene body methylation (gbM) is an epigenetic mark where gene exons are methylated in the CG context only, as opposed to CHG and CHH contexts (where H stands for A, C, or T). CG methylation is transmitted transgenerationally in plants, opening the possibility that gbM may be shaped by adaptation. This presupposes, however, that gbM has a function that affects phenotype, which has been a topic of debate in the literature. Here, we review our current knowledge of gbM in plants. We start by presenting the well-elucidated mechanisms of plant gbM establishment and maintenance. We then review more controversial topics: the evolution of gbM and the potential selective pressures that act on it. Finally, we discuss the potential functions of gbM that may affect organismal phenotypes: gene expression stabilization and upregulation, inhibition of aberrant transcription (reverse and internal), prevention of aberrant intron retention, and protection against TE insertions. To bolster the review of these topics, we include novel analyses to assess the effect of gbM on transcripts. Overall, a growing body of literature finds that gbM correlates with levels and patterns of gene expression. It is not clear, however, if this is a causal relationship. Altogether, functional work suggests that the effects of gbM, if any, must be relatively small, but there is nonetheless evidence that it is shaped by natural selection. We conclude by discussing the potential adaptive character of gbM and its implications for an updated view of the mechanisms of adaptation in plants.

CHH Methylation Islands: A Nonconserved Feature of Grass Genomes That Is Positively Associated with Transposable Elements but Negatively Associated with Gene-Body Methylation.

Methylated CHH (mCHH) islands are peaks of CHH methylation that occur primarily upstream to genes. These regions are actively targeted by the methylation machinery, occur at boundaries between heterochromatin and euchromatin, and tend to be near highly expressed genes. Here we took an evolutionary perspective by studying upstream mCHH islands across a sample of eight grass species. Using a statistical approach to define mCHH islands as regions that differ from genome-wide background CHH methylation levels, we demonstrated that mCHH islands are common and associate with 39% of genes, on average. We hypothesized that islands should be more frequent in genomes of large size, because they have more heterochromatin and hence more need for defined boundaries. We found, however, that smaller genomes tended to have a higher proportion of genes associated with 5' mCHH islands. Consistent with previous work suggesting that islands reflect the silencing of the edge of transposable elements (TEs), genes with nearby TEs were more likely to have mCHH islands. However, the presence of mCHH islands was not a function solely of TEs, both because the underlying sequences of islands were often not homologous to TEs and because genic properties also predicted the presence of 5' mCHH islands. These genic properties included length and gene-body methylation (gbM); in fact, in three of eight species, the absence of gbM was a stronger predictor of a 5' mCHH island than TE proximity. In contrast, gene expression level was a positive but weak predictor of the presence of an island. Finally, we assessed whether mCHH islands were evolutionarily conserved by focusing on a set of 2,720 orthologs across the eight species. They were generally not conserved across evolutionary time. Overall, our data establish additional genic properties that are associated with mCHH islands and suggest that they are not just a consequence of the TE silencing machinery.

Gene body methylation is under selection in Arabidopsis thaliana.

In plants, mammals and insects, some genes are methylated in the CG dinucleotide context, a phenomenon called gene body methylation (gbM). It has been controversial whether this phenomenon has any functional role. Here, we took advantage of the availability of 876 leaf methylomes in Arabidopsis thaliana to characterize the population frequency of methylation at the gene level and to estimate the site-frequency spectrum of allelic states. Using a population genetics model specifically designed for epigenetic data, we found that genes with ancestral gbM are under significant selection to remain methylated. Conversely, ancestrally unmethylated genes were under selection to remain unmethylated. Repeating the analyses at the level of individual cytosines confirmed these results. Estimated selection coefficients were small, on the order of 4 Nes = 1.4, which is similar to the magnitude of selection acting on codon usage. We also estimated that A. thaliana is losing gbM threefold more rapidly than gaining it, which could be due to a recent reduction in the efficacy of selection after a switch to selfing. Finally, we investigated the potential function of gbM through its link with gene expression. Across genes with polymorphic methylation states, the expression of gene body methylated alleles was consistently and significantly higher than unmethylated alleles. Although it is difficult to disentangle genetic from epigenetic effects, our work suggests that gbM has a small but measurable effect on fitness, perhaps due to its association to a phenotype-like gene expression.

Phylogenetic Shifts in Gene Body Methylation Correlate with Gene Expression and Reflect Trait Conservation.

A subset of genes in plant genomes are labeled with DNA methylation specifically at CG residues. These genes, known as gene-body methylated (gbM), have a number of associated characteristics. They tend to have longer sequences, to be enriched for intermediate expression levels, and to be associated with slower rates of molecular evolution. Most importantly, gbM genes tend to maintain their level of DNA methylation between species, suggesting that this trait is under evolutionary constraint. Given the degree of conservation in gbM, we still know surprisingly little about its function in plant genomes or whether gbM is itself a target of selection. To address these questions, we surveyed DNA methylation across eight grass (Poaceae) species that span a gradient of genome sizes. We first established that genome size correlates with genome-wide DNA methylation levels, but less so for genic levels. We then leveraged genomic data to identify a set of 2,982 putative orthologs among the eight species and examined shifts of methylation status for each ortholog in a phylogenetic context. A total of 55% of orthologs exhibited a shift in gbM, but these shifts occurred predominantly on terminal branches, indicating that shifts in gbM are rarely conveyed over time. Finally, we found that the degree of conservation of gbM across species is associated with increased gene length, reduced rates of molecular evolution, and increased gene expression level, but reduced gene expression variation across species. Overall, these observations suggest a basis for evolutionary pressure to maintain gbM status over evolutionary time.

The genome-wide dynamics of purging during selfing in maize.

Self-fertilization (also known as selfing) is an important reproductive strategy in plants and a widely applied tool for plant genetics and plant breeding. Selfing can lead to inbreeding depression by uncovering recessive deleterious variants, unless these variants are purged by selection. Here we investigated the dynamics of purging in a set of eleven maize lines that were selfed for six generations. We show that heterozygous, putatively deleterious single nucleotide polymorphisms are preferentially lost from the genome during selfing. Deleterious single nucleotide polymorphisms were lost more rapidly in regions of high recombination, presumably because recombination increases the efficacy of selection by uncoupling linked variants. Overall, heterozygosity decreased more slowly than expected, by an estimated 35% to 40% per generation instead of the expected 50%, perhaps reflecting pervasive associative overdominance. Finally, three lines exhibited marked decreases in genome size due to the purging of transposable elements. Genome loss was more likely to occur for lineages that began with larger genomes with more transposable elements and chromosomal knobs. These three lines purged an average of 398 Mb from their genomes, an amount equivalent to three Arabidopsis thaliana genomes per lineage, in only a few generations.

Living with Two Genomes: Grafting and Its Implications for Plant Genome-to-Genome Interactions, Phenotypic Variation, and Evolution.

Plant genomes interact when genetically distinct individuals join, or are joined, together. Individuals can fuse in three contexts: artificial grafts, natural grafts, and host-parasite interactions. Artificial grafts have been studied for decades and are important platforms for studying the movement of RNA, DNA, and protein. Yet several mysteries about artificial grafts remain, including the factors that contribute to graft incompatibility, the prevalence of genetic and epigenetic modifications caused by exchanges between graft partners, and the long-term effects of these modifications on phenotype. Host-parasite interactions also lead to the exchange of materials, and RNA exchange actively contributes to an ongoing arms race between parasite virulence and host resistance. Little is known about natural grafts except that they can be frequent and may provide opportunities for evolutionary innovation through genome exchange. In this review, we survey our current understanding about these three mechanisms of contact, the genomic interactions that result, and the potential evolutionary implications.

Transmission ratio distortion is frequent in Arabidopsis thaliana controlled crosses.

The equal probability of transmission of alleles from either parent during sexual reproduction is a central tenet of genetics and evolutionary biology. Yet, there are many cases where this rule is violated. The preferential transmission of alleles or genotypes is termed transmission ratio distortion (TRD). Examples of TRD have been identified in many species, implying that they are universal, but the resolution of species-wide studies of TRD are limited. We have performed a species-wide screen for TRD in over 500 segregating F2 populations of Arabidopsis thaliana using pooled reduced-representation genome sequencing. TRD was evident in up to a quarter of surveyed populations. Most populations exhibited distortion at only one genomic region, with some regions being repeatedly affected in multiple populations. Our results begin to elucidate the species-level architecture of biased transmission of genetic material in A. thaliana, and serve as a springboard for future studies into the biological basis of TRD in this species.

Demography and its effects on genomic variation in crop domestication.

Over two thousand plant species have been modified morphologically through cultivation and human use. Here, we review three aspects of crop domestication that are currently undergoing marked revisions, due to analytical advancements and their application to whole genome resequencing (WGS) data. We begin by discussing the duration and demographic history of domestication. There has been debate as to whether domestication occurred quickly or slowly. The latter is tentatively supported both by fossil data and application of WGS data to sequentially Markovian coalescent methods that infer the history of effective population size. This history suggests the possibility of extended human impacts on domesticated lineages prior to their purposeful cultivation. We also make the point that demographic history matters, because it shapes patterns and levels of extant genetic diversity. We illustrate this point by discussing the evolutionary processes that contribute to the empirical observation that most crops examined to date have more putatively deleterious alleles than their wild relatives. These deleterious alleles may contribute to genetic load within crops and may be fitting targets for crop improvement. Finally, the same demographic factors are likely to shape the spectrum of structural variants (SVs) within crops. SVs are known to underlie many of the phenotypic changes associated with domestication and crop improvement, but we currently lack sufficient knowledge about the mechanisms that create SVs, their rates of origin, their population frequencies and their phenotypic effects.

The Generation of Doubled Haploid Lines for QTL Mapping.

Recombinant inbred lines (RILs) are an essential tool for quantitative trait locus (QTL) mapping in Arabidopsis thaliana. Conventionally, the development of these lines is a time-consuming and tedious process requiring six to eight generations of selfing. Here, we describe an alternative approach: the rapid generation of RILs in A. thaliana via the creation of doubled haploids. In this method, F1 plants are crossed to an engineered haploid inducer to produce haploid plants. The chromosomes of these haploids then spontaneously double, generating immortalized homozygous F2 lines called doubled haploid RILs (DH RILs). Finally, DH RILs are genotyped using low-coverage whole-genome sequencing and are ready to be used for QTL mapping.

The Evolutionary Dynamics of Orthologs That Shift in Gene Body Methylation between Arabidopsis Species.

DNA methylation labels a specific subset of genes in plant genomes. Recent work has shown that this gene-body methylation (gbM) is a conserved feature of orthologs, because highly methylated genes in one species tend to be highly methylated in another. In this study, we examined the exceptions to that rule by identifying genes that differ in gbM status between two plant species-Arabidopsis thaliana and Arabidopsis lyrata. Using Capsella grandiflora as an outgroup, we polarized the loss and gain of gbM for orthologs in the Arabidopsis lineage. Our survey identified a few hundred genes that differed between ingroup species, out of ∼9,000 orthologs. The estimated rate of gbM gain was ∼2 × 10-9 per gene per year for both ingroup taxa and was similar to the loss rate in A. lyrata. In contrast, A. thaliana had a ∼3-fold higher estimated rate of gbM loss per gene, consistent with a recent diminishment of genome size. As in previous studies, we found that body-methylated genes were expressed broadly across tissues and were also longer than other genic sets. Genes that differed in gbM status exhibited higher variance in expression between species than genes that were body-methylated in both species. Moreover, the gain of gbM in one lineage tended to be associated with an increase of expression in that lineage. The genes that varied in gbM status between species varied more significantly in length between species than other sets of genes; we hypothesize that length is a key feature in the transition between body-methylated and nonmethylated genes.

The causes and consequences of DNA methylome variation in plants.

Epigenetic variation - polymorphisms at the level of DNA methylation or histone modifications - modulates chromatin accessibility, which can perturb transcriptional activity and spur phenotypic variation. Determining the origin, frequency spectrum, and consequences of epigenetic variants is key to understanding the role of this variation in generating stable phenotypic variation in plants. Here we review recent literature on DNA methylation variation in both model and crop plant species with a focus on the link between genotype, epigenotype, and transcription. We highlight population epigenomics studies that explore the relationship between epigenetic variants and genetic diversity. Moreover, we provide an overview of relevant studies that together advocate a minor, albeit significant role for epigenetic variation in directing specific transcriptional changes.

Methods for Genotyping-by-Sequencing.

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.

Genetic architecture of nonadditive inheritance in Arabidopsis thaliana hybrids.

The ubiquity of nonparental hybrid phenotypes, such as hybrid vigor and hybrid inferiority, has interested biologists for over a century and is of considerable agricultural importance. Although examples of both phenomena have been subject to intense investigation, no general model for the molecular basis of nonadditive genetic variance has emerged, and prediction of hybrid phenotypes from parental information continues to be a challenge. Here we explore the genetics of hybrid phenotype in 435 Arabidopsis thaliana individuals derived from intercrosses of 30 parents in a half diallel mating scheme. We find that nonadditive genetic effects are a major component of genetic variation in this population and that the genetic basis of hybrid phenotype can be mapped using genome-wide association (GWA) techniques. Significant loci together can explain as much as 20% of phenotypic variation in the surveyed population and include examples that have both classical dominant and overdominant effects. One candidate region inherited dominantly in the half diallel contains the gene for the MADS-box transcription factor AGAMOUS-LIKE 50 (AGL50), which we show directly to alter flowering time in the predicted manner. Our study not only illustrates the promise of GWA approaches to dissect the genetic architecture underpinning hybrid performance but also demonstrates the contribution of classical dominance to genetic variance.

KH domain protein RCF3 is a tissue-biased regulator of the plant miRNA biogenesis cofactor HYL1.

The biogenesis of microRNAs (miRNAs), which regulate mRNA abundance through posttranscriptional silencing, comprises multiple well-orchestrated processing steps. We have identified the Arabidopsis thaliana K homology (KH) domain protein REGULATOR OF CBF GENE EXPRESSION 3 (RCF3) as a cofactor affecting miRNA biogenesis in specific plant tissues. MiRNA and miRNA-target levels were reduced in apex-enriched samples of rcf3 mutants, but not in other tissues. Mechanistically, RCF3 affects miRNA biogenesis through nuclear interactions with the phosphatases C-TERMINAL DOMAIN PHOSPHATASE-LIKE1 and 2 (CPL1 and CPL2). These interactions are essential to regulate the phosphorylation status, and thus the activity, of the double-stranded RNA binding protein and DICER-LIKE1 (DCL1) cofactor HYPONASTIC LEAVES1 (HYL1).

Improving the Annotation of Arabidopsis lyrata Using RNA-Seq Data.

Gene model annotations are important community resources that ensure comparability and reproducibility of analyses and are typically the first step for functional annotation of genomic regions. Without up-to-date genome annotations, genome sequences cannot be used to maximum advantage. It is therefore essential to regularly update gene annotations by integrating the latest information to guarantee that reference annotations can remain a common basis for various types of analyses. Here, we report an improvement of the Arabidopsis lyrata gene annotation using extensive RNA-seq data. This new annotation consists of 31,132 protein coding gene models in addition to 2,089 genes with high similarity to transposable elements. Overall, ~87% of the gene models are corroborated by evidence of expression and 2,235 of these models feature multiple transcripts. Our updated gene annotation corrects hundreds of incorrectly split or merged gene models in the original annotation, and as a result the identification of alternative splicing events and differential isoform usage are vastly improved.

Evolution of DNA methylation patterns in the Brassicaceae is driven by differences in genome organization.

DNA methylation is an ancient molecular modification found in most eukaryotes. In plants, DNA methylation is not only critical for transcriptionally silencing transposons, but can also affect phenotype by altering expression of protein coding genes. The extent of its contribution to phenotypic diversity over evolutionary time is, however, unclear, because of limited stability of epialleles that are not linked to DNA mutations. To dissect the relative contribution of DNA methylation to transposon surveillance and host gene regulation, we leveraged information from three species in the Brassicaceae that vary in genome architecture, Capsella rubella, Arabidopsis lyrata, and Arabidopsis thaliana. We found that the lineage-specific expansion and contraction of transposon and repeat sequences is the main driver of interspecific differences in DNA methylation. The most heavily methylated portions of the genome are thus not conserved at the sequence level. Outside of repeat-associated methylation, there is a surprising degree of conservation in methylation at single nucleotides located in gene bodies. Finally, dynamic DNA methylation is affected more by tissue type than by environmental differences in all species, but these responses are not conserved. The majority of DNA methylation variation between species resides in hypervariable genomic regions, and thus, in the context of macroevolution, is of limited phenotypic consequence.