Abstract Salt tolerance is a complex trait with much of the underlying genetic variation and integrated response strategies yet to be discovered from stress adapted plants. Schrenkiella parvula and Eutrema salsugineum are salt-tolerant extremophytes related to Arabidopsis thaliana in Brassicaceae. We investigated their response strategies contrasted against the salt-sensitive model, A. thaliana to cope with salt stresses via transcriptomic, metabolomic, and ionomic adjustments. The extremophytes exemplified divergent routes to achieve nutrient balance, build osmotolerance, boost antioxidant capacity, and extend transcriptomic support for modified ion transport and stress signaling. Those led to similar molecular phenotypes adapted to salt stress in the extremophytes, absent in A. thaliana . The predominant transcriptomic signals in all three species were associated with salt stress. However, root architecture modulation mediated by negative regulation of auxin and ABA signaling supported minimally-affected root growth unique to each extremophyte during salt treatments. Overall, E. salsugineum exhibited pre-adapted responses at the metabolome level, whereas S. parvula showed dynamic metabolomic responses coupled to a pre-adapted transcriptome to survive salt stress. Our work shows that the two extremophytes share common salt tolerance features, but differ substantially in pathways leading to the convergent, adaptive traits.

Abstract Boron toxicity is a worldwide problem for crop production, yet we have only a limited understanding of the genetic responses and adaptive mechanisms to this environmental stress in plants. Here we identified responses to excess boron in boron stress-sensitive Arabidopsis thaliana and its boron stress-tolerant extremophyte relative Schrenkiella parvula using comparative genomics, transcriptomics, metabolomics, and ionomics . S. parvula maintains a lower level of total boron and free boric acid in its roots and shoots and sustains growth for longer durations than A. thaliana when grown with excess boron. S. parvula likely excludes boron more efficiently than A. thaliana , which we propose is partly driven by BOR5, a boron transporter that we functionally characterized in the current study. Both species allocate significant transcriptomic and metabolomic resources to enable their cell walls to serve as a partial sink for excess boron, particularly discernable in A. thaliana shoots. We provide evidence that the S. parvula transcriptome is pre-adapted to boron toxicity, exhibiting substantial overlap with the boron-stressed transcriptome of A. thaliana . Our transcriptomic and metabolomics data also suggest that RNA metabolism is a primary target of boron toxicity. Cytoplasmic boric acid likely forms complexes with ribose and ribose-containing compounds critical to RNA and other primary metabolic functions. A model depicting some of the cellular responses that enable a plant to grow in the presence of normally toxic levels of boron is presented.

Abstract The rapid invasion of the non-native Phragmites australis (Poaceae, subfamily Arundinoideae) is a major threat to native ecosystems in North America. We describe a 1.14 Gbp reference genome for P. australis and compare invasive (ssp. australis ) and native (ssp. americanus ) genotypes collected across the Laurentian Great Lakes to deduce genomic bases driving its invasive success. We report novel genomic features including a lineage-specific whole genome duplication, followed by gene loss and preferential retention of genes associated with transcription factors and regulatory functions in the remaining duplicates. The comparative transcriptomic analyses revealed that genes associated with biotic stress and defense responses were expressed at a higher basal level in invasive genotypes, but the native genotypes showed a stronger induction of defense responses following fungal inoculation. The reference genome and transcriptomes, combined with previous ecological and environmental data, support the development of novel, genomics-assisted management approaches for invasive Phragmites .

We report a new visualization tool for analysis of whole genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by NCBI using assembly alignment algorithms developed by us and others. Researchers can examine the alignments between the two assemblies using two alternate views: a chromosome ideogram-based view or a 2D genome dotplot. Whole genome alignment views expose large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed in the ideogram view where available. CGV currently provides approximately 700 alignments from over 300 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

Phosphorus (P) is an essential plant macronutrient vital to fundamental metabolic processes. Plant-available P is low in most soils, making it a frequent limiter of growth. Declining P reserves for fertilizer production exasperates this agricultural challenge. Plants modulate complex responses to fluctuating P levels via global transcriptional regulatory networks. Although chromatin structure plays a substantial role in controlling gene expression, the chromatin dynamics involved in regulating P homeostasis have not been determined. Here we define distinct chromatin states across the rice genome by integrating multiple aspects of chromatin structure, including the H2A.Z histone variant, H3K4me3 modification, and nucleosome positioning. In response to P starvation, 40% of all protein-coding genes exhibit a transition from one chromatin state to another at their transcription start site. Several of these transitions are enriched in subsets of genes differentially expressed by P deficiency. The most prominent subset supports the presence of a coordinated signaling network that targets cell wall structure and is regulated in part via a decrease of H3K4me3 at the transcription start site. The P-starvation induced chromatin dynamics and correlated genes identified here will aid in enhancing P-use efficiency in crop plants, benefitting global agriculture.

Local adaptation is often a product of environmental variations in the geographical space and has implications for biodiversity conservation. We investigated the role of latitudinal heterogeneity in climate on the organization of genetic and phenotypic variation in the dominant coastal tree, Avicennia schaueriana . In a common garden experiment, samples from an equatorial region, marked by rainy/dry seasons, accumulated less biomass, showed lower stomatal conductance and transpiration, narrower xylem vessels, smaller leaves and higher reflectance of long wavelengths (red light) on the stem epidermis, than samples from a subtropical region, marked by warm/cold seasons. Transcriptome differences identified between trees sampled under field conditions at equatorial and subtropical sites, were enriched in functional categories as responses to temperature, solar radiation, water deficit, photosynthesis and cell wall biosynthesis. The diversity based on thousands of SNP loci revealed a north-south genetic structure. Remarkably, signatures of selection were identified in loci associated with photosynthesis, anthocyanin accumulation and the responses to osmotic and hypoxia stresses. Our results suggest the existence of divergence in key resource-use characteristics, likely driven by climate seasonality, based on water-deficit and solar radiation. These findings provide a basis for conservation plans and for predictions for coastal plant responses to climate change.

We developed the CLfinder-OrthNet pipeline that detects co-linearity in gene arrangement among multiple closely related genomes; find ortholog groups; and encodes the evolutionary history of each ortholog group into a representative network (OrthNet). Using a search based on network topology, out of a total of 17,432 OrthNets in six Brassicaceae genomes, we identified 1,394 that included gene transposition-duplication (tr-d) events in one or more genomes. Occurrences of tr-d shared by subsets of Brassicaceae genomes mirrored the divergence times between the genomes and their repeat contents. The majority of tr-d events resulted in truncated open reading frames (ORFs) in the duplicated loci. However, the duplicates with complete ORFs were significantly more frequent than expected from random events. They also had a higher chance of being expressed and derived from older tr-d events. We also found an enrichment, compared to random chance, of tr-d events with complete loss of intergenic sequence conservation between the original and duplicated loci. Finally, we identified tr-d events uniquely found in two extremophytes among the six Brassicaceae genomes, including tr-d of SALT TOLERANCE 32 and ZINC TRANSPORTER 3. The CLfinder-OrthNet pipeline provides a flexible and a modular toolkit to compare gene order, encode and visualize evolutionary paths among orthologs as networks, and identify all gene loci that share the same evolutionary history using network topology searches.

Environmental variation along the geographical space can shape populations by natural selection. In the context of global warming, accompanied by substantial changes in precipitation regimes, it is crucial to understand the role of environmental heterogeneity in tropical trees adaptation, given their disproportional contribution to water and carbon biogeochemical cycles. Here we investigated how heterogeneity in freshwater availability along tropical wetlands has influenced molecular variations of the Black-Mangrove ( Avicennia germinans ). Fifty-seven trees were sampled in seven sites differing markedly on precipitation regime and riverine freshwater inputs. Using 2,297 genome-wide single nucleotide polymorphic markers, we found signatures of natural selection by the genotype association with the precipitation of the warmest quarter and the annual precipitation. We also found candidate loci for selection, based on statistical deviations from neutral expectations of interpopulation differentiation. Most candidate loci present within coding sequences were functionally associated with central aspects of drought-tolerance or plant response to drought. Complementarily, our results suggest the occurrence of rapid evolution of a population, likely in response to sudden and persistent limitation in plant access to soil water, following a road construction in 1974. Observations supporting rapid evolution included reduction in tree size and changes in allele frequencies and in transcripts expression levels associated with increased drought-tolerance, through accumulation of osmoprotectants and antioxidants, biosynthesis of plant cuticle, protection against stress-induced proteins degradation, stomatal closure, photorespiration and photosynthesis. We describe a major role of spatial heterogeneity in freshwater availability in the specialization of this typically tropical tree.

Abstract Background Autographa californica Multiple Nucleopolyhedrovirus (AcMNPV) is a baculovirus with a high potential for its use as a biopesticide against arthropod pests. The budded form of the virus causes a systemic infection when it escapes the midgut to enter the hemolymph of susceptible hosts. Yet, the specific molecular processes underlying the biocidal activity of AcMNPV on its insect hosts are largely unknown. Results In this study, we describe the transcriptional responses in two major pests, Spodoptera frugiperda and Trichoplusia ni , to determine the host-pathogen responses during AcMNPV infection, concurrently with the viral response to the host. We assembled species-specific de novo reference transcriptomes of the hemolymph to identify key transcripts that respond during pathogenesis in these arthropod models where genomic resources are sparse. We found that the suppression of transcriptional processes related to chitin, a metabolite critical for basement membrane stability and tracheal development are central in establishing a systemic infection. Synergistic transcriptional support was observed to suggest suppression of immune responses and induction of oxidative stress indicating disease progression in the host. The entire AcMNPV core genome was expressed in the host hemolymph and viral genes predominantly associated with the budded virus replication, structure, and movement were more abundant than those associated with the occlusion-derived virus. Genes known to directly arrest host cell cycle and development were among the most abundant AcMNPV transcripts in infected hosts. Interestingly, several of the host genes (e.g. Chitin synthase ) that were targeted by the pathogen as revealed by our study are also targets of several chemical insecticides currently used commercially to control arthropod pests. Conclusions Our results reveal an extensive overlap between biological processes represented by genes differently expressed in both hosts, as well as convergence on highly abundant viral genes expressed in the two hosts, providing an overview of the host-pathogen transcriptomic landscape during systemic infection. Given the diversity of AcMNPV strains that infect a wide range of insect hosts, our study provides a framework where pathogen strains could be selected to target specific host genes that facilitates modulation of the infection strength and specificity of the susceptible hosts.