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.

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.

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.

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.

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.

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.

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.

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.

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).

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.

LNK genes integrate light and clock signaling networks at the core of the Arabidopsis oscillator.

Light signaling pathways and the circadian clock interact to help organisms synchronize physiological and developmental processes with periodic environmental cycles. The plant photoreceptors responsible for clock resetting have been characterized, but signaling components that link the photoreceptors to the clock remain to be identified. Here we describe a family of night light-inducible and clock-regulated genes (LNK) that play a key role linking light regulation of gene expression to the control of daily and seasonal rhythms in Arabidopsis thaliana. A genomewide transcriptome analysis revealed that most light-induced genes respond more strongly to light during the subjective day, which is consistent with the diurnal nature of most physiological processes in plants. However, a handful of genes, including the homologous genes LNK1 and LNK2, are more strongly induced by light in the middle of the night, when the clock is most responsive to this signal. Further analysis revealed that the morning phased LNK1 and LNK2 genes control circadian rhythms, photomorphogenic responses, and photoperiodic dependent flowering, most likely by regulating a subset of clock and flowering time genes in the afternoon. LNK1 and LNK2 themselves are directly repressed by members of the TIMING OF CAB1 EXPRESSION/PSEUDO RESPONSE REGULATOR family of core-clock genes in the afternoon and early night. Thus, LNK1 and LNK2 integrate early light signals with temporal information provided by core oscillator components to control the expression of afternoon genes, allowing plants to keep track of seasonal changes in day length.

Rapid creation of Arabidopsis doubled haploid lines for quantitative trait locus mapping.

Quantitative trait loci (QTL) mapping is a powerful tool for investigating the genetic basis of natural variation. QTL can be mapped using a number of different population designs, but recombinant inbred lines (RILs) are among the most effective. Unfortunately, homozygous RIL populations are time consuming to construct, typically requiring at least six generations of selfing starting from a heterozygous F(1). Haploid plants produced from an F(1) combine the two parental genomes and have only one allele at every locus. Converting these sterile haploids into fertile diploids (termed "doubled haploids," DHs) produces immortal homozygous lines in only two steps. Here we describe a unique technique for rapidly creating recombinant doubled haploid populations in Arabidopsis thaliana: centromere-mediated genome elimination. We generated a population of 238 doubled haploid lines that combine two parental genomes and genotyped them by reduced representation Illumina sequencing. The recombination rate and parental allele frequencies in our population are similar to those found in existing RIL sets. We phenotyped this population for traits related to flowering time and for petiole length and successfully mapped QTL controlling each trait. Our work demonstrates that doubled haploid populations offer a rapid, easy alternative to RILs for Arabidopsis genetic analysis.

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.

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.

Comparative mapping in watermelon [Citrullus lanatus (Thunb.) Matsum. et Nakai].

The first single-nucleotide polymorphism (SNP) maps for watermelon [Citrullus lanatus (Thunb.) Matsum. et Nakai] were constructed and compared. Three populations were developed from crosses between two elite cultivars, Klondike Black Seeded × New Hampshire Midget (KBS × NHM), an elite cultivar and wild egusi accession, Strain II × PI 560023 (SII × Egusi) and an elite cultivar and a wild citron accession, ZWRM50 × PI 244019 (ZWRM × Citroides). The SII × Egusi and ZWRM × Citroides F(2) populations consisted of 187 and 182 individuals respectively while the KBS × NHM recombinant inbred line (RIL) population consisted of 164 lines. The length of the genetic maps were 1,438, 1,514 and 1,144 cM with average marker distances of 3.8, 4.2, and 3.4 cM for the KBS × NHM, SII × Egusi and ZWRM × Citroides populations, respectively. Shared markers were used to align the three maps so that the linkage groups (LGs) represented the 11 chromosomes of the species. Marker segregation distortion were observed in all three populations, but was highest (12.7 %) in the ZWRM × Citroides population, where Citroides alleles were favored. The three maps were used to construct a consensus map containing 378 SNP markers with an average distance of 5.1 cM between markers. Phenotypic data was collected for fruit weight (FWT), fruit length (FL), fruit width (FWD), fruit shape index (FSI), rind thickness (RTH) and Brix (BRX) and analyzed for quantitative trait loci (QTL) associated with these traits. A total of 40 QTL were identified in the three populations, including major QTL for fruit size and shape that were stable across genetic backgrounds and environments. The present study reports the first SNP maps for Citrullus and the first map constructed using two elite parents. We also report the first stable QTL associated with fruit size and shape in Citrullus lanatus. These maps, QTL and SNPs should be useful for the watermelon community and represent a significant step towards the potential use of molecular tools in watermelon breeding.

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.

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.