skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 24, 2026

Title: The spectrum of diversity of nucleotide-binding leucine-rich repeat (NLR) genes in citrus and its relatives
Abstract Genomic clusters of immune genes, including those encoding nucleotide-binding leucine-rich repeat (NLR) proteins, are a model for exploring the dynamics of genomic regions in flux. Rapid sequence evolution of immune genes, including NLRs, and variation in their gene content, may enable long-lived plants, which lack adaptive immune systems, to keep pace with the fast evolution of pathogens. To explore the patterns and processes shaping the evolution of NLR gene content in a genus of long-lived tree species, we unified the annotation of NLR genes across 11 accessions (or 15 haplotypes) from the genusCitrusand its relatives, including three new diploid genome assemblies. A majority of NLRs were arranged in genomic clusters composed of paralogous genes, typically from a single gene family. Even larger clusters, with 10 or more NLRs, were limited to genes derived from one or few gene families. These patterns suggested that genomic clustering of NLRs arose through local expansion of phylogenetically related NLRs, but the mechanistic processes driving these patterns are not clear. Local gene duplication can be mediated by multiple processes, including transposon-mediated gene capture and subsequent proliferation, and non-allelic repair of double stranded breaks, including unequal recombination. Examples of retrotransposon-mediated duplication of NLRs were identified, but these were not sufficient to explain massive regional expansions. Signatures of unequal recombination are challenging to identify. Focusing on recent lineage-specific sequence duplications, at least one case of unequal recombination was identified, supporting a role for unequal recombination in shaping genomic variation in these regions.  more » « less
Award ID(s):
2215705
PAR ID:
10656454
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
bioRxiv
Date Published:
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Plant innate immunity relies on nucleotide binding leucine-rich repeat receptors (NLRs) that recognize pathogen-derived molecules and activate downstream signaling pathways. We analyzed the variation in NLR gene copy number and identified plants with a low number of NLR genes relative to sister species. We specifically focused on four plants from two distinct lineages, one monocot lineage (Alismatales) and one eudicot lineage (Lentibulariaceae). In these lineages, the loss of NLR genes coincides with loss of the well-known downstream immune signaling complex ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1)/PHYTOALEXIN DEFICIENT 4 (PAD4). We expanded our analysis across whole proteomes and found that other characterized immune genes were absent only in Lentibulariaceae and Alismatales. Additionally, we identified genes of unknown function that were convergently lost together with EDS1/PAD4 in five plant species. Gene expression analyses in Arabidopsis (Arabidopsis thaliana) and Oryza sativa revealed that several homologs of the candidates are differentially expressed during pathogen infection, drought, and abscisic acid treatment. Our analysis provides evolutionary evidence for the rewiring of plant immunity in some plant lineages, as well as the coevolution of the EDS1/PAD4 pathway and drought responses. 
    more » « less
  2. Genomic data can provide valuable insights into the evolutionary history of rapidly diversifying groups and the genetic basis of phenotypic differences among lineages. We used whole-genome sequencing of the warbler genus Myioborus to investigate dynamics of its recent diversification in Neotropical mountains. We found that mitochondrial and UCE phylogenies are mostly, but not fully, concordant, and we found phylogenetic support for a pattern of north-to-south and low-to-high elevation colonization in the genus. Within the ornatus-melanocephalus complex, which showed topological incongruence between our phylogenies, we found that genetic structure generally coincides with geographic variation in plumage, although three subspecies with striking plumage differences exhibit low mitochondrial divergence. The hybridizing taxa M. o. chrysops and M. m. bairdi show very shallow genomic differentiation, with marked peaks of divergence. Most of these are shared with other parulid warbler pairs, pointing to broad genomic features, like recombination rate, as the processes shaping these regions. However, other highly differentiated regions were unique to Myioborus, including one containing the gene CCDC91, which is associated with melanin-based plumage differences in several other birds. Lastly, we found higher levels of differentiation on the Z chromosome relative to autosomes, including two putative chromosomal inversions. Together, these results highlight the interplay of deep ancestral divergence, recent hybridization, and shared genomic architecture in shaping the evolution of phenotypic and genomic diversity within Myioborus. 
    more » « less
  3. Genes involved in disease resistance are some of the fastest evolving and most diverse components of genomes. Large numbers of nucleotide-binding, leucine-rich repeat (NLR) genes are found in plant genomes and are required for disease resistance. However, NLRs can trigger autoimmunity, disrupt beneficial microbiota or reduce fitness. It is therefore crucial to understand how NLRs are controlled. Here, we show that the RNA-binding protein FPA mediates widespread premature cleavage and polyadenylation of NLR transcripts, thereby controlling their functional expression and impacting immunity. Using long-read Nanopore direct RNA sequencing, we resolved the complexity of NLR transcript processing and gene annotation. Our results uncover a co-transcriptional layer of NLR control with implications for understanding the regulatory and evolutionary dynamics of NLRs in the immune responses of plants. 
    more » « less
  4. The genes that encode the α- and β-chain subunits of vertebrate hemoglobin have served as a model system for elucidating general principles of gene family evolution, but little is known about patterns of evolution in amniotes other than mammals and birds. Here, we report a comparative genomic analysis of the α- and β-globin gene clusters in sauropsids (archosaurs and nonavian reptiles). The objectives were to characterize changes in the size and membership composition of the α- and β-globin gene families within and among the major sauropsid lineages, to reconstruct the evolutionary history of the sauropsid α- and β-globin genes, to resolve orthologous relationships, and to reconstruct evolutionary changes in the developmental regulation of gene expression. Our comparisons revealed contrasting patterns of evolution in the unlinked α- and β-globin gene clusters. In the α-globin gene cluster, which has remained in the ancestral chromosomal location, evolutionary changes in gene content are attributable to the differential retention of paralogous gene copies that were present in the common ancestor of tetrapods. In the β-globin gene cluster, which was translocated to a new chromosomal location, evolutionary changes in gene content are attributable to differential gene gains (via lineage-specific duplication events) and gene losses (via lineage-specific deletions and inactivations). Consequently, all major groups of amniotes possess unique repertoires of embryonic and postnatally expressed β-type globin genes that diversified independently in each lineage. These independently derived β-type globins descend from a pair of tandemly linked paralogs in the most recent common ancestor of sauropsids. 
    more » « less
  5. Rokas, A (Ed.)
    Abstract Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation. 
    more » « less