Gene functional-state covariance reveals modular structure in the Arabidopsis thaliana dispensable genome
- 1Royal Botanic Gardens, Kew
Abstract
Classical multilocus theory predicts that selection acting on combinations of loci can maintain non-random associations among allelic states, causing selection to act on multilocus genotypes rather than individual alleles. Across eukaryotic species, pangenome studies have documented thousands of genes that vary in functional status between individuals, yet this variation is analyzed as if each gene evolves independently. Here we test for non-random association within the dispensable pangenome across 1,012 natural genotypes of Arabidopsis thaliana to better understand the structure and evolution of eukaryotic pangenomes. Genome-wide pairwise association testing of loss-of-function alleles across genotypes identified 150,033 significantly covarying gene pairs among 5,679 unique genes after controlling for kinship and population structure. Most associations were between genes on different chromosomes, and covarying gene pairs showed greater co-expression and functional similarity than expected by chance. Network analysis revealed a non-random higher-order architecture comprising stable functional communities with distinct biological functions and structured connectivity among communities. The enriched functions and connectivity of communities revealed a tight coupling of modules involving plant immunity, environmental response, and specialized metabolism, highlighting how these systems integrate and evolve in response to dynamic ecological pressures. Taken together, these results show that the dispensable genome has a complex higher-order structure driven by functionally interacting genes, and that variation in gene content provides a population-scale resource for investigating the organization and evolution of biological systems.
Related
- Overexpression of SiLEA_2-6 from Saussurea involucrata Enhances Salt and Osmotic Tolerance in Arabidopsis thaliana
- Hypomagnetic field inhibits seed germination via ROS signaling and reprogramming transcriptome in Arabidopsis thaliana
- TOC159 receptors use a galactolipid-specific bipartite signal for chloroplast targeting in Arabidopsis thaliana
- A miniaturized Arabidopsis thaliana - Botrytis cinerea seedling pathosystem and the GreenDayz workflow for longitudinal image-based phenotyping
- Loss of ALN Function in Arabidopsis thaliana Induces Allantoin Accumulation and Enhances Antioxidant Enzyme Responses to Copper Toxicity
BibTeX
@misc{bird2026gene,
title = {{Gene functional-state covariance reveals modular structure in the Arabidopsis thaliana dispensable genome}},
author = {K. A. Bird and K. Zhao and D. Kliebenstein and J. G. Monroe},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.08.757691},
url = {https://www.biorxiv.org/content/10.64898/2026.10.08.757691v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Gene functional-state covariance reveals modular structure in the Arabidopsis thaliana dispensable genome
AU - K. A. Bird
AU - K. Zhao
AU - D. Kliebenstein
AU - J. G. Monroe
PY - 2026
DA - 2026/10/08
PB - bioRxiv
DO - 10.64898/2026.10.08.757691
UR - https://www.biorxiv.org/content/10.64898/2026.10.08.757691v1
AB - Classical multilocus theory predicts that selection acting on combinations of loci can maintain non-random associations among allelic states, causing selection to act on multilocus genotypes rather than individual alleles. Across eukaryotic species, pangenome studies have documented thousands of genes that vary in functional status between individuals, yet this variation is analyzed as if each gene evolves independently. Here we test for non-random association within the dispensable pangenome across 1,012 natural genotypes of Arabidopsis thaliana to better understand the structure and evolution of eukaryotic pangenomes. Genome-wide pairwise association testing of loss-of-function alleles across genotypes identified 150,033 significantly covarying gene pairs among 5,679 unique genes after controlling for kinship and population structure. Most associations were between genes on different chromosomes, and covarying gene pairs showed greater co-expression and functional similarity than expected by chance. Network analysis revealed a non-random higher-order architecture comprising stable functional communities with distinct biological functions and structured connectivity among communities. The enriched functions and connectivity of communities revealed a tight coupling of modules involving plant immunity, environmental response, and specialized metabolism, highlighting how these systems integrate and evolve in response to dynamic ecological pressures. Taken together, these results show that the dispensable genome has a complex higher-order structure driven by functionally interacting genes, and that variation in gene content provides a population-scale resource for investigating the organization and evolution of biological systems.
ER -