Transgenic Expression of Rubisco Accumulation Factor2 and Rubisco Subunits Increases Photosynthesis and Growth in Maize.

Carbon assimilation by Rubisco is often a limitation to photosynthesis and therefore plant productivity. We have previously shown that transgenic co-expression of the Rubisco large (LS) and small (SS) subunits along with an essential Rubisco assembly factor, Raf1, leads to faster growth, increased photosynthesis, and enhanced chilling tolerance in maize (Zea mays). Maize also requires Rubisco Accumulation Factor2 (Raf2) for full accumulation of Rubisco. Here we have analyzed transgenic maize lines with increased expression of Raf2 or Raf2 plus LS and SS. We show that increasing Raf2 expression alone had minor effects on photosynthesis, whereas expressing RAF2 with Rubisco subunits led to increased Rubisco content, more rapid carbon assimilation, and greater plant height, most notably in plants at least six weeks of age. The magnitude of the effects was similar to what was observed previously for expression of Raf1 together with Rubisco subunits. Taken together, this suggests that increasing the amount of either assembly factor with Rubisco subunits can independently enhance Rubisco abundance and some aspects of plant performance. These results could also imply either synergy, or a degree of functional redundancy for Raf1 and Raf2, the latter of whose precise role in Rubisco assembly is currently unknown.

Delivery mode impacts gut bacteriophage colonization during infancy.

Background: Cesarean section delivery is associated with altered early-life bacterial colonization and later adverse inflammatory and immune health outcomes. Although gut bacteriophages can alter gut microbiome composition and impact host immune responses, little is known about how delivery mode impacts bacteriophage colonization over time. To begin to address this we examined how delivery mode affected bacteriophage colonization over the first two years of life. Results: Shotgun metagenomic sequencing was conducted on 272 serial stool samples from 55 infants, collected at 1-2 days of life and 2, 6, 12 and 24 months. 33/55 (60%) infants were born by vaginal delivery. DNA viruses were identified, and by host inference, 94% of the viral sequences were found to be bacteriophages. Alpha diversity of the virome was increased in vaginally delivered infants compared to cesarean section delivered infants at 2 months (Shannon index, p=0.022). Beta diversity significantly differed by delivery mode at 2, 6, and 12 months when stratified by peripartum antibiotic use (Bray-Curtis dissimilarity, all p<0.05). Significant differentially abundant predicted bacteriophage hosts by delivery mode were seen at all time points. Moreover, there were differences in predicted bacteriophage functional gene abundances up to 24 months by delivery mode. Many of the functions considered to play a role in host response were increased in vaginal delivery. Conclusions: Clear differences in bacteriophage composition and function were seen by delivery mode over the first two years of life. Given that phages are known to affect host immune response, our results suggest that future investigation into how delivery mode may lead to adverse inflammatory outcomes should not only include bacterial microbial colonization but also the potential role of bacteriophages and transkingdom interactions.

Amphibian myelopoiesis.

Macrophage-lineage cells are indispensable to immunity and physiology of all vertebrates. Amongst these, amphibians represent a key stage in vertebrate evolution and are facing decimating population declines and extinctions, in large part due to emerging infectious agents. While recent studies indicate that macrophages and related innate immune cells are critically involved during these infections, much remains unknown regarding the ontogeny and functional differentiation of these cell types in amphibians. Accordingly, in this review we coalesce what has been established to date about amphibian blood cell development (hematopoiesis), the development of key amphibian innate immune cells (myelopoiesis) and the differentiation of amphibian macrophage subsets (monopoiesis). We explore the current understanding of designated sites of larval and adult hematopoiesis across distinct amphibian species and consider what mechanisms may lend to these species-specific adaptations. We discern the identified molecular mechanisms governing the functional differentiation of disparate amphibian (chiefly Xenopus laevis) macrophage subsets and describe what is known about the roles of these subsets during amphibian infections with intracellular pathogens. Macrophage lineage cells are at the heart of so many vertebrate physiological processes. Thus, garnering greater understanding of the mechanisms responsible for the ontogeny and functionality of these cells in amphibians will lend to a more comprehensive view of vertebrate evolution.

A Phylogenetic Framework to Simulate Synthetic Interspecies RNA-Seq Data.

Interspecies RNA-Seq datasets are increasingly common, and have the potential to answer new questions about the evolution of gene expression. Single-species differential expression analysis is now a well-studied problem that benefits from sound statistical methods. Extensive reviews on biological or synthetic datasets have provided the community with a clear picture on the relative performances of the available methods in various settings. However, synthetic dataset simulation tools are still missing in the interspecies gene expression context. In this work, we develop and implement a new simulation framework. This tool builds on both the RNA-Seq and the phylogenetic comparative methods literatures to generate realistic count datasets, while taking into account the phylogenetic relationships between the samples. We illustrate the usefulness of this new framework through a targeted simulation study, that reproduces the features of a recently published dataset, containing gene expression data in adult eye tissue across blind and sighted freshwater crayfish species. Using our simulated datasets, we perform a fair comparison of several approaches used for differential expression analysis. This benchmark reveals some of the strengths and weaknesses of both the classical and phylogenetic approaches for interspecies differential expression analysis, and allows for a reanalysis of the crayfish dataset. The tool has been integrated in the R package compcodeR, freely available on Bioconductor.

Removing systemic barriers to equity, diversity, and inclusion: Report of the 2019 Plant Science Research Network workshop "Inclusivity in the Plant Sciences".

A future in which scientific discoveries are valued and trusted by the general public cannot be achieved without greater inclusion and participation of diverse communities. To envision a path towards this future, in January 2019 a diverse group of researchers, educators, students, and administrators gathered to hear and share personal perspectives on equity, diversity, and inclusion (EDI) in the plant sciences. From these broad perspectives, the group developed strategies and identified tactics to facilitate and support EDI within and beyond the plant science community. The workshop leveraged scenario planning and the richness of its participants to develop recommendations aimed at promoting systemic change at the institutional level through the actions of scientific societies, universities, and individuals and through new funding models to support research and training. While these initiatives were formulated specifically for the plant science community, they can also serve as a model to advance EDI in other disciplines. The proposed actions are thematically broad, integrating into discovery, applied and translational science, requiring and embracing multidisciplinarity, and giving voice to previously unheard perspectives. We offer a vision of barrier-free access to participation in science, and a plant science community that reflects the diversity of our rapidly changing nation, and supports and invests in the training and well-being of all its members. The relevance and robustness of our recommendations has been tested by dramatic and global events since the workshop. The time to act upon them is now.

Genome-wide signatures of synergistic epistasis during parallel adaptation in a Baltic Sea copepod.

The role of epistasis in driving adaptation has remained an unresolved problem dating back to the Evolutionary Synthesis. In particular, whether epistatic interactions among genes could promote parallel evolution remains unexplored. To address this problem, we employ an Evolve and Resequence (E&R) experiment, using the copepod Eurytemora affinis, to elucidate the evolutionary genomic response to rapid salinity decline. Rapid declines in coastal salinity at high latitudes are a predicted consequence of global climate change. Based on time-resolved pooled whole-genome sequencing, we uncover a remarkably parallel, polygenic response across ten replicate selection lines, with 79.4% of selected alleles shared between lines by the tenth generation of natural selection. Using extensive computer simulations of our experiment conditions, we find that this polygenic parallelism is consistent with positive synergistic epistasis among alleles, far more so than other mechanisms tested. Our study provides experimental and theoretical support for a novel mechanism promoting repeatable polygenic adaptation, a phenomenon that may be common for selection on complex physiological traits.

The genus Creaserinus Hobbs, 1973 (Decapoda: Cambaridae) in Texas.

Presented is an account of the crayfish genus Creaserinus Hobbs, 1973 for Texas, based on materials gathered during a 13-year survey of the state. Home to Texas are six members of the genus, including C. hedgpethi (Hobbs, 1948) stat. rev., n. comb., which is resurrected from the synonymy of C. fodiens; and five species new to science described herein, including C. brevistylus n. sp., C. clausus n. sp., C. crenastylus n. sp., C. limulus n. sp., and C. trinensis n. sp. Collections of these species except for C. trinensis n. sp. were previously known and studied but ascribed to C. fodiens (Cottle, 1863), which is removed from the fauna of the state. Support for the taxonomic acts comes from genetics, morphology, distribution, life history, habitat, and syntopy. Accounts are provided for each species and include illustrations and information on distribution, color pattern, relationships, life history, ecology, size, variations, and crayfish associates. A key to the species in the state based on form I males is provided. Creaserinus limulus n. sp. is extraordinary in that a majority of its populations sampled have been composed mostly or entirely of females. Additions to the faunas of Texass neighboring states include C. clausus n. sp. (Louisiana), C. crenastylus n. sp. (Louisiana), and C. limulus n. sp. (Arkansas and Oklahoma).

Rubisco production in maize mesophyll cells through ectopic expression of subunits and chaperones.

C4 plants, such as maize, strictly compartmentalize Rubisco to bundle sheath chloroplasts. The molecular basis for the restriction of Rubisco from the more abundant mesophyll chloroplasts is not fully understood. Mesophyll chloroplasts transcribe the Rubisco large subunit gene and, when normally quiescent transcription of the nuclear Rubisco small subunit gene family is overcome by ectopic expression, mesophyll chloroplasts still do not accumulate measurable Rubisco. Here we show that a combination of five ubiquitin promoter-driven nuclear transgenes expressed in maize leads to mesophyll accumulation of assembled Rubisco. These encode the Rubisco large and small subunits, Rubisco assembly factors 1 and 2, and the assembly factor Bundle sheath defective 2. In these plants, Rubisco large subunit accumulates in mesophyll cells, and appears to be assembled into a holoenzyme capable of binding the substrate analog CABP (carboxyarabinitol bisphosphate). Isotope discrimination assays suggest, however, that mesophyll Rubisco is not participating in carbon assimilation in these plants, most probably due to a lack of the substrate ribulose 1,5-bisphosphate and/or Rubisco activase. Overall, this work defines a minimal set of Rubisco assembly factors in planta and may help lead to methods of regulating the C4 pathway.

CFM1, a member of the CRM-domain protein family, functions in chloroplast group II intron splicing in Setaria viridis.

The chloroplast RNA splicing and ribosome maturation (CRM) domain is a RNA-binding domain found in a plant-specific protein family whose characterized members play essential roles in splicing group I and group II introns in mitochondria and chloroplasts. Together, these proteins are required for splicing of the majority of the approximately 20 chloroplast introns in land plants. Here, we provide evidence from Setaria viridis and maize that an uncharacterized member of this family, CRM Family Member1 (CFM1), promotes the splicing of most of the introns that had not previously been shown to require a CRM domain protein. A Setaria mutant expressing mutated CFM1 was strongly disrupted in the splicing of three chloroplast tRNAs: trnI, trnV and trnA. Analyses by RNA gel blot and polysome association suggest that the tRNA deficiencies lead to compromised chloroplast protein synthesis and the observed whole-plant chlorotic phenotypes. Co-immunoprecipitation data demonstrate that the maize CFM1 ortholog is bound to introns whose splicing is disrupted in the cfm1 mutant. With these results, CRM domain proteins have been shown to promote the splicing of all but two of the introns found in angiosperm chloroplast genomes.

Setaria viridis chlorotic and seedling-lethal mutants define critical functions for chloroplast gene expression.

Deep insights into chloroplast biogenesis have been obtained by mutant analysis; however, in C4 plants a relevant mutant collection has only been developed and exploited for maize. Here, we report the initial characterization of an ethyl methyl sulfonate-induced mutant population for the C4 model Setaria viridis. Approximately 1000 M2 families were screened for the segregation of pale-green seedlings in the M3 generation, and a subset of these was identified to be deficient in post-transcriptional steps of chloroplast gene expression. Causative mutations were identified for three lines using deep sequencing-based bulked segregant analysis, and in one case confirmed by transgenic complementation. Using chloroplast RNA-sequencing and other molecular assays, we describe phenotypes of mutants deficient in PSRP7, a plastid-specific ribosomal protein, OTP86, an RNA editing factor, and cpPNP, the chloroplast isozyme of polynucleotide phosphorylase. The psrp mutant is globally defective in chloroplast translation, and has varying deficiencies in the accumulation of chloroplast-encoded proteins. The otp86 mutant, like its Arabidopsis counterpart, is specifically defective in editing of the rps14 mRNA; however, the conditional pale-green mutant phenotype contrasts with the normal growth of the Arabidopsis mutant. The pnp mutant exhibited multiple defects in 3' end maturation as well as other qualitative changes in the chloroplast RNA population. Overall, our collection opens the door to global analysis of photosynthesis and early seedling development in an emerging C4 model.

Hemiptera phylogenomic resources: Tree-based orthology prediction and conserved exon identification.

High-throughput sequencing of transcriptomes and targeted genomic regions are advancing our knowledge of The Tree of Life. Building phylogenies with regions of the genome requires 1-to-1 orthologue resources of genes and noncoding loci. One organismal group that has received little attention in this area is the Hemiptera, the fifth largest insect order represented by ~103,590 named species. Here, we present a set of 3,872 Hemiptera 1-to-1 orthogroups based on tree-based orthology inference of eight Hemiptera species with publicly available genome sequences. We also estimate a set of 406 orthologous exons with similar mRNA splice sites that can be used for Sanger sequencing and develop enrichment probes for targeted genome sequencing for phylogenomic inference. We show this novel set of orthologues is informative at the protein, coding sequence and exon molecular levels and provides robust branch support in both gene tree-species tree methods and concatenated sequence phylogenies. In addition, we demonstrate the utility of these loci to resolve relationships in whiteflies, Bemisia tabaci, a large species complex with few phylogenomic resources. Last, we compare our Hemiptera phylogeny with previously published phylogenies and other orthologue databases, while providing suggestions on further improvement to this phylogenomic resource.

Consumption of Diet Soda Sweetened with Sucralose and Acesulfame-Potassium Alters Inflammatory Transcriptome Pathways in Females with Overweight and Obesity.

SCOPE: Low-calorie sweetener (LCS) consumption is associated with metabolic disease in observational studies. However, physiologic mechanisms underlying LCS-induced metabolic impairments in humans are unclear. This study is aimed at identifying molecular pathways in adipose impacted by LCSs. METHODS AND RESULTS: Seven females with overweight or obesity, who did not report LCS use, consumed 12 ounces of diet soda containing sucralose and acesulfame-potassium (Ace-K) three times daily for 8 weeks. A subcutaneous adipose biopsy from the left abdomen and a fasting blood sample were collected at baseline and post-intervention. Global gene expression were assessed using RNA-sequencing followed by functional pathway analysis. No differences in circulating metabolic or inflammatory biomarkers were observed. However, ANOVA detected 828 differentially expressed annotated genes after diet soda consumption (p < 0.05), including transcripts for inflammatory cytokines. Fifty-eight of 140 canonical pathways represented in pathway analyses regulated inflammation, and several key upstream regulators of inflammation (e.g., TNF-alpha) were also represented. CONCLUSION: Consumption of diet soda with sucralose and Ace-K alters inflammatory transcriptomic pathways (e.g., NF-κB signaling) in subcutaneous adipose tissue but does not significantly alter circulating biomarkers. Findings highlight the need to examine molecular and metabolic effects of LCS exposure in a larger randomized control trial for a longer duration.

Plant Ribonuclease J: An Essential Player in Maintaining Chloroplast RNA Quality Control for Gene Expression.

RNA quality control is an indispensable but poorly understood process that enables organisms to distinguish functional RNAs from nonfunctional or inhibitory ones. In chloroplasts, whose gene expression activities are required for photosynthesis, retrograde signaling, and plant development, RNA quality control is of paramount importance, as transcription is relatively unregulated. The functional RNA population is distilled from this initial transcriptome by a combination of RNA-binding proteins and ribonucleases. One of the key enzymes is RNase J, a 5'→3' exoribonuclease and an endoribonuclease that has been shown to trim 5' RNA termini and eliminate deleterious antisense RNA. In the absence of RNase J, embryo development cannot be completed. Land plant RNase J contains a highly conserved C-terminal domain that is found in GT-1 DNA-binding transcription factors and is not present in its bacterial, archaeal, and algal counterparts. The GT-1 domain may confer specificity through DNA and/or RNA binding and/or protein-protein interactions and thus be an element in the mechanisms that identify target transcripts among diverse RNA populations. Further understanding of chloroplast RNA quality control relies on discovering how RNase J is regulated and how its specificity is imparted.

Increased Rubisco content in maize mitigates chilling stress and speeds recovery.

Many C4 plants, including maize, perform poorly under chilling conditions. This phenomenon has been linked in part to decreased Rubisco abundance at lower temperatures. An exception to this is chilling-tolerant Miscanthus, which is able to maintain Rubisco protein content under such conditions. The goal of this study was to investigate whether increasing Rubisco content in maize could improve performance during or following chilling stress. Here, we demonstrate that transgenic lines overexpressing Rubisco large and small subunits and the Rubisco assembly factor RAF1 (RAF1-LSSS), which have increased Rubisco content and growth under control conditions, maintain increased Rubisco content and growth during chilling stress. RAF1-LSSS plants exhibited 12% higher CO2 assimilation relative to nontransgenic controls under control growth conditions, and a 17% differential after 2 weeks of chilling stress, although assimilation rates of all genotypes were ~50% lower in chilling conditions. Chlorophyll fluorescence measurements showed RAF1-LSSS and WT plants had similar rates of photochemical quenching during chilling, suggesting Rubisco may not be the primary limiting factor that leads to poor performance in maize under chilling conditions. In contrast, RAF1-LSSS had improved photochemical quenching before and after chilling stress, suggesting that increased Rubisco may help plants recover faster from chilling conditions. Relatively increased leaf area, dry weight and plant height observed before chilling in RAF1-LSSS were also maintained during chilling. Together, these results demonstrate that an increase in Rubisco content allows maize plants to better cope with chilling stress and also improves their subsequent recovery, yet additional modifications are required to engineer chilling tolerance in maize.

Systematic sequencing of chloroplast transcript termini from Arabidopsis thaliana reveals >200 transcription initiation sites and the extensive imprints of RNA-binding proteins and secondary structures.

Chloroplast transcription requires numerous quality control steps to generate the complex but selective mixture of accumulating RNAs. To gain insight into how this RNA diversity is achieved and regulated, we systematically mapped transcript ends by developing a protocol called Terminome-seq. Using Arabidopsis thaliana as a model, we catalogued >215 primary 5' ends corresponding to transcription start sites (TSS), as well as 1628 processed 5' ends and 1299 3' ends. While most termini were found in intergenic regions, numerous abundant termini were also found within coding regions and introns, including several major TSS at unexpected locations. A consistent feature was the clustering of both 5' and 3' ends, contrasting with the prevailing description of discrete 5' termini, suggesting an imprecision of the transcription and/or RNA processing machinery. Numerous termini correlated with the extremities of small RNA footprints or predicted stem-loop structures, in agreement with the model of passive RNA protection. Terminome-seq was also implemented for pnp1-1, a mutant lacking the processing enzyme polynucleotide phosphorylase. Nearly 2000 termini were altered in pnp1-1, revealing a dominant role in shaping the transcriptome. In summary, Terminome-seq permits precise delineation of the roles and regulation of the many factors involved in organellar transcriptome quality control.

Overexpression of Rubisco subunits with RAF1 increases Rubisco content in maize.

Rubisco catalyses a rate-limiting step in photosynthesis and has long been a target for improvement due to its slow turnover rate. An alternative to modifying catalytic properties of Rubisco is to increase its abundance within C4 plant chloroplasts, which might increase activity and confer a higher carbon assimilation rate. Here, we overexpress the Rubisco large (LS) and small (SS) subunits with the Rubisco assembly chaperone RUBISCO ASSEMBLY FACTOR 1 (RAF1). While overexpression of LS and/or SS had no discernable impact on Rubisco content, addition of RAF1 overexpression resulted in a >30% increase in Rubisco content. Gas exchange showed a 15% increase in CO2 assimilation (ASAT) in UBI-LSSS-RAF1 transgenic plants, which correlated with increased fresh weight and in vitro Vcmax calculations. The divergence of Rubisco content and assimilation could be accounted for by the Rubisco activation state, which decreased up to 23%, suggesting that Rubisco activase may be limiting Vcmax, and impinging on the realization of photosynthetic potential from increased Rubisco content.

A Guide to the Chloroplast Transcriptome Analysis Using RNA-Seq.

Since its first use in plants in 2007, high-throughput RNA sequencing (RNA-Seq) has generated a vast amount of data for both model and nonmodel species. Organellar transcriptomes, however, are virtually always overlooked at the data analysis step. We therefore developed ChloroSeq, a bioinformatic pipeline aimed at facilitating the systematic analysis of chloroplast RNA metabolism, and we provide here a step-by-step user's manual. Following the alignment of quality-controlled data to the genome of interest, ChloroSeq measures genome expression level along with splicing and RNA editing efficiencies. When used in combination with the Tuxedo suite (TopHat and Cufflinks), ChloroSeq allows the simultaneous analysis of organellar and nuclear transcriptomes, opening the way to a better understanding of nucleus-organelle cross talk. We also describe the use of R commands to produce publication-quality figures based on ChloroSeq outputs. The effectiveness of the pipeline is illustrated through analysis of an RNA-Seq dataset covering the transition from growth to maturation to senescence of Arabidopsis thaliana leaves.

The Evolution of Gene Expression Underlying Vision Loss in Cave Animals.

Dissecting the evolutionary genetic processes underlying eye reduction and vision loss in obligate cave-dwelling organisms has been a long-standing challenge in evolutionary biology. Independent vision loss events in related subterranean organisms can provide critical insight into these processes as well as into the nature of convergent loss of complex traits. Advances in evolutionary developmental biology have illuminated the significant role of heritable gene expression variation in the evolution of new forms. Here, we analyze gene expression variation in adult eye tissue across the freshwater crayfish, representing four independent vision-loss events in caves. Species and individual expression patterns cluster by eye function rather than phylogeny, suggesting convergence in transcriptome evolution in independently blind animals. However, this clustering is not greater than what is observed in surface species with conserved eye function after accounting for phylogenetic expectations. Modeling expression evolution suggests that there is a common increase in evolutionary rates in the blind lineages, consistent with a relaxation of selective constraint maintaining optimal expression levels. This is evidence for a repeated loss of expression constraint in the transcriptomes of blind animals and that convergence occurs via a similar trajectory through genetic drift.

Phototransduction Gene Expression and Evolution in Cave and Surface Crayfishes.

In the absence of light in caves, animals have repeatedly evolved reduced eyes and visual systems. Whether the underlying genetic components remain intact in blind species remains unanswered across taxa. The freshwater crayfish have evolved to live in caves multiple times throughout their history; therefore, this system provides an opportunity to probe the genetic patterns and processes underlying repeated vision loss. Using transcriptomic data from the eyes of 14 species of cave and surface crayfishes, we identify the expression of 17 genes putatively related to visual phototransduction. We find a similarly complete repertoire of phototransduction gene families expressed in cave and surface species, but that the expression levels of those transcripts are consistently lower in cave species. We find statistical support for episodic positive selection, increased and decreased selection strength in caves, depending on the gene family. Analyses of gene expression evolution suggest convergent and possibly adaptive downregulation of these genes across eye-reduction events. Our results reveal a combination of evolutionary processes acting on the sequences and gene expression levels of vision-related genes underlying the loss of vision in caves.