Erratum: Citrus genomic resources unravel putative genetic determinants of Huanglongbing pathogenicity.

[This corrects the article DOI: 10.1016/j.isci.2023.106024.].

Nitrogen Journey in Plants: From Uptake to Metabolism, Stress Response, and Microbe Interaction.

Plants uptake and assimilate nitrogen from the soil in the form of nitrate, ammonium ions, and available amino acids from organic sources. Plant nitrate and ammonium transporters are responsible for nitrate and ammonium translocation from the soil into the roots. The unique structure of these transporters determines the specificity of each transporter, and structural analyses reveal the mechanisms by which these transporters function. Following absorption, the nitrogen metabolism pathway incorporates the nitrogen into organic compounds via glutamine synthetase and glutamate synthase that convert ammonium ions into glutamine and glutamate. Different isoforms of glutamine synthetase and glutamate synthase exist, enabling plants to fine-tune nitrogen metabolism based on environmental cues. Under stressful conditions, nitric oxide has been found to enhance plant survival under drought stress. Furthermore, the interaction between salinity stress and nitrogen availability in plants has been studied, with nitric oxide identified as a potential mediator of responses to salt stress. Conversely, excessive use of nitrate fertilizers can lead to health and environmental issues. Therefore, alternative strategies, such as establishing nitrogen fixation in plants through diazotrophic microbiota, have been explored to reduce reliance on synthetic fertilizers. Ultimately, genomics can identify new genes related to nitrogen fixation, which could be harnessed to improve plant productivity.

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.

Citrus genomic resources unravel putative genetic determinants of Huanglongbing pathogenicity.

Citrus HLB caused by Candidatus Liberibacter asiaticus is a pathogen-triggered immune disease. Here, we identified putative genetic determinants of HLB pathogenicity by integrating citrus genomic resources to characterize the pan-genome of accessions that differ in their response to HLB. Genome-wide association mapping and analysis of allele-specific expression between susceptible, tolerant, and resistant accessions further refined candidates underlying the response to HLB. We first developed a phased diploid assembly of Citrus sinensis 'Newhall' genome and produced resequencing data for 91 citrus accessions that differ in their response to HLB. These data were combined with previous resequencing data from 356 accessions for genome-wide association mapping of the HLB response. Genes determinants for HLB pathogenicity were associated with host immune response, ROS production, and antioxidants. Overall, this study has provided a significant resource of citrus genomic data and identified candidate genes to be further explored to understand the genetic determinants of HLB pathogenicity.

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.

Gene capture by transposable elements leads to epigenetic conflict in maize.

Transposable elements (TEs) regularly capture fragments of genes. When the host silences these TEs, siRNAs homologous to the captured regions may also target the genes. This epigenetic crosstalk establishes an intragenomic conflict: silencing the TEs has the cost of silencing the genes. If genes are important, however, natural selection may maintain function by moderating the silencing response, which may also advantage the TEs. In this study, we examined this model by focusing on Helitrons, Pack-MULEs, and Sirevirus LTR retrotransposons in the maize genome. We documented 1263 TEs containing exon fragments from 1629 donor genes. Consistent with epigenetic conflict, donor genes mapped more siRNAs and were more methylated than genes with no evidence of capture. However, these patterns differed between syntelog versus translocated donor genes. Syntelogs appeared to maintain function, as measured by gene expression, consistent with moderation of silencing for functionally important genes. Epigenetic marks did not spread beyond their captured regions and 24nt crosstalk siRNAs were linked with CHH methylation. Translocated genes, in contrast, bore the signature of silencing. They were highly methylated and less expressed, but also overrepresented among donor genes and located away from chromosomal arms, which suggests a link between capture and gene movement. Splitting genes into potential functional categories based on evolutionary constraint supported the synteny-based findings. TE families captured genes in different ways, but the evidence for their advantage was generally less obvious; nevertheless, TEs with captured fragments were older, mapped fewer siRNAs, and were slightly less methylated than TEs without captured fragments. Collectively, our results argue that TE capture triggers an intragenomic conflict that may not affect the function of important genes but may lead to the pseudogenization of less-constrained genes.

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 genomic diversification of grapevine clones.

BACKGROUND: Vegetatively propagated clones accumulate somatic mutations. The purpose of this study was to better appreciate clone diversity and involved defining the nature of somatic mutations throughout the genome. Fifteen Zinfandel winegrape clone genomes were sequenced and compared to one another using a highly contiguous genome reference produced from one of the clones, Zinfandel 03. RESULTS: Though most heterozygous variants were shared, somatic mutations accumulated in individual and subsets of clones. Overall, heterozygous mutations were most frequent in intergenic space and more frequent in introns than exons. A significantly larger percentage of CpG, CHG, and CHH sites in repetitive intergenic space experienced transition mutations than in genic and non-repetitive intergenic spaces, likely because of higher levels of methylation in the region and because methylated cytosines often spontaneously deaminate. Of the minority of mutations that occurred in exons, larger proportions of these were putatively deleterious when they occurred in relatively few clones. CONCLUSIONS: These data support three major conclusions. First, repetitive intergenic space is a major driver of clone genome diversification. Second, clones accumulate putatively deleterious mutations. Third, the data suggest selection against deleterious variants in coding regions or some mechanism by which mutations are less frequent in coding than noncoding regions of the genome.

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.