Interspecific transfer of genetic information through polyploid bridges.

Hybridization blurs species boundaries and leads to intertwined lineages resulting in reticulate evolution. Polyploidy, the outcome of whole genome duplication (WGD), has more recently been implicated in promoting and facilitating hybridization between polyploid species, potentially leading to adaptive introgression. However, because polyploid lineages are usually ephemeral states in the evolutionary history of life it is unclear whether WGD-potentiated hybridization has any appreciable effect on their diploid counterparts. Here, we develop a model of cytotype dynamics within mixed-ploidy populations to demonstrate that polyploidy can in fact serve as a bridge for gene flow between diploid lineages, where introgression is fully or partially hampered by the species barrier. Polyploid bridges emerge in the presence of triploid organisms, which despite critically low levels of fitness, can still allow the transfer of alleles between diploid states of independently evolving mixed-ploidy species. Notably, while marked genetic divergence prevents polyploid-mediated interspecific gene flow, we show that increased recombination rates can offset these evolutionary constraints, allowing a more efficient sorting of alleles at higher-ploidy levels before introgression into diploid gene pools. Additionally, we derive an analytical approximation for the rate of gene flow at the tetraploid level necessary to supersede introgression between diploids with nonzero introgression rates, which is especially relevant for plant species complexes, where interspecific gene flow is ubiquitous. Altogether, our results illustrate the potential impact of polyploid bridges on the (re)distribution of genetic material across ecological communities during evolution, representing a potential force behind reticulation.

The duplication of genomes and genetic networks and its potential for evolutionary adaptation and survival during environmental turmoil.

The importance of whole-genome duplication (WGD) for evolution is controversial. Whereas some view WGD mainly as detrimental and an evolutionary dead end, there is growing evidence that polyploidization can help overcome environmental change, stressful conditions, or periods of extinction. However, despite much research, the mechanistic underpinnings of why and how polyploids might be able to outcompete or outlive nonpolyploids at times of environmental upheaval remain elusive, especially for autopolyploids, in which heterosis effects are limited. On the longer term, WGD might increase both mutational and environmental robustness due to redundancy and increased genetic variation, but on the short-or even immediate-term, selective advantages of WGDs are harder to explain. Here, by duplicating artificially generated Gene Regulatory Networks (GRNs), we show that duplicated GRNs-and thus duplicated genomes-show higher signal output variation than nonduplicated GRNs. This increased variation leads to niche expansion and can provide polyploid populations with substantial advantages to survive environmental turmoil. In contrast, under stable environments, GRNs might be maladaptive to changes, a phenomenon that is exacerbated in duplicated GRNs. We believe that these results provide insights into how genome duplication and (auto)polyploidy might help organisms to adapt quickly to novel conditions and to survive ecological uproar or even cataclysmic events.

Rhizogenic Agrobacterium protein RolB interacts with the TOPLESS repressor proteins to reprogram plant immunity and development.

Rhizogenic Agrobacterium strains comprise biotrophic pathogens that cause hairy root disease (HRD) on hydroponically grown Solanaceae and Cucurbitaceae crops, besides being widely explored agents for the creation of hairy root cultures for the sustainable production of plant-specialized metabolites. Hairy root formation is mediated through the expression of genes encoded on the T-DNA of the root-inducing (Ri) plasmid, of which several, including root oncogenic locus B (rolB), play a major role in hairy root development. Despite decades of research, the exact molecular function of the proteins encoded by the rol genes remains enigmatic. Here, by means of TurboID-mediated proximity labeling in tomato (Solanum lycopersicum) hairy roots, we identified the repressor proteins TOPLESS (TPL) and Novel Interactor of JAZ (NINJA) as direct interactors of RolB. Although these interactions allow RolB to act as a transcriptional repressor, our data hint at another in planta function of the RolB oncoprotein. Hence, by a series of plant bioassays, transcriptomic and DNA-binding site enrichment analyses, we conclude that RolB can mitigate the TPL functioning so that it leads to a specific and partial reprogramming of phytohormone signaling, immunity, growth, and developmental processes. Our data support a model in which RolB manipulates host transcription, at least in part, through interaction with TPL, to facilitate hairy root development. Thereby, we provide important mechanistic insights into this renowned oncoprotein in HRD.

Structure-function study of a Ca2+-independent metacaspase involved in lateral root emergence.

Metacaspases are part of an evolutionarily broad family of multifunctional cysteine proteases, involved in disease and normal development. As the structure-function relationship of metacaspases remains poorly understood, we solved the X-ray crystal structure of an Arabidopsis thaliana type II metacaspase (AtMCA-IIf) belonging to a particular subgroup not requiring calcium ions for activation. To study metacaspase activity in plants, we developed an in vitro chemical screen to identify small molecule metacaspase inhibitors and found several hits with a minimal thioxodihydropyrimidine-dione structure, of which some are specific AtMCA-IIf inhibitors. We provide mechanistic insight into the basis of inhibition by the TDP-containing compounds through molecular docking onto the AtMCA-IIf crystal structure. Finally, a TDP-containing compound (TDP6) effectively hampered lateral root emergence in vivo, probably through inhibition of metacaspases specifically expressed in the endodermal cells overlying developing lateral root primordia. In the future, the small compound inhibitors and crystal structure of AtMCA-IIf can be used to study metacaspases in other species, such as important human pathogens, including those causing neglected diseases.

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)6, Ge] and Gossypium stephensii [(AD)7, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)1, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates-including phenotypic differentiation, genetic isolation, and genetic convergence-that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

Efficient CRISPR-mediated base editing in Agrobacterium spp.

Agrobacterium spp. are important plant pathogens that are the causative agents of crown gall or hairy root disease. Their unique infection strategy depends on the delivery of part of their DNA to plant cells. Thanks to this capacity, these phytopathogens became a powerful and indispensable tool for plant genetic engineering and agricultural biotechnology. Although Agrobacterium spp. are standard tools for plant molecular biologists, current laboratory strains have remained unchanged for decades and functional gene analysis of Agrobacterium has been hampered by time-consuming mutation strategies. Here, we developed clustered regularly interspaced short palindromic repeats (CRISPR)-mediated base editing to enable the efficient introduction of targeted point mutations into the genomes of both Agrobacterium tumefaciens and Agrobacterium rhizogenes As an example, we generated EHA105 strains with loss-of-function mutations in recA, which were fully functional for maize (Zea mays) transformation and confirmed the importance of RolB and RolC for hairy root development by A. rhizogenes K599. Our method is highly effective in 9 of 10 colonies after transformation, with edits in at least 80% of the cells. The genomes of EHA105 and K599 were resequenced, and genome-wide off-target analysis was applied to investigate the edited strains after curing of the base editor plasmid. The off-targets present were characteristic of Cas9-independent off-targeting and point to TC motifs as activity hotspots of the cytidine deaminase used. We anticipate that CRISPR-mediated base editing is the start of "engineering the engineer," leading to improved Agrobacterium strains for more efficient plant transformation and gene editing.

Evolutionary history and pan-genome dynamics of strawberry (Fragaria spp.).

Strawberry (Fragaria spp.) has emerged as a model system for various fundamental and applied research in recent years. In total, the genomes of five different species have been sequenced over the past 10 y. Here, we report chromosome-scale reference genomes for five strawberry species, including three newly sequenced species' genomes, and genome resequencing data for 128 additional accessions to estimate the genetic diversity, structure, and demographic history of key Fragaria species. Our analyses obtained fully resolved and strongly supported phylogenies and divergence times for most diploid strawberry species. These analyses also uncovered a new diploid species (Fragaria emeiensis Jia J. Lei). Finally, we constructed a pan-genome for Fragaria and examined the evolutionary dynamics of gene families. Notably, we identified multiple independent single base mutations of the MYB10 gene associated with white pigmented fruit shared by different strawberry species. These reference genomes and datasets, combined with our phylogenetic estimates, should serve as a powerful comparative genomic platform and resource for future studies in strawberry.

The evolutionary origin and domestication history of goldfish (Carassius auratus).

Goldfish have been subjected to over 1,000 y of intensive domestication and selective breeding. In this report, we describe a high-quality goldfish genome (2n = 100), anchoring 95.75% of contigs into 50 pseudochromosomes. Comparative genomics enabled us to disentangle the two subgenomes that resulted from an ancient hybridization event. Resequencing 185 representative goldfish variants and 16 wild crucian carp revealed the origin of goldfish and identified genomic regions that have been shaped by selective sweeps linked to its domestication. Our comprehensive collection of goldfish varieties enabled us to associate genetic variations with a number of well-known anatomical features, including features that distinguish traditional goldfish clades. Additionally, we identified a tyrosine-protein kinase receptor as a candidate causal gene for the first well-known case of Mendelian inheritance in goldfish-the transparent mutant. The goldfish genome and diversity data offer unique resources to make goldfish a promising model for functional genomics, as well as domestication.

Mining for protein S-sulfenylation in Arabidopsis uncovers redox-sensitive sites.

Hydrogen peroxide (H2O2) is an important messenger molecule for diverse cellular processes. H2O2 oxidizes proteinaceous cysteinyl thiols to sulfenic acid, also known as S-sulfenylation, thereby affecting the protein conformation and functionality. Although many proteins have been identified as S-sulfenylation targets in plants, site-specific mapping and quantification remain largely unexplored. By means of a peptide-centric chemoproteomics approach, we mapped 1,537 S-sulfenylated sites on more than 1,000 proteins in Arabidopsis thaliana cells. Proteins involved in RNA homeostasis and metabolism were identified as hotspots for S-sulfenylation. Moreover, S-sulfenylation frequently occurred on cysteines located at catalytic sites of enzymes or on cysteines involved in metal binding, hinting at a direct mode of action for redox regulation. Comparison of human and Arabidopsis S-sulfenylation datasets provided 155 conserved S-sulfenylated cysteines, including Cys181 of the Arabidopsis MITOGEN-ACTIVATED PROTEIN KINASE4 (AtMAPK4) that corresponds to Cys161 in the human MAPK1, which has been identified previously as being S-sulfenylated. We show that, by replacing Cys181 of recombinant AtMAPK4 by a redox-insensitive serine residue, the kinase activity decreased, indicating the importance of this noncatalytic cysteine for the kinase mechanism. Altogether, we quantitatively mapped the S-sulfenylated cysteines in Arabidopsis cells under H2O2 stress and thereby generated a comprehensive view on the S-sulfenylation landscape that will facilitate downstream plant redox studies.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude.

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Histone 2B monoubiquitination complex integrates transcript elongation with RNA processing at circadian clock and flowering regulators.

HISTONE MONOUBIQUITINATION1 (HUB1) and its paralog HUB2 act in a conserved heterotetrameric complex in the chromatin-mediated transcriptional modulation of developmental programs, such as flowering time, dormancy, and the circadian clock. The KHD1 and SPEN3 proteins were identified as interactors of the HUB1 and HUB2 proteins with in vitro RNA-binding activity. Mutants in SPEN3 and KHD1 had reduced rosette and leaf areas. Strikingly, in spen3 mutants, the flowering time was slightly, but significantly, delayed, as opposed to the early flowering time in the hub1-4 mutant. The mutant phenotypes in biomass and flowering time suggested a deregulation of their respective regulatory genes CIRCADIAN CLOCK-ASSOCIATED1 (CCA1) and FLOWERING LOCUS C (FLC) that are known targets of the HUB1-mediated histone H2B monoubiquitination (H2Bub). Indeed, in the spen3-1 and hub1-4 mutants, the circadian clock period was shortened as observed by luciferase reporter assays, the levels of the CCA1α and CCA1ß splice forms were altered, and the CCA1 expression and H2Bub levels were reduced. In the spen3-1 mutant, the delay in flowering time was correlated with an enhanced FLC expression, possibly due to an increased distal versus proximal ratio of its antisense COOLAIR transcript. Together with transcriptomic and double-mutant analyses, our data revealed that the HUB1 interaction with SPEN3 links H2Bub during transcript elongation with pre-mRNA processing at CCA1 Furthermore, the presence of an intact HUB1 at the FLC is required for SPEN3 function in the formation of the FLC-derived antisense COOLAIR transcripts.

TNF-α inhibits glucocorticoid receptor-induced gene expression by reshaping the GR nuclear cofactor profile.

Glucocorticoid resistance (GCR) is defined as an unresponsiveness to the therapeutic effects, including the antiinflammatory ones of glucocorticoids (GCs) and their receptor, the glucocorticoid receptor (GR). It is a problem in the management of inflammatory diseases and can be congenital as well as acquired. The strong proinflammatory cytokine TNF-alpha (TNF) induces an acute form of GCR, not only in mice, but also in several cell lines: e.g., in the hepatoma cell line BWTG3, as evidenced by impaired Dexamethasone (Dex)-stimulated direct GR-dependent gene up- and down-regulation. We report that TNF has a significant and broad impact on this transcriptional performance of GR, but no impact on nuclear translocation, dimerization, or DNA binding capacity of GR. Proteome-wide proximity-mapping (BioID), however, revealed that the GR interactome was strongly modulated by TNF. One GR cofactor that interacted significantly less with the receptor under GCR conditions is p300. NFκB activation and p300 knockdown both reduced direct transcriptional output of GR whereas p300 overexpression and NFκB inhibition reverted TNF-induced GCR, which is in support of a cofactor reshuffle model. This hypothesis was supported by FRET studies. This mechanism of GCR opens avenues for therapeutic interventions in GCR diseases.

Genetic mapping of species differences via in vitro crosses in mouse embryonic stem cells.

Discovering the genetic changes underlying species differences is a central goal in evolutionary genetics. However, hybrid crosses between species in mammals often suffer from hybrid sterility, greatly complicating genetic mapping of trait variation across species. Here, we describe a simple, robust, and transgene-free technique to generate "in vitro crosses" in hybrid mouse embryonic stem (ES) cells by inducing random mitotic cross-overs with the drug ML216, which inhibits the DNA helicase Bloom syndrome (BLM). Starting with an interspecific F1 hybrid ES cell line between the Mus musculus laboratory mouse and Mus spretus (∼1.5 million years of divergence), we mapped the genetic basis of drug resistance to the antimetabolite tioguanine to a single region containing hypoxanthine-guanine phosphoribosyltransferase (Hprt) in as few as 21 d through "flow mapping" by coupling in vitro crosses with fluorescence-activated cell sorting (FACS). We also show how our platform can enable direct study of developmental variation by rederiving embryos with contribution from the recombinant ES cell lines. We demonstrate how in vitro crosses can overcome major bottlenecks in mouse complex trait genetics and address fundamental questions in evolutionary biology that are otherwise intractable through traditional breeding due to high cost, small litter sizes, and/or hybrid sterility. In doing so, we describe an experimental platform toward studying evolutionary systems biology in mouse and potentially in human and other mammals, including cross-species hybrids.

Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen.

The Populus genus is one of the major plant model systems, but genomic resources have thus far primarily been available for poplar species, and primarily Populus trichocarpa (Torr. & Gray), which was the first tree with a whole-genome assembly. To further advance evolutionary and functional genomic analyses in Populus, we produced genome assemblies and population genetics resources of two aspen species, Populus tremula L. and Populus tremuloides Michx. The two aspen species have distributions spanning the Northern Hemisphere, where they are keystone species supporting a wide variety of dependent communities and produce a diverse array of secondary metabolites. Our analyses show that the two aspens share a similar genome structure and a highly conserved gene content with P. trichocarpa but display substantially higher levels of heterozygosity. Based on population resequencing data, we observed widespread positive and negative selection acting on both coding and noncoding regions. Furthermore, patterns of genetic diversity and molecular evolution in aspen are influenced by a number of features, such as expression level, coexpression network connectivity, and regulatory variation. To maximize the community utility of these resources, we have integrated all presented data within the PopGenIE web resource (PopGenIE.org).

Identification of genes from the general phenylpropanoid and monolignol-specific metabolism in two sugarcane lignin-contrasting genotypes.

The phenylpropanoid pathway is an important route of secondary metabolism involved in the synthesis of different phenolic compounds such as phenylpropenes, anthocyanins, stilbenoids, flavonoids, and monolignols. The flux toward monolignol biosynthesis through the phenylpropanoid pathway is controlled by specific genes from at least ten families. Lignin polymer is one of the major components of the plant cell wall and is mainly responsible for recalcitrance to saccharification in ethanol production from lignocellulosic biomass. Here, we identified and characterized sugarcane candidate genes from the general phenylpropanoid and monolignol-specific metabolism through a search of the sugarcane EST databases, phylogenetic analysis, a search for conserved amino acid residues important for enzymatic function, and analysis of expression patterns during culm development in two lignin-contrasting genotypes. Of these genes, 15 were cloned and, when available, their loci were identified using the recently released sugarcane genomes from Saccharum hybrid R570 and Saccharum spontaneum cultivars. Our analysis points out that ShPAL1, ShPAL2, ShC4H4, Sh4CL1, ShHCT1, ShC3H1, ShC3H2, ShCCoAOMT1, ShCOMT1, ShF5H1, ShCCR1, ShCAD2, and ShCAD7 are strong candidates to be bona fide lignin biosynthesis genes. Together, the results provide information about the candidate genes involved in monolignol biosynthesis in sugarcane and may provide useful information for further molecular genetic studies in sugarcane.

Complete overview of protein-inactivating sequence variations in 36 sequenced mouse inbred strains.

Mouse inbred strains remain essential in science. We have analyzed the publicly available genome sequences of 36 popular inbred strains and provide lists for each strain of protein-coding genes that acquired sequence variations that cause premature STOP codons, loss of STOP codons and single nucleotide polymorphisms, and short in-frame insertions and deletions. Our data give an overview of predicted defective proteins, including predicted impact scores, of all these strains compared with the reference mouse genome of C57BL/6J. These data can also be retrieved via a searchable website (mousepost.be) and allow a global, better interpretation of genetic background effects and a source of naturally defective alleles in these 36 sequenced classical and high-priority mouse inbred strains.

Network-based integration of systems genetics data reveals pathways associated with lignocellulosic biomass accumulation and processing.

As a consequence of their remarkable adaptability, fast growth, and superior wood properties, eucalypt tree plantations have emerged as key renewable feedstocks (over 20 million ha globally) for the production of pulp, paper, bioenergy, and other lignocellulosic products. However, most biomass properties such as growth, wood density, and wood chemistry are complex traits that are hard to improve in long-lived perennials. Systems genetics, a process of harnessing multiple levels of component trait information (e.g., transcript, protein, and metabolite variation) in populations that vary in complex traits, has proven effective for dissecting the genetics and biology of such traits. We have applied a network-based data integration (NBDI) method for a systems-level analysis of genes, processes and pathways underlying biomass and bioenergy-related traits using a segregating Eucalyptus hybrid population. We show that the integrative approach can link biologically meaningful sets of genes to complex traits and at the same time reveal the molecular basis of trait variation. Gene sets identified for related woody biomass traits were found to share regulatory loci, cluster in network neighborhoods, and exhibit enrichment for molecular functions such as xylan metabolism and cell wall development. These findings offer a framework for identifying the molecular underpinnings of complex biomass and bioprocessing-related traits. A more thorough understanding of the molecular basis of plant biomass traits should provide additional opportunities for the establishment of a sustainable bio-based economy.

Genome of wild olive and the evolution of oil biosynthesis.

Here we present the genome sequence and annotation of the wild olive tree (Olea europaea var. sylvestris), called oleaster, which is considered an ancestor of cultivated olive trees. More than 50,000 protein-coding genes were predicted, a majority of which could be anchored to 23 pseudochromosomes obtained through a newly constructed genetic map. The oleaster genome contains signatures of two Oleaceae lineage-specific paleopolyploidy events, dated at ∼28 and ∼59 Mya. These events contributed to the expansion and neofunctionalization of genes and gene families that play important roles in oil biosynthesis. The functional divergence of oil biosynthesis pathway genes, such as FAD2, SACPD, EAR, and ACPTE, following duplication, has been responsible for the differential accumulation of oleic and linoleic acids produced in olive compared with sesame, a closely related oil crop. Duplicated oleaster FAD2 genes are regulated by an siRNA derived from a transposable element-rich region, leading to suppressed levels of FAD2 gene expression. Additionally, neofunctionalization of members of the SACPD gene family has led to increased expression of SACPD2, 3, 5, and 7, consequently resulting in an increased desaturation of steric acid. Taken together, decreased FAD2 expression and increased SACPD expression likely explain the accumulation of exceptionally high levels of oleic acid in olive. The oleaster genome thus provides important insights into the evolution of oil biosynthesis and will be a valuable resource for oil crop genomics.

Grass pea (Lathyrus sativus L.): orphan crop, nutraceutical or just plain food?

MAIN CONCLUSION: Although grass pea is an environmentally successful robust legume with major traits of interest for food and nutrition security, the genetic potential of this orphan crop has long been neglected. Grass pea (Lathyrus sativus L.) is a Neolithic plant that has survived millennia of cultivation and has spread over three continents. It is a robust legume crop that is considered one of the most resilient to climate changes and to be survival food during drought-triggered famines. The hardy penetrating root system allows the cultivation of grass pea in various soil types, including marginal ones. As an efficient nitrogen fixer, it meets its own nitrogen requirements and positively benefits subsequent crops. However, already in ancient India and Greece, overconsumption of the seeds and a crippling neurological disorder, later coined neurolathyrism, had been linked. Overemphasis of their suspected toxic properties has led to disregard the plant's exceptionally positive agronomic properties and dietary advantages. In normal socio-economic and environmental situations, in which grass pea is part of a balanced diet, neurolathyrism is virtually non-existent. The etiology of neurolathyrism has been oversimplified and the deficiency in methionine in the diet has been overlooked. In view of the global climate change, this very adaptable and nutritious orphan crop deserves more attention. Grass pea can become a wonder crop if the double stigma on its reputation as a toxic plant and as food of the poor can be disregarded. Additionally, recent research has exposed the potential of grass pea as a health-promoting nutraceutical. Development of varieties with an improved balance in essential amino acids and diet may be relevant to enhance the nutritional value without jeopardizing the multiple stress tolerance of this promising crop.

Efficient analysis of mouse genome sequences reveal many nonsense variants.

Genetic polymorphisms in coding genes play an important role when using mouse inbred strains as research models. They have been shown to influence research results, explain phenotypical differences between inbred strains, and increase the amount of interesting gene variants present in the many available inbred lines. SPRET/Ei is an inbred strain derived from Mus spretus that has ∼1% sequence difference with the C57BL/6J reference genome. We obtained a listing of all SNPs and insertions/deletions (indels) present in SPRET/Ei from the Mouse Genomes Project (Wellcome Trust Sanger Institute) and processed these data to obtain an overview of all transcripts having nonsynonymous coding sequence variants. We identified 8,883 unique variants affecting 10,096 different transcripts from 6,328 protein-coding genes, which is about 28% of all coding genes. Because only a subset of these variants results in drastic changes in proteins, we focused on variations that are nonsense mutations that ultimately resulted in a gain of a stop codon. These genes were identified by in silico changing the C57BL/6J coding sequences to the SPRET/Ei sequences, converting them to amino acid (AA) sequences, and comparing the AA sequences. All variants and transcripts affected were also stored in a database, which can be browsed using a SPRET/Ei M. spretus variants web tool (www.spretus.org), including a manual. We validated the tool by demonstrating the loss of function of three proteins predicted to be severely truncated, namely Fas, IRAK2, and IFNγR1.