Gene Assembly Research Articles

Full-length transcriptome revealed the accumulation of polyunsaturated fatty acids in developing seeds of Plukenetia volubilis.

BackgroundPlukenetia volubilis is cultivated as a valuable oilseed crop, and its mature seeds are rich in polyunsaturated fatty acids (FAs), which are widely used in food and pharmaceutical industries. Recently, next-generation sequencing (NGS) transcriptome studies in P. volubilis indicated that some candidate genes were involved in oil biosynthesis. The NGS were inaccuracies in assembly of some candidate genes, leading to unknown errors in date analyses. However, single molecular real-time (SMRT) sequencing can overcome these assembled errors. Unfortunately, this technique has not been reported in P. volubilis.MethodsThe total oil content of P. volubilis seed (PVS) was determined using Soxhlet extraction system. The FA composition were analyzed by gas chromatography. Combining PacBio SMRT and Illumina technologies, the transcriptome analysis of developing PVS was performed. Functional annotation and differential expression were performed by BLAST software (version 2.2.26) and RSEM software (version 1.2.31), respectively. The lncRNA-targeted transcripts were predicted in developing PVS using LncTar tool.ResultsBy Soxhlet extraction system, the oil content of superior plant-type (SPT) was 13.47% higher than that of inferior plant-type (IPT) at mature PVS. The most abundant FAs were C18:2 and C18:3, among which C18:3 content of SPT was 1.11-fold higher than that of IPT. Combined with PacBio and Illumina platform, 68,971 non-redundant genes were obtained, among which 7,823 long non-coding RNAs (lncRNAs) and 7,798 lncRNA-targeted genes were predicted. In developing seed, the expressions of 57 TFs showed a significantly positive correlation with oil contents, including WRI1-like1, LEC1-like1, and MYB44-like. Comparative analysis of expression profiles between SPT and IPT implied that orthologs of FAD3, PDCT, PDAT, and DAGT2 were possibly important for the accumulation of polyunsaturated FAs. Together, these results provide a reference for oil biosynthesis of P. volubilis and genetic improvement of oil plants.

A Cache Efficient One Hashing Blocked Bloom Filter (OHBB) for Random Strings and the K-mer Strings in DNA Sequence

Bloom filters are widely used in genome assembly, IoT applications and several network applications such as symmetric encryption algorithms, and blockchain applications owing to their advantages of fast querying, despite some false positives in querying the input elements. There are many research works carried out to improve both the insertion and querying speed or reduce the false-positive or reduce the storage requirements separately. However, the optimization of all the aforementioned parameters is quite challenging with the existing reported systems. This work proposes to simultaneously improve the insertion and querying speeds by introducing a Cache-efficient One-Hashing Blocked Bloom filter. The proposed method aims to reduce the number of memory accesses required for querying elements into one by splitting the memory into blocks where the block size is equal to the cache line size of the memory. In the proposed filter, each block has further been split into partitions where the size of each partition is the prime number. For insertion and query, one hash value is required, which yields different values when modulo divided with prime numbers. The speed is accelerated using simple hash functions where the hash function is called only once. The proposed method has been implemented and validated using random strings and symmetric K-mer datasets used in the gene assembly. The simulation results show that the proposed filter outperforms the Standard Bloom Filter in terms of the insertion and querying speed.

De novo transcriptome assembly and analysis of genes involved in desiccation tolerance in Grimmia pilifera

The anatomically simple plants that transition from the aquatic to the terrestrial have a certain mechanism to deal with the damage caused by the changing living environment. Grimmia pilifera is a type of resurrection plants that can survive a long period of desiccation. To understand the molecular mechanisms of desiccation tolerance, nine cDNA libraries were constructed in triplicate from mRNA obtained from G. pilifera for the 0 h, 6 h and 18 h desiccation treatment. RNA-Seq generated ∼ 666 million reads which were assembled into 135,850 unigenes. The differentially expressed genes (DEGs) were identified between different period of time of desiccation. Gene ontology analysis showed that a significant number of genes involved in water deprivation, chloroplast organization, xyloglucan metabolic process, regulation of signaling pathway. In addition, genes involved in osmotic stress, oxidative stress, accumulation of soluble matter, such as gene encoding antioxidant enzymes, trehalose synthase and channel protein, were up-regulated in response to drought stress. These results will be helpful for further studying on the molecular mechanisms of desiccation responses in G. pilifera.

The RAG1 Ubiquitin Ligase Domain Stimulates Recombination of TCRβ and TCRα Genes and Influences Development of αβ T Cell Lineages.

RAG1/RAG2 (RAG) endonuclease-mediated assembly of diverse lymphocyte Ag receptor genes by V(D)J recombination is critical for the development and immune function of T and B cells. The RAG1 protein contains a ubiquitin ligase domain that stabilizes RAG1 and stimulates RAG endonuclease activity in vitro. We report in this study that mice with a mutation that inactivates the Rag1 ubiquitin ligase in vitro exhibit decreased rearrangements and altered repertoires of TCRβ and TCRα genes in thymocytes and impaired thymocyte developmental transitions that require the assembly and selection of functional TCRβ and/or TCRα genes. These Rag1 mutant mice present diminished positive selection and superantigen-mediated negative selection of conventional αβ T cells, decreased genesis of invariant NK T lineage αβ T cells, and mature CD4+ αβ T cells with elevated autoimmune potential. Our findings reveal that the Rag1 ubiquitin ligase domain functions in vivo to stimulate TCRβ and TCRα gene recombination and influence differentiation of αβ T lineage cells, thereby establishing replete diversity of αβ TCRs and populations of αβ T cells while restraining generation of potentially autoreactive conventional αβ T cells.

A consensus protocol for the recovery of mercury methylation genes from metagenomes.

Mercury (Hg) methylation genes (hgcAB) mediate the formation of the toxic methylmercury and have been identified from diverse environments, including freshwater and marine ecosystems, Arctic permafrost, forest and paddy soils, coal-ash amended sediments, chlor-alkali plants discharges and geothermal springs. Here we present the first attempt at a standardized protocol for the detection, identification and quantification of hgc genes from metagenomes. Our Hg-cycling microorganisms in aquatic and terrestrial ecosystems (Hg-MATE) database, a catalogue of hgc genes, provides the most accurate information to date on the taxonomic identity and functional/metabolic attributes of microorganisms responsible for Hg methylation in the environment. Furthermore, we introduce "marky-coco", a ready-to-use bioinformatic pipeline based on de novo single-metagenome assembly, for easy and accurate characterization of hgc genes from environmental samples. We compared the recovery of hgc genes from environmental metagenomes using the marky-coco pipeline with an approach based on coassembly of multiple metagenomes. Our data show similar efficiency in both approaches for most environments except those with high diversity (i.e., paddy soils) for which a coassembly approach was preferred. Finally, we discuss the definition of true hgc genes and methods to normalize hgc gene counts from metagenomes.

Transcriptional and functional characterizations of multiple flagellin genes in spirochetes.

The flagellar filament is a helical propeller for bacterial locomotion. In external flagellates, the filaments are mostly homopolymers of a single flagellin protein. By contrast, the flagellar filaments of spirochetes are mostly heteropolymers of multiple flagellin proteins. This report seeks to investigate the role of multiple flagellin proteins using the oral spirochete Treponema denticola as a model. First, biochemical and genetic studies uncover that the flagellar filaments of T. denticola mainly comprise four proteins, FlaA, FlaB1, FlaB2, and FlaB3, in a defined stoichiometry. Second, transcriptional analyses reveal that the genes encoding these four proteins are regulated by two different transcriptional factors, sigma28 and sigma70. Third, loss‐of‐function studies demonstrate that each individual flagellin protein contributes to spirochete motility, but none of them is absolutely required. Last, we provide genetic and structural evidence that FlaA forms a “seam”‐like structure around the core and that deletion of individual flagellin protein alters the flagellar homeostasis. Collectively, these results demonstrate that T. denticola has evolved a unique mechanism to finely regulate its flagellar filament gene expression and assembly which renders the organelle with the right number, shape, strength, and structure for its distinct motility.

Monogenic TCRβ Assembly and Expression Are Paramount for Uniform Antigen Receptor Specificity of Individual αβ T Lymphocytes.

The ability of individual T and B cells to display Ag receptors of unique uniform specificity is the molecular basis of adaptive immunity. Most αβ T cells achieve uniform specificity by assembling in-frame genes on only one allelic copy of TCRβ and TCRα loci, while others prevent incorporation of TCRα protein from both alleles into TCRs. Analysis of mice expressing TCR proteins from a restricted combination of transgenes showed that TCR protein pairing restrictions achieve uniform specificity of cells expressing two types of TCRβ protein. However, whether this mechanism operates in the physiological context where each dual-TCRβ cell expresses one set of a vast number of different TCRβ proteins remains an open question, largely because there is a low, but significant, portion of cells carrying two in-frame TCRβ genes. To resolve this issue, we inactivated one allelic copy of the TCRα locus in a new mouse strain that assembles two in-frame TCRβ genes in an elevated fraction of cells. This genetic manipulation has no effect on the frequency of cells that display multiple types of αβ TCR, yet increases the representation of cells displaying TCRβ proteins that generate more highly expressed TCRs. Our data demonstrate that some TCRβ proteins exhibit differential functional pairing with TCRα proteins, but these restrictions have negligible contribution for ensuring uniform specificity of cells that express two types of TCRβ protein. Therefore, we conclude that mechanisms governing monogenic assembly and expression of TCRβ genes in individual cells are paramount for uniform specificity of αβ T lymphocytes.

DINOS: Data INspired Oligo Synthesis for DNA Data Storage

As interest in DNA-based information storage grows, the costs of synthesis have been identified as a key bottleneck. A potential direction is to tune synthesis for data. Data strands tend to be composed of a small set of recurring code word sequences, and they contain longer sequences of repeated data. To exploit these properties, we propose a new framework called DINOS. DINOS consists of three key parts: (i) The first is a hierarchical strand assembly algorithm, inspired by gene assembly techniques that can assemble arbitrary data strands from a small set of primitive blocks. (ii) The assembly algorithm relies on our novel formulation for how to construct primitive blocks, spanning a variety of useful configurations from a set of code words and overhangs. Each primitive block is a code word flanked by a pair of overhangs that are created by a cyclic pairing process that keeps the number of primitive blocks small. Using these primitive blocks, any data strand of arbitrary length can be assembled, theoretically. We show a minimal system for a binary code with as few as six primitive blocks, and we generalize our processes to support an arbitrary set of overhangs and code words. (iii) We exploit our hierarchical assembly approach to identify redundant sequences and coalesce the reactions that create them to make assembly more efficient. We evaluate DINOS and describe its key characteristics. For example, the number of reactions needed to make a strand can be reduced by increasing the number of overhangs or the number of code words, but increasing the number of overhangs offers a small advantage over increasing code words while requiring substantially fewer primitive blocks. However, density is improved more by increasing the number of code words. We also find that a simple redundancy coalescing technique is able to reduce reactions by 90.6% and 41.2% on average for decompressed and compressed data, respectively, even when the smallest data fragments being assembled are 16 bits. With a simple padding heuristic that finds even more redundancy, we can further decrease reactions for the same operating point up to 91.1% and 59% for decompressed and compressed data, respectively, on average. Our approach offers greater density by up to 80% over a prior general purpose gene assembly technique. Finally, in an analysis of synthesis costs in which we make 1 GB volume using de novo synthesis versus making only the primitive blocks with de novo synthesis and otherwise assembling using DINOS, we estimate DINOS as 10 5 × cheaper than de novo synthesis.

Engineering Saccharomyces cerevisiae for production of the capsaicinoid nonivamide

BackgroundCapsaicinoids are produced by plants in the Capsicum genus and are the main reason for the pungency of chili pepper fruits. They are strong agonists of TRPV1 (the transient receptor potential cation channel subfamily V member 1) and used as active ingredients in pharmaceuticals for the treatment of pain. The use of bioengineered microorganisms in a fermentation process may be an efficient route for their preparation, as well as for the discovery of (bio-)synthetic capsaicinoids with improved or novel bioactivities.ResultsSaccharomyces cerevisiae was engineered to over-express a selection of amide-forming N-acyltransferase and CoA-ligase enzyme cascades using a combinatorial gene assembly method, and was screened for nonivamide production from supplemented vanillylamine and nonanoic acid. Data from this work demonstrate that Tyramine N-hydroxycinnamoyl transferase from Capsicum annuum (CaAT) was most efficient for nonivamide formation in yeast, outcompeting the other candidates including AT3 (Pun1) from Capsicum spp. The CoA-ligase partner with highest activity from the ones evaluated here were from Petunia hybrida (PhCL) and Spingomonas sp. Ibu-2 (IpfF). A yeast strain expressing CaAT and IpfF produced 10.6 mg L−1 nonivamide in a controlled bioreactor setup, demonstrating nonivamide biosynthesis by S. cerevisiae for the first time.ConclusionsBaker’s yeast was engineered for production of nonivamide as a model capsaicinoid, by expressing N-acyltransferases and CoA-ligases of plant and bacterial origin. The constructed yeast platform holds potential for in vivo biocatalytic formation of capsaicinoids and could be a useful tool for the discovery of novel drugs.

Leaf and rhizome transcriptome assembly and expression analysis of genes involved in terpene biosynthesis in Dioscorea opposita

Leaf and rhizome transcriptome assembly and expression analysis of genes involved in terpene biosynthesis in Dioscorea opposita

Modular Cloning of MUC1 Recombinant Antibodies by Assembly of Synthetic Domain Genes

Antibody immunotherapy has emerged as a constructive way to fight cancer. New immunotherapies require the identification of a tumor‐specific antigen. The protein mucin 1 (MUC1) has long been recognized as an important specific antigen and a potential immunotherapy target. In healthy human tissues, MUC1 is responsible for protecting the epithelial lining. However, in cancer cells, MUC1 displays truncated glycosylation, revealing tumor‐specific epitopes. Recombinant antibodies are preferred for therapeutic development, and in some cases, it is desirable to produce different antibody formats, including intact IgG or Fab fragments. Using the MUC1 specific antibody 4H5 as a model, we have developed a modular cloning approach that can rapidly produce different antibody modalities to facilitate the rapid production of both IgG and Fab antibody formats. The 4H5 variable domains were PCR amplified to contain 20 base‐pair (bp) overlaps with the pcDNA3.1 vector on the 5’ end, and 20 bp overlaps to either the human CH1 or Ck domains on the 3’end of the genes. To produce the Fab, a His‐tag and a 20 base‐pair (bp) overlap to the vector was added to the 3’ end of a human CH1 domain gene sequence. For IgG production, human CH1‐CH2‐CH3 domains were PCR amplified, and the same vector overlap sequence was added to the 3’ end of the sequence. Finally, a 20 bp overlap to the vector was added to the 3’ end of a human Ck domain. The PCR products were purified, and the constructs were assembled using the NEBuilder®HiFi DNA Assembly Cloning Kit. DNA sequencing confirmed the successful assembly of the coding sequence, and the 4H5 Fab fragments and IgG were produced by transient transfection in CHO cells. The antibodies (Fab and IgG) were purified by affinity chromatography and binding to recombinant MUC1 antigens confirmed by ELISA. The cloning method presented here simultaneously produces multiple antibody formats. This approach will permit rapid assembly and screening of different antibody formats and isotypes for therapeutic antibody discovery and analysis.

XAF1 destabilizes estrogen receptor α through the assembly of a BRCA1-mediated destruction complex and promotes estrogen-induced apoptosis.

X-linked inhibitor of apoptosis-associated factor 1 (XAF1) is a pro-apoptotic tumor suppressor that is frequently inactivated in multiple human cancers. However, its candidacy as a suppressor in the pathogenesis of breast cancer remains undefined. Here, we report that XAF1 acts as a molecular switch in estrogen (E2)-mediated cell-fate decisions favoring apoptosis over cell proliferation. XAF1 promoter hypermethylation is observed predominantly in estrogen receptor α (ERα)-positive versus ERα-negative tumor cells and associated with attenuated apoptotic response to E2. XAF1 is activated by E2 through a G protein-coupled estrogen receptor-mediated non-genomic pathway and induces ERα degradation and apoptosis while it is repressed by ERα for E2 stimulation of cell proliferation. The XAF1-ERα mutual antagonism dictates the outcomes of E2 signaling and its alteration is linked to the development of E2-resistant tumors. Mechanistically, XAF1 destabilizes ERα through the assembly of breast cancer-associated gene 1 (BRCA1)-mediated destruction complex. XAF1 interacts with ERα and BRCA1 via the zinc finger (ZF) domains 5/6 and 4, respectively, and the mutants lacking either of these domains fail to drive ERα ubiquitination and apoptosis. E2-induced regression of XAF1+/+ tumors is abolished by XAF1 depletion while XAF1-/- tumors recover E2 response by XAF1 restoration. XAF1 and ERα expression show an inverse correlation in primary breast tumors, and XAF1 expression is associated with the overall survival of patients with ERα-positive but not ERα-negative cancer. Together, this study uncovers an important role for the XAF1-ERα antagonism as a linchpin to govern E2-mediated cell-fate decisions, illuminating the mechanistic consequence of XAF1 alteration in breast tumorigenesis.

Investigating the Prognostic Relevance of Tumor Immune Microenvironment and Immune Gene Assembly in Breast Carcinoma Subtypes.

Simple SummaryThe tumor immune microenvironment of different breast carcinoma (BC) subtypes and immune gene assembly of metastasizing and non-metastasizing HER2-negative BCs was analyzed. Examination of 309 Hematoxylin and Eosin-stained slides highlighted that the distribution of tumor infiltrating lymphocytes (TILs) from peritumoral, stromal and intratumoral regions varied greatly within all subtypes, with most tumors (66.01%) belonging to the immunologically “cold” group. Hormone receptor (HR) negative subtypes generally showed higher immune activity in all analyzed regions. No survival benefit was detected based on the spatial distribution of TILs. A lower CD4+/CD8+ ratio at the stromal internal tumor region indicated longer distant metastasis-free survival. When assessing immune gene expression between metastatic and non-metastatic BCs, the list of differentially expressed genes were non-identical across luminal and TNBCs, suggesting that these subtypes may use different mechanisms to bypass the immunological surveillance. Understanding these differences in the immune gene assembly may pave the way to the development of new immune-modulation therapies.We hypothesized that different BC subtypes are characterized by spatially distinct tumor immune microenvironment (TIME) and that immune gene assembly of metastatic (Met) and non-metastatic (Ctrl) BCs vary across subtypes. Peritumoral, stromal and intratumoral TIL was assessed on 309 BC cases. Hot, cold and immune-excluded groups were defined, and the prognostic role of this classification was assessed. CD4+/CD8+ positivity was analyzed in 75 cases in four systematically predefined tumor regions. Immune gene expression of Met and Ctrl HER2-negative BCs was compared by using NanoString nCounter technology. The amount of TIL infiltration varied greatly within all BC subtypes. Two-third of the cases were cold tumors with no significant survival difference compared to hot tumors. A lower CD4+/CD8+ ratio at the stromal internal tumor region was significantly associated with longer distant metastasis-free survival. The differentially expressed immune genes between Met and Ctrl varied across the studied BC subtypes with TNBC showing distinct features from the luminal subtypes. The TIME is characterized by a considerable heterogeneity; however, low level of TILs does not equate to disease progression. The differences in immune gene expression observed between Met and Ctrl breast carcinomas call attention to the important role of altered immune function in BC progression.

RoboMoClo: A Robotics-Assisted Modular Cloning Framework for Multiple Gene Assembly in Biofoundry.

Efficient and versatile DNA assembly frameworks have had an impact on promoting synthetic biology to build complex biological systems. To accelerate system development, laboratory automation (or biofoundry) provides an opportunity to construct organisms and DNA assemblies via computer-aided design. However, a modular cloning (MoClo) system for multiple DNA assemblies limits the biofoundry workflow in terms of simplicity and feasibility by preparing the number of cloning materials such as destination vectors prior to the automation process. Herein, we propose robot-assisted MoClo (RoboMoClo) to accelerate a synthetic biology project with multiple gene expressions at the biofoundry. The architecture of the RoboMoClo framework provides a hybrid strategy of hierarchical gene assembly and iterative gene assembly, and fewer destination vectors compared with other MoClo systems. An industrial bacterium, Corynebacterium glutamicum, was used as a model host for RoboMoClo. After building a biopart library (promoter and terminator; level 0) and evaluating its features (level 1), various transcriptional directions in multiple gene assemblies (level 2) were studied using the RoboMoClo vectors. Among the constructs, the convergent construct exhibited potential transcriptional interference through the collision of RNA polymerases. To study design of experiment-guided lycopene biosynthesis in C. glutamicum (levels 1, 2, and 3), the biofoundry-assisted multiple gene assembly was demonstrated as a proof-of-concept by constructing various sub-pathway units (level 2) and pathway units (level 3) for C. glutamicum. The RoboMoClo framework provides an improved MoClo toolkit for laboratory automation in a synthetic biology application.

S-Nitrosation of E3 Ubiquitin Ligase Complex Components Regulates Hormonal Signalings in Arabidopsis.

E3 ubiquitin ligases mediate the last step of the ubiquitination pathway in the ubiquitin-proteasome system (UPS). By targeting transcriptional regulators for their turnover, E3s play a crucial role in every aspect of plant biology. In plants, SKP1/CULLIN1/F-BOX PROTEIN (SCF)-type E3 ubiquitin ligases are essential for the perception and signaling of several key hormones including auxins and jasmonates (JAs). F-box proteins, TRANSPORT INHIBITOR RESPONSE 1 (TIR1) and CORONATINE INSENSITIVE 1 (COI1), bind directly transcriptional repressors AUXIN/INDOLE-3-ACETIC ACID (AUX/IAA) and JASMONATE ZIM-DOMAIN (JAZ) in auxin- and JAs-depending manner, respectively, which permits the perception of the hormones and transcriptional activation of signaling pathways. Redox modification of proteins mainly by S-nitrosation of cysteines (Cys) residues via nitric oxide (NO) has emerged as a valued regulatory mechanism in physiological processes requiring its rapid and versatile integration. Previously, we demonstrated that TIR1 and Arabidopsis thaliana SKP1 (ASK1) are targets of S-nitrosation, and these NO-dependent posttranslational modifications enhance protein-protein interactions and positively regulate SCFTIR1 complex assembly and expression of auxin response genes. In this work, we confirmed S-nitrosation of Cys140 in TIR1, which was associated in planta to auxin-dependent developmental and stress-associated responses. In addition, we provide evidence on the modulation of the SCFCOI1 complex by different S-nitrosation events. We demonstrated that S-nitrosation of ASK1 Cys118 enhanced ASK1-COI1 protein-protein interaction. Overexpression of non-nitrosable ask1 mutant protein impaired the activation of JA-responsive genes mediated by SCFCOI1 illustrating the functional relevance of this redox-mediated regulation in planta. In silico analysis positions COI1 as a promising S-nitrosation target, and demonstrated that plants treated with methyl JA (MeJA) or S-nitrosocysteine (NO-Cys, S-nitrosation agent) develop shared responses at a genome-wide level. The regulation of SCF components involved in hormonal perception by S-nitrosation may represent a key strategy to determine the precise time and site-dependent activation of each hormonal signaling pathway and highlights NO as a pivotal molecular player in these scenarios.

Slag: A program for seeded local assembly of genes in complex genomes.

Although finished genomes have become more common, there is still a need for assemblies of individual genes or chromosomal regions when only unassembled reads are available. slag (Seeded Local Assembly of Genes) fulfils this need by performing iterative local assembly based on cycles of matching‐read retrieval with blast and assembly with cap3, phrap, spades, canu or unicycler. The target sequence can be nucleotide or protein. Read fragmentation allows slag to use phrap or cap3 to assemble long reads at lower coverage (e.g., 5×) than is possible with canu or unicycler. In simple, nonrepetitive genomes, a slag assembly can cover a whole chromosome, but in complex genomes the growth of target‐matching contigs is limited as additional reads are consumed by consensus contigs consisting of repetitive elements. Apart from genomic complexity, contig length and correctness depend on read length and accuracy. With pyrosequencing or Illumina reads, slag‐assembled contigs are accurate enough to allow design of PCR primers, while contigs assembled from Oxford Nanopore or pre‐HiFi Pacific Biosciences long reads are generally only accurate enough to design baiting sequences for further targeted sequencing. In an application with real reads, slag successfully extended sequences for four wheat genes, which were verified by cloning and Sanger sequencing of overlapping amplicons. slag is a robust alternative to atram 2 for local assemblies, especially for read sets with less than 20× coverage. slag is freely available at https://github.com/cfcrane/SLAG.

Metagenomic Sequencing Reveals that the Assembly of Functional Genes and Taxa Varied Highly and Lacked Redundancy in the Earthworm Gut Compared with Soil under Vanadium Stress.

ABSTRACTExploring the ecological mechanism of microbial community assembly in soil and the earthworm gut in a vanadium polluted environment could help us understand the effects of vanadium stress on microbial diversity maintenance and function, as well as the mechanism of microbial mitigation of vanadium stress. Combining metagenomic sequencing and abundance distribution models, we explored the assembly of earthworm intestinal bacteria and native soil bacteria after 21 days of earthworm exposure to a gradient level of vanadate (0 to 300 mg kg−1) in soil. Stochastic processes dominated the assembly of both genes and taxa in earthworm gut and soil. Both the composition of taxa and functional genes in earthworm gut varied highly with the vanadium concentration, while in soil, only the taxa changed significantly, whereas the functional genes were relatively stable. The functional redundancy in soil, but not in the earthworm gut, was confirmed by a Mantel test and analysis of similarities (ANOSIM) test. In addition, vanadium detoxifying gene (VDG)-carrying taxa were more diverse but less abundant in soil than in the worm gut; and VDGs were more abundant in soil than in the worm gut. Their wider niche breadth indicated that VDG-carrying taxa were generalists in soil, in contrast to their role in the worm gut. These results suggested that earthworm intestinal and soil microbes adopted different strategies to counteract vanadium stress. The results provide new insights into the effects of soil vanadium stress on the assembly of earthworm gut and soil microbiota from both bacterial taxa and genetic function perspectives.IMPORTANCE Metagenomic sequencing revealed the variation of functional genes in the microbial community in soil and earthworm gut with increasing vanadium concentrations, which provided a new insight to explore the effect of vanadium stress on microbial community assembly from the perspective of functional genes. Our results reinforced the view that functional genes and taxa do not appear to have a simple corresponding relationship. Taxa are more sensitive compared with functional genes, suggesting the existence of bacterial functional redundancy in soil, but not in the earthworm gut. These observations indicate different assembly patterns of earthworm intestinal and soil bacteria under vanadium stress. Thus, it is important and necessary to include genetic functions to comprehensively understand microbial community assembly.

A new combinatorial megaplasmid library assembly method designed to screen for minimal pathways by using SCRaMbLE.

Human proteins expressed in yeast are common to enhance protein production while the expression of functional human pathways remain challenging. Here, we propose a simple and economical high-throughput gene assembly method to create a yeast megaplasmid library from human cDNA to screen for minimal human functional pathways. We introduced artificial promoters followed by symmetric loxP sites into the megaplasmids using Golden Gate assembly coupled with streptavidin-bead-based purification. The isolated high molecular weight, randomly assembled cDNA megaplasmid library may be useful for high-throughput directed evolution experiments and may be adapted for use in other model organisms.

Function and Molecular Mechanism of the DNA Damage Response in Immunity and Cancer Immunotherapy.

The DNA damage response (DDR) is an organized network of multiple interwoven components evolved to repair damaged DNA and maintain genome fidelity. Conceptually the DDR includes damage sensors, transducer kinases, and effectors to maintain genomic stability and accurate transmission of genetic information. We have recently gained a substantially improved molecular and mechanistic understanding of how DDR components are interconnected to inflammatory and immune responses to stress. DDR shapes both innate and adaptive immune pathways: (i) in the context of innate immunity, DDR components mainly enhance cytosolic DNA sensing and its downstream STimulator of INterferon Genes (STING)-dependent signaling; (ii) in the context of adaptive immunity, the DDR is needed for the assembly and diversification of antigen receptor genes that is requisite for T and B lymphocyte development. Imbalances between DNA damage and repair impair tissue homeostasis and lead to replication and transcription stress, mutation accumulation, and even cell death. These impacts from DDR defects can then drive tumorigenesis, secretion of inflammatory cytokines, and aberrant immune responses. Yet, DDR deficiency or inhibition can also directly enhance innate immune responses. Furthermore, DDR defects plus the higher mutation load in tumor cells synergistically produce primarily tumor-specific neoantigens, which are powerfully targeted in cancer immunotherapy by employing immune checkpoint inhibitors to amplify immune responses. Thus, elucidating DDR-immune response interplay may provide critical connections for harnessing immunomodulatory effects plus targeted inhibition to improve efficacy of radiation and chemotherapies, of immune checkpoint blockade, and of combined therapeutic strategies.

Resolving the Microalgal Gene Landscape at the Strain Level: a Novel Hybrid Transcriptome of Emiliania huxleyi CCMP3266.

Microalgae are key ecological players with a complex evolutionary history. Genomic diversity, in addition to limited availability of high-quality genomes, challenge studies that aim to elucidate molecular mechanisms underlying microalgal ecophysiology. Here, we present a novel and comprehensive transcriptomic hybrid approach to generate a reference for genetic analyses and resolve the microalgal gene landscape at the strain level. The approach is demonstrated for a strain of the coccolithophore microalga Emiliania huxleyi, which is a species complex with considerable genome variability. The investigated strain is commonly studied as a model for algal-bacterial interactions and was therefore sequenced in the presence of bacteria to elicit the expression of interaction-relevant genes. We applied complementary PacBio Iso-Seq full-length cDNA and poly(A)-independent Illumina total RNA sequencing, which resulted in a de novo-assembled, near-complete hybrid transcriptome. In particular, hybrid sequencing improved the reconstruction of long transcripts and increased the recovery of full-length transcript isoforms. To use the resulting hybrid transcriptome as a reference for genetic analyses, we demonstrate a method that collapses the transcriptome into a genome-like data set, termed "synthetic genome" (sGenome). We used the sGenome as a reference to visually confirm the robustness of the CCMP3266 gene assembly, to conduct differential gene expression analysis, and to characterize novel E. huxleyi genes. The newly identified genes contribute to our understanding of E. huxleyi genome diversification and are predicted to play a role in microbial interactions. Our transcriptomic toolkit can be implemented in various microalgae to facilitate mechanistic studies on microalgal diversity and ecology. IMPORTANCE Microalgae are key players in the ecology and biogeochemistry of our oceans. Efforts to implement genomic and transcriptomic tools in laboratory studies involving microalgae suffer from the lack of published genomes. In the case of coccolithophore microalgae, the problem has long been recognized; the model species Emiliania huxleyi is a species complex with genomes composed of a core and a large variable portion. To study the role of the variable portion in niche adaptation, and specifically in microbial interactions, strain-specific genetic information is required. Here, we present a novel transcriptomic hybrid approach, and generated strain-specific genome-like information. We demonstrate our approach on an E. huxleyi strain that is cocultivated with bacteria. By constructing a "synthetic genome," we generated comprehensive gene annotations that enabled accurate analyses of gene expression patterns. Importantly, we unveiled novel genes in the variable portion of E. huxleyi that play putative roles in microbial interactions.