Search

Genome editing tools, especially CRISPR/Cas9-based strategies, have transformed biomedical research and opened opportunities for developing curative treatments for genetic diseases. Despite rapid progress, low efficiency of targeted DNA integration and generation of undesired mutations represent major limitations for genome editing applications. Both issues arise from the interplay between the main DNA Double-Strand Break (DSB) repair pathways, Homology-Directed Repair (HDR), Non-Homologous End Joining (NHEJ), and Microhomology-Mediated End Joining (MMEJ). To improve efficiencies of targeted CRISPR-Cas9 genome editing, we screened a large compound library. This led to the discovery of AZD7648, a DNA-dependent protein kinase (DNA-PK) inhibitor and potent enhancer of CRISPR-Cas9-mediated integration. We demonstrated that AZD7648 increased HDR and decreased mutagenic NHEJ repair, thus resulting in improved performance of precise gene editing. Furthermore, we observed additional improvement of integration efficiency by impairing MMEJ repair through DNA polymerase ⊖ (Pol⊖) inhibition. Combined treatment with AZD7648 and Pol⊖ inhibitors (which we named 2iHDR) substantially increased precision of templated insertions, with efficiencies of up to 80%, and nearly no formation of undesired Insertion-Deletions (InDels). Importantly, 2iHDR also decreased Cas9-associated off-target activity, dramatically improving the performance and fidelity of CRISPR-Cas9 gene editing.

SUMMARYMost bacteria cells are protected by a peptidoglycan cell wall. Defining the chemical structure of the peptidoglycan has been instrumental to characterize cell wall associated proteins and to illuminate the mode of action of cell wall-acting antibiotics. However, a major roadblock for a comprehensive understanding of peptidoglycan homeostasis has been the lack of methods to conduct large-scale, systematic studies. Here we have developed and applied an innovative high throughput peptidoglycan analytical pipeline to analyze the entire non-essential, arrayed mutant library of Vibrio cholerae. The unprecedented breadth of these analyses revealed that peptidoglycan homeostasis is preserved by a large percentage of the genome organized in complex networks that functionally link peptidoglycan features with genetic determinants. As an example, we discovered a novel bifunctional penicillin-binding protein in V. cholerae. Collectively, genome-wide peptidoglycan profiling provides a fast, easy, and unbiased method for systematic identification of the genetic determinants of peptidoglycan synthesis and remodeling.

Leaf shape is a defining feature of how we recognize and classify plant species. Although there is extensive variation in leaf shape within many species, few studies have disentangled the underlying genetic architecture. We characterized the genetic architecture of leaf shape variation in Eurasian aspen (Populus tremula L.) by performing genome‐wide association study (GWAS) for physiognomy traits. To ascertain the roles of identified GWAS candidate genes within the leaf development transcriptional program, we generated RNA‐Seq data that we used to perform gene co‐expression network analyses from a developmental series, which is publicly available within the PlantGenIE resource. We additionally used existing gene expression measurements across the population to analyze GWAS candidate genes in the context of a population‐wide co‐expression network and to identify genes that were differentially expressed between groups of individuals with contrasting leaf shapes. These data were integrated with expression GWAS (eQTL) results to define a set of candidate genes associated with leaf shape variation. Our results identified no clear adaptive link to leaf shape variation and indicate that leaf shape traits are genetically complex, likely determined by numerous small‐effect variations in gene expression. Genes associated with shape variation were peripheral within the population‐wide co‐expression network, were not highly connected within the leaf development co‐expression network, and exhibited signatures of relaxed selection. As such, our results are consistent with the omnigenic model., We characterized the genetic architecture of leaf shape variation in Eurasian aspen (Populus tremula L.) by performing a genome‐wide association studies (GWAS) for physiognomy traits. Our results identified no clear adaptive link to leaf shape variation and indicate that leaf shape traits are genetically complex, likely determined by numerous small‐effect variations in gene expression. Genes associated with shape variation were peripheral within the population‐wide co‐expression network, were not highly connected within the leaf development co‐expression network, and exhibited signatures of relaxed selection.

Plant mitogenomes can be difficult to assemble because they are structurally dynamic and prone to intergenomic DNA transfers, leading to the unusual situation where an organelle genome is far outnumbered by its nuclear counterparts. As a result, comparative mitogenome studies are in their infancy and some key aspects of genome evolution are still known mainly from pre-genome, qualitative methods. To help address these limitations, we combined machine learning andin silicoenrichment of mitochondrial-like long reads to assemble the bacterial-sized mitogenome of Norway spruce (Pinaceae:Picea abies). We conducted comparative analyses of repeat abundance, intergenomic transfers, substitution and rearrangement rates, and estimated repeat-by-repeat homologous recombination rates. Prompted by our discovery of highly recombinogenic small repeats inP. abies, we assessed the genomic support for the prevailing hypothesis that intramolecular recombination is predominantly driven by repeat length, with larger repeats facilitating DNA exchange more readily. Overall, we found mixed support for this view: recombination dynamics were heterogeneous across vascular plants and highly active small repeats (ca. 200 bp) were present in about a third of studied mitogenomes. As in previous studies, we did not observe any robust relationships among commonly-studied genome attributes, but we identify variation in recombination rates as a underinvestigated source of plant mitogenome diversity.

Seasonal cues influence several aspects of the secondary growth of tree stems, including cambial activity, wood chemistry, and transition to latewood formation. We investigated seasonal changes in cambial activity, secondary cell wall formation, and tracheid cell death in woody tissues of Norway spruce (Picea abies) throughout one seasonal cycle. RNA sequencing was performed simultaneously in both the xylem and cambium/phloem tissues of the stem. Principal component analysis revealed gradual shifts in the transcriptomes that followed a chronological order throughout the season. A notable remodeling of the transcriptome was observed in the winter, with many genes having maximal expression during the coldest months of the year. A highly coexpressed set of monolignol biosynthesis genes showed high expression during the period of secondary cell wall formation as well as a second peak in midwinter. This midwinter peak in expression did not trigger lignin deposition, as determined by pyrolysis-gas chromatography/mass spectrometry. Coexpression consensus network analyses suggested the involvement of transcription factors belonging to the ASYMMETRIC LEAVES2/LATERAL ORGAN BOUNDARIES and MYELOBLASTOSIS-HOMEOBOX families in the seasonal control of secondary cell wall formation of tracheids. Interestingly, the lifetime of the latewood tracheids stretched beyond the winter dormancy period, correlating with a lack of cell death-related gene expression. Our transcriptomic analyses combined with phylogenetic and microscopic analyses also identified the cellulose and lignin biosynthetic genes and putative regulators for latewood formation and tracheid cell death in Norway spruce, providing a toolbox for further physiological and functional assays of these important phase transitions.

The rate of meiotic recombination is one of the central factors determining genome-wide levels of linkage disequilibrium which has important consequences for the efficiency of natural selection and for the dissection of quantitative traits. Here we present a new, high-resolution linkage map for Populus tremula that we use to anchor approximately two thirds of the P. tremula draft genome assembly on to the expected 19 chromosomes, providing us with the first chromosome-scale assembly for P. tremula (Table 2). We then use this resource to estimate variation in recombination rates across the P. tremula genome and compare these results to recombination rates based on linkage disequilibrium in a large number of unrelated individuals. We also assess how variation in recombination rates is associated with a number of genomic features, such as gene density, repeat density and methylation levels. We find that recombination rates obtained from the two methods largely agree, although the LD-based method identifies a number of genomic regions with very high recombination rates that the map-based method fails to detect. Linkage map and LD-based estimates of recombination rates are positively correlated and show similar correlations with other genomic features, showing that both methods can accurately infer recombination rate variation across the genome. Recombination rates are positively correlated with gene density and negatively correlated with repeat density and methylation levels, suggesting that recombination is largely directed toward gene regions in P. tremula.

The genome assembly of the European aspen Populus tremula proved difficult for a short-read based strategy due to high genomic variation. As a consequence, the fragmented sequence is impeding studies that benefit from highly contiguous data, particularly genome-wide association studies (GWAS) and comparative genomics. Here we present an updated assembly based on long-read sequences, optical mapping and genetic mapping. This assembly - henceforth referred to as Potra V2 - is assembled into 19 contiguous chromosomes which provides a powerful tool for future association studies. The genome sequence and any feature files are available from the PopGenIE resource.

The rate of meiotic recombination is one of the central factors determining levels of linkage disequilibrium and the efficiency of natural selection, and many organisms show a positive correlation between local rates of recombination and levels of nucleotide diversity indicating that linked selection is an important factor determining genome-wide levels of nucleotide diversity. Several methods for estimating recombination rates from segregating polymorphisms in natural populations have recently been developed. These methods have been extensively used in part because they are relatively simple to implement even in many non-model organisms, but also because they potentially offer higher resolution than traditional map-based methods. However, thorough comparisons of LD and map-based estimates of recombination are not readily available in plants. Here we present a new, high-resolution linkage map forPopulus tremulaand use this to estimate variation in recombination rates across theP. tremulagenome. We compare these results to recombination rates estimated based on linkage disequilibrium in a large number of unrelated individuals. We also assess how variation in recombination rates is associated with genomic features, such as gene density, repeat density and methylation levels. We find that recombination rates obtained from the two methods largely agree, although the LD-based method identify a number of genomic regions with very high recombination rates that the map-based method fail to detect. Linkage map and LD-based estimates of recombination rates are positively correlated and show similar correlations with other genomic features, showing that both methods can accurately infer recombination rate variation across the genome.

Current extinction rates are comparable to five prior mass extinctions in the earth’s history, and are strongly affected by human activities that have modified more than half of the earth’s terrestrial surface. Increasing human activity restricts animal movements and isolates formerly connected populations, a particular concern for the conservation of large carnivores, but no prior research has used high throughput sequencing in a standardized manner to examine genetic connectivity for multiple species of large carnivores and multiple ecosystems. Here, we used RAD SNP genotypes to test for differences in connectivity between multiple ecosystems for African wild dogs (Lycaon pictus) and lions (Panthera leo), and to test correlations between genetic distance, geographic distance and landscape resistance due to human activity. We found weaker connectivity, a stronger correlation between genetic distance and geographic distance, and a stronger correlation between genetic distance and landscape resistance for lions than for wild dogs, and propose a new hypothesis that adaptations to interspecific competition may help to explain differences in vulnerability to isolation by humans.

Significance We performed de novo, full-genome sequence analysis of two Populus species, North American quaking and Eurasian trembling aspen, that contain striking levels of genetic variation. Our results showed that positive and negative selection broadly affects patterns of genomic variation, but to varying degrees across coding and noncoding regions. The strength of selection and rates of sequence divergence were strongly related to differences in gene expression and coexpression network connectivity. These results highlight the importance of both positive and negative selection in shaping genome-wide levels of genetic variation in an obligately outcrossing, perennial plant. The resources we present establish aspens as a powerful study system enabling future studies for understanding the genomic determinants of adaptive evolution., 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).

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).

Norway spruce (Picea abies(L.) Karst.) is a conifer species of substanital economic and ecological importance. In common with most conifers, theP. abiesgenome is very large (∼20 Gbp) and contains a high fraction of repetitive DNA. The currentP. abiesgenome assembly (v1.0) covers approximately 60% of the total genome size but is highly fragmented, consisting of >10 million scaffolds. The genome annotation contains 66,632 gene models that are at least partially validated (www.congenie.org), however, the fragmented nature of the assembly means that there is currently little information available on how these genes are physically distributed over the 12P. abieschromosomes. By creating an ultra-dense genetic linkage map, we anchored and ordered scaffolds into linkage groups, which complements the fine-scale information available in assembly contigs. Our ultra-dense haploid consensus genetic map consists of 21,056 markers derived from 14,336 scaffolds that contain 17,079 gene models (25.6% of the validated gene models) that we have anchored to the 12 linkage groups. We used data from three independent component maps, as well as comparisons with previously publishedPiceamaps to evaluate the accuracy and marker ordering of the linkage groups. We demonstrate that approximately 3.8% of the anchored scaffolds and 1.6% of the gene models covered by the consensus map have likely assembly errors as they contain genetic markers that map to different regions within or between linkage groups. We further evaluate the utility of the genetic map for the conifer research community by using an independent data set of unrelated individuals to assess genome-wide variation in genetic diversity using the genomic regions anchored to linkage groups. The results show that our map is sufficiently dense to enable detailed evolutionary analyses across theP. abiesgenome.

Xylan is one of the main compounds determining wood properties in hardwood species. The xylan backbone is thought to be synthesized by a synthase complex comprising two members of the GT43 family. We downregulated all GT43 genes in hybrid aspen (Populus tremula × tremuloides) to understand their involvement in xylan biosynthesis. All three clades of the GT43 family were targeted for downregulation using RNA interference individually or in different combinations, either constitutively or specifically in developing wood. Simultaneous downregulation in developing wood of the B (IRX9) and C (IRX14) clades resulted in reduced xylan Xyl content relative to reducing end sequence, supporting their role in xylan backbone biosynthesis. This was accompanied by a higher lignocellulose saccharification efficiency. Unexpectedly, GT43 suppression in developing wood led to an overall growth stimulation, xylem cell wall thinning and a shift in cellulose orientation. Transcriptome profiling of these transgenic lines indicated that cell cycling was stimulated and secondary wall biosynthesis was repressed. We suggest that the reduced xylan elongation is sensed by the cell wall integrity surveying mechanism in developing wood. Our results show that wood-specific suppression of xylan-biosynthetic GT43 genes activates signaling responses, leading to increased growth and improved lignocellulose saccharification.

Apoplastic events such as monolignol oxidation and lignin polymerization are difficult to study in intact trees. To investigate the role of apoplastic hydrogen peroxide (H2O2) in gymnosperm phenolic metabolism, an extracellular lignin-forming cell culture of Norway spruce (Picea abies) was used as a research model. Scavenging of apoplastic H2O2 by potassium iodide repressed lignin formation, in line with peroxidases activating monolignols for lignin polymerization. Time-course analyses coupled to candidate substrate-product pair network propagation revealed differential accumulation of low-molecular-weight phenolics, including (glycosylated) oligolignols, (glycosylated) flavonoids, and proanthocyanidins, in lignin-forming and H2O2-scavenging cultures and supported that monolignols are oxidatively coupled not only in the cell wall but also in the cytoplasm, where they are coupled to other monolignols and proanthocyanidins. Dilignol glycoconjugates with reduced structures were found in the culture medium, suggesting that cells are able to transport glycosylated dilignols to the apoplast. Transcriptomic analyses revealed that scavenging of apoplastic H2O2 resulted in remodulation of the transcriptome, with reduced carbon flux into the shikimate pathway propagating down to monolignol biosynthesis. Aggregated coexpression network analysis identified candidate enzymes and transcription factors for monolignol oxidation and apoplastic H2O2 production in addition to potential H2O2 receptors. The results presented indicate that the redox state of the apoplast has a profound influence on cellular metabolism.

Plastid sequences are a cornerstone in plant systematic studies and key aspects of their evolution, such as uniparental inheritance and absent recombination, are often treated as axioms. While exceptions to these assumptions can profoundly influence evolutionary inference, detecting them can require extensive sampling, abundant sequence data, and detailed testing. Using advancements in high-throughput sequencing, we analyzed the whole plastomes of 65 accessions ofPicea,a genus of ~35 coniferous forest tree species, to test for deviations from canonical plastome evolution. Using complementary hypothesis and data-driven tests, we found evidence for chimeric plastomes generated by interspecific hybridization and recombination in the clade comprising Norway spruce (P. abies) and ten other species. Support for interspecific recombination remained after controlling for sequence saturation, positive selection, and potential alignment artifacts. These results reconcile previous conflicting plastid-based phylogenies and strengthen the mounting evidence of reticulate evolution inPicea.Given the relatively high frequency of hybridization and biparental plastid inheritance in plants, we suggest interspecific plastome recombination may be more widespread than currently appreciated and could underlie reported cases of discordant plastid phylogenies.

The advancement of high-throughput sequencing (HTS) technologies and the rapid development of numerous analysis algorithms and pipelines in this field has resulted in an unprecedentedly high demand for training scientists in HTS data analysis. Embarking on developing new training materials is challenging for many reasons. Trainers often do not have prior experience in preparing or delivering such materials and struggle to keep them up to date. A repository of curated HTS training materials would support trainers in materials preparation, reduce the duplication of effort by increasing the usage of existing materials, and allow for the sharing of teaching experience among the HTS trainers’ community. To achieve this, we have developed a strategy for materials’ curation and dissemination. Standards for describing training materials have been proposed and applied to the curation of existing materials. A Git repository has been set up for sharing annotated materials that can now be reused, modified, or incorporated into new courses. This repository uses Git; hence, it is decentralized and self-managed by the community and can be forked/built-upon by all users. The repository is accessible at http://bioinformatics.upsc.se/htmr., Author Summary In recent years, the advancement of high-throughput sequencing (HTS) and the rapid development of numerous analysis algorithms and pipelines in this field have resulted in an unprecedentedly high demand for training scientists in HTS data analysis. Generating effective training materials is time-consuming, and a large body of training materials on HTS data analysis has already been generated but is rarely shared among trainers. In this paper we provide guidelines to trainers for describing training materials to increase their reusability. The best practices standards proposed here have been used to annotate a collection of HTS training materials, which is now available to the trainers’ community in Git and discoverable through the ELIXIR and GOBLET portals. Efforts are now underway to utilize the strategy presented in this paper to annotate a wider collection of training materials and define a generic approach for the curation and dissemination of materials that should be adopted by existing training portals and new emerging initiatives.

With the rapid advancement of high throughput sequencing, large numbers of genetic markers can be readily and cheaply acquired, but most current software packages for genetic map construction cannot handle such dense input. Modern computer architectures and server farms represent untapped resources that can be used to enable higher marker densities to be processed in tractable time. Here we present a pipeline using a modified version of OneMap that parallelizes over bottleneck functions and achieves substantial speedups for producing a high density linkage map (N = 20,000). Using simulated data we show that the outcome is as accurate as the traditional pipeline. We further demonstrate that there is a direct relationship between the number of markers used and the level of deviation between true and estimated order, which in turn impacts the final size of a genetic map.

Background Evolutionary theory suggests that males and females may evolve sexually dimorphic phenotypic and biochemical traits concordant with each sex having different optimal strategies of resource investment to maximise reproductive success and fitness. Such sexual dimorphism would result in sex biased gene expression patterns in non-floral organs for autosomal genes associated with the control and development of such phenotypic traits. Results We examined morphological, biochemical and herbivory traits to test for sexually dimorphic resource allocation strategies within collections of sexually mature and immature Populus tremula (European aspen) trees. In addition we profiled gene expression in mature leaves of sexually mature wild trees using whole-genome oligonucleotide microarrays and RNA-Sequencing. Conclusions We found no evidence of sexual dimorphism or differential resource investment strategies between males and females in either sexually immature or mature trees. Similarly, single-gene differential expression and machine learning approaches revealed no evidence of large-scale sex biased gene expression. However, two significantly differentially expressed genes were identified from the RNA-Seq data, one of which is a robust diagnostic marker of sex in P. tremula. Electronic supplementary material The online version of this article (doi:10.1186/s12870-014-0276-5) contains supplementary material, which is available to authorized users.