WRKY6 transcription factor modulates root potassium acquisition through promoting expression of AKT1 in Arabidopsis.

Potassium (K+), being an essential macronutrient in plants, plays a central role in many aspects. Root growth is highly plastic and is affected by many different abiotic stresses including nutrient deficiency. The Shaker-type K+ channel Arabidopsis (Arabidopsis thaliana) K+ Transporter 1 (AKT1) is responsible for K+ uptake under both low and high external K+ conditions. However, the upstream transcription factor of AKT1 is not clear. Here, we demonstrated that the WRKY6 transcription factor modulates root growth to low potassium (LK) stress in Arabidopsis. WRKY6 showed a quick response to LK stress and also to many other abiotic stress treatments. The two wrky6 T-DNA insertion mutants were highly sensitive to LK treatment, whose primary root lengths were much shorter, less biomass and lower K+ content in roots than those of wild-type plants, while WRKY6-overexpression lines showed opposite phenotypes. A further investigation showed that WRKY6 regulated the expression of the AKT1 gene via directly binding to the W-box elements in its promoter through EMSA and ChIP-qPCR assays. A dual luciferase reporter analysis further demonstrated that WRKY6 enhanced the transcription of AKT1. Genetic analysis further revealed that the overexpression of AKT1 greatly rescued the short root phenotype of the wrky6 mutant under LK stress, suggesting AKT1 is epistatic to WRKY6 in the control of LK response. Further transcriptome profiling suggested that WRKY6 modulates LK response through a complex regulatory network. Thus, this study unveils a transcription factor that modulates root growth under potassium deficiency conditions by affecting the potassium channel gene AKT1 expression.

Hybrid Genome Assembly of Berkeleyomyces rouxiae, an Emerging Cannabis Fungal Pathogen Causing Black Root Rot in an Aeroponic Facility.

The resurged interest in cultivation of Cannabis sativa has presented an array of new challenges. Among them are the difficult-to-control pests and pathogens that infect cannabis plants. The limited methods for disease control available to cannabis growers necessitates early detection of plant pathogens, something that molecular techniques such as DNA sequencing has greatly improved. This study reports for the first time the fungal plant pathogen Berkeleyomyces rouxiae causing black root rot in high THC-containing cannabis. Aeroponically grown cannabis plants at a licenced production facility in Cranbrook BC, Canada, rapidly displayed root discoloration and rot symptoms despite testing negative for all commercially available pathogen tests. Developing sequencing-based disease diagnostics requires genomic information, so this study presents the first whole genome sequence of the multihost, widespread black root rot pathogen B. rouxiae. Hybrid genome assembly using Oxford Nanopore long-reads and Illumina short-reads yielded a genome size of 28.2 Mb represented over 404 contigs with an N50 of 267 kb. Genome annotation predicted 6,960 protein-coding genes with 59,477 functional annotations. The availability of this genome will assist in sequence-based diagnostic development, comparative genomics, and taxonomic resolution of this globally important plant pathogen.

Evaluation of Differentially Expressed Genes in Leaves vs. Roots Subjected to Drought Stress in Flax (Linum usitatissimum L.).

Drought stress is a common environmental challenge that plants face, severely constraining plant growth and reducing crop yield and quality. Several studies have highlighted distinct responses between monocotyledonous and dicotyledonous plants. However, the mechanisms underlying flax tolerance to abiotic stress, such as drought, remain unclear. In this study, we investigated the morphological, physiological, and biochemical characteristics and the genome-wide gene expression of oil flax and fiber flax in response to drought stress. The results revealed that drought stress caused significant wilting of flax leaves. Within the first 24 h of stress, various physiological and biochemical characteristics exhibited rapid responses. These included fresh weight, relative water content (RWC), proline, soluble protein, soluble sugar, superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) in the leaves or roots of flax. Additionally, drought stress led to a significant rise in lignin content in fiber flax. In addition, the transcriptome analysis demonstrated genome-wide variations in gene expression induced by drought stress. Specifically, genes associated with photosynthesis, proline biosynthesis, and phytohormone metabolism exhibited significant differences in expression levels under stress conditions in flax. These findings highlight the rapid response of flax to drought stress within a short-term period. Our experiment also revealed that, although there were variations in the levels of small compound content or gene expression between Longya10 and Fany under drought stress, most stress-resistance responses were similar. Furthermore, the results provide additional evidence supporting the existence of mechanisms underlying the response to drought stress in plants.

Genome-Wide Analysis of Flax (Linum usitatissimum L.) Growth-Regulating Factor (GRF) Transcription Factors.

Flax is an important cash crop globally with a variety of commercial uses. It has been widely used for fiber, oil, nutrition, feed and in composite materials. Growth regulatory factor (GRF) is a transcription factor family unique to plants, and is involved in regulating many processes of growth and development. Bioinformatics analysis of the GRF family in flax predicted 17 LuGRF genes, which all contained the characteristic QLQ and WRC domains. Equally, 15 of 17 LuGRFs (88%) are predicted to be regulated by lus-miR396 miRNA. Phylogenetic analysis of GRFs from flax and several other well-characterized species defined five clades; LuGRF genes were found in four clades. Most LuGRF gene promoters contained cis-regulatory elements known to be responsive to hormones and stress. The chromosomal locations and collinearity of LuGRF genes were also analyzed. The three-dimensional structure of LuGRF proteins was predicted using homology modeling. The transcript expression data indicated that most LuGRF family members were highly expressed in flax fruit and embryos, whereas LuGRF3, LuGRF12 and LuGRF16 were enriched in response to salt stress. Real-time quantitative fluorescent PCR (qRT-PCR) showed that both LuGRF1 and LuGRF11 were up-regulated under ABA and MeJA stimuli, indicating that these genes were involved in defense. LuGRF1 was demonstrated to be localized to the nucleus as expected for a transcription factor. These results provide a basis for further exploration of the molecular mechanism of LuGRF gene function and obtaining improved flax breeding lines.

Ectopic overexpression of a membrane-tethered transcription factor gene NAC60 from oilseed rape positively modulates programmed cell death and age-triggered leaf senescence.

Senescence is an integrative final stage of plant development that is governed by internal and external cues. The NAM, ATAF1/2, CUC2 (NAC) transcription factor (TF) family is specific to plants and membrane-tethered NAC TFs (MTTFs) constitute a unique and sophisticated mechanism in stress responses and development. However, the function of MTTFs in oilseed rape (Brassica napus L.) remains unknown. Here, we report that BnaNAC60 is an MTTF associated with the endoplasmic reticulum (ER) membrane. Expression of BnaNAC60 was induced during the progression of leaf senescence. Translocation of BnaNAC60 into nuclei was induced by ER stress and oxidative stress treatments. It binds to the NTLBS motif, rather than the canonical NAC recognition site. Overexpression of BnaNAC60 devoid of the transmembrane domain, but not the full-length BnaNAC60, induces significant reactive oxygen species (ROS) accumulation and hypersensitive response-like cell death in both tobacco (Nicotiana benthamiana) and oilseed rape protoplasts. Moreover, ectopic overexpression of BnaNAC60 devoid of the transmembrane domain, but not the full-length BnaNAC60, in Arabidopsis also induces precocious leaf senescence. Furthermore, screening and expression profiling identified an array of functional genes that are significantly induced by BnaNAC60 expression. Further it was found that BnaNAC60 can activate the promoter activities of BnaNYC1, BnaRbohD, BnaBFN1, BnaZAT12, and multiple BnaVPEs in a dual-luciferase reporter assay. Electrophoretic mobility shift assay and chromatin immunoprecipitation coupled to quantitative PCR assays revealed that BnaNAC60 directly binds to the promoter regions of these downstream target genes. To summarize, our data show that BnaNAC60 is an MTTF that modulates cell death, ROS accumulation, and leaf senescence.

WRKY42 transcription factor positively regulates leaf senescence through modulating SA and ROS synthesis in Arabidopsis thaliana.

Leaf senescence represents the final stage of leaf growth and development, and its onset and progression are strictly regulated; however, the underlying regulatory mechanisms remain largely unknown. In this study we found that WRKY42 was highly induced during leaf senescence. Loss-of-function wrky42 mutants showed delayed leaf senescence whereas the overexpression of WRKY42 accelerated senescence. Transcriptome analysis revealed 2721 differentially expressed genes between wild-type and WRKY42-overexpressing plants, including genes involved in salicylic acid (SA) and reactive oxygen species (ROS) synthesis as well as several senescence-associated genes (SAGs). Moreover, WRKY42 activated the transcription of isochorismate synthase 1 (ICS1), respiratory burst oxidase homolog F (RbohF) and a few SAG genes. Consistently, the expression of these genes was reduced in wrky42 mutants but was markedly increased in transgenic Arabidopsis overexpressing WRKY42. Both in vitro electrophoretic mobility shift assays (EMSAs) and in vivo chromatin immunoprecipitation and dual luciferase assays demonstrated that WRKY42 directly bound to the promoters of ICS1 and RbohF, as well as a few SAGs, to activate their expression. Genetic analysis further showed that mutations of ICS1 and RbohF suppressed the early senescence phenotype evoked by WRKY42 overexpression. Thus, we have identified WRKY42 as a novel transcription factor positively regulating leaf senescence by directly activating the transcription of ICS1, RbohF and SAGs, without any seed yield penalty.

Organellomic data sets confirm a cryptic consensus on (unrooted) land-plant relationships and provide new insights into bryophyte molecular evolution.

PREMISE: Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution. METHODS: We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements. RESULTS: Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses. CONCLUSIONS: A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure.

Chromosome-scale pseudomolecules refined by optical, physical and genetic maps in flax.

Genomes of varying sizes have been sequenced with next-generation sequencing platforms. However, most reference sequences include draft unordered scaffolds containing chimeras caused by mis-scaffolding. A BioNano genome (BNG) optical map was constructed to improve the previously sequenced flax genome (Linum usitatissimum L., 2n = 30, about 373 Mb), which consisted of 3852 scaffolds larger than 1 kb and totalling 300.6 Mb. The high-resolution BNG map of cv. CDC Bethune totalled 317 Mb and consisted of 251 BNG contigs with an N50 of 2.15 Mb. A total of 622 scaffolds (286.6 Mb, 94.9%) aligned to 211 BNG contigs (298.6 Mb, 94.2%). Of those, 99 scaffolds, diagnosed to contain assembly errors, were refined into 225 new scaffolds. Using the newly refined scaffold sequences and the validated bacterial artificial chromosome-based physical map of CDC Bethune, the 211 BNG contigs were scaffolded into 94 super-BNG contigs (N50 of 6.64 Mb) that were further assigned to the 15 flax chromosomes using the genetic map. The pseudomolecules total about 316 Mb, with individual chromosomes of 15.6 to 29.4 Mb, and cover 97% of the annotated genes. Evidence from the chromosome-scale pseudomolecules suggests that flax has undergone palaeopolyploidization and mesopolyploidization events, followed by rearrangements and deletions or fusion of chromosome arms from an ancient progenitor with a haploid chromosome number of eight.

A Novel NAC-Type Transcription Factor, NAC87, from Oilseed Rape Modulates Reactive Oxygen Species Accumulation and Cell Death.

Reactive oxygen species (ROS) are thought to play a dual role in plants by functioning as signaling molecules and toxic by-products of aerobic metabolism. The hypersensitive response (HR) is a typical feature of immune responses in plants and also a type of programmed cell death (PCD). How these two processes are regulated in oilseed rape (Brassica napus L.) at the transcriptional level remains largely unknown. In this study, we report that an oilseed rape (Brassica napus L.) NAM-ATAF-CUC (NAC)-type transcription factor NAC87 modulates ROS and cell death accompanied by typical changes at the morphological and cellular levels. The BnaNAC87 gene was induced by multiple stress and hormone treatments and was highly expressed in senescent leaves by quantitative reverse transcription-PCR (qRT-PCR). BnaNAC87 is located in nuclei and has transcriptional activation activity. Expression of BnaNAC87 promoted significant ROS production, cell death as well as death of protoplasts, as indicated by histological staining. In addition, putative downstream target genes of NAC87 were identified through both qRT-PCR and dual luciferase reporter assays. We found that genes implicated in ROS generation (RbohB), cell death (VPE1a, ZEN1), leaf senescence (WRKY6, ZAT12) and defense (PR2, PR5 and HIN1) were significantly induced. Through an electrophoretic mobility shift assay (EMSA), we confirmed that BnaNAC87 directly binds to the NACRS-containing promoter fragments of ZEN1, ZAT12, HIN1 and PR5 genes. From these results, we conclude that oilseed rape NAC87 is a novel NAC transcription factor that acts as a positive regulator of ROS metabolism and cell death.

Insights into the Evolution of Hydroxyproline-Rich Glycoproteins from 1000 Plant Transcriptomes.

The carbohydrate-rich cell walls of land plants and algae have been the focus of much interest given the value of cell wall-based products to our current and future economies. Hydroxyproline-rich glycoproteins (HRGPs), a major group of wall glycoproteins, play important roles in plant growth and development, yet little is known about how they have evolved in parallel with the polysaccharide components of walls. We investigate the origins and evolution of the HRGP superfamily, which is commonly divided into three major multigene families: the arabinogalactan proteins (AGPs), extensins (EXTs), and proline-rich proteins. Using motif and amino acid bias, a newly developed bioinformatics pipeline, we identified HRGPs in sequences from the 1000 Plants transcriptome project (www.onekp.com). Our analyses provide new insights into the evolution of HRGPs across major evolutionary milestones, including the transition to land and the early radiation of angiosperms. Significantly, data mining reveals the origin of glycosylphosphatidylinositol (GPI)-anchored AGPs in green algae and a 3- to 4-fold increase in GPI-AGPs in liverworts and mosses. The first detection of cross-linking (CL)-EXTs is observed in bryophytes, which suggests that CL-EXTs arose though the juxtaposition of preexisting SPn EXT glycomotifs with refined Y-based motifs. We also detected the loss of CL-EXT in a few lineages, including the grass family (Poaceae), that have a cell wall composition distinct from other monocots and eudicots. A key challenge in HRGP research is tracking individual HRGPs throughout evolution. Using the 1000 Plants output, we were able to find putative orthologs of Arabidopsis pollen-specific GPI-AGPs in basal eudicots.

Impact of Low-Intensity Pulsed Ultrasound on Transcript and Metabolite Abundance in Saccharomyces cerevisiae.

The interactions of ultrasound with biological materials are exploited for diagnostic, interventional, and therapeutic applications in humans and can improve productivity in industrial-scale generation of organic molecules such as biofuels, vaccines, and antibodies. Accordingly, there is great interest in better understanding the biological effects of ultrasound. We studied the impact of low-intensity pulsed ultrasound (LIPUS) on RNA expression and metabolism of S. cerevisiae. Although the transcript expression signature of LIPUS-treated cells does not differ significantly from that of untreated cells after 5 days, metabolomic profiling by chemical-isotopic-labeling-liquid-chromatography-mass-spectrometry suggests that LIPUS has an impact on the pathways of pyrimidine, proline, alanine, aspartate, glutamate, and arginine metabolism. Therefore, LIPUS triggers metabolic effects beyond reprogramming of the core pathways of carbon metabolism. Further characterization of metabolism will likely be important for elucidation of the biological effects of LIPUS.

Phylotranscriptomic analysis of the origin and early diversification of land plants.

Reconstructing the origin and evolution of land plants and their algal relatives is a fundamental problem in plant phylogenetics, and is essential for understanding how critical adaptations arose, including the embryo, vascular tissue, seeds, and flowers. Despite advances in molecular systematics, some hypotheses of relationships remain weakly resolved. Inferring deep phylogenies with bouts of rapid diversification can be problematic; however, genome-scale data should significantly increase the number of informative characters for analyses. Recent phylogenomic reconstructions focused on the major divergences of plants have resulted in promising but inconsistent results. One limitation is sparse taxon sampling, likely resulting from the difficulty and cost of data generation. To address this limitation, transcriptome data for 92 streptophyte taxa were generated and analyzed along with 11 published plant genome sequences. Phylogenetic reconstructions were conducted using up to 852 nuclear genes and 1,701,170 aligned sites. Sixty-nine analyses were performed to test the robustness of phylogenetic inferences to permutations of the data matrix or to phylogenetic method, including supermatrix, supertree, and coalescent-based approaches, maximum-likelihood and Bayesian methods, partitioned and unpartitioned analyses, and amino acid versus DNA alignments. Among other results, we find robust support for a sister-group relationship between land plants and one group of streptophyte green algae, the Zygnematophyceae. Strong and robust support for a clade comprising liverworts and mosses is inconsistent with a widely accepted view of early land plant evolution, and suggests that phylogenetic hypotheses used to understand the evolution of fundamental plant traits should be reevaluated.

Horizontal transfer of an adaptive chimeric photoreceptor from bryophytes to ferns.

Ferns are well known for their shade-dwelling habits. Their ability to thrive under low-light conditions has been linked to the evolution of a novel chimeric photoreceptor--neochrome--that fuses red-sensing phytochrome and blue-sensing phototropin modules into a single gene, thereby optimizing phototropic responses. Despite being implicated in facilitating the diversification of modern ferns, the origin of neochrome has remained a mystery. We present evidence for neochrome in hornworts (a bryophyte lineage) and demonstrate that ferns acquired neochrome from hornworts via horizontal gene transfer (HGT). Fern neochromes are nested within hornwort neochromes in our large-scale phylogenetic reconstructions of phototropin and phytochrome gene families. Divergence date estimates further support the HGT hypothesis, with fern and hornwort neochromes diverging 179 Mya, long after the split between the two plant lineages (at least 400 Mya). By analyzing the draft genome of the hornwort Anthoceros punctatus, we also discovered a previously unidentified phototropin gene that likely represents the ancestral lineage of the neochrome phototropin module. Thus, a neochrome originating in hornworts was transferred horizontally to ferns, where it may have played a significant role in the diversification of modern ferns.

Elucidating steroid alkaloid biosynthesis in Veratrum californicum: production of verazine in Sf9 cells.

Steroid alkaloids have been shown to elicit a wide range of pharmacological effects that include anticancer and antifungal activities. Understanding the biosynthesis of these molecules is essential to bioengineering for sustainable production. Herein, we investigate the biosynthetic pathway to cyclopamine, a steroid alkaloid that shows promising antineoplastic activities. Supply of cyclopamine is limited, as the current source is solely derived from wild collection of the plant Veratrum californicum. To elucidate the early stages of the pathway to cyclopamine, we interrogated a V. californicum RNA-seq dataset using the cyclopamine accumulation profile as a predefined model for gene expression with the pattern-matching algorithm Haystack. Refactoring candidate genes in Sf9 insect cells led to discovery of four enzymes that catalyze the first six steps in steroid alkaloid biosynthesis to produce verazine, a predicted precursor to cyclopamine. Three of the enzymes are cytochromes P450 while the fourth is a γ-aminobutyrate transaminase; together they produce verazine from cholesterol.

Ty1-copia elements reveal diverse insertion sites linked to polymorphisms among flax (Linum usitatissimum L.) accessions.

BACKGROUND: Initial characterization of the flax genome showed that Ty1-copia retrotransposons are abundant, with several members being recently inserted, and in close association with genes. Recent insertions indicate a potential for ongoing transpositional activity that can create genomic diversity among accessions, cultivars or varieties. The polymorphisms generated constitute a good source of molecular markers that may be associated with phenotype if the insertions alter gene activity. Flax, where accessions are bred mainly for seed nutritional properties or for fibers, constitutes a good model for studying the relationship of transpositional activity with diversification and breeding. In this study, we estimated copy number and used a type of transposon display known as Sequence-Specific Amplification Polymorphisms (SSAPs), to characterize six families of Ty1-copia elements across 14 flax accessions. Polymorphic insertion sites were sequenced to find insertions that could potentially alter gene expression, and a preliminary test was performed with selected genes bearing transposable element (TE) insertions. RESULTS: Quantification of six families of Ty1-copia elements indicated different abundances among TE families and between flax accessions, which suggested diverse transpositional histories. SSAPs showed a high level of polymorphism in most of the evaluated retrotransposon families, with a trend towards higher levels of polymorphism in low-copy number families. Ty1-copia insertion polymorphisms among cultivars allowed a general distinction between oil and fiber types, and between spring and winter types, demonstrating their utility in diversity studies. Characterization of polymorphic insertions revealed an overwhelming association with genes, with insertions disrupting exons, introns or within 1 kb of coding regions. A preliminary test on the potential transcriptional disruption by TEs of four selected genes evaluated in three different tissues, showed one case of significant impact of the insertion on gene expression. CONCLUSIONS: We demonstrated that specific Ty1-copia families have been active since breeding commenced in flax. The retrotransposon-derived polymorphism can be used to separate flax types, and the close association of many insertions with genes defines a good source of potential mutations that could be associated with phenotypic changes, resulting in diversification processes.

Functional characterization of NAC55 transcription factor from oilseed rape (Brassica napus L.) as a novel transcriptional activator modulating reactive oxygen species accumulation and cell death.

NAC transcription factors (TFs) are plant-specific and play important roles in development, responses to biotic and abiotic cues and hormone signaling. So far, only a few NAC genes have been reported to regulate cell death. In this study, we identified and characterized a NAC55 gene isolated from oilseed rape (Brassica napus L.). BnaNAC55 responds to multiple stresses, including cold, heat, abscisic acid (ABA), jasmonic acid (JA) and a necrotrophic fungal pathogen Sclerotinia sclerotiorum. BnaNAC55 has transactivation activity and is located in the nucleus. BnaNAC55 is able to form homodimers in planta. Unlike ANAC055, full-length BnaNAC55, but not either the N-terminal NAC domain or C-terminal regulatory domain, induces ROS accumulation and hypersensitive response (HR)-like cell death when expressed both in oilseed rape protoplasts and Nicotiana benthamiana. Furthermore, BnaNAC55 expression causes obvious nuclear DNA fragmentation. Moreover, quantitative reverse transcription PCR (qRT-PCR) analysis identified that the expression levels of multiple genes regulating ROS production and scavenging, defense response as well as senescence are significantly induced. Using a dual luciferase reporter assay, we further confirm that BnaNAC55 could activate the expression of a few ROS and defense-related gene expression. Taken together, our work has identified a novel NAC TF from oilseed rape that modulates ROS accumulation and cell death.

Multiple polyploidy events in the early radiation of nodulating and nonnodulating legumes.

Unresolved questions about evolution of the large and diverse legume family include the timing of polyploidy (whole-genome duplication; WGDs) relative to the origin of the major lineages within the Fabaceae and to the origin of symbiotic nitrogen fixation. Previous work has established that a WGD affects most lineages in the Papilionoideae and occurred sometime after the divergence of the papilionoid and mimosoid clades, but the exact timing has been unknown. The history of WGD has also not been established for legume lineages outside the Papilionoideae. We investigated the presence and timing of WGDs in the legumes by querying thousands of phylogenetic trees constructed from transcriptome and genome data from 20 diverse legumes and 17 outgroup species. The timing of duplications in the gene trees indicates that the papilionoid WGD occurred in the common ancestor of all papilionoids. The earliest diverging lineages of the Papilionoideae include both nodulating taxa, such as the genistoids (e.g., lupin), dalbergioids (e.g., peanut), phaseoloids (e.g., beans), and galegoids (=Hologalegina, e.g., clovers), and clades with nonnodulating taxa including Xanthocercis and Cladrastis (evaluated in this study). We also found evidence for several independent WGDs near the base of other major legume lineages, including the Mimosoideae-Cassiinae-Caesalpinieae (MCC), Detarieae, and Cercideae clades. Nodulation is found in the MCC and papilionoid clades, both of which experienced ancestral WGDs. However, there are numerous nonnodulating lineages in both clades, making it unclear whether the phylogenetic distribution of nodulation is due to independent gains or a single origin followed by multiple losses.

Diversity of ABC transporter genes across the plant kingdom and their potential utility in biotechnology.

BACKGROUND: The ATP-binding cassette (ABC) transporter gene superfamily is ubiquitous among extant organisms and prominently represented in plants. ABC transporters act to transport compounds across cellular membranes and are involved in a diverse range of biological processes. Thus, the applicability to biotechnology is vast, including cancer resistance in humans, drug resistance among vertebrates, and herbicide and other xenobiotic resistance in plants. In addition, plants appear to harbor the highest diversity of ABC transporter genes compared with any other group of organisms. This study applied transcriptome analysis to survey the kingdom-wide ABC transporter diversity in plants and suggest biotechnology applications of this diversity. RESULTS: We utilized sequence similarity-based informatics techniques to infer the identity of ABC transporter gene candidates from 1295 phylogenetically-diverse plant transcriptomes. A total of 97,149 putative (approximately 25 % were full-length) ABC transporter gene members were identified; each RNA-Seq library (plant sample) had 88 ± 30 gene members. As expected, simpler organisms, such as algae, had fewer unique members than vascular land plants. Differences were also noted in the richness of certain ABC transporter subfamilies. Land plants had more unique ABCB, ABCC, and ABCG transporter gene members on average (p < 0.005), and green algae, red algae, and bryophytes had significantly more ABCF transporter gene members (p < 0.005). Ferns had significantly fewer ABCA transporter gene members than all other plant groups (p < 0.005). CONCLUSIONS: We present a transcriptomic overview of ABC transporter gene members across all major plant groups. An increase in the number of gene family members present in the ABCB, ABCC, and ABCD transporter subfamilies may indicate an expansion of the ABC transporter superfamily among green land plants, which include all crop species. The striking difference between the number of ABCA subfamily transporter gene members between ferns and other plant taxa is surprising and merits further investigation. Discussed is the potential exploitation of ABC transporters in plant biotechnology, with an emphasis on crops.

Wild soybean roots depend on specific transcription factors and oxidation reduction related genesin response to alkaline stress.

Soil alkalinity is an important environmental problem limiting agricultural productivity. Wild soybean (Glycine soja) shows strong alkaline stress tolerance, so it is an ideal plant candidate for studying the molecular mechanisms of alkaline tolerance and identifying alkaline stress-responsive genes. However, limited information is available about G. soja responses to alkaline stress on a genomic scale. Therefore, in the present study, we used RNA sequencing to compare transcript profiles of G. soja root responses to sodium bicarbonate (NaHCO3) at six time points, and a total of 68,138,478 pairs of clean reads were obtained using the Illumina GAIIX. Expression patterns of 46,404 G. soja genes were profiled in all six samples based on RNA-seq data using Cufflinks software. Then, t12 transcription factors from MYB, WRKY, NAC, bZIP, C2H2, HB, and TIFY families and 12 oxidation reduction related genes were chosen and verified to be induced in response to alkaline stress by using quantitative real-time polymerase chain reaction (qRT-PCR). The GO functional annotation analysis showed that besides "transcriptional regulation" and "oxidation reduction," these genes were involved in a variety of processes, such as "binding" and "response to stress." This is the first comprehensive transcriptome profiling analysis of wild soybean root under alkaline stress by RNA sequencing. Our results highlight changes in the gene expression patterns and identify a set of genes induced by NaHCO3 stress. These findings provide a base for the global analyses of G. soja alkaline stress tolerance mechanisms.

The evolutionary history of ferns inferred from 25 low-copy nuclear genes.

UNLABELLED: • PREMISE OF THE STUDY: Understanding fern (monilophyte) phylogeny and its evolutionary timescale is critical for broad investigations of the evolution of land plants, and for providing the point of comparison necessary for studying the evolution of the fern sister group, seed plants. Molecular phylogenetic investigations have revolutionized our understanding of fern phylogeny, however, to date, these studies have relied almost exclusively on plastid data.• METHODS: Here we take a curated phylogenomics approach to infer the first broad fern phylogeny from multiple nuclear loci, by combining broad taxon sampling (73 ferns and 12 outgroup species) with focused character sampling (25 loci comprising 35877 bp), along with rigorous alignment, orthology inference and model selection.• KEY RESULTS: Our phylogeny corroborates some earlier inferences and provides novel insights; in particular, we find strong support for Equisetales as sister to the rest of ferns, Marattiales as sister to leptosporangiate ferns, and Dennstaedtiaceae as sister to the eupolypods. Our divergence-time analyses reveal that divergences among the extant fern orders all occurred prior to ∼200 MYA. Finally, our species-tree inferences are congruent with analyses of concatenated data, but generally with lower support. Those cases where species-tree support values are higher than expected involve relationships that have been supported by smaller plastid datasets, suggesting that deep coalescence may be reducing support from the concatenated nuclear data.• CONCLUSIONS: Our study demonstrates the utility of a curated phylogenomics approach to inferring fern phylogeny, and highlights the need to consider underlying data characteristics, along with data quantity, in phylogenetic studies.