Multiple mechanisms explain loss of anthocyanins from betalain-pigmented Caryophyllales, including repeated wholesale loss of a key anthocyanidin synthesis enzyme.

In this study, we investigate the genetic mechanisms responsible for the loss of anthocyanins in betalain-pigmented Caryophyllales, considering our hypothesis of multiple transitions to betalain pigmentation. Utilizing transcriptomic and genomic datasets across 357 species and 31 families, we scrutinize 18 flavonoid pathway genes and six regulatory genes spanning four transitions to betalain pigmentation. We examined evidence for hypotheses of wholesale gene loss, modified gene function, altered gene expression, and degeneration of the MBW (MYB-bHLH-WD40) trasnscription factor complex, within betalain-pigmented lineages. Our analyses reveal that most flavonoid synthesis genes remain conserved in betalain-pigmented lineages, with the notable exception of TT19 orthologs, essential for the final step in anthocyanidin synthesis, which appear to have been repeatedly and entirely lost. Additional late-stage flavonoid pathway genes upstream of TT19 also manifest strikingly reduced expression in betalain-pigmented species. Additionally, we find repeated loss and alteration in the MBW transcription complex essential for canonical anthocyanin synthesis. Consequently, the loss and exclusion of anthocyanins in betalain-pigmented species appear to be orchestrated through several mechanisms: loss of a key enzyme, downregulation of synthesis genes, and degeneration of regulatory complexes. These changes have occurred iteratively in Caryophyllales, often coinciding with evolutionary transitions to betalain pigmentation.

Transcriptional regulation of wound suberin deposition in potato cultivars with differential wound healing capacity.

Wounding during mechanical harvesting and post-harvest handling results in tuber desiccation and provides an entry point for pathogens resulting in substantial post​-harvest crop losses. Poor wound healing is a major culprit of these losses. Wound tissue in potato (Solanum tuberosum) tubers, and all higher plants, is composed of a large proportion of suberin that is deposited in a specialized tissue called the wound periderm. However, the genetic regulatory pathway controlling wound-induced suberization remains unknown. Here, we implicate two potato transcription factors, StMYB102 (PGSC0003DMG400011250) and StMYB74 (PGSC0003DMG400022399), as regulators of wound suberin biosynthesis and deposition. Using targeted metabolomics and transcript profiling from the wound healing tissues of two commercial potato cultivars, as well as heterologous expression, we provide evidence for the molecular-genetic basis of the differential wound suberization capacities of different potato cultivars. Our results suggest that (i) the export of suberin from the cytosol to the apoplast and ligno-suberin deposition may be limiting factors for wound suberization, (ii) StMYB74 and StMYB102 are important regulators of the wound suberization process in tubers, and (iii) polymorphisms in StMYB102 may influence cultivar-specific wound suberization capacity. These results represent an important step in understanding the regulated biosynthesis and deposition of wound suberin and provide a practical foundation for targeted breeding approaches aimed at improving potato tuber storage life.

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l.

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations. [Amaranthaceae; gene tree discordance; hybridization; incomplete lineage sorting; phylogenomics; species network; species tree; transcriptomics.].

Membrane Profiling by Free Flow Electrophoresis and SWATH-MS to Characterize Subcellular Compartment Proteomes in Mesembryanthemum crystallinum.

The study of subcellular membrane structure and function facilitates investigations into how biological processes are divided within the cell. However, work in this area has been hampered by the limited techniques available to fractionate the different membranes. Free Flow Electrophoresis (FFE) allows for the fractionation of membranes based on their different surface charges, a property made up primarily of their varied lipid and protein compositions. In this study, high-resolution plant membrane fractionation by FFE, combined with mass spectrometry-based proteomics, allowed the simultaneous profiling of multiple cellular membranes from the leaf tissue of the plant Mesembryanthemum crystallinum. Comparisons of the fractionated membranes' protein profile to that of known markers for specific cellular compartments sheds light on the functions of proteins, as well as provides new evidence for multiple subcellular localization of several proteins, including those involved in lipid metabolism.

Evolution of l-DOPA 4,5-dioxygenase activity allows for recurrent specialisation to betalain pigmentation in Caryophyllales.

The evolution of l-DOPA 4,5-dioxygenase activity, encoded by the gene DODA, was a key step in the origin of betalain biosynthesis in Caryophyllales. We previously proposed that l-DOPA 4,5-dioxygenase activity evolved via a single Caryophyllales-specific neofunctionalisation event within the DODA gene lineage. However, this neofunctionalisation event has not been confirmed and the DODA gene lineage exhibits numerous gene duplication events, whose evolutionary significance is unclear. To address this, we functionally characterised 23 distinct DODA proteins for l-DOPA 4,5-dioxygenase activity, from four betalain-pigmented and five anthocyanin-pigmented species, representing key evolutionary transitions across Caryophyllales. By mapping these functional data to an updated DODA phylogeny, we then explored the evolution of l-DOPA 4,5-dioxygenase activity. We find that low l-DOPA 4,5-dioxygenase activity is distributed across the DODA gene lineage. In this context, repeated gene duplication events within the DODA gene lineage give rise to polyphyletic occurrences of elevated l-DOPA 4,5-dioxygenase activity, accompanied by convergent shifts in key functional residues and distinct genomic patterns of micro-synteny. In the context of an updated organismal phylogeny and newly inferred pigment reconstructions, we argue that repeated convergent acquisition of elevated l-DOPA 4,5-dioxygenase activity is consistent with recurrent specialisation to betalain synthesis in Caryophyllales.

Temporal and spatial transcriptomic and microRNA dynamics of CAM photosynthesis in pineapple.

The altered carbon assimilation pathway of crassulacean acid metabolism (CAM) photosynthesis results in an up to 80% higher water-use efficiency than C3 photosynthesis in plants making it a potentially useful pathway for engineering crop plants with improved drought tolerance. Here we surveyed detailed temporal (diel time course) and spatial (across a leaf gradient) gene and microRNA (miRNA) expression patterns in the obligate CAM plant pineapple [Ananas comosus (L.) Merr.]. The high-resolution transcriptome atlas allowed us to distinguish between CAM-related and non-CAM gene copies. A differential gene co-expression network across green and white leaf diel datasets identified genes with circadian oscillation, CAM-related functions, and source-sink relations. Gene co-expression clusters containing CAM pathway genes are enriched with clock-associated cis-elements, suggesting circadian regulation of CAM. About 20% of pineapple microRNAs have diel expression patterns, with several that target key CAM-related genes. Expression and physiology data provide a model for CAM-specific carbohydrate flux and long-distance hexose transport. Together these resources provide a list of candidate genes for targeted engineering of CAM into C3 photosynthesis crop species.

Sporobolus stapfianus: Insights into desiccation tolerance in the resurrection grasses from linking transcriptomics to metabolomics.

BACKGROUND: Understanding the response of resurrection angiosperms to dehydration and rehydration is critical for deciphering the mechanisms of how plants cope with the rigors of water loss from their vegetative tissues. We have focused our studies on the C4 resurrection grass, Sporobolus stapfianus Gandoger, as a member of a group of important forage grasses. METHODS: We have combined non-targeted metabolomics with transcriptomics, via a NimbleGen array platform, to develop an understanding of how gene expression and metabolite profiles can be linked to generate a more detailed mechanistic appreciation of the cellular response to both desiccation and rehydration. RESULTS: The rehydration transcriptome and metabolome are primarily geared towards the rapid return of photosynthesis, energy metabolism, protein turnover, and protein synthesis during the rehydration phase. However, there are some metabolites associated with ROS protection that remain elevated during rehydration, most notably the tocopherols. The analysis of the dehydration transcriptome reveals a strong concordance between transcript abundance and the associated metabolite abundance reported earlier, but only in responses that are directly related to cellular protection during dehydration: carbohydrate metabolism and redox homeostasis. The transcriptome response also provides strong support for the involvement of cellular protection processes as exemplified by the increases in the abundance of transcripts encoding late embryogenesis abundant (LEA) proteins, anti-oxidant enzymes, early light-induced proteins (ELIP) proteins, and cell-wall modification enzymes. There is little concordance between transcript and metabolite abundance for processes such as amino acid metabolism that do not appear to contribute directly to cellular protection, but are nonetheless important for the desiccation tolerant phenotype of S. stapfianus. CONCLUSIONS: The transcriptomes of both dehydration and rehydration offer insight into the complexity of the regulation of responses to these processes that involve complex signaling pathways and associated transcription factors. ABA appears to be important in the control of gene expression in both the latter stages of the dehydration and the early stages of rehydration. These findings add to the growing body of information detailing how plants tolerate and survive the severe cellular perturbations of dehydration, desiccation, and rehydration.

An rbcL mRNA-binding protein is associated with C3 to C4 evolution and light-induced production of Rubisco in Flaveria.

Nuclear-encoded RLSB protein binds chloroplastic rbcL mRNA encoding the Rubisco large subunit. RLSB is highly conserved across all groups of land plants and is associated with positive post-transcriptional regulation of rbcL expression. In C3 leaves, RLSB and Rubisco occur in all chlorenchyma cell chloroplasts, while in C4 leaves these accumulate only within bundle sheath (BS) chloroplasts. RLSB's role in rbcL expression makes modification of its localization a likely prerequisite for the evolutionary restriction of Rubisco to BS cells. Taking advantage of evolutionarily conserved RLSB orthologs in several C3, C3-C4, C4-like, and C4 photosynthetic types within the genus Flaveria, we show that low level RLSB sequence divergence and modification to BS specificity coincided with ontogeny of Rubisco specificity and Kranz anatomy during C3 to C4 evolution. In both C3 and C4 species, Rubisco production reflected RLSB production in all cell types, tissues, and conditions examined. Co-localization occurred only in photosynthetic tissues, and both proteins were co-ordinately induced by light at post-transcriptional levels. RLSB is currently the only mRNA-binding protein to be associated with rbcL gene regulation in any plant, with variations in sequence and acquisition of cell type specificity reflecting the progression of C4 evolution within the genus Flaveria.

The highly improved genome of Ixodes scapularis with X and Y pseudochromosomes.

Ixodes scapularis, the black-legged tick, is the principal vector of the Lyme disease spirochete, Borrelia burgdorferi, and is responsible for most of the ∼470,000 estimated Lyme disease cases annually in the USA. Ixodes scapularis can transmit six additional pathogens of human health significance. Because of its medical importance, I. scapularis was the first tick genome to be sequenced and annotated. However, the first assembly, I. scapularis Wikel (IscaW), was highly fragmented because of the technical challenges posed by the long, repetitive genome sequences characteristic of arthropod genomes and the lack of long-read sequencing techniques. Although I. scapularis has emerged as a model for tick research because of the availability of new tools such as embryo injection and CRISPR-Cas9-mediated gene editing yet the lack of chromosome-scale scaffolds has slowed progress in tick biology and the development of tools for their control. Here we combine diverse technologies to produce the I. scapularis Gulia-Nuss (IscGN) genome assembly and gene set. We used DNA from eggs and male and female adult ticks and took advantage of Hi-C, PacBio HiFi sequencing, and Illumina short-read sequencing technologies to produce a chromosome-level assembly. In this work, we present the predicted pseudochromosomes consisting of 13 autosomes and the sex pseudochromosomes: X and Y, and a markedly improved genome annotation compared with the existing assemblies and annotations.

Cas9-mediated gene editing in the black-legged tick, Ixodes scapularis, by embryo injection and ReMOT Control.

Despite their capacity to acquire and pass on an array of debilitating pathogens, research on ticks has lagged behind other arthropod vectors, such as mosquitoes, largely because of challenges in applying available genetic and molecular tools. CRISPR-Cas9 is transforming non-model organism research; however, successful gene editing has not yet been reported in ticks. Technical challenges for injecting tick embryos to attempt gene editing have further slowed research progress. Currently, no embryo injection protocol exists for any chelicerate species, including ticks. Herein, we report a successful embryo injection protocol for the black-legged tick, Ixodes scapularis, and the use of this protocol for genome editing with CRISPR-Cas9. We also demonstrate that the ReMOT Control technique could be successfully used to generate genome mutations outside Insecta. Our results provide innovative tools to the tick research community that are essential for advancing our understanding of the molecular mechanisms governing pathogen transmission by tick vectors and the underlying biology of host-vector-pathogen interactions.

Mitochondrial haplotypes are not associated with mice selectively bred for high voluntary wheel running.

Mitochondrial haplotypes have been associated with human and rodent phenotypes, including nonshivering thermogenesis capacity, learning capability, and disease risk. Although the mammalian mitochondrial D-loop is highly polymorphic, D-loops in laboratory mice are identical, and variation occurs elsewhere mainly between nucleotides 9820 and 9830. Part of this region codes for the tRNAArg gene and is associated with mitochondrial densities and number of mtDNA copies. We hypothesized that the capacity for high levels of voluntary wheel-running behavior would be associated with mitochondrial haplotype. Here, we analyzed the mtDNA polymorphic region in mice from each of four replicate lines selectively bred for 54 generations for high voluntary wheel running (HR) and from four control lines (Control) randomly bred for 54 generations. Sequencing the polymorphic region revealed a variable number of adenine repeats. Single nucleotide polymorphisms (SNPs) varied from 2 to 3 adenine insertions, resulting in three haplotypes. We found significant genetic differentiations between the HR and Control groups (Fst = 0.779, p ≤ 0.0001), as well as among the replicate lines of mice within groups (Fsc = 0.757, p ≤ 0.0001). Haplotypes, however, were not strongly associated with voluntary wheel running (revolutions run per day), nor with either body mass or litter size. This system provides a useful experimental model to dissect the physiological processes linking mitochondrial, genomic SNPs, epigenetics, or nuclear-mitochondrial cross-talk to exercise activity.

Crassulacean Acid Metabolism Abiotic Stress-Responsive Transcription Factors: a Potential Genetic Engineering Approach for Improving Crop Tolerance to Abiotic Stress.

This perspective paper explores the utilization of abiotic stress-responsive transcription factors (TFs) from crassulacean acid metabolism (CAM) plants to improve abiotic stress tolerance in crop plants. CAM is a specialized type of photosynthetic adaptation that enhances water-use efficiency (WUE) by shifting CO2 uptake to all or part of the nighttime when evaporative water losses are minimal. Recent studies have shown that TF-based genetic engineering could be a useful approach for improving plant abiotic stress tolerance because of the role of TFs as master regulators of clusters of stress-responsive genes. Here, we explore the use of abiotic stress-responsive TFs from CAM plants to improve abiotic stress tolerance and WUE in crops by controlling the expression of gene cohorts that mediate drought-responsive adaptations. Recent research has revealed several TF families including AP2/ERF, MYB, WRKY, NAC, NF-Y, and bZIP that might regulate water-deficit stress responses and CAM in the inducible CAM plant Mesembryanthemum crystallinum under water-deficit stress-induced CAM and in the obligate CAM plant Kalanchoe fedtschenkoi. Overexpression of genes from these families in Arabidopsis thaliana can improve abiotic stress tolerance in A. thaliana in some instances. Therefore, we propose that TF-based genetic engineering with a small number of CAM abiotic stress-responsive TFs will be a promising strategy for improving abiotic stress tolerance and WUE in crop plants in a projected hotter and drier landscape in the 21st-century and beyond.

The bracteatus pineapple genome and domestication of clonally propagated crops.

Domestication of clonally propagated crops such as pineapple from South America was hypothesized to be a 'one-step operation'. We sequenced the genome of Ananas comosus var. bracteatus CB5 and assembled 513 Mb into 25 chromosomes with 29,412 genes. Comparison of the genomes of CB5, F153 and MD2 elucidated the genomic basis of fiber production, color formation, sugar accumulation and fruit maturation. We also resequenced 89 Ananas genomes. Cultivars 'Smooth Cayenne' and 'Queen' exhibited ancient and recent admixture, while 'Singapore Spanish' supported a one-step operation of domestication. We identified 25 selective sweeps, including a strong sweep containing a pair of tandemly duplicated bromelain inhibitors. Four candidate genes for self-incompatibility were linked in F153, but were not functional in self-compatible CB5. Our findings support the coexistence of sexual recombination and a one-step operation in the domestication of clonally propagated crops. This work guides the exploration of sexual and asexual domestication trajectories in other clonally propagated crops.

Draft Nuclear Genome Sequence of the Halophilic and Beta-Carotene-Accumulating Green Alga Dunaliella salina Strain CCAP19/18.

The halotolerant alga Dunaliella salina is a model for stress tolerance and is used commercially for production of beta-carotene (=pro-vitamin A). The presented draft genome of the genuine strain CCAP19/18 will allow investigations into metabolic processes involved in regulation of stress responses, including carotenogenesis and adaptations to life in high-salinity environments.

Identification of Ice Plant (Mesembryanthemum crystallinum L.) MicroRNAs Using RNA-Seq and Their Putative Roles in High Salinity Responses in Seedlings.

The halophyte Mesembryanthemum crystallinum (common or crystalline ice plant) is a useful model for studying molecular mechanisms of salt tolerance. The morphology, physiology, metabolism, and gene expression of ice plant have been studied and large-scale analyses of gene expression profiling have drawn an outline of salt tolerance in ice plant. A rapid root growth to a sudden increase in salinity was observed in ice plant seedlings. Using a fluorescent dye to detect Na(+), we found that ice plant roots respond to an increased flux of Na(+) by either secreting or storing Na(+) in specialized cells. High-throughput sequencing was used to identify small RNA profiles in 3-day-old seedlings treated with or without 200 mM NaCl. In total, 135 conserved miRNAs belonging to 21 families were found. The hairpin precursor of 19 conserved mcr-miRNAs and 12 novel mcr-miRNAs were identified. After 6 h of salt stress, the expression of most mcr-miRNAs showed decreased relative abundance, whereas the expression of their corresponding target genes showed increased mRNA relative abundance. The cognate target genes are involved in a broad range of biological processes: transcription factors that regulate growth and development, enzymes that catalyze miRNA biogenesis for the most conserved mcr-miRNA, and proteins that are involved in ion homeostasis and drought-stress responses for some novel mcr-miRNAs. Analyses of the functions of target genes revealed that cellular processes, including growth and development, metabolism, and ion transport activity are likely to be enhanced in roots under salt stress. The expression of eleven conserved miRNAs and two novel miRNAs were correlated reciprocally with predicted targets within hours after salt stress exposure. Several conserved miRNAs have been known to regulate root elongation, root apical meristem activity, and lateral root formation. Based upon the expression pattern of miRNA and target genes in combination with the observation of Na(+) distribution, ice plant likely responds to increased salinity by using Na(+) as an osmoticum for cell expansion and guard cell opening. Excessive Na(+) could either be secreted through the root epidermis or stored in specialized leaf epidermal cells. These responses are regulated in part at the miRNA-mediated post-transcriptional level.

The pineapple genome and the evolution of CAM photosynthesis.

Pineapple (Ananas comosus (L.) Merr.) is the most economically valuable crop possessing crassulacean acid metabolism (CAM), a photosynthetic carbon assimilation pathway with high water-use efficiency, and the second most important tropical fruit. We sequenced the genomes of pineapple varieties F153 and MD2 and a wild pineapple relative, Ananas bracteatus accession CB5. The pineapple genome has one fewer ancient whole-genome duplication event than sequenced grass genomes and a conserved karyotype with seven chromosomes from before the ρ duplication event. The pineapple lineage has transitioned from C3 photosynthesis to CAM, with CAM-related genes exhibiting a diel expression pattern in photosynthetic tissues. CAM pathway genes were enriched with cis-regulatory elements associated with the regulation of circadian clock genes, providing the first cis-regulatory link between CAM and circadian clock regulation. Pineapple CAM photosynthesis evolved by the reconfiguration of pathways in C3 plants, through the regulatory neofunctionalization of preexisting genes and not through the acquisition of neofunctionalized genes via whole-genome or tandem gene duplication.