Comparative analyses suggest a link between mRNA splicing, stability, and RNA covalent modifications in flowering plants.

BACKGROUND: In recent years, covalent modifications on RNA nucleotides have emerged as pivotal moieties influencing the structure, function, and regulatory processes of RNA Polymerase II transcripts such as mRNAs and lncRNAs. However, our understanding of their biological roles and whether these roles are conserved across eukaryotes remains limited. RESULTS: In this study, we leveraged standard polyadenylation-enriched RNA-sequencing data to identify and characterize RNA modifications that introduce base-pairing errors into cDNA reads. Our investigation incorporated data from three Poaceae (Zea mays, Sorghum bicolor, and Setaria italica), as well as publicly available data from a range of stress and genetic contexts in Sorghum and Arabidopsis thaliana. We uncovered a strong enrichment of RNA covalent modifications (RCMs) deposited on a conserved core set of nuclear mRNAs involved in photosynthesis and translation across these species. However, the cohort of modified transcripts changed based on environmental context and developmental program, a pattern that was also conserved across flowering plants. We determined that RCMs can partly explain accession-level differences in drought tolerance in Sorghum, with stress-associated genes receiving a higher level of RCMs in a drought tolerant accession. To address function, we determined that RCMs are significantly enriched near exon junctions within coding regions, suggesting an association with splicing. Intriguingly, we found that these base-pair disrupting RCMs are associated with stable mRNAs, are highly correlated with protein abundance, and thus likely associated with facilitating translation. CONCLUSIONS: Our data point to a conserved role for RCMs in mRNA stability and translation across the flowering plant lineage.

Regulation of a single inositol 1-phosphate synthase homeologue by HSFA6B contributes to fibre yield maintenance under drought conditions in upland cotton.

Drought stress substantially impacts crop physiology resulting in alteration of growth and productivity. Understanding the genetic and molecular crosstalk between stress responses and agronomically important traits such as fibre yield is particularly complicated in the allopolyploid species, upland cotton (Gossypium hirsutum), due to reduced sequence variability between A and D subgenomes. To better understand how drought stress impacts yield, the transcriptomes of 22 genetically and phenotypically diverse upland cotton accessions grown under well-watered and water-limited conditions in the Arizona low desert were sequenced. Gene co-expression analyses were performed, uncovering a group of stress response genes, in particular transcription factors GhDREB2A-A and GhHSFA6B-D, associated with improved yield under water-limited conditions in an ABA-independent manner. DNA affinity purification sequencing (DAP-seq), as well as public cistrome data from Arabidopsis, were used to identify targets of these two TFs. Among these targets were two lint yield-associated genes previously identified through genome-wide association studies (GWAS)-based approaches, GhABP-D and GhIPS1-A. Biochemical and phylogenetic approaches were used to determine that GhIPS1-A is positively regulated by GhHSFA6B-D, and that this regulatory mechanism is specific to Gossypium spp. containing the A (old world) genome. Finally, an SNP was identified within the GhHSFA6B-D binding site in GhIPS1-A that is positively associated with yield under water-limiting conditions. These data lay out a regulatory connection between abiotic stress and fibre yield in cotton that appears conserved in other systems such as Arabidopsis.

Development of a mobile, high-throughput, and low-cost image-based plant growth phenotyping system.

Nondestructive plant phenotyping forms a key technique for unraveling molecular processes underlying plant development and response to the environment. While the emergence of high-throughput phenotyping facilities can further our understanding of plant development and stress responses, their high costs greatly hinder scientific progress. To democratize high-throughput plant phenotyping, we developed sets of low-cost image- and weight-based devices to monitor plant shoot growth and evapotranspiration. We paired these devices to a suite of computational pipelines for integrated and straightforward data analysis. The developed tools were validated for their suitability for large genetic screens by evaluating a cowpea (Vigna unguiculata) diversity panel for responses to drought stress. The observed natural variation was used as an input for a genome-wide association study, from which we identified nine genetic loci that might contribute to cowpea drought resilience during early vegetative development. The homologs of the candidate genes were identified in Arabidopsis (Arabidopsis thaliana) and subsequently evaluated for their involvement in drought stress by using available T-DNA insertion mutant lines. These results demonstrate the varied applicability of this low-cost phenotyping system. In the future, we foresee these setups facilitating the identification of genetic components of growth, plant architecture, and stress tolerance across a wide variety of plant species.

Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants.

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

Composition and function of stress granules and P-bodies in plants.

Stress Granules (SGs) and Processing-bodies (P-bodies) are biomolecular condensates formed in the cell with the highly conserved purpose of maintaining balance between storage, translation, and degradation of mRNA. This balance is particularly important when cells are exposed to different environmental conditions and adjustments have to be made in order for plants to respond to and tolerate stressful conditions. While P-bodies are constitutively present in the cell, SG formation is a stress-induced event. Typically thought of as protein-RNA aggregates, SGs and P-bodies are formed by a process called liquid-liquid phase separation (LLPS), and both their function and composition are very dynamic. Both foci are known to contain proteins involved in translation, protein folding, and ATPase activity, alluding to their roles in regulating mRNA and protein expression levels. From an RNA perspective, SGs and P-bodies primarily consist of mRNAs, though long non-coding RNAs (lncRNAs) have also been observed, and more focus is now being placed on the specific RNAs associated with these aggregates. Recently, metabolites such as nucleotides and amino acids have been reported in purified plant SGs with implications for the energetic dynamics of these condensates. Thus, even though the field of plant SGs and P-bodies is relatively nascent, significant progress has been made in understanding their composition and biological role in stress responses. In this review, we discuss the most recent discoveries centered around SG and P-body function and composition in plants.

Multiple-biomarkers show the importance of blue carbon to commercially important fishery species.

Coastal blue carbon ecosystems (BCE) support nearshore food webs and provide habitat for many commercially important fish and crustacean species. However, the complex links between catchment vegetation and the carbon food-base of estuarine systems are difficult to disern. We employed a multi-biomarker approach (stable isotope ratios - δ13C and δ15N, fatty acid trophic markers - FATMs and metabolomics - central carbon metabolism metabolites) to test links between estuarine vegetation and the food sources available to commercially important crabs and fish occurring within the river systems of the near-pristine eastern coastline of the Gulf of Carpentaria, Australia. Stable isotope analysis confirmed the dietary importance of fringing macrophytes to consumer diet, but showed that this is modulated by their dominance along the riverbank. FATMs indicative of specific food sources further confirmed the differences among upper intertidal macrophytes (driven by concentrations of 16: 1ω7, 18:1ω9, 18:2ω6, 18:3ω3 & 22.0) and seagrass (driven by 18:2ω6, 18:3ω3). These dietary patterns were also reflected in the concentration of central carbon metabolism metabolites. Overall, our study demonstrates the congruence of different biomarker approaches to resolve biochemical links between blue carbon ecosystems and important nekton species, and provides fresh insights into the pristine tropical estuaries of northern Australia.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs.

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Detection of marine oil-like features in Sentinel-1 SAR images by supplementary use of deep learning and empirical methods: Performance assessment for the Great Barrier Reef marine park.

Continuous monitoring of oil discharges in coastal and open ocean waters using Earth Observation (EO) has undeniably contributed to diminishing their occurrence wherever a detection system was in place, such as in Europe (EMSA's CleanSeaNet) or in the United States (NOAA's OR&R). This study describes the development and testing of a semi-automated oil slick detection system tailored to the Great Barrier Reef (GBR) marine park solely based on EO data as no such service was routinely available in Australia until recently. In this study, a large, curated, historical global dataset of SAR imagery acquired by Sentinel-1 SAR, now publicly available, is used to assess classification techniques, namely an empirical approach and a deep learning model, to discriminate between oil-like features and look-alikes in the scenes acquired over the marine park. An evaluation of this detection system on 10 Sentinel-1 SAR images of the GBR using two performance metrics - the detection accuracy and the false-positive rate (FPR) - shows that the classifiers perform best when combined (accuracy >98 %; FPR 0.01) rather than when used separately. This study demonstrates the benefit of sequentially combining classifiers to improve the detection and monitoring of unreported oil discharge events in SAR imagery. The workflow has also been tested outside the GBR, demonstrating its robustness when applied to other regions such as Australia's Northwest Shelf, Southeast Asia and the Pacific.

Telomerase RNA in Hymenoptera (Insecta) switched to plant/ciliate-like biogenesis.

In contrast to the catalytic subunit of telomerase, its RNA subunit (TR) is highly divergent in size, sequence and biogenesis pathways across eukaryotes. Current views on TR evolution assume a common origin of TRs transcribed with RNA polymerase II in Opisthokonta (the supergroup including Animalia and Fungi) and Trypanosomida on one hand, and TRs transcribed with RNA polymerase III under the control of type 3 promoter, found in TSAR and Archaeplastida supergroups (including e.g. ciliates and Viridiplantae taxa, respectively). Here, we focus on unknown TRs in one of the largest Animalia order - Hymenoptera (Arthropoda) with more than 300 available representative genomes. Using a combination of bioinformatic and experimental approaches, we identify their TRs. In contrast to the presumed type of TRs (H/ACA box snoRNAs transcribed with RNA Polymerase II) corresponding to their phylogenetic position, we find here short TRs of the snRNA type, likely transcribed with RNA polymerase III under the control of the type 3 promoter. The newly described insect TRs thus question the hitherto assumed monophyletic origin of TRs across Animalia and point to an evolutionary switch in TR type and biogenesis that was associated with the divergence of Arthropods.

Evolutionary analysis of the LORELEI gene family in plants reveals regulatory subfunctionalization.

A signaling complex comprising members of the LORELEI (LRE)-LIKE GPI-anchored protein (LLG) and Catharanthus roseus RECEPTOR-LIKE KINASE 1-LIKE (CrRLK1L) families perceive RAPID ALKALINIZATION FACTOR (RALF) peptides and regulate growth, reproduction, immunity, and stress responses in Arabidopsis (Arabidopsis thaliana). Genes encoding these proteins are members of multigene families in most angiosperms and could generate thousands of signaling complex variants. However, the links between expansion of these gene families and the functional diversification of this critical signaling complex as well as the evolutionary factors underlying the maintenance of gene duplicates remain unknown. Here, we investigated LLG gene family evolution by sampling land plant genomes and explored the function and expression of angiosperm LLGs. We found that LLG diversity within major land plant lineages is primarily due to lineage-specific duplication events, and that these duplications occurred both early in the history of these lineages and more recently. Our complementation and expression analyses showed that expression divergence (i.e. regulatory subfunctionalization), rather than functional divergence, explains the retention of LLG paralogs. Interestingly, all but one monocot and all eudicot species examined had an LLG copy with preferential expression in male reproductive tissues, while the other duplicate copies showed highest levels of expression in female or vegetative tissues. The single LLG copy in Amborella trichopoda is expressed vastly higher in male compared to in female reproductive or vegetative tissues. We propose that expression divergence plays an important role in retention of LLG duplicates in angiosperms.

High-Throughput Evolutionary Comparative Analysis of Long Intergenic Noncoding RNAs in Multiple Organisms.

Comparative genomic and transcriptomic analyses can help prioritize and facilitate the functional analysis of long noncoding RNAs (lncRNAs). Evolinc-II is a bioinformatic pipeline that automates comparative analyses, searching for sequence and structural conservation for thousands of lncRNAs at once. In addition, Evolinc-II takes a phylogenetic approach to infer key evolutionary events that may have occurred during the emergence of each query lncRNA. Here, we describe how to use command line or GUI (CyVerse's Discovery Environment) versions of Evolinc-II to identify lncRNA homologs and prioritize them for functional analysis.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae.

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

A Conserved Long Intergenic Non-coding RNA Containing snoRNA Sequences, lncCOBRA1, Affects Arabidopsis Germination and Development.

Long non-coding RNAs (lncRNAs) are an increasingly studied group of non-protein coding transcripts with a wide variety of molecular functions gaining attention for their roles in numerous biological processes. Nearly 6,000 lncRNAs have been identified in Arabidopsis thaliana but many have yet to be studied. Here, we examine a class of previously uncharacterized lncRNAs termed CONSERVED IN BRASSICA RAPA (lncCOBRA) transcripts that were previously identified for their high level of sequence conservation in the related crop species Brassica rapa, their nuclear-localization and protein-bound nature. In particular, we focus on lncCOBRA1 and demonstrate that its abundance is highly tissue and developmental specific, with particularly high levels early in germination. lncCOBRA1 contains two snoRNAs domains within it, making it the first sno-lincRNA example in a non-mammalian system. However, we find that it is processed differently than its mammalian counterparts. We further show that plants lacking lncCOBRA1 display patterns of delayed germination and are overall smaller than wild-type plants. Lastly, we identify the proteins that interact with lncCOBRA1 and propose a novel mechanism of lincRNA action in which it may act as a scaffold with the RACK1A protein to regulate germination and development, possibly through a role in ribosome biogenesis.

2',3'-cAMP treatment mimics the stress molecular response in Arabidopsis thaliana.

The role of the RNA degradation product 2',3'-cyclic adenosine monophosphate (2',3'-cAMP) is poorly understood. Recent studies have identified 2',3'-cAMP in plant material and determined its role in stress signaling. The level of 2',3'-cAMP increases upon wounding, in the dark, and under heat, and 2',3'-cAMP binding to an RNA-binding protein, Rbp47b, promotes stress granule (SG) assembly. To gain further mechanistic insights into the function of 2',3'-cAMP, we used a multi-omics approach by combining transcriptomics, metabolomics, and proteomics to dissect the response of Arabidopsis (Arabidopsis thaliana) to 2',3'-cAMP treatment. We demonstrated that 2',3'-cAMP is metabolized into adenosine, suggesting that the well-known cyclic nucleotide-adenosine pathway of human cells might also exist in plants. Transcriptomics analysis revealed only minor overlap between 2',3'-cAMP- and adenosine-treated plants, suggesting that these molecules act through independent mechanisms. Treatment with 2',3'-cAMP changed the levels of hundreds of transcripts, proteins, and metabolites, many previously associated with plant stress responses, including protein and RNA degradation products, glucosinolates, chaperones, and SG components. Finally, we demonstrated that 2',3'-cAMP treatment influences the movement of processing bodies, confirming the role of 2',3'-cAMP in the formation and motility of membraneless organelles.

Omics-based ecosurveillance uncovers the influence of estuarine macrophytes on sediment microbial function and metabolic redundancy in a tropical ecosystem.

Vertical zonation within estuarine ecosystems can strongly influence microbial diversity and function by regulating competition, predation, and environmental stability. The degree to which microbial communities exhibit horizontal patterns through an estuary has received comparatively less attention. Here, we take a multi-omics ecosurveillance approach to study environmental gradients created by the transition between dominant vegetation types along a near pristine tropical river system (Wenlock River, Far North Queensland, Australia). The study sites included intertidal mudflats fringed by saltmarsh, mangrove or mixed soft substrata habitats. Collected sediments were analyzed for eukaryotes and prokaryotes using small sub-unit (SSU) rRNA gene amplicons to profile the relative taxonomic composition. Central carbon metabolism metabolites and other associated organic polar metabolites were analyzed using established metabolomics-based approaches, coupled with total heavy metals analysis. Eukaryotic taxonomic information was found to be more informative of habitat type. Bacterial taxonomy and community composition also showed habitat-specificity, with phyla Proteobacteria and Cyanobacteria strongly linked to mangroves and saltmarshes, respectively. In contrast, metabolite profiling was critical for understanding the biochemical pathways and expressed functional outputs in these systems that were tied to predicted microbial gene function (16S rRNA). A high degree of metabolic redundancy was observed in the bacterial communities, with the metabolomics data suggesting varying degrees of metabolic criticality based on habitat type. The predicted functions of the bacterial taxa combined with annotated metabolites accounted for the conservative perspective of microbial community redundancy against the putative metabolic pathway impacts in the metabolomics data. Coupling these data demonstrates that habitat-mediated estuarine gradients drive patterns of community diversity and metabolic function and highlights the real redundancy potential of habitat microbiomes. This information is useful as a point of comparison for these sensitive ecosystems and provides a framework for identifying potentially vulnerable or at-risk systems before they are significantly degraded.

Evolution of plant telomerase RNAs: farther to the past, deeper to the roots.

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats.

Chloroplast quality control pathways are dependent on plastid DNA synthesis and nucleotides provided by cytidine triphosphate synthase two.

Reactive oxygen species (ROS) produced in chloroplasts cause oxidative damage, but also signal to initiate chloroplast quality control pathways, cell death, and gene expression. The Arabidopsis thaliana plastid ferrochelatase two (fc2) mutant produces the ROS singlet oxygen in chloroplasts that activates such signaling pathways, but the mechanisms are largely unknown. Here we characterize one fc2 suppressor mutation and map it to CYTIDINE TRIPHOSPHATE SYNTHASE TWO (CTPS2), which encodes one of five enzymes in Arabidopsis necessary for de novo cytoplasmic CTP (and dCTP) synthesis. The ctps2 mutation reduces chloroplast transcripts and DNA content without similarly affecting mitochondria. Chloroplast nucleic acid content and singlet oxygen signaling are restored by exogenous feeding of the dCTP precursor deoxycytidine, suggesting ctps2 blocks signaling by limiting nucleotides for chloroplast genome maintenance. An investigation of CTPS orthologs in Brassicaceae showed CTPS2 is a member of an ancient lineage distinct from CTPS3. Complementation studies confirmed this analysis; CTPS3 was unable to compensate for CTPS2 function in providing nucleotides for chloroplast DNA and signaling. Our studies link cytoplasmic nucleotide metabolism with chloroplast quality control pathways. Such a connection is achieved by a conserved clade of CTPS enzymes that provide nucleotides for chloroplast function, thereby allowing stress signaling to occur.

Impact of catchment-derived nutrients and sediments on marine water quality on the Great Barrier Reef: An application of the eReefs marine modelling system.

Water quality of the Great Barrier Reef (GBR) is determined by a range of natural and anthropogenic drivers that are resolved in the eReefs coupled hydrodynamic - biogeochemical marine model forced by a process-based catchment model, GBR Dynamic SedNet. Model simulations presented here quantify the impact of anthropogenic catchment loads of sediments and nutrients on a range of marine water quality variables. Simulations of 2011-2018 show that reduction of anthropogenic catchment loads results in improved water quality, especially within river plumes. Within the 16 resolved river plumes, anthropogenic loads increased chlorophyll concentration by 0.10 (0.02-0.25) mg Chl m-3. Reductions of anthropogenic loads following proposed Reef 2050 Water Quality Improvement Plan targets reduced chlorophyll concentration in the plumes by 0.04 (0.01-0.10) mg Chl m-3. Our simulations demonstrate the impact of anthropogenic loads on GBR water quality and quantify the benefits of improved catchment management.

Identification and Characterization of the Heat-Induced Plastidial Stress Granules Reveal New Insight Into Arabidopsis Stress Response.

Plants exhibit different physiological and molecular responses to adverse changes in their environment. One such molecular response is the sequestration of proteins, RNAs, and metabolites into cytoplasmic bodies called stress granules (cSGs). Here we report that, in addition to cSGs, heat stress also induces the formation of SG-like foci (cGs) in the chloroplasts of the model plant Arabidopsis thaliana. Similarly to the cSGs, (i) cpSG assemble rapidly in response to stress and disappear when the stress ceases, (ii) cpSG formation is inhibited by treatment with a translation inhibitor (lincomycin), and (iii) cpSG are composed of a stable core and a fluid outer shell. A previously published protocol for cSG extraction was successfully adapted to isolate cpSG, followed by protein, metabolite, and RNA analysis. Analogously to the cSGs, cpSG sequester proteins essential for SG formation, dynamics, and function, also including RNA-binding proteins with prion-like domain, ATPases and chaperones, and the amino acids proline and glutamic acid. However, the most intriguing observation relates to the cpSG localization of proteins, such as a complete magnesium chelatase complex, which is involved in photosynthetic acclimation to stress. These data suggest that cpSG have a role in plant stress tolerance.

Biased Gene Retention in the Face of Introgression Obscures Species Relationships.

Phylogenomic analyses are recovering previously hidden histories of hybridization, revealing the genomic consequences of these events on the architecture of extant genomes. We applied phylogenomic techniques and several complementary statistical tests to show that introgressive hybridization appears to have occurred between close relatives of Arabidopsis, resulting in cytonuclear discordance and impacting our understanding of species relationships in the group. The composition of introgressed and retained genes indicates that selection against incompatible cytonuclear and nuclear-nuclear interactions likely acted during introgression, whereas linkage also contributed to genome composition through the retention of ancient haplotype blocks. We also applied divergence-based tests to determine the species branching order and distinguish donor from recipient lineages. Surprisingly, these analyses suggest that cytonuclear discordance arose via extensive nuclear, rather than cytoplasmic, introgression. If true, this would mean that most of the nuclear genome was displaced during introgression whereas only a small proportion of native alleles were retained.