Upstream Sequence Research Articles

Genome-Wide Identification and Comparison of Cysteine Proteases in the Pollen Coat and Other Tissues in Maize.

Cysteine proteases, belonging to the C1-papain family, play a major role in plant growth and development, senescence, and immunity. There is evidence to suggest that pollen cysteine protease (CP) (ZmCP03) is involved in regulating the anther development and pollen formation in maize. However, there is no report on the genome-wide identification and comparison of CPs in the pollen coat and other tissues in maize. In this study, a total of 38 homologous genes of ZmCP03 in maize were identified. Subsequently, protein motifs, conserved domains, gene structures, and duplication patterns of 39 CPs are analyzed to explore their evolutionary relationship and potential functions. The cis-elements were identified in the upstream sequence of 39 CPs, especially those that are related to regulating growth and development and responding to environmental stresses and hormones. The expression patterns of these genes displayed remarked difference at a tissue or organ level in maize based on the available transcriptome data in the public database. Quantitative reverse transcription PCR (RT-qPCR) analysis showed that ZmCP03 was preferably expressed at a high level in maize pollen. Analyses by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and immunoblot, immunofluorescence and immunogold electron microscopy all validated the cellular localization of ZmCP03 in both the pollen coat and pollen cytoplasm. In addition, 142 CP genes from Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa) and cotton (Gossypium hirsutum), together with 39 maize CPs, were retrieved to analyze their evolution by comparing with orthologous genes. The results suggested that ZmCP03 was relatively conservative and stable during evolution. This study may provide a referential evidence on the function of ZmCP03 in pollen development and germination in maize.

Genetic Diversity in Casein Gene Cluster in a Dromedary Camel (C. dromedarius) Population from the United Arab Emirates.

Genetic polymorphisms, causing variation in casein genes (CSN1S1, CSN1S2, CSN2, and CSN3), have been extensively studied in goats and cows, but there are only few studies reported in camels. Therefore, we aimed to identify alleles with functional roles in the United Arab Emirates dromedary camel (Camelus dromedarius) population to complement previous studies conducted on the same species. Using targeted next-generation sequencing, we sequenced all genes in the casein gene cluster in 93 female camels to identify and characterize novel gene variants. Most variants were found in noncoding introns and upstream sequences, but a few variants showed the possibility of functional impact. CSN2 was found to be most polymorphic, with total 91 different variants, followed by CSN1S1, CSN3 and CSN1S2. CSN1S1, CSN1S2 and CSN2 each had at least two variants while CSN3 had only one functional allele. In future research, the functional impact of these variants should be investigated further.

GmPGL2, Encoding a Pentatricopeptide Repeat Protein, Is Essential for Chloroplast RNA Editing and Biogenesis in Soybean.

Chloroplast biogenesis and development are highly complex processes requiring interactions between plastids and nuclear genomic products. Pentatricopeptide repeat (PPR) proteins play an essential role in the development of chloroplasts; however, it remains unclear how RNA editing factors influence soybean development. In this study, a Glycine max pale green leaf 2 mutant (Gmpgl2) was identified with decreased chlorophyll contents. Genetic mapping revealed that a single-nucleotide deletion at position 1949 bp in the Glyma.05g132700 gene in the Gmpgl2 mutant, resulting in a truncated GmPGL2 protein. The nuclear-encoded GmPGL2 is a PLS-type PPR protein that localizes to the chloroplasts. The C-to-U editing efficiencies of rps16, rps18, ndhB, ndhD, ndhE, and ndhF were reduced in the Gmpgl2 mutant. RNA electrophoresis mobility shift assay (REMSA) analysis further revealed that GmPGL2 binds to the immediate upstream sequences at RNA editing sites of rps16 and ndhB in vitro, respectively. In addition, GmPGL2 was found to interact with GmMORF8, GmMORF9, and GmORRM6. These results suggest that GmPGL2 participates in C-to-U RNA editing via the formation of a complex RNA editosome in soybean chloroplasts.

Utr.annotation: a tool for annotating genomic variants that could influence post-transcriptional regulation.

Whole genome sequencing of patient populations is identifying thousands of new variants in untranslated regions (UTRs). While the consequences of UTR mutations are not as easily predicted from primary sequence as coding mutations are, there are some known features of UTRs that modulate their function. utr.annotation is an R package that can be used to annotate potential deleterious variants in the UTR regions for both human and mouse species. Given a CSV or VCF format variant file, utr.annotation provides information of each variant on whether and how it alters known translational regulators including upstream open reading frames, upstream Kozak sequences, polyA signals, Kozak sequences at the annotated translation start site, start codons and stop codons, conservation scores in the variant position, and whether and how it changes ribosome loading based on a model derived from empirical data. utr.annotation is freely available on Bitbucket (https://bitbucket.org/jdlabteam/utr.annotation/src/master/) and CRAN (https://cran.r-project.org/web/packages/utr.annotation/index.html). Supplementary data are available at Bioinformatics online.

Pitx controls amphioxus asymmetric morphogenesis by promoting left-side development and repressing right-side formation

BackgroundLeft-right (LR) asymmetry is an essential feature of bilateral animals. Studies in vertebrates show that LR asymmetry formation comprises three major steps: symmetry breaking, asymmetric gene expression, and LR morphogenesis. Although much progress has been made in the first two events, mechanisms underlying asymmetric morphogenesis remain largely unknown due to the complex developmental processes deployed by vertebrate organs.ResultsWe here addressed this question by studying Pitx gene function in the basal chordate amphioxus whose asymmetric organogenesis, unlike that in vertebrates, occurs essentially in situ and does not rely on cell migration. Pitx null mutation in amphioxus causes loss of all left-sided organs and incomplete ectopic formation of all right-sided organs on the left side, whereas Pitx partial loss-of-function leads to milder phenotypes with only some LR organs lost or ectopically formed. At the N1 to N3 stages, Pitx expression is gradually expanded from the dorsal anterior domain to surrounding regions. This leads to activation of genes like Lhx3 and/or Prop1 and Pit, which are essential for left-side organs, and downregulation of genes like Hex and/or Nkx2.1 and FoxE4, which are required for right-side organs to form ectopically on the left side. In Pitx mutants, the left-side expressed genes are not activated, while the right-side genes fail to decrease expression on the left side. In contrast, in embryos overexpressing Pitx genes, the left-side genes are induced ectopically on the right side, and the right-side genes are inhibited. Several Pitx binding sites are identified in the upstream sequences of the left-side and right-side genes which are essential for activation of the former and repression of the latter by Pitx.ConclusionsOur results demonstrate that (1) Pitx is a major (although not the only) determinant of asymmetric morphogenesis in amphioxus, (2) the development of different LR organs have distinct requirements for Pitx activity, and (3) Pitx controls amphioxus LR morphogenesis probably through inducing left-side organs and inhibiting right-side organs directly. These findings show much more dependence of LR organogenesis on Pitx in amphioxus than in vertebrates. They also provide insight into the molecular developmental mechanism of some vertebrate LR organs like the lungs and atria, since they show a right-isomerism phenotype in Pitx2 knockout mice like right-sided organs in Pitx mutant amphioxus. Our results also explain why some organs like the adenohypophysis are asymmetrically located in amphioxus but symmetrically positioned in vertebrates.

Genome-wide identification and analysis of the heat shock transcription factor family in moso bamboo (Phyllostachys edulis)

Heat shock transcription factors (HSFs) are central elements in the regulatory network that controls plant heat stress response. They are involved in multiple transcriptional regulatory pathways and play important roles in heat stress signaling and responses to a variety of other stresses. We identified 41 members of the HSF gene family in moso bamboo, which were distributed non-uniformly across its 19 chromosomes. Phylogenetic analysis showed that the moso bamboo HSF genes could be divided into three major subfamilies; HSFs from the same subfamily shared relatively conserved gene structures and sequences and encoded similar amino acids. All HSF genes contained HSF signature domains. Subcellular localization prediction indicated that about 80% of the HSF proteins were located in the nucleus, consistent with the results of GO enrichment analysis. A large number of stress response–associated cis-regulatory elements were identified in the HSF upstream promoter sequences. Synteny analysis indicated that the HSFs in the moso bamboo genome had greater collinearity with those of rice and maize than with those of Arabidopsis and pepper. Numerous segmental duplicates were found in the moso bamboo HSF gene family. Transcriptome data indicated that the expression of a number of PeHsfs differed in response to exogenous gibberellin (GA) and naphthalene acetic acid (NAA). A number of HSF genes were highly expressed in the panicles and in young shoots, suggesting that they may have functions in reproductive growth and the early development of rapidly-growing shoots. This study provides fundamental information on members of the bamboo HSF gene family and lays a foundation for further study of their biological functions in the regulation of plant responses to adversity.

Efficient integration of transmembrane domains depends on the folding properties of the upstream sequences

The topology of most membrane proteins is defined by the successive integration of α-helical transmembrane domains at the Sec61 translocon. The translocon provides a pore for the transfer of polypeptide segments across the membrane while giving them lateral access to the lipid. For each polypeptide segment of ∼20 residues, the combined hydrophobicities of its constituent amino acids were previously shown to define the extent of membrane integration. Here, we discovered that different sequences preceding a potential transmembrane domain substantially affect its hydrophobicity requirement for integration. Rapidly folding domains, sequences that are intrinsically disordered or very short or capable of binding chaperones with high affinity, allow for efficient transmembrane integration with low-hydrophobicity thresholds for both orientations in the membrane. In contrast, long protein fragments, folding-deficient mutant domains, and artificial sequences not binding chaperones interfered with membrane integration, requiring higher hydrophobicity. We propose that the latter sequences, as they compact on their hydrophobic residues, partially folded but unable to reach a native state, expose hydrophobic surfaces that compete with the translocon for the emerging transmembrane segment, reducing integration efficiency. The results suggest that rapid folding or strong chaperone binding is required for efficient transmembrane integration.

Cryptic promoter activation occurs by at least two different mechanisms in the Arabidopsis genome.

In gene-trap screening of plant genomes, promoterless reporter constructs are often expressed without trapping of annotated gene promoters. The molecular basis of this phenomenon, which has been interpreted as the trapping of cryptic promoters, is poorly understood. Here, we found that cryptic promoter activation occurs by at least two different mechanisms using Arabidopsis gene-trap lines in which a firefly luciferase (LUC) open reading frame (ORF) without an apparent promoter sequence was expressed from intergenic regions: one mechanism is 'cryptic promoter capturing', in which the LUC ORF captured pre-existing promoter-like chromatin marked by H3K4me3 and H2A.Z, and the other is 'promoter de novo origination', in which the promoter chromatin was newly formed near the 5' end of the inserted LUC ORF. The latter finding raises a question as to how the inserted LUC ORF sequence is involved in this phenomenon. To examine this, we performed a model experiment with chimeric LUC genes in transgenic plants. Using Arabidopsis psaH1 promoter-LUC constructs, we found that the functional core promoter region, where transcription start sites (TSSs) occur, cannot simply be determined by the upstream nor core promoter sequences; rather, its positioning proximal to the inserted LUC ORF sequence was more critical. This result suggests that the insertion of the coding sequence alters the local distribution of TSSs in the plant genome. The possible impact of the two types of cryptic promoter activation mechanisms on plant genome evolution and endosymbiotic gene transfer is discussed.

Antibody upstream sequence diversity and its biological implications revealed by repertoire sequencing

The sequence upstream of the antibody variable region (antibody upstream sequence [AUS]) consists of a 5′ untranslated region (5′ UTR) and a preceding leader region. The sequence variations in AUS affect antibody engineering and PCR based antibody quantification and may also be implicated in mRNA transcription and translation. However, the diversity of AUSs remains elusive. Using 5′ rapid amplification of cDNA ends and high-throughput antibody repertoire sequencing technique, we acquired full-length AUSs for human, rhesus macaque, cynomolgus macaque, mouse, and rat. We designed a bioinformatics pipeline and identified 3307 unique AUSs, corresponding to 3026 and 1457 unique sequences for 5′ UTR and leader region, respectively. Comparative analysis indicated that 928 (63.69%) leader sequences are novel relative to those recorded in the international ImMunoGeneTics information system. Evolutionarily, leader sequences are more conserved than 5′ UTR and seem to coevolve with their downstream V genes. Besides, single-nucleotide polymorphisms are position dependent for leader regions and may contribute to the functional reversal of the downstream V genes. Finally, the AUGs in AUSs were found to have little impact on gene expression. Taken together, our findings can facilitate primer design for capturing antibodies efficiently and provide a valuable resource for antibody engineering and molecule-level antibody studies.

QuantifyPoly(A): reshaping alternative polyadenylation landscapes of eukaryotes with weighted density peak clustering.

The dynamic choice of different polyadenylation sites in a gene is referred to as alternative polyadenylation, which functions in many important biological processes. Large-scale messenger RNA 3' end sequencing has revealed that cleavage sites for polyadenylation are presented with microheterogeneity. To date, the conventional determination of polyadenylation site clusters is subjective and arbitrary, leading to inaccurate annotations. Here, we present a weighted density peak clustering method, QuantifyPoly(A), to accurately quantify genome-wide polyadenylation choices. Applying QuantifyPoly(A) on published 3' end sequencing datasets from both animals and plants, their polyadenylation profiles are reshaped into myriads of novel polyadenylation site clusters. Most of these novel polyadenylation site clusters show significantly dynamic usage across different biological samples or associate with binding sites of trans-acting factors. Upstream sequences of these clusters are enriched with polyadenylation signals UGUA, UAAA and/or AAUAAA in a species-dependent manner. Polyadenylation site clusters also exhibit species specificity, while plants ones generally show higher microheterogeneity than that of animals. QuantifyPoly(A) is broadly applicable to any types of 3' end sequencing data and species for accurate quantification and construction of the complex and dynamic polyadenylation landscape and enables us to decode alternative polyadenylation events invisible to conventional methods at a much higher resolution.

Select Porcine Elongation Factor 1α Sequences Mediate Stable High-Level and Upregulated Expression of Heterologous Genes in Porcine Cells in Response to Primate Serum.

Genetically engineered (GE) pigs with various combinations of genetic profiles have been developed using heterologous promoters. This study aimed to identify autologous promoters for high and ubiquitous expression of xenotransplantation relevant genes in GE pigs. A 1.4 kb upstream regulatory sequence of porcine elongation factor 1α (pEF1α) gene was selected and isolated for use as a promoter. Activity of the pEF1α promoter was subsequently compared with that of the cytomegalovirus (CMV) promoter, CMV enhancer/chicken β-actin (CAG) promoter, and human EF1α (hEF1α) promoter in different types of pig-derived cells. Comparative analysis of luciferase and mutant human leukocyte antigen class E-F2A-β-2 microglobulin (HLA-E) expression driven by pEF1α, CMV, CAG, and hEF1α promoters revealed the pEF1α promoter mediated comparable expression levels with those of the CAG promoter in porcine ear skin fibroblasts (PEFs) and porcine kidney-15 (PK-15) cells, but lower than those of the CAG promoter in porcine aortic endothelial cells (PAECs). The pEF1α promoter provided long-term stable HLA-E expression in PEFs, but the CAG promoter failed to sustain those levels of expression. For xenogeneic serum-induced cytotoxicity assays, the cells were cultured for several hours in growth medium supplemented with primate serum. Notably, the pEF1α promoter induced significant increases in luciferase and HLA-E expression in response to primate serum in PAECs compared with those driven by the CAG promoter, suggesting the pEF1α promoter could regulate temporal expression of heterologous genes under xenogeneic-cytotoxic conditions. These results suggest the pEF1α promoter may be valuable for development of GE pigs spatiotemporally and stably expressing immunomodulatory genes for xenotransplantation.

Hybrid promoter engineering strategies in Yarrowia lipolytica: isoamyl alcohol production as a test study

BackgroundIn biological cells, promoters drive gene expression by specific binding of RNA polymerase. They determine the starting position, timing and level of gene expression. Therefore, rational fine-tuning of promoters to regulate the expression levels of target genes for optimizing biosynthetic pathways in metabolic engineering has recently become an active area of research.ResultsIn this study, we systematically detected and characterized the common promoter elements in the unconventional yeast Yarrowia lipolytica, and constructed an artificial hybrid promoter library that covers a wide range of promoter strength. The results indicate that the hybrid promoter strength can be fine-tuned by promoter elements, namely, upstream activation sequences (UAS), TATA box and core promoter. Notably, the UASs of Saccharomyces cerevisiae promoters were reported for the first time to be functionally transferred to Y. lipolytica. Subsequently, using the production of a versatile platform chemical isoamyl alcohol as a test study, the hybrid promoter library was applied to optimize the biosynthesis pathway expression in Y. lipolytica. By expressing the key pathway gene, ScARO10, with the promoter library, 1.1–30.3 folds increase in the isoamyl alcohol titer over that of the control strain Y. lipolytica Po1g KU70∆ was achieved. Interestingly, the highest titer increase was attained with a weak promoter PUAS1B4-EXPm to express ScARO10. These results suggest that our hybrid promoter library can be a powerful toolkit for identifying optimum promoters for expressing metabolic pathways in Y. lipolytica.ConclusionWe envision that this promoter engineering strategy and the rationally engineered promoters constructed in this study could also be extended to other non-model fungi for strain improvement.

Polymorphisms in coding and non-coding regions of rabbit (Oryctolagus cuniculus) myogenin (MyoG) gene

In animal breeding, selection based on growth is very often used, as this trait affects the profitability of animal production. Identification of polymorphisms within the genes affecting the growth process seems to be very important. Therefore, we decided to analyse rabbit myogenin (<em>MyoG</em> gene) for potential polymorphic sites and their association with growth and carcass traits in Termond White (TER), Belgian Giant Grey (BGG) and crossbred New Zealand White×Belgian Giant Grey (NZW×BGG) rabbits. We found three single nucleotide polymorphisms (SNPs) – in 5’ upstream sequence g.68679476 C>T, in exon 1 – silent mutation g.68680096 T>C and g.68680097 G>A resulting in change of GTG triplet (valine) into ATG triplet (methionine). Association analysis showed that GG genotype weaning weight was statistically higher compared to GA in TER population (<em>P</em>=0.005), and that the hind parts for GG genotypes were heavier compared to those of GA (<em>P</em>=0.024), but association analysis of dissectible parts showed this was caused by higher bone weight (<em>P</em>=0.015). For g.68679476 C>T in NZW×BGG population, the CC genotypes for fore (678±35) and hind part (615±29) weights were heavier compared to CT (588±16 and 549±13, respectively); moreover, association analysis of dissectible parts showed that weight of dissectible meat in hind part. Unfortunately, we did not find similar associations for other analysed breeds. For g.68679476 C>T in NZWxBGG musculus longissimus lumborum pH leg after 24 h chilling (pH24L) were statistically lower for CC genotypes compared to CT (<em>P</em>=0.027). For g.68680097 G>A in Termond White population L* value on the hind leg after 24 h chilling (L*24H) was higher for GA genotypes compared to GG (<em>P</em>=0.03), while for g.68679476 C>T for musculus longissimus lumborum L* value after 24 h (L*24L) CC genotypes had higher value compared to CT (<em>P</em>=0.016) in BGG population. Moreover, in BGG population CT genotypes had higher weaning weight compared to CC (<em>P</em>=0.018). Our results show that SNPs within the <em>MyoG </em>gene may influence growth traits in some rabbit breeds, but the evolutionary conserved sequence may not be favourable for changes within coding sequences. For a better understanding thereof, additional analysis is required.

A chromodomain protein mediates heterochromatin-directed piRNA expression

PIWI-interacting RNAs (piRNAs) play significant roles in suppressing transposons, maintaining genome integrity, and defending against viral infections. How piRNA source loci are efficiently transcribed is poorly understood. Here, we show that in Caenorhabditis elegans, transcription of piRNA clusters depends on the chromatin microenvironment and a chromodomain-containing protein, UAD-2. piRNA clusters form distinct focus in germline nuclei. We conducted a forward genetic screening and identified UAD-2 that is required for piRNA focus formation. In the absence of histone 3 lysine 27 methylation or proper chromatin-remodeling status, UAD-2 is depleted from the piRNA focus. UAD-2 recruits the upstream sequence transcription complex (USTC), which binds the Ruby motif to piRNA promoters and promotes piRNA generation. Vice versa, the USTC complex is required for UAD-2 to associate with the piRNA focus. Thus, transcription of heterochromatic small RNA source loci relies on coordinated recruitment of both the readers of histone marks and the core transcriptional machinery to DNA.

Genome-wide identification and expression analysis of LBD transcription factor genes in Moso bamboo (Phyllostachys edulis)

BackgroundMoso bamboo, the fastest growing plant on earth, is an important source for income in large areas of Asia, mainly cultivated in China. Lateral organ boundaries domain (LBD) proteins, a family of transcription factors unique to plants, are involved in multiple transcriptional regulatory pathways and play important roles in lateral organ development, pathogen response, secondary growth, and hormone response. The LBD gene family has not previously been characterized in moso bamboo (Phyllostachys edulis).ResultsIn this study, we identified 55 members of the LBD gene family from moso bamboo and found that they were distributed non-uniformly across its 18 chromosomes. Phylogenetic analysis showed that the moso bamboo LBD genes could be divided into two classes. LBDs from the same class share relatively conserved gene structures and sequences encoding similar amino acids. A large number of hormone response–associated cis-regulatory elements were identified in the LBD upstream promoter sequences. Synteny analysis indicated that LBDs in the moso bamboo genome showed greater collinearity with those of O. sativa (rice) and Zea mays (maize) than with those of Arabidopsis and Capsicum annuum (pepper). Numerous segmental duplicates were found in the moso bamboo LBD gene family. Gene expression profiles in four tissues showed that the LBD genes had different spatial expression patterns. qRT–PCR assays with the Short Time-series Expression Miner (STEM) temporal expression analysis demonstrated that six genes (PeLBD20, PeLBD29, PeLBD46, PeLBD10, PeLBD38, and PeLBD06) were consistently up-regulated during the rapid growth and development of bamboo shoots. In addition, 248 candidate target genes that function in a variety of pathways were identified based on consensus LBD binding motifs.ConclusionsIn the current study, we identified 55 members of the moso bamboo transcription factor LBD and characterized for the first time. Based on the short-time sequence expression software and RNA-seq data, the PeLBD gene expression was analyzed. We also investigated the functional annotation of all PeLBDs, including PPI network, GO, and KEGG enrichment based on String database. These results provide a theoretical basis and candidate genes for studying the molecular breeding mechanism of rapid growth of moso bamboo.

Evolution of plant telomerase RNAs: farther to the past, deeper to the roots.

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats.

Fusion Transcripts without Corresponding Cytogenetic Abnormalities in Acute Myeloid Leukemia: Implications for AML Pathogenesis

Chromosomal translocations and inversions that create fusion genes are important causal events in AML. 7 recurrent fusion genes define AML entities in the 2016 revision of the WHO classification of hematologic malignancies. We have developed a comprehensive diagnostic AML platform based on next-generation mRNA sequencing (RNAseq). The adaptable integrated bioinformatic pipeline HAMLET (Human AML Expedited Transcriptomics) analyses the raw data simultaneously for fusion transcripts and for sequence variants and expression levels of genes with relevance to AML pathogenesis, prognosis, or actionability.Poly(A)+ RNA from 100 cryopreserved AML samples with a blast count of 10-99% (median: 75%) from 96 patients (3 diagnosis-relapse pairs, 1 de novo AML and subsequent tAML) was sequenced at a depth of 59x106 paired-end reads/sample with median read length of 126 bp in an ISO17025-accredited Illumina HiSeq 2500 pipeline. Fusion transcripts were called by STAR-Fusion and FusionCatcher. HAMLET detected the recurrent fusion transcripts RUNX1-RUNX1T1, CBFB-MYH11, PML-RARA, ERG-FUS, DEK-NUP214, and KMT2A-MLLT3 plus the variants KMT2A-MLLT4, KMT2A-MLLT1, and KMT2A-MLLT6 with 100% sensitivity and 100% specificity compared to metaphase cytogenetics, FISH, and RT-PCR (Griffioen, ASH 2016).In addition, HAMLET identified fusion transcripts from 54 combinations of genes in 57 AML samples that lacked corresponding chromosomal aberrations by metaphase cytogenetics. 14 “unusual” fusions involved genes located on different chromosomes; 40 fusions represented intrachromosomal events. 8 fusions were detected recurrently, i.e. in at least 2 cases (Fig. A). 32 AML carried at least 2 unusual fusions, with a maximum of 9 different fusions in an AML with a complex karyotype (Fig. B). Unusual fusions were not correlated with classical recurrent AML fusion genes. The 5 most prevalent fusions and an ETV6-LYN fusion were successfully validated in a total of 18 AML by RT-PCR. The only fusion that could not be validated was a YWHAZ-KLF10 fusion (2 cases).The most prevalent RT-PCR-validated fusion (seen in 30 AML cases) has been described in ALL (Roberts, 2012; Safavi, 2015) and was an intrachromosomal TPM4-KLF2 fusion between the tropomyosin 4 and an alternative reading frame (ARF) of the transcription factor Kruppel-like factor 2 genes. Unexpectedly, we detected this fusion in all 3 tested normal bone marrows (BM) but not in EBV-transformed B cells. The hitherto unvalidated C15orf57-CBX3 fusion may be generated by retrotransposition (Schrider, 2013) and was present in 18 AML as well as in 2/3 normal BM and 2/3 EBV B-cells with CBX3 in ARF. The intrachromosomal NCOR2-UBC fusion (9 AML) is listed in the database ChimerDB without further characterization and involves the genes of Nuclear corepressor 2 and Ubiquitin C in the normal reading frame (NRF). The NUP98-NSD1 fusion (3 AML) with NSD1 in NRF is frequently found in hematological malignancies including AML and originates from a cryptic t(5;11) due to telomeric positions of both genes. The PIM3-SCO2 fusion between the serine/threonine kinase PIM3 and an ARF of Cytochrome C oxidase assembly protein can originate from a 0.6 Mb inv(22) in chronic neutrophilic leukemia (Menezes, 2013) and occurs frequently in pediatric AML (Kuhn, 2013). The well-characterized ETV6-LYN fusion (1 AML) fusing the TEL1 oncogene to the SRC family kinase LYN in NRF originates from a complex translocation with inversion involving chromosomes 8 and 12 (Telford, 2016) and has oncogenic activity by activation of STAT5 (Takeda, 2011).Fusion transcripts without recognizable chromosomal abnormalities on metaphase cytogenetics occur frequently in AML and involve actionable targets. The upstream sequence can be fused to the normal or alternative reading frame of the downstream coding sequences and involve potential or validated proto-oncogenes. Apart from cryptic translocations or inversions, trans-splicing is suggested as an alternative mechanism by the presence of some fusion genes in non-malignant hematopoietic cells. Trans-splicing without DNA rearrangement can occur physiologically and alter cellular development towards potential premalignant differentiation (Yuan, 2013). Increasing use of RNAseq technology in routine AML diagnostics, e.g. by HAMLET, will permit to unravel the prognostic significance of unexpected fusion transcripts and to guide novel targeted therapies. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

Defining the Transcriptional Code That Specifies Sinusoidal Endothelial Cells in the HSPC Niche

Hematopoietic stem and progenitor cell (HSPC) niches are specialized microenvironments that support the maintenance and self-renewal of HSPCs. In the fetal liver, bone marrow and spleen, vascular endothelial cells are thought to function as core components of the HSPC niche. In the bone marrow, sinusoidal endothelial cells are known to secrete soluble factors and express cell surface adhesion molecules that promote the recruitment, maintenance and function of HSPCs. Using tomo-seq (RNA tomography) - a technique that involves extracting and barcoding the RNA from sequentially collected cryosections and then performing RNA-seq - we identified ~300 genes enriched in the zebrafish caudal hematopoietic tissue (CHT; the hematopoietic equivalent of the mammalian fetal liver). We also used FACS with tissue-specific fluorescent reporter transgenes to isolate different cell populations from whole embryos and then performed RNA-seq. This included endothelial cells, macrophages, neutrophils and erythrocytes. By overlaying these datasets with the tomo-seq results, we determined the cell types in which many of the CHT-enriched genes were expressed. Unexpectedly, we found a set of ~20 genes that appeared to be enriched specifically in endothelial cells within the CHT. Whole embryo in situ hybridization revealed that these genes had similar spatial and temporal expression patterns, being specifically enriched in sinusoidal endothelial cells in the CHT but absent in other endothelial cells throughout the body. We cloned the upstream promoter regions for several of these genes, fused them to a GFP reporter and then injected them into embryos. Of note, we found that a 5.3 kb sequence upstream of selectin-e (sele ; a gene known to be expressed by bone marrow endothelial cells where it promotes HSPC homing) drove GFP expression specifically in sinusoidal endothelial cells in the CHT. By crossing this sele:GFP transgene to a pan-endothelial marker (kdrl:mCherry), we were able to use FACS to isolate CHT endothelial cells, as well as endothelial cells from tissues outside the CHT. We performed ATAC-seq - a transposon-based method used to detect regions of open chromatin - on these two cell populations and identified 6,710 peaks (regions of open chromatin) that were unique to sele:GFP+; kdrl:mCherry+ CHT niche endothelial cells. Several of these peaks were associated with CHT endothelial genes we had initially identified by tomo-seq, including mrc1a, gpr182, stab1, stab2 and exoc3l2a . Informatic analysis of the sequences under the 6,710 CHT endothelial peaks identified Ets, Sox, Nuclear Hormone and Gata sites as the most enriched transcription factor binding motifs compared to peaks from the non-CHT endothelial cells. These results suggested a transcription factor code involving members of these families might specify sinusoidal endothelial cells in the HSPC niche. Upon closer examination of the sele locus itself, we observed a narrow ATAC-seq peak within the 5.3 kb upstream sequence that was uniquely open in the sele:GFP+; kdrl:mCherry+ CHT niche endothelial cells. Promoter truncation experiments identified a short 158 bp sequence that perfectly aligned with this unique ATAC-seq peak and which was sufficient to drive GFP expression in CHT endothelial cells. This 158 bp sequence contained three Ets, three Sox and one RORA (nuclear hormone receptor) binding sites. Synthetic variants were made in which the different transcription factor binding sites within the 158 bp sequence were disrupted. Injection of these variants driving GFP demonstrated that the Sox and RORA sites were required for expression in CHT endothelial cells. Disruption of the Ets sites, in contrast, led to dramatically enhanced GFP expression, suggesting these sites could be repressive. Together these studies demonstrate that the transcriptional regulation of niche endothelial cells include transcription factors that bind to Ets, Sox, and RORA motifs. The results of this work has important implications for culturing HSPCs in vitro with endothelial cells or for modulating the HSPC niche in the context of human disease or transplantation. DisclosuresZon:Fate, Inc.: Consultancy, Equity Ownership; Marauder, Inc.: Consultancy, Equity Ownership; Scholar Rock, Inc: Consultancy, Equity Ownership; Stemgent: Consultancy.

Distribution and phasing of sequence motifs that facilitate CRISPR adaptation

CRISPR-associated proteins (Cas1 and Cas2) integrate foreign DNA at the "leader" end of CRISPR loci. Several CRISPR leader sequences are reported to contain a binding site for a DNA-bending protein called integration host factor (IHF). IHF-induced DNA bending kinks the leader of type I-E CRISPRs, recruiting an upstream sequence motif that helps dock Cas1-2 onto the first repeat of the CRISPR locus. To determine the prevalence of IHF-directed CRISPR adaptation, we analyzed 15,274 bacterial and archaeal CRISPR leaders. These experiments reveal multiple IHF binding sites and diverse upstream sequence motifs in a subset of the I-C, I-E, I-F, and II-C CRISPR leaders. We identify subtype-specific motifs and show that the phase of these motifs is critical for CRISPR adaptation. Collectively, this work clarifies the prevalence and mechanism(s) of IHF-dependent CRISPR adaptation and suggests that leader sequences and adaptation proteins may coevolve under the selective pressures of foreign genetic elements like plasmids or phages.

EDLm6APred: ensemble deep learning approach for mRNA m6A site prediction

BackgroundAs a common and abundant RNA methylation modification, N6-methyladenosine (m6A) is widely spread in various species' transcriptomes, and it is closely related to the occurrence and development of various life processes and diseases. Thus, accurate identification of m6A methylation sites has become a hot topic. Most biological methods rely on high-throughput sequencing technology, which places great demands on the sequencing library preparation and data analysis. Thus, various machine learning methods have been proposed to extract various types of features based on sequences, then occupied conventional classifiers, such as SVM, RF, etc., for m6A methylation site identification. However, the identification performance relies heavily on the extracted features, which still need to be improved.ResultsThis paper mainly studies feature extraction and classification of m6A methylation sites in a natural language processing way, which manages to organically integrate the feature extraction and classification simultaneously, with consideration of upstream and downstream information of m6A sites. One-hot, RNA word embedding, and Word2vec are adopted to depict sites from the perspectives of the base as well as its upstream and downstream sequence. The BiLSTM model, a well-known sequence model, was then constructed to discriminate the sequences with potential m6A sites. Since the above-mentioned three feature extraction methods focus on different perspectives of m6A sites, an ensemble deep learning predictor (EDLm6APred) was finally constructed for m6A site prediction. Experimental results on human and mouse data sets show that EDLm6APred outperforms the other single ones, indicating that base, upstream, and downstream information are all essential for m6A site detection. Compared with the existing m6A methylation site prediction models without genomic features, EDLm6APred obtains 86.6% of the area under receiver operating curve on the human data sets, indicating the effectiveness of sequential modeling on RNA. To maximize user convenience, a webserver was developed as an implementation of EDLm6APred and made publicly available at www.xjtlu.edu.cn/biologicalsciences/EDLm6APred.ConclusionsOur proposed EDLm6APred method is a reliable predictor for m6A methylation sites.