Rapid amplification of cDNA ends (RACE)

B-Cell Receptor Repertoire: Recent Advances in Autoimmune Diseases

In the field of contemporary medicine, autoimmune diseases (AIDs) are a prevalent and debilitating group of illnesses. However, they present extensive and profound challenges in terms of etiology, pathogenesis, and treatment. A major reason for this is the elusive pathophysiological mechanisms driving disease onset. Increasing evidence suggests the indispensable role of B cells in the pathogenesis of autoimmune diseases. Interestingly, B-cell receptor (BCR) repertoires in autoimmune diseases display a distinct skewing that can provide insights into disease pathogenesis. Over the past few years, advances in high-throughput sequencing have provided powerful tools for analyzing B-cell repertoire to understand the mechanisms during the period of B-cell immune response. In this paper, we have provided an overview of the mechanisms and analytical methods for generating BCR repertoire diversity and summarize the latest research progress on BCR repertoire in autoimmune diseases, including systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), primary Sjögren's syndrome (pSS), multiple sclerosis (MS), and type 1 diabetes (T1D). Overall, B-cell repertoire analysis is a potent tool to understand the involvement of B cells in autoimmune diseases, facilitating the creation of innovative therapeutic strategies targeting specific B-cell clones or subsets.

The immunoglobulin A isotype of the Arabian camel (Camelus dromedarius) preserves the dualistic structure of unconventional single-domain and canonical heavy chains

Introduction

The evolution of adaptive immunity in Camelidae resulted in the concurrent expression of classic heterotetrameric and unconventional homodimeric heavy chain-only IgG antibodies. Heavy chain-only IgG bears a single variable domain and lacks the constant heavy (CH) γ1 domain required for pairing with the light chain. It has not been reported whether this distinctive feature of IgG is also observed in the IgA isotype.

Methods

Gene-specific primers were used to generate an IgA heavy chain cDNA library derived from RNA extracted from the dromedary's third eyelid where isolated lymphoid follicles and plasma cells abound at inductive and effector sites, respectively.

Results

Majority of the cDNA clones revealed hallmarks of heavy chain-only antibodies, i.e. camelid-specific amino acid substitutions in framework region 1 and 2, broad length distribution of complementarity determining region 3, and the absence of the CHα1 domain. In a few clones, however, the cDNA of the canonical IgA heavy chain was amplified which included the CHα1 domain, analogous to CHγ1 domain in IgG1 subclass. Moreover, we noticed a short, proline-rich hinge, and, at the N-terminal end of the CHα3 domain, a unique, camelid-specific pentapeptide of undetermined function, designated as the inter-α region. Immunoblots using rabbit anti-camel IgA antibodies raised against CHα2 and CHα3 domains as well as the inter-α region revealed the expression of a ~52 kDa and a ~60 kDa IgA species, corresponding to unconventional and canonical IgA heavy chain, respectively, in the third eyelid, trachea, small and large intestine. In contrast, the leporine anti-CHα1 antibody detected canonical, but not unconventional IgA heavy chain, in all the examined tissues, milk, and serum, in addition to another hitherto unexplored species of ~45 kDa in milk and serum. Immunohistology using anti-CHα domain antibodies confirmed the expression of both variants of IgA heavy chains in plasma cells in the third eyelid's lacrimal gland, conjunctiva, tracheal and intestinal mucosa.

Conclusion

We found that in the dromedary, the IgA isotype has expanded the immunoglobulin repertoire by co-expressing unconventional and canonical IgA heavy chains, comparable to the IgG class, thus underscoring the crucial role of heavy chain-only antibodies not only in circulation but also at the mucosal frontiers.

Human Satellite 1A analysis provides evidence of pericentromeric transcription

BACKGROUND

Pericentromeric regions of human chromosomes are composed of tandem-repeated and highly organized sequences named satellite DNAs. Human classical satellite DNAs are classified into three families named HSat1, HSat2, and HSat3, which have historically posed a challenge for the assembly of the human reference genome where they are misrepresented due to their repetitive nature. Although being known for a long time as the most AT-rich fraction of the human genome, classical satellite HSat1A has been disregarded in genomic and transcriptional studies, falling behind other human satellites in terms of functional knowledge. Here, we aim to characterize and provide an understanding on the biological relevance of HSat1A.

RESULTS

The path followed herein trails with HSat1A isolation and cloning, followed by in silico analysis. Monomer copy number and expression data was obtained in a wide variety of human cell lines, with greatly varying profiles in tumoral/non-tumoral samples. HSat1A was mapped in human chromosomes and applied in in situ transcriptional assays. Additionally, it was possible to observe the nuclear organization of HSat1A transcripts and further characterize them by 3' RACE-Seq. Size-varying polyadenylated HSat1A transcripts were detected, which possibly accounts for the intricate regulation of alternative polyadenylation.

CONCLUSION

As far as we know, this work pioneers HSat1A transcription studies. With the emergence of new human genome assemblies, acrocentric pericentromeres are becoming relevant characters in disease and other biological contexts. HSat1A sequences and associated noncoding RNAs will most certainly prove significant in the future of HSat research.

High-throughput optofluidic screening of single B cells identifies novel cross-reactive antibodies as inhibitors of uPAR with antibody-dependent effector functions

The urokinase-type plasminogen activator receptor (uPAR) is an essential regulator for cell signaling in tumor cell proliferation, adhesion, and metastasis. The ubiquitous nature of uPAR in many aggressive cancer types makes uPAR an attractive target for immunotherapy. Here, we present a rapid and successful workflow for developing cross-reactive anti-uPAR recombinant antibodies (rAbs) using high-throughput optofluidic screening of single B-cells from human uPAR-immunized mice. A total of 80 human and cynomolgus uPAR cross-reactive plasma cells were identified, and selected mouse VH/VL domains were linked to the trastuzumab (Herceptin®) constant domains for the expression of mouse-human chimeric antibodies. The resulting rAbs were characterized by their tumor-cell recognition, binding activity, and cell adhesion inhibition on triple-negative breast cancer cells. In addition, the rAbs were shown to enact antibody-dependent cellular cytotoxicity (ADCC) in the presence of either human natural killer cells or peripheral blood mononuclear cells, and were evaluated for the potential use of uPAR-targeting antibody-drug conjugates (ADCs). Three lead antibodies (11857, 8163, and 3159) were evaluated for their therapeutic efficacy in vivo and were shown to suppress tumor growth. Finally, the binding epitopes of the lead antibodies were characterized, providing information on their unique binding modes to uPAR. Altogether, the strategy identified unique cross-reactive antibodies with ADCC, ADC, and functional inhibitory effects by targeting cell-surface uPAR, that can be tested in safety studies and serve as potential immunotherapeutics.

Engineered bispecific antibodies targeting the interleukin-6 and -8 receptors potently inhibit cancer cell migration and tumor metastasis

Simultaneous inhibition of interleukin-6 (IL-6) and interleukin-8 (IL-8) signaling diminishes cancer cell migration, and combination therapy has recently been shown to synergistically reduce metastatic burden in a preclinical model of triple-negative breast cancer. Here, we have engineered two novel bispecific antibodies that target the IL-6 and IL-8 receptors to concurrently block the signaling activity of both ligands. We demonstrate that a first-in-class bispecific antibody design has promising therapeutic potential, with enhanced selectivity and potency compared with monoclonal antibody and small-molecule drug combinations in both cellular and animal models of metastatic triple-negative breast cancer. Mechanistic characterization revealed that our engineered bispecific antibodies have no impact on cell viability, but profoundly reduce the migratory potential of cancer cells; hence they constitute a true anti-metastatic treatment. Moreover, we demonstrate that our antibodies can be readily combined with standard-of-care anti-proliferative drugs to develop effective anti-cancer regimens. Collectively, our work establishes an innovative metastasis-focused direction for cancer drug development.

Advances in antibody discovery from human BCR repertoires

Antibodies make up an important and growing class of compounds used for the diagnosis or treatment of disease. While traditional antibody discovery utilized immunization of animals to generate lead compounds, technological innovations have made it possible to search for antibodies targeting a given antigen within the repertoires of B cells in humans. Here we group these innovations into four broad categories: cell sorting allows the collection of cells enriched in specificity to one or more antigens; BCR sequencing can be performed on bulk mRNA, genomic DNA or on paired (heavy-light) mRNA; BCR repertoire analysis generally involves clustering BCRs into specificity groups or more in-depth modeling of antibody-antigen interactions, such as antibody-specific epitope predictions; validation of antibody-antigen interactions requires expression of antibodies, followed by antigen binding assays or epitope mapping. Together with innovations in Deep learning these technologies will contribute to the future discovery of diagnostic and therapeutic antibodies directly from humans.

Analysis of TCR Repertoire by High-Throughput Sequencing Indicates the Feature of T Cell Immune Response after SARS-CoV-2 Infection

Coronavirus disease 2019 (COVID-19) is a global infectious disease caused by the SARS-CoV-2 coronavirus. T cells play an essential role in the body’s fighting against the virus invasion, and the T cell receptor (TCR) is crucial in T cell-mediated virus recognition and clearance. However, little has been known about the features of T cell response in convalescent COVID-19 patients. In this study, using 5′RACE technology and PacBio sequencing, we analyzed the TCR repertoire of COVID-19 patients after recovery for 2 weeks and 6 months compared with the healthy donors. The TCR clustering and CDR3 annotation were exploited to discover groups of patient-specific TCR clonotypes with potential SARS-CoV-2 antigen specificities. We first identified CD4⁺ and CD8⁺ T cell clones with certain clonal expansion after infection, and then observed the preferential recombination usage of V(D) J gene segments in CD4⁺ and CD8⁺ T cells of COVID-19 patients with different convalescent stages. More important, the TRBV6-5-TRBD2-TRBJ2-7 combination with high frequency was shared between CD4⁺ T and CD8⁺ T cells of different COVID-19 patients. Finally, we found the dominant characteristic motifs of the CDR3 sequence between recovered COVID-19 and healthy control. Our study provides novel insights on TCR in COVID-19 with different convalescent phases, contributing to our understanding of the immune response induced by SARS-CoV-2.

A New Type of Satellite Associated with Cassava Mosaic Begomoviruses

Cassava mosaic disease (CMD), which is caused by single-stranded DNA begomoviruses, severely limits cassava production across Africa. A previous study showed that CMD symptom severity and viral DNA accumulation increase in cassava in the presence of a DNA sequence designated SEGS-2 (sequence enhancing geminivirus symptoms). We report here that when SEGS-2 is coinoculated with African cassava mosaic virus (ACMV) onto Arabidopsis thaliana, viral symptoms increase. Transgenic Arabidopsis with an integrated copy of SEGS-2 inoculated with ACMV also display increased symptom severity and viral DNA levels. Moreover, SEGS-2 enables Cabbage leaf curl virus (CaLCuV) to infect a geminivirus-resistant Arabidopsis thaliana accession. Although SEGS-2 is related to cassava genomic sequences, an earlier study showed that it occurs as episomes and is packaged into virions in CMD-infected cassava and viruliferous whiteflies. We identified SEGS-2 episomes in SEGS-2 transgenic Arabidopsis. The episomes occur as both double-stranded and single-stranded DNA, with the single-stranded form packaged into virions. In addition, SEGS-2 episomes replicate in tobacco protoplasts in the presence, but not the absence, of ACMV DNA-A. SEGS-2 episomes contain a SEGS-2 derived promoter and an open reading frame with the potential to encode a 75-amino acid protein. An ATG mutation at the beginning of the SEGS-2 coding region does not enhance ACMV infection in A. thaliana. Together, the results established that SEGS-2 is a new type of begomovirus satellite that enhances viral disease through the action of an SEGS-2-encoded protein that may also be encoded by the cassava genome. IMPORTANCE Cassava is an important root crop in the developing world and a food and income crop for more than 300 million African farmers. Cassava is rising in global importance and trade as the demands for biofuels and commercial starch increase. More than half of the world's cassava is produced in Africa, where it is primarily grown by smallholder farmers, many of whom are from the poorest villages. Although cassava can grow under high temperature, drought, and poor soil conditions, its production is severely limited by viral diseases. Cassava mosaic disease (CMD) is one of the most important viral diseases of cassava and can cause up to 100% yield losses. We provide evidence that SEGS-2, which was originally isolated from cassava crops displaying severe and atypical CMD symptoms in Tanzanian fields, is a novel begomovirus satellite that can compromise the development of durable CMD resistance.

A CpG island promoter drives the CXXC5 gene expression

CXXC5 is a member of the zinc-finger CXXC family that binds to unmethylated CpG dinucleotides. CXXC5 modulates gene expressions resulting in diverse cellular events mediated by distinct signaling pathways. However, the mechanism responsible for CXXC5 expression remains largely unknown. We found here that of the 14 annotated CXXC5 transcripts with distinct 5' untranslated regions encoding the same protein, transcript variant 2 with the highest expression level among variants represents the main transcript in cell models. The DNA segment in and at the immediate 5'-sequences of the first exon of variant 2 contains a core promoter within which multiple transcription start sites are present. Residing in a region with high G-C nucleotide content and CpG repeats, the core promoter is unmethylated, deficient in nucleosomes, and associated with active RNA polymerase-II. These findings suggest that a CpG island promoter drives CXXC5 expression. Promoter pull-down revealed the association of various transcription factors (TFs) and transcription co-regulatory proteins, as well as proteins involved in histone/chromatin, DNA, and RNA processing with the core promoter. Of the TFs, we verified that ELF1 and MAZ contribute to CXXC5 expression. Moreover, the first exon of variant 2 may contain a G-quadruplex forming region that could modulate CXXC5 expression.

Gene Validation and Remodelling Using Proteogenomics of Phytophthora cinnamomi, the Causal Agent of Dieback

Phytophthora cinnamomi is a pathogenic oomycete that causes plant dieback disease across a range of natural ecosystems and in many agriculturally important crops on a global scale. An annotated draft genome sequence is publicly available (JGI Mycocosm) and suggests 26,131 gene models. In this study, soluble mycelial, extracellular (secretome), and zoospore proteins of P. cinnamomi were exploited to refine the genome by correcting gene annotations and discovering novel genes. By implementing the diverse set of sub-proteomes into a generated proteogenomics pipeline, we were able to improve the P. cinnamomi genome annotation. Liquid chromatography mass spectrometry was used to obtain high confidence peptides with spectral matching to both the annotated genome and a generated 6-frame translation. Two thousand seven hundred sixty-four annotations from the draft genome were confirmed by spectral matching. Using a proteogenomic pipeline, mass spectra were used to edit the P. cinnamomi genome and allowed identification of 23 new gene models and 60 edited gene features using high confidence peptides obtained by mass spectrometry, suggesting a rate of incorrect annotations of 3% of the detectable proteome. The novel features were further validated by total peptide support, alongside functional analysis including the use of Gene Ontology and functional domain identification. We demonstrated the use of spectral data in combination with our proteogenomics pipeline can be used to improve the genome annotation of important plant diseases and identify missed genes. This study presents the first use of spectral data to edit and manually annotate an oomycete pathogen.

The Identity Card of T Cells-Clinical Utility of T-cell Receptor Repertoire Analysis in Transplantation

There is a clear medical need to change the current strategy of "one-size-fits-all" immunosuppression for controlling transplant rejection to precision medicine and targeted immune intervention. As T cells play a key role in both undesired graft rejection and protection, a better understanding of the fate and function of both alloreactive graft-deteriorating T cells and those protecting to infections is required. The T-cell receptor (TCR) is the individual identity card of each T cell clone and can help to follow single specificities. In this context, tracking of lymphocytes with certain specificity in blood and tissue in clinical follow up is of especial importance. After overcoming technical limitations of the past, novel molecular technologies opened new avenues of diagnostics. Using advantages of next generation sequencing, a method was established for T-cell tracing by detection of variable TCR region as identifiers of individual lymphocyte clones. The current review describes principles of laboratory and computational methods of TCR repertoire analysis, and gives an overview on applications for the basic understanding of transplant biology and immune monitoring. The review also delineates methodological pitfalls and challenges. With the outlook on prediction of antigens in immune-mediated processes including those of unknown causative pathogens, monitoring the fate and function of individual T cell clones, and the adoptive transfer of protective effector or regulatory T cells, this review highlights the current and future capability of TCR repertoire analysis.

The role of microRNAs in the legume-Rhizobium nitrogen-fixing symbiosis

Under nitrogen starvation, most legume plants form a nitrogen-fixing symbiosis with Rhizobium bacteria. The bacteria induce the formation of a novel organ called the nodule in which rhizobia reside as intracellular symbionts and convert atmospheric nitrogen into ammonia. During this symbiosis, miRNAs are essential for coordinating the various plant processes required for nodule formation and function. miRNAs are non-coding, endogenous RNA molecules, typically 20-24 nucleotides long, that negatively regulate the expression of their target mRNAs. Some miRNAs can move systemically within plant tissues through the vascular system, which mediates, for example, communication between the stem/leaf tissues and the roots. In this review, we summarize the growing number of miRNAs that function during legume nodulation focusing on two model legumes, Lotus japonicus and Medicago truncatula, and two important legume crops, soybean (Glycine max) and common bean (Phaseolus vulgaris). This regulation impacts a variety of physiological processes including hormone signaling and spatial regulation of gene expression. The role of mobile miRNAs in regulating legume nodule number is also highlighted.

Unveiling ncRNA regulatory axes in atherosclerosis progression

Completion of the human genome sequencing project highlighted the richness of the cellular RNA world, and opened the door to the discovery of a plethora of short and long non-coding RNAs (the dark transcriptome) with regulatory or structural potential, which shifted the balance of pathological gene alterations from coding to non-coding RNAs. Thus, disease risk assessment currently has to also evaluate the expression of new RNAs such as small micro RNAs (miRNAs), long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), competing endogenous RNAs (ceRNAs), retrogressed elements, 3'UTRs of mRNAs, etc. We are interested in the pathogenic mechanisms of atherosclerosis (ATH) progression in patients suffering Chronic Kidney Disease, and in this review, we will focus in the role of the dark transcriptome (non-coding RNAs) in ATH progression. We will focus in miRNAs and in the formation of regulatory axes or networks with their mRNA targets and with the lncRNAs that function as miRNA sponges or competitive inhibitors of miRNA activity. In this sense, we will pay special attention to retrogressed genomic elements, such as processed pseudogenes and Alu repeated elements, that have been recently seen to also function as miRNA sponges, as well as to the use or miRNA derivatives in gene silencing, anti-ATH therapies. Along the review, we will discuss technical developments associated to research in lncRNAs, from sequencing technologies to databases, repositories and algorithms to predict miRNA targets, as well as new approaches to miRNA function, such as integrative or enrichment analysis and their potential to unveil RNA regulatory networks.

Assembly and functionality of the ribosome with tethered subunits

Ribo-T is an engineered ribosome whose small and large subunits are tethered together by linking 16S rRNA and 23S rRNA in a single molecule. Although Ribo-T can support cell proliferation in the absence of wild type ribosomes, Ribo-T cells grow slower than those with wild type ribosomes. Here, we show that cell growth defect is likely explained primarily by slow Ribo-T assembly rather than its imperfect functionality. Ribo-T maturation is stalled at a late assembly stage. Several post-transcriptional rRNA modifications and some ribosomal proteins are underrepresented in the accumulated assembly intermediates and rRNA ends are incompletely trimmed. Ribosome profiling of Ribo-T cells shows no defects in translation elongation but reveals somewhat higher occupancy by Ribo-T of the start codons and to a lesser extent stop codons, suggesting that subunit tethering mildly affects the initiation and termination stages of translation. Understanding limitations of Ribo-T system offers ways for its future development.

The tethered ribosome system Ribo-T supports cell proliferation though at a reduced rate. Here the authors show this is due to slower ribosome assembly instead of reduced functionality.

Genome organization and molecular characterization of the three Formica exsecta viruses-FeV1, FeV2 and FeV4

We present the genome organization and molecular characterization of the three Formica exsecta viruses, along with ORF predictions, and functional annotation of genes. The Formica exsecta virus-4 (FeV4; GenBank ID: MF287670) is a newly discovered negative-sense single-stranded RNA virus representing the first identified member of order Mononegavirales in ants, whereas the Formica exsecta virus-1 (FeV1; GenBank ID: KF500001), and the Formica exsecta virus-2 (FeV2; GenBank ID: KF500002) are positive single-stranded RNA viruses initially identified (but not characterized) in our earlier study. The new virus FeV4 was found by re-analyzing data from a study published earlier. The Formica exsecta virus-4 genome is 9,866 bp in size, with an overall G + C content of 44.92%, and containing five predicted open reading frames (ORFs). Our bioinformatics analysis indicates that gaps are absent and the ORFs are complete, which based on our comparative genomics analysis suggests that the genomes are complete. Following the characterization, we validate virus infection for FeV1, FeV2 and FeV4 for the first time in field-collected worker ants. Some colonies were infected by multiple viruses, and the viruses were observed to infect all castes, and multiple life stages of workers and queens. Finally, highly similar viruses were expressed in adult workers and queens of six other Formica species: F. fusca, F. pressilabris, F. pratensis, F. aquilonia, F. truncorum and F. cinerea. This research indicates that viruses can be shared between ant species, but further studies on viral transmission are needed to understand viral infection pathways.

New fusion sarcomas: histopathology and clinical significance of selected entities

Many sarcomas contain gene fusions that can be pathogenetic mechanisms and diagnostic markers. In this article we review selected fusion sarcomas and techniques for their detection. CIC-DUX4 fusion sarcoma is a round cell tumor now considered an entity separate from Ewing sarcoma with a more aggressive clinical course, occurrence in older age, and predilection to soft tissues. It is composed of larger cells than Ewing sarcoma and often has prominent necrosis. Nuclear DUX4 expression is a promising immuno histochemical marker. BCOR-CCNB3 fusion sarcoma is cyclin B3-positive, usually occurs in bone or soft tissue of children, and may mimic a poorly differentiated synovial sarcoma. EWSR1-NFATC2 sarcoma may present in bone or soft tissue. It is typically composed of small round cells in a trabecular pattern in a myxoid matrix resembling myoepithelioma. ACTB-GLI1 fusion sarcoma may mimic a skin adnexal carcinoma, showing focal expression of epithelial markers and S100 protein. NTRK-fusion sarcomas include, in addition to infantile fibrosarcoma with ETV6-NTRK3 fusion, LMNA-NTRK1 fusion sarcoma, a low-grade spindle cell sarcoma seen in peripheral soft tissues in children and young adults. Methods to detect gene fusions include next-generation sequencing panels, anchored multiplex polymerase chain reaction systems to detect partner for a known fusion gene, and comprehensive RNA sequencing to detect virtually all gene fusions. In situ hybridization testing using probes for both fusion partners can be used as an alternative confirmation technique, especially in the absence of satisfactory RNA yield. In addition, fusion protein-related and other immunohistochemical markers can have a high specificity for fusion sarcomas.

A comprehensive catalogue of polyketide synthase gene clusters in lichenizing fungi

Lichens are fungi that form symbiotic partnerships with algae. Although lichens produce diverse polyketides, difficulties in establishing and maintaining lichen cultures have prohibited detailed studies of their biosynthetic pathways. Creative, albeit non-definitive, methods have been developed to assign function to biosynthetic gene clusters in lieu of techniques such as gene knockout and heterologous expressions that are commonly applied to easily cultivatable organisms. We review a total of 81 completely sequenced polyketide synthase (PKS) genes from lichenizing fungi, comprising to our best efforts all complete and reported PKS genes in lichenizing fungi to date. This review provides an overview of the approaches used to locate and sequence PKS genes in lichen genomes, current approaches to assign function to lichen PKS gene clusters, and what polyketides are proposed to be biosynthesized by these PKS. We conclude with remarks on prospects for genomics-based natural products discovery in lichens. We hope that this review will serve as a guide to ongoing research efforts on polyketide biosynthesis in lichenizing fungi.

The complete genomic sequence of a novel botybirnavirus isolated from a phytopathogenic Bipolaris maydis

Bipolaris maydis is the causal agent of corn southern leaf blight. Here, we report a novel double-stranded RNA (dsRNA) mycovirus designated Bipolaris maydis botybirnavirus 1 (BmBRV1) from B. maydis strain JZ11 in Jingzhou, Hubei province of China. BmBRV1 has a genome consisting of two dsRNAs (dsRNA1 and dsRNA2) with a size of 6435 and 5987 bp, respectively, each of which contains a single open reading frame (ORF). The two polyproteins encoded by dsRNA1 and dsRNA2 share the highest amino acid identities of 81.8 and 75.3%, respectively, with the RdRp and coat protein of Sclerotinia sclerotiorum botybirnavirus 1 (SsBRV1), a tentative species of the genus Botybirnavirus. Phylogenetic analysis based on the amino acid sequences of RdRp indicated that BmBRV1 belongs to a distinct species of the newly proposed family Botybirnaviridae.

Analyzing Immunoglobulin Repertoires

Somatic assembly of T cell receptor and B cell receptor (BCR) genes produces a vast diversity of lymphocyte antigen recognition capacity. The advent of efficient high-throughput sequencing of lymphocyte antigen receptor genes has recently generated unprecedented opportunities for exploration of adaptive immune responses. With these opportunities have come significant challenges in understanding the analysis techniques that most accurately reflect underlying biological phenomena. In this regard, sample preparation and sequence analysis techniques, which have largely been borrowed and adapted from other fields, continue to evolve. Here, we review current methods and challenges of library preparation, sequencing and statistical analysis of lymphocyte receptor repertoire studies. We discuss the general steps in the process of immune repertoire generation including sample preparation, platforms available for sequencing, processing of sequencing data, measurable features of the immune repertoire, and the statistical tools that can be used for analysis and interpretation of the data. Because BCR analysis harbors additional complexities, such as immunoglobulin (Ig) (i.e., antibody) gene somatic hypermutation and class switch recombination, the emphasis of this review is on Ig/BCR sequence analysis.

Integrative Analysis of Three RNA Sequencing Methods Identifies Mutually Exclusive Exons of MADS-Box Isoforms During Early Bud Development in Picea abies

Recent efforts to sequence the genomes and transcriptomes of several gymnosperm species have revealed an increased complexity in certain gene families in gymnosperms as compared to angiosperms. One example of this is the gymnosperm sister clade to angiosperm TM3-like MADS-box genes, which at least in the conifer lineage has expanded in number of genes. We have previously identified a member of this sub-clade, the conifer gene DEFICIENS AGAMOUS LIKE 19 (DAL19), as being specifically upregulated in cone-setting shoots. Here, we show through Sanger sequencing of mRNA-derived cDNA and mapping to assembled conifer genomic sequences that DAL19 produces six mature mRNA splice variants in Picea abies. These splice variants use alternate first and last exons, while their four central exons constitute a core region present in all six transcripts. Thus, they are likely to be transcript isoforms. Quantitative Real-Time PCR revealed that two mutually exclusive first DAL19 exons are differentially expressed across meristems that will form either male or female cones, or vegetative shoots. Furthermore, mRNA in situ hybridization revealed that two mutually exclusive last DAL19 exons were expressed in a cell-specific pattern within bud meristems. Based on these findings in DAL19, we developed a sensitive approach to transcript isoform assembly from short-read sequencing of mRNA. We applied this method to 42 putative MADS-box core regions in P. abies, from which we assembled 1084 putative transcripts. We manually curated these transcripts to arrive at 933 assembled transcript isoforms of 38 putative MADS-box genes. 152 of these isoforms, which we assign to 28 putative MADS-box genes, were differentially expressed across eight female, male, and vegetative buds. We further provide evidence of the expression of 16 out of the 38 putative MADS-box genes by mapping PacBio Iso-Seq circular consensus reads derived from pooled sample sequencing to assembled transcripts. In summary, our analyses reveal the use of mutually exclusive exons of MADS-box gene isoforms during early bud development in P. abies, and we find that the large number of identified MADS-box transcripts in P. abies results not only from expansion of the gene family through gene duplication events but also from the generation of numerous splice variants.

Overexpression of the Prunus sogdiana NBS-LRR Subgroup Gene PsoRPM2 Promotes Resistance to the Root-Knot Nematode Meloidogyne incognita in Tobacco

Root-knot nematodes (RKNs), particularly Meloidogyne incognita, are the most devastating soil-borne pathogens that significantly affect the production of Prunus spp. fruit. RKN infection is difficult to control and consequently causes massive yield losses each year. However, several germplasms of wild Prunus spp. have been shown to display resistance to M. incognita. Consequently, both the isolation of novel plant resistance (R) genes and the characterization of their resistance mechanisms are important strategies for future disease control. R proteins require the co-chaperone protein HSP90-SGT1-RAR1 to achieve correct folding, maturation, and stabilization. Here, we used homologous cloning to isolate the R gene PsoRPM2 from the RKN-resistant species Prunus sogdiana. PsoRPM2 was found to encode a TIR-NB-LRR-type protein and react with significantly elevated PsoRPM2 expression levels in response to RKN infection. Transient expression assays indicated PsoRPM2 to be located in both the cytoplasm and the nucleus. Four transgenic tobacco lines that heterologously expressed PsoRPM2 showed enhanced resistance to M. incognita. Yeast two-hybrid analysis and bimolecular fluorescence complementation analysis demonstrated that both PsoRAR1 and PsoRPM2 interacted with PsoHSP90-1 and PsoSGT1, but not with one another. These results indicate that the observed PsoRPM2-mediated RKN resistance requires both PsoHSP90-1 and PsoSGT1, further suggesting that PsoRAR1 plays a functionally redundant role in the HSP90-SGT1-RAR1 co-chaperone.

Expression of the heterologous Dunaliella tertiolecta fatty acyl-ACP thioesterase leads to increased lipid production in Chlamydomonas reinhardtii

Biofuel production from genetically-engineered microalgae is currently among the most widely studied strategies in generating renewable energy. However, microalgae currently suffer from low oil yields which limit the commercial feasibility of industrial-scale production. A major bottleneck in cost-efficient biofuel production from microalgae is the dilemma between biomass productivity and lipid accumulation. When grown under stressful culture conditions such as nitrogen depletion, microalgae accumulate large amounts of neutral lipids, but it comes at the expense of growth which negatively impacts overall lipid productivity. Overexpression of acyl-ACP thioesterases (TE) had been successful in increasing the production of fatty acids (FA) in prokaryotes such as E. coli and cyanobacteria, but has not been effectively tested in microalgae. In this study, we introduced a TE from D. tertiolecta (DtTE) into C. reinhardtii to investigate its effects on FA production without compromising growth. The results indicate that C. reinhardtii transformants were able to produce 63 and 94% more neutral lipids than the wild-type, which translates to an approximately 56% improvement in total lipids, without compromising growth. These findings demonstrate the cross-species functionality of TE, and provide a platform for further studies into using TE as a strategy to increase biofuel production from microalgae.

Event-specific detection of transgenic potato AV43-6-G7 using real-time and digital PCR methods

Background

The isolation of unknown DNA sequences flanked by known sequences is an important task in the event-specific detection of GMOs. None of event-specific detection method was developed based on the junction sequence of an exogenous integrant in the transgenic potato AV43-6-G7.

Results

The flanking sequence between the exogenous fragment and recombinant chromosome of this potato was successfully acquired through exogenous gene 5′-RACE. The event-specific primers and Taqman probe were designed to amplify fragments spanning the exogenous DNA and potato genomic DNA. The specific real-time PCR and digital PCR detection methods for AV43-6-G7 potato were established based on primers designed according to the flanking sequences. The detection limit of the qualitative PCR assay was 0.01 % for AV43-6-G7 potato in 100 ng of potato genomic DNA, corresponding to approximately 11.6 copies of the potato haploid genome. The ddPCR assays for Potato AV43-6-G7 achieved a limit of quantification of approximately 58 target copies, with RSD ≤ 25 %. The aLOQ of this system was approximately 1.2 copies.

Conclusions

These results indicated that these event-specific methods would be useful for the identification of potato AV43-6-G7.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-016-0303-8) contains supplementary material, which is available to authorized users.

Extension of human lncRNA transcripts by RACE coupled with long-read high-throughput sequencing (RACE-Seq)

Long non-coding RNAs (lncRNAs) constitute a large, yet mostly uncharacterized fraction of the mammalian transcriptome. Such characterization requires a comprehensive, high-quality annotation of their gene structure and boundaries, which is currently lacking. Here we describe RACE-Seq, an experimental workflow designed to address this based on RACE (rapid amplification of cDNA ends) and long-read RNA sequencing. We apply RACE-Seq to 398 human lncRNA genes in seven tissues, leading to the discovery of 2,556 on-target, novel transcripts. About 60% of the targeted loci are extended in either 5′ or 3′, often reaching genomic hallmarks of gene boundaries. Analysis of the novel transcripts suggests that lncRNAs are as long, have as many exons and undergo as much alternative splicing as protein-coding genes, contrary to current assumptions. Overall, we show that RACE-Seq is an effective tool to annotate an organism's deep transcriptome, and compares favourably to other targeted sequencing techniques.

Long non-coding RNAs are increasingly recognised to be important factors in regulating cellular processes and comprise a large faction of the transcriptome, however most are uncharacterised. Here the authors present RACE-Seq, a tool to improve and extend the annotation of low-expression transcripts.

Identification of 6-Hydroxymellein Synthase and Accessory Genes in the Lichen Cladonia uncialis

A transcribed polyketide synthase (PKS) gene has been identified in the lichen Cladonia uncialis. The complete nucleotide sequence of this PKS was determined from the amplified cDNA, and an assignment of individual domains was accomplished by homology searching using AntiSMASH. A scan of the complete genome sequence of C. uncialis revealed the accessory genes associated with this PKS gene. A homology search has identified that several genes in this cluster are similar to genes responsible for the biosynthesis of terrein in Aspergillus terreus. This permitted assignment of putative function to each of the genes in this new C. uncialis cluster. It is proposed that this gene cluster is responsible for the biosynthesis of a halogenated iscoumarin. This is the first report linking a gene cluster to a halogenated metabolite in lichen.

Transcriptome-enabled discovery and functional characterization of enzymes related to (2S)-pinocembrin biosynthesis from Ornithogalum caudatum and their application for metabolic engineering

Background

(2S)-Pinocembrin is a chiral flavanone with versatile pharmacological and biological activities. Its health-promoting effects have spurred on research effects on the microbial production of (2S)-pinocembrin. However, an often-overlooked salient feature in the analysis of microbial (2S)-pinocembrin is its chirality.

Results

Here, we presented a full characterization of absolute configuration of microbial (2S)-pinocembrin from engineered Escherichia coli. Specifically, a transcriptome-wide search for genes related to (2S)-pinocembrin biosynthesis from Ornithogalum caudatum, a plant rich in flavonoids, was first performed in the present study. A total of 104,180 unigenes were finally generated with an average length of 520 bp. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway mapping assigned 26 unigenes, representing three enzyme families of 4-coumarate:coenzyme A ligase (4CL), chalcone synthase (CHS) and chalcone isomerase(CHI), onto (2S)-pinocembrin biosynthetic pathway. A total of seven, three and one full-length candidates encoding 4CL, CHS and CHI were then verified by reverse transcription polymerase chain reaction, respectively. These candidates were screened by functional expression in E. coli individual or coupled multienzyme reaction systems based on metabolic engineering processes. Oc4CL1, OcCHS2 and OcCHI were identified to be bona fide genes encoding respective pathway enzymes of (2S)-pinocembrin biosynthesis. Then Oc4CL1, OcCHS2 and MsCHI from Medicago sativa, assembled as artificial gene clusters in different organizations, were used for fermentation production of (2S)-pinocembrin in E. coli. The absolute configuration of the resulting microbial pinocembrin at C-2 was assigned to be 2S-configured by combination of retention time, UV spectrum, LC–MS, NMR, optical rotation and circular dichroism spectroscopy. Improvement of (2S)-pinocembrin titres was then achieved by optimization of gene organizations, using of codon-optimized pathway enzymes and addition of cerulenin for increasing intracellular malonyl CoA pools. Overall, the optimized strain can produce (2S)-pinocembrin of 36.92 ± 4.1 mg/L.

Conclusions

High titre of (2S)-pinocembrin can be obtained from engineered E. coli by an efficient method. The fermentative production of microbial (2S)-pinocembrin in E. coli paved the way for yield improvement and further pharmacological testing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0424-8) contains supplementary material, which is available to authorized users.

Identification and Validation of Putative Nesprin Variants

Nesprins are a family of multi-isomeric scaffolding proteins that were originally identified at the nuclear envelope (NE), where they bind to lamin A/C, emerin, and SUN-domain containing proteins, to form the LInker of Nucleoskeleton-and-Cytoskeleton (LINC) complex that connects the NE to the actin cytoskeleton. However, nesprin genes also give rise to a variety of tissue-specific variants of different sizes with potential roles beyond the NE. These variants are generated through alternative initiation, termination, and splicing, which makes nesprin biology very complex to study due to the difficulty in generating specific antibodies and/or short interfering RNAs (siRNA) to particular isoforms. In order to distinguish genuine nesprin variants and eliminate confusion with degradation products of larger nesprin isoforms, in this chapter we discuss methods including 5' and 3' Rapid Amplification of cDNA Ends (RACE) and RT-PCR in combination with EST database searching, for identifying and validating putative nesprin isoforms. This information is essential to allow a better understanding of nesprin functions in different cell types.

Discovery of a Novel Human Pegivirus in Blood Associated with Hepatitis C Virus Co-Infection

Hepatitis C virus (HCV) and human pegivirus (HPgV), formerly GBV-C, are the only known human viruses in the Hepacivirus and Pegivirus genera, respectively, of the family Flaviviridae. We present the discovery of a second pegivirus, provisionally designated human pegivirus 2 (HPgV-2), by next-generation sequencing of plasma from an HCV-infected patient with multiple bloodborne exposures who died from sepsis of unknown etiology. HPgV-2 is highly divergent, situated on a deep phylogenetic branch in a clade that includes rodent and bat pegiviruses, with which it shares <32% amino acid identity. Molecular and serological tools were developed and validated for high-throughput screening of plasma samples, and a panel of 3 independent serological markers strongly correlated antibody responses with viral RNA positivity (99.9% negative predictive value). Discovery of 11 additional RNA-positive samples from a total of 2440 screened (0.45%) revealed 93–94% nucleotide identity between HPgV-2 strains. All 12 HPgV-2 RNA-positive cases were identified in individuals also testing positive for HCV RNA (12 of 983; 1.22%), including 2 samples co-infected with HIV, but HPgV-2 RNA was not detected in non-HCV-infected individuals (p<0.0001), including those singly infected by HIV (p = 0.0075) or HBV (p = 0.0077), nor in volunteer blood donors (p = 0.0082). Nine of the 12 (75%) HPgV-2 RNA positive samples were reactive for antibodies to viral serologic markers, whereas only 28 of 2,429 (1.15%) HPgV-2 RNA negative samples were seropositive. Longitudinal sampling in two individuals revealed that active HPgV-2 infection can persist in blood for at least 7 weeks, despite the presence of virus-specific antibodies. One individual harboring both HPgV-2 and HCV RNA was found to be seronegative for both viruses, suggesting a high likelihood of simultaneous acquisition of HCV and HPgV-2 infection from an acute co-transmission event. Taken together, our results indicate that HPgV-2 is a novel bloodborne infectious virus of humans and likely transmitted via the parenteral route.

To date, only one human hepacivirus (HCV) and one human pegivirus (HPgV-1/GBV-C) in the family Flaviviridae are known to exist. Using unbiased metagenomic next-generation sequencing, we discovered and assembled the genome of a novel pegivirus from plasma corresponding to an HCV-infected patient who died from unknown sepsis. This virus, provisionally named human pegivirus 2 (HPgV-2), is highly divergent, sharing <32% amino acid identity with its nearest relatives, a bat and rodent pegivirus. Identification and sequencing of 11 additional HPgV-2 viruses, revealing 93–94% identity between strains, as well as documented antibody responses using multiple markers, confirm that HPgV-2 is a bona fide novel infectious virus of humans. Several lines of evidence, including (1) a documented history of multiple bloodborne exposures in the index patient, (2) parallel detection of HPgV-2 and HCV RNA in an individual during the "window period" between infection and the appearance of detectable antibody, and (3) a tight observed association between HPgV-2 and HCV co-infection, suggest that HPgV-2 is an infectious agent capable of bloodborne transmission.

Isolation and Analysis of the Cppsy Gene and Promoter from Chlorella protothecoides CS-41

Phytoene synthase (PSY) catalyzes the condensation of two molecules of geranylgeranyl pyrophosphate to form phytoene, the first colorless carotene in the carotenoid biosynthesis pathway. So it is regarded as the crucial enzyme for carotenoid production, and has unsurprisingly been involved in genetic engineering studies of carotenoid production. In this study, the psy gene from Chlorella protothecoides CS-41, designated Cppsy, was cloned using rapid amplification of cDNA ends. The full-length DNA was 2488 bp, and the corresponding cDNA was 1143 bp, which encoded 380 amino acids. Computational analysis suggested that this protein belongs to the Isoprenoid_Biosyn_C1 superfamily. It contained the consensus sequence, including three predicted substrate-Mg²⁺ binding sites. The Cppsy gene promoter was also cloned and characterized. Analysis revealed several candidate motifs for the promoter, which exhibited light- and methyl jasmonate (MeJA)-responsive characteristics, as well as some typical domains universally discovered in promoter sequences, such as the TATA-box and CAAT-box. Light- and MeJA treatment showed that the Cppsy expression level was significantly enhanced by light and MeJA. These results provide a basis for genetically modifying the carotenoid biosynthesis pathway in C. protothecoides.

Cloning of expansin genes in ramie (Boehmeria nivea L.) based on universal fast walking

Gene cloning is the first step to study the expression profiles and functions of a particular gene; considerable cloning methods have been developed. Expansin, thought to involve in the cell-wall modification events, was not cloned in ramie (Boehmeria nivea L.), which is one of the most important bast fiber crops with little conducted molecular research, especially on its fiber development. Studying the expansin gene family will uncover its possible relationship with ramie fiber development and other growth events. As a result, five expansin genes were cloned with full-length and their sequence information was investigated. Additionally, the phylogenetic analysis was conducted, which suggested that the cloned genes belong to the α-subfamily, and these genes expressed differently during ramie fiber developmental process. In this study, we aimed to apply a strategy for cloning novel full-length genes from genomic DNA of ramie, based on using degenerate primers, touchdown polymerase chain reaction and universal fast walking protocols. By cloning five full-length expansin genes, we believe the polymerase chain reaction-based gene cloning strategy could be applied to general gene studies in ramie and other crops.

Systematic dissection of the sequence determinants of gene 3' end mediated expression control

The 3'end genomic region encodes a wide range of regulatory process including mRNA stability, 3' end processing and translation. Here, we systematically investigate the sequence determinants of 3' end mediated expression control by measuring the effect of 13,000 designed 3' end sequence variants on constitutive expression levels in yeast. By including a high resolution scanning mutagenesis of more than 200 native 3' end sequences in this designed set, we found that most mutations had only a mild effect on expression, and that the vast majority (~90%) of strongly effecting mutations localized to a single positive TA-rich element, similar to a previously described 3' end processing efficiency element, and resulted in up to ten-fold decrease in expression. Measurements of 3' UTR lengths revealed that these mutations result in mRNAs with aberrantly long 3'UTRs, confirming the role for this element in 3' end processing. Interestingly, we found that other sequence elements that were previously described in the literature to be part of the polyadenylation signal had a minor effect on expression. We further characterize the sequence specificities of the TA-rich element using additional synthetic 3' end sequences and show that its activity is sensitive to single base pair mutations and strongly depends on the A/T content of the surrounding sequences. Finally, using a computational model, we show that the strength of this element in native 3' end sequences can explain some of their measured expression variability (R = 0.41). Together, our results emphasize the importance of efficient 3' end processing for endogenous protein levels and contribute to an improved understanding of the sequence elements involved in this process.

Rapid mapping of RNA 3' and 5' ends

In recent years, an enormous progress has been made in applied genomics leading to identification and isolation of novel cDNAs. However, most attempts result in the acquisition of transcribed sequences that represent only a part of the mRNA's complete sequence. Rapid Amplification of cDNA Ends (RACE) is a technique used in molecular biology to obtain the full length sequence of an RNA transcript found within a cell. Since the first report of this technique, many significant improvements have been made on the basic approach. This chapter describes the most recent update of the relatively simple and versatile classic RACE protocol.

Construction, characterization and expression of full length cDNA clone of sheep YAP1 gene

RT-PCR, 5'RACE, 3'RACE were used to clone sheep full length cDNA sequence of YAP1 (Yes-associated protein 1), eukaryotic expression plasmid and a mutant that cannot be phosphorylated at Ser42 was successfully constructed. The amino acid sequence analysis revealed that sheep YAP1 gene encoded water-soluble protein and its relative molecular weight and isoelectric point was 44,079.0 Da and 4.91, respectively. Sub-cellular localization of YAP1 was in the nucleus, it is hydrophilic non-transmembrane and non-secreted protein. YAP1 protein contained 33 phosphorylation sites, seven glycosylation sites and two WW domains. The secondary structure of YAP1 was mainly composed of random coil, while the tertiary structure of domain area showed a forniciform helix structure. YAP1 gene was expressed in different tissues, the highest expression was in kidney and the lowest was in hypothalamus. The CDS of sheep YAP1was amplified by RT-PCR from healthy sheep longissimus dorsi muscle, cloned into pMD19-T simple vector by T/A ligation. YAP1 coding region was further sub-cloned into pEGFP-C1 vector by T4 Ligase to construct a eukaryotic expression plasmid and then make the eukaryotic expression vector as the template to construct the phosphorylation site mutant. PCR, restriction enzyme and sequencing were used to confirm the recombinant plasmid. The sheep full-length YAP1 cDNA sequence is 1712 in length encoding 403 amino acids. It was confirmed that the sheep YAP1 CDS was correctly inserted into eukaryotic expression vector and serine had been mutated to alanine by PCR, restriction digestion and sequencing. The result showed that the recombinant plasmid pEGFP-C1-YAP1 and pEGFP-C1-YAP1 S42A was constructed correctly, this will help for further studies on the YAP1 protein expression and its biological activities.

Molecular cloning and characterization of a novel P450 gene encoding CYP6BK18 from Dastarcus helophoroides (Coleoptera: Bothrideridae)

A novel cDNA clone encoding a cytochrome P450 gene, named CYP6BK18 (GenBank KC683905), was isolated by reverse transcription PCR from Dastarcus helophoroides (Fairmaire) (Coleoptera: Bothrideridae), a natural enemy of beetles. The full-length cDNA sequence is 1,659 bp, containing a 1,533 bp open reading frame predicting a 510-amino acid protein possessing a transmembrane domain with a calculated molecular weight of 59.4 kDa and a theoretical pI of 8.94. The deduced amino acid sequence of CYP6BK18 showed a 59% identity with CYP6BK17 (GenBank XP_970481.1) from Tribolium castaneum. Phylogenetic analysis indicated that CYP6BK18 was most closely related to CYP6BK17 and CYP6BK14 (GenBank EFA05731.1) from T. castaneum. Expression patterns of CYP6BK18 in different tissues (head, oviduct, midgut, fat bodies, and Malpighian tubules), developmental stages (first- to sixth-instar larvae and adult) and 10 age groups of adult were analyzed by real-time quantitative PCR (RT-qPCR). The results showed that CYP6BK18 was highly expressed in adulthood. Also, RT-qPCR analysis among different age groups of adult showed that CYP6BK18 transcripts were abundant in the spawning period and peaked at the early stage of the adult development. Moreover, the tissue-specific expression levels of CYP6BK18 were in the order of midgut, Malpighian tubes and fatty body from high to low. These results suggest that cytochrome P450 CYP6BK18 may play a role in regulating the development and aging of D. helophoroides.

Measurements of the impact of 3' end sequences on gene expression reveal wide range and sequence dependent effects

A full understanding of gene regulation requires an understanding of the contributions that the various regulatory regions have on gene expression. Although it is well established that sequences downstream of the main promoter can affect expression, our understanding of the scale of this effect and how it is encoded in the DNA is limited. Here, to measure the effect of native S. cerevisiae 3′ end sequences on expression, we constructed a library of 85 fluorescent reporter strains that differ only in their 3′ end region. Notably, despite being driven by the same strong promoter, our library spans a continuous twelve-fold range of expression values. These measurements correlate with endogenous mRNA levels, suggesting that the 3′ end contributes to constitutive differences in mRNA levels. We used deep sequencing to map the 3′UTR ends of our strains and show that determination of polyadenylation sites is intrinsic to the local 3′ end sequence. Polyadenylation mapping was followed by sequence analysis, we found that increased A/T content upstream of the main polyadenylation site correlates with higher expression, both in the library and genome-wide, suggesting that native genes differ by the encoded efficiency of 3′ end processing. Finally, we use single cells fluorescence measurements, in different promoter activation levels, to show that 3′ end sequences modulate protein expression dynamics differently than promoters, by predominantly affecting the size of protein production bursts as opposed to the frequency at which these bursts occur. Altogether, our results lead to a more complete understanding of gene regulation by demonstrating that 3′ end regions have a unique and sequence dependent effect on gene expression.

A basic question in gene expression is the relative contribution of different regulatory layers and genomic regions to the differences in protein levels. In this work we concentrated on the effect of 3′ end sequences. For this, we constructed a library of yeast strains that differ only by a native 3′ end region integrated downstream to a reported gene driven by a constant inducible promoter. Thus we could attribute all differences in reporter expression between the strains to the different 3′ end sequences. Interestingly, we found that despite being driven by the same strong, inducible promoter, our library spanned a wide and continuous range of expression levels of more than twelve-fold. As these measurements represent the sole effect of the 3′ end region, we quantify the contribution of these sequences to the variance in mRNA levels by comparing our measurements to endogenous mRNA levels. We follow by sequence analysis to find a simple sequence signature that correlates with expression. In addition, single cell analysis reveals distinct noise dynamics of 3′ end mediated differences in expression compared to different levels of promoter activation leading to a more complete understanding of gene expression which also incorporates the effect of these regions.

A hybridization-based approach for quantitative and low-bias single-stranded DNA ligation

Single-stranded DNA (ssDNA) ligation is a crucial step in many biochemical assays. Efficient ways of carrying out this reaction, however, are lacking. We show here that existing ssDNA ligation methods suffer from slow kinetics, poor yield, and severe nucleotide preference. To resolve these issues, we introduce a hybridization-based strategy that provides efficient and low-bias ligation of ssDNA. Our method uses a hairpin DNA to hybridize to any incoming acceptor ssDNA with low bias, with ligation of these strands mediated by T4 DNA ligase. This technique potentially can be applied in protocols that require ligation of ssDNA, including ligation-mediated polymerase chain reaction (LMPCR) and complementary DNA (cDNA) library construction.

Linking genome-scale metabolic modeling and genome annotation

Genome-scale metabolic network reconstructions, assembled from annotated genomes, serve as a platform for integrating data from heterogeneous sources and generating hypotheses for further experimental validation. Implementing constraint-based modeling techniques such as flux balance analysis (FBA) on network reconstructions allows for interrogating metabolism at a systems level, which aids in identifying and rectifying gaps in knowledge. With genome sequences for various organisms from prokaryotes to eukaryotes becoming increasingly available, a significant bottleneck lies in the structural and functional annotation of these sequences. Using topologically based and biologically inspired metabolic network refinement, we can better characterize enzymatic functions present in an organism and link annotation of these functions to candidate transcripts; both steps can be experimentally validated.