ZBED evolution: repeated utilization of DNA transposons as regulators of diverse host functions

Association of Polymorphism in Locus of rs274503 (ZBED5/GALNT18) with the Risk of Idiopathic Clubfoot in Chinese Children: An 11-Center Case-Control Study

Background: Idiopathic clubfoot (IC) can be corrected initially using the Ponseti method, but still there is a high recurrence rate. The etiology of IC may include many undetermined genetic and environmental factors. Single nucleotide polymorphism of rs274503 in ZBED5/GALNT18 has been found to be associated with IC in Caucasian children. Therefore, we decided to investigate the association between this polymorphism and the risk of IC in the Chinese population. Methods: We conducted an 11-center case-control study of 516 patients with IC and 661 IC-free children. The rs274503 (A>G) polymorphism was genotyped using TaqMan. Odds ratios (ORs) and adjusted ORs, as well as 95% confidence intervals (CIs) and adjusted 95% CIs, were calculated to explore the association between rs274503 polymorphism and IC risk. Results: G of rs274503 was found to be associated with increased IC risk (AG vs. AA: adjusted OR = 1.40, 95% CI = 1.03-1.92, p = 0.0327; and GG/AG vs. AA: adjusted OR = 1.38, 95% CI = 1.02-1.87, p = 0.0357) after adjusting for age and sex. Furthermore, the risk effect of rs274503 GG/AG with IC was observed in patients with bilateral feet (adjusted OR = 1.68, 95% CI = 1.12-2.54, p = 0.0133), while AA in nonrelapsed groups (OR = 0.70, 95% CI = 0.53-0.92, p = 0.0095) in the stratified analysis. However, the association was not significant in the recessive model of G (GG vs. AA/AG: adjusted OR = 1.06, 95% CI = 0.44-2.58, p = 0.8906). Conclusions: The rs274503 polymorphism is associated with the risk of clubfoot occurrence. G of rs274503 appeared to be a risk factor of IC as it may increase the bilateral case rate. However, further studies are required to confirm these findings.

Perspectives on Codebook: sequence specificity of uncharacterized human transcription factors

We describe an effort ("Codebook") to determine the sequence specificity of 332 putative and largely uncharacterized human transcription factors (TFs), as well as 61 control TFs. Nearly 5,000 independent experiments across multiple in vitro and in vivo assays produced motifs for just over half of the putative TFs analyzed (177, or 53%), of which most are unique to a single TF. The data highlight the extensive contribution of transposable elements to TF evolution, both in cis and trans, and identify tens of thousands of conserved, base-level binding sites in the human genome. The use of multiple assays provides an unprecedented opportunity to benchmark and analyze TF sequence specificity, function, and evolution, as further explored in accompanying manuscripts. 1,421 human TFs are now associated with a DNA binding motif. Extrapolation from the Codebook benchmarking, however, suggests that many of the currently known binding motifs for well-studied TFs may inaccurately describe the TF's true sequence preferences.

An Aegilops longissima NLR protein with integrated CC-BED module mediates resistance to wheat powdery mildew

Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), reduces wheat yields and grain quality, thus posing a significant threat to global food security. Wild relatives of wheat serve as valuable resources for resistance to powdery mildew. Here, the powdery mildew resistance gene Pm6Sl is cloned from the wild wheat species Aegilops longissima. It encodes a nucleotide-binding leucine-rich repeat (NLR) protein featuring a CC-BED module formed by a zinc finger BED (Znf-BED) domain integrated into the coiled-coil (CC) domain. The function of Pm6Sl is validated via mutagenesis, gene silencing, and transgenic assays. In addition, we develop a resistant germplasm harbouring Pm6Sl in a very small segment with no linkage drag along with the diagnostic gene marker pm6sl-1 to facilitate Pm6Sl deployment in wheat breeding programs. The cloning of Pm6Sl, a resistance gene with BED-NLR architecture, will increase our understanding of the molecular mechanisms underlying BED-NLR-mediated resistance to various pathogens.

Genome-wide identification and characterization of SLEEPER, a transposon-derived gene family and their expression pattern in Brassica napus L

BACKGROUND

The transposons of the hAT superfamily are the most widespread transposons ever known. SLEEPER genes encode domesticated transposases from the hAT superfamily, which may have lost their transposable functions during long-term evolution and transformed into host proteins that regulate plant growth and development.

RESULTS

This study identified 162 members of the SLEEPER gene family from Brassica napus. These members are widely distributed on 19 chromosomes, mainly in the Cn subgenome, and have promoters with various cis-acting elements related to hormone regulation, abiotic stress, and growth and development regulation. Most of the genes in this family contain similar conserved domains and motifs, and the closer the genes are distributed on evolutionary branches, the more similar their structures are. Transcriptome sequencing performed on tissues at different growth stages from B. napus line 3529 indicated that these genes had different expression patterns, and nearly half of the genes were not detectably expressed in all samples.

CONCLUSIONS

This study investigated the gene structure, expression patterns, evolutionary features, and gene localization of the SLEEPER family members to confirm the significance of these genes in the growth of B. napus, providing a reference for the study of transposon domestication and outstanding genetic resources for the genetic improvement of B. napus.

The hAT family hopper transposon exists as highly similar yet discontinuous elements in the Bactrocera tephritid fly genus

The hAT family transposable element, hopper, was originally discovered as a defective 3120-bp full-length element in a wild-type strain of the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), and subsequently a functional 3131-bp element, hopperBdwe, was isolated from a white eye mutant strain. The latter study showed that closely related elements exist in melonfly, Zeugodacus cucurbitae (Coquillett) (Diptera: Tephritidae), a closely related subgenus, suggesting that hopper could have a widespread presence in the Bactrocera genus. To further understand the distribution of hopper within and beyond the B. dorsalis species complex, primer pairs from hopperBdwe and its adjacent genomic insertion site were used to survey the presence and relatedness of hopper in five species within the complex and four species beyond the complex. Based on sequence identity of a 1.94 kb internal nucleotide sequence, the closest relationships were with mutated elements from B. dorsalis s.s. and species synonymized with B. dorsalis including B. papayae, B. philippinensis and B. invadens, ranging in identity between 88.4% and 99.5%. Notably, Bactrocera carambolae (Drew & Hancock) (Diptera: Tephritidae), which is most closely related to B. dorsalis beyond the synonymized species, shared hopper identities of 97.3%-99.5%. Beyond the B. dorsalis complex, Z. cucurbitae, Bactrocera tryoni (Froggatt) (Diptera: Tephritidae) and Bactrocera zonata (Saunders) (Diptera: Tephritidae) shared identities of 83.1%-97.1%, while hopper was absent from the Bactrocera oleae (Gmelin) (Diptera: Tephritidae) strain tested. While the functional autonomous hopperBdwe element was not detected in these species, another closely related hopper element isolated from a B. dorsalis genetic sexing strain has an uninterrupted transposase open reading frame. The discontinuous presence of hopper in the Bactrocera genus has implications for its use for genomic manipulation and understanding the phylogenetic relationship of these species.

Parallel Evolution of Transcription Factors in Basal Metazoans

Transcription factors (TFs) play a pivotal role as regulators of gene expression, orchestrating the formation and maintenance of diverse animal body plans and innovations. However, the precise contributions of TFs and the underlying mechanisms driving the origin of basal metazoan body plans, particularly in ctenophores, remain elusive. Here, we present a comprehensive catalog of TFs in 2 ctenophore species, Pleurobrachia bachei and Mnemiopsis leidyi, revealing 428 and 418 TFs in their respective genomes. In contrast, morphologically simpler metazoans have a reduced TF representation compared to ctenophores, cnidarians, and bilaterians: the sponge Amphimedon encodes 277 TFs, and the placozoan Trichoplax adhaerens encodes 274 TFs. The emergence of complex ctenophore tissues and organs coincides with significant lineage-specific diversification of the zinc finger C2H2 (ZF-C2H2) and homeobox superfamilies of TFs. Notable, the lineages leading to Amphimedon and Trichoplax exhibit independent expansions of leucine zipper (BZIP) TFs. Some lineage-specific TFs may have evolved through the domestication of mobile elements, thereby supporting alternative mechanisms of parallel TF evolution and body plan diversification across the Metazoa.

PiggyBac Transposon Mining in the Small Genomes of Animals

TEs, including DNA transposons, are major contributors of genome expansions, and have played a very significant role in shaping the evolution of animal genomes, due to their capacity to jump from one genomic position to the other. In this study, we investigated the evolution landscapes of PB transposons, including their distribution, diversity, activity and structure organization in 79 species of small (compact) genomes of animals comprising both vertebrate and invertebrates. Overall, 212 PB transposon types were detected from almost half (37) of the total number of the small genome species (79) investigated. The detected PB transposon types, which were unevenly distributed in various genera and phyla, have been classified into seven distinct clades or families with good bootstrap support (>80%). The PB transposon types that were identified have a length ranging from 1.23 kb to 9.51 kb. They encode transposases of approximately ≥500 amino acids in length, and possess terminal inverted repeats (TIRs) ranging from 4 bp to 24 bp. Though some of the transposon types have long TIRs (528 bp), they still maintain the consistent and reliable 4 bp target site duplication (TSD) of TTAA. However, PiggyBac-2_Cvir transposon originating from the Crassostrea virginica species exhibits a unique TSD of TATG. The TIRs of the transposons in all the seven families display high divergence, with a highly conserved 5' end motif. The core transposase domains (DDD) were better conserved among the seven different families compared to the other protein domains, which were less prevalent in the vertebrate genome. The divergent evolution dynamics analysis also indicated that the majority of the PB transposon types identified in this study are either relatively young or old, with some being active. Additionally, numerous invasions of PB transposons were found in the genomes of both vertebrate and invertebrate animals. The data reveals that the PB superfamily is widely distributed in these species. PB transposons exhibit high diversity and activity in the small genomes of animals, and might play a crucial role in shaping the evolution of these small genomes of animals.

Zinc-finger BED domains drive the formation of the active Hermes transpososome by asymmetric DNA binding

The Hermes DNA transposon is a member of the eukaryotic hAT superfamily, and its transposase forms a ring-shaped tetramer of dimers. Our investigation, combining biochemical, crystallography and cryo-electron microscopy, and in-cell assays, shows that the full-length Hermes octamer extensively interacts with its transposon left-end through multiple BED domains of three Hermes protomers contributed by three dimers explaining the role of the unusual higher-order assembly. By contrast, the right-end is bound to no BED domains at all. Thus, this work supports a model in which Hermes multimerizes to gather enough BED domains to find its left-end among the abundant genomic DNA, facilitating the subsequent interaction with the right-end.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Genomic Structure, Protein Character, Phylogenic Implication, and Embryonic Expression Pattern of a Zebrafish New Member of Zinc Finger BED-Type Gene Family

We reported a new member of the C₂H₂-zinc-finger BED-type (ZBED) protein family found in zebrafish (Danio rerio). It was previously assigned as an uncharacterized protein LOC569044 encoded by the Zgc:161969 gene, the transcripts of which were highly expressed in the CNS after the spinal cord injury of zebrafish. As such, this novel gene deserves a more detailed investigation. The 2.79-kb Zgc:161969 gene contains one intron located on Chromosome 6 at 16,468,776-16,475,879 in the zebrafish genome encoding a 630-aa protein LOC569044. This protein is composed of a DNA-binding BED domain, which is highly conserved among the ZBED protein family, and a catalytic domain consisting of an α-helix structure and an hAT dimerization region. Phylogenetic analysis revealed the LOC569044 protein to be clustered into the monophyletic clade of the ZBED protein family of golden fish. Specifically, the LOC569044 protein was classified as closely related to the monophyletic clades of zebrafish ZBED4-like isoforms and ZBED isoform 2. Furthermore, Zgc:161969 transcripts represented maternal inheritance, expressed in the brain and eyes at early developmental stages and in the telencephalon ventricular zone at late developmental stages. After characterizing the LOC569044 protein encoded by the Zgc:161969 gene, it was identified as a new member of the zebrafish ZBED protein family, named the ZBEDX protein.

Transposon-derived transcription factors across metazoans

Transposable elements (TE) could serve as sources of new transcription factors (TFs) in plants and some other model species, but such evidence is lacking for most animal lineages. Here, we discovered multiple independent co-options of TEs to generate 788 TFs across Metazoa, including all early-branching animal lineages. Six of ten superfamilies of DNA transposon-derived conserved TF families (ZBED, CENPB, FHY3, HTH-Psq, THAP, and FLYWCH) were identified across nine phyla encompassing the entire metazoan phylogeny. The most extensive convergent domestication of potentially TE-derived TFs occurred in the hydroid polyps, polychaete worms, cephalopods, oysters, and sea slugs. Phylogenetic reconstructions showed species-specific clustering and lineage-specific expansion; none of the identified TE-derived TFs revealed homologs in their closest neighbors. Together, our study established a framework for categorizing TE-derived TFs and informing the origins of novel genes across phyla.

TGF-β generates a population of cancer cells residing in G1 phase with high motility and metastatic potential via KRTAP2-3

Transforming growth factor β (TGF-β) increases epithelial cancer cell migration and metastasis by inducing epithelial-mesenchymal transition (EMT). TGF-β also inhibits cell proliferation by inducing G1 phase cell-cycle arrest. However, the correlation between these tumor-promoting and -suppressing effects remains unclear. Here, we show that TGF-β confers higher motility and metastatic ability to oral cancer cells in G1 phase. Mechanistically, keratin-associated protein 2-3 (KRTAP2-3) is a regulator of these dual effects of TGF-β, and its expression is correlated with tumor progression in patients with head and neck cancer and migratory and metastatic potentials of oral cancer cells. Furthermore, single-cell RNA sequencing reveals that TGF-β generates two populations of mesenchymal cancer cells with differential cell-cycle status through two distinctive EMT pathways mediated by Slug/HMGA2 and KRTAP2-3. Thus, TGF-β-induced KRTAP2-3 orchestrates cancer cell proliferation and migration by inducing EMT, suggesting motile cancer cells arrested in G1 phase as a target to suppress metastasis.

Taming, Domestication and Exaptation: Trajectories of Transposable Elements in Genomes

During evolution, several types of sequences pass through genomes. Along with mutations and internal genetic tinkering, they are a useful source of genetic variability for adaptation and evolution. Most of these sequences are acquired by horizontal transfers (HT), but some of them may come from the genomes themselves. If they are not lost or eliminated quickly, they can be tamed, domesticated, or even exapted. Each of these processes results from a series of events, depending on the interactions between these sequences and the host genomes, but also on environmental constraints, through their impact on individuals or population fitness. After a brief reminder of the characteristics of each of these states (taming, domestication, exaptation), the evolutionary trajectories of these new or acquired sequences will be presented and discussed, emphasizing that they are not totally independent insofar as the first can constitute a step towards the second, and the second is another step towards the third.

Diverse Eukaryotic CGG-Binding Proteins Produced by Independent Domestications of hAT Transposons

The human transcription factor (TF) CGGBP1 (CGG-binding protein) is conserved only in amniotes and is believed to derive from the zf-BED and Hermes transposase DNA-binding domains (DBDs) of a hAT DNA transposon. Here, we show that sequence-specific DNA-binding proteins with this bipartite domain structure have resulted from dozens of independent hAT domestications in different eukaryotic lineages. CGGBPs display a wide range of sequence specificity, usually including preferences for CGG or CGC trinucleotides, whereas some bind AT-rich motifs. The CGGBPs are almost entirely nonsyntenic, and their protein sequences, DNA-binding motifs, and patterns of presence or absence in genomes are uncharacteristic of ancestry via speciation. At least eight CGGBPs in the coelacanth Latimeria chalumnae bind distinct motifs, and the expression of the corresponding genes varies considerably across tissues, suggesting tissue-restricted function.

A novel genomic model for predicting the likelihood of delayed graft function in DCD kidney transplantation

Background

The high incidence of delayed graft function (DGF) following kidney transplantation with donation after cardiac death allografts (DCD-KT) poses great challenges to transplant clinicians. This study aimed to explore the DGF-related biomarkers and establish a genomic model for DGF prediction specific to DCD KT.

Methods

By data mining a public dataset (GSE43974), the key DGF-related genes in DCD kidney biopsies taken after short-time reperfusion (45–60 min) were identified by differential expression analysis and a LASSO-penalized logistic regression model. Their coefficients for modeling were calculated by multivariate logistic regression. Receiver operating characteristic curves and a nomogram were generated to evaluate its predictive ability for DGF occurrence. Gene set enrichment analysis (GSEA) was performed to explore biological pathways underlying DGF in DCD KT.

Results

Five key DGF-related genes (CHST3, GOLPH3, ZBED5, AKR1C4, and ERRFI1) were first identified, all of which displayed good discrimination for DGF occurrence after DCD KT (all P<0.05). A five-mRNA-based risk score was further established and showed excellent predictive ability (AUC =0.9708, P<0.0001), which was obviously higher than that of the five genes alone. Eight DGF-related biological pathways in DCD kidneys, such as “arachidonic acid metabolism”, “lysosome”, “proximal tubule bicarbonate reclamation”, “glutathione metabolism”, were identified by GSEA (all P<0.05). Moreover, a convenient and visual nomogram based on the genomic risk score was also constructed and displayed high accuracy for DGF prediction specific to DCD KT.

Conclusions

The novel genomic model may effectively predict the likelihood of DGF immediately after DCD KT or even prior to transplantation in the context of normothermic machine perfusion in the future.

Recurrent evolution of vertebrate transcription factors by transposase capture

Genes with novel cellular functions may evolve through exon shuffling, which can assemble novel protein architectures. Here, we show that DNA transposons provide a recurrent supply of materials to assemble protein-coding genes through exon shuffling. We find that transposase domains have been captured-primarily via alternative splicing-to form fusion proteins at least 94 times independently over the course of ~350 million years of tetrapod evolution. We find an excess of transposase DNA binding domains fused to host regulatory domains, especially the Krüppel-associated box (KRAB) domain, and identify four independently evolved KRAB-transposase fusion proteins repressing gene expression in a sequence-specific fashion. The bat-specific KRABINER fusion protein binds its cognate transposons genome-wide and controls a network of genes and cis-regulatory elements. These results illustrate how a transcription factor and its binding sites can emerge.

De Novo Transcriptome Profiling of Brain Tissue from the Annual Killifish Nothobranchius guentheri

Nothobranchius is a genus of small annual killifish found in Africa. Due to the relatively short lifespan, as well as easy breeding and care, Nothobranchius fish are becoming widely used as a vertebrate model system. Studying the genome and transcriptome of these fish is essential for advancing the field. In this study, we performed de novo transcriptome assembly of brain tissues from Nothobranchius guentheri using Trinity. Annotation of 104,271 potential genes (with transcripts longer than 500 bp) was carried out; for 24,967 genes (53,654 transcripts), in which at least one GO annotation was derived. We also analyzed the effect of a long-term food supplement with Torin 2, second-generation ATP-competitive inhibitor of mTOR, on the gene expression changes in brain tissue of adult N. guentheri. Overall, 1491 genes in females and 249 genes in males were differently expressed under Torin 2-supplemented diet. According to the Gene Set Enrichment Analysis (GSEA), the majority of identified genes were predominantly involved in the regulation of metabolic process, dendritic spine maintenance, circadian rhythms, retrotransposition, and immune response. Thus, we have provided the first transcriptome assembly and assessed the differential gene expression in response to exposure to Torin 2, which allow a better understanding of molecular changes in the brain tissues of adult fish in the mTOR pathway inhibition.

Transposable element-derived sequences in vertebrate development

Transposable elements (TEs) are major components of all vertebrate genomes that can cause deleterious insertions and genomic instability. However, depending on the specific genomic context of their insertion site, TE sequences can sometimes get positively selected, leading to what are called "exaptation" events. TE sequence exaptation constitutes an important source of novelties for gene, genome and organism evolution, giving rise to new regulatory sequences, protein-coding exons/genes and non-coding RNAs, which can play various roles beneficial to the host. In this review, we focus on the development of vertebrates, which present many derived traits such as bones, adaptive immunity and a complex brain. We illustrate how TE-derived sequences have given rise to developmental innovations in vertebrates and how they thereby contributed to the evolutionary success of this lineage.

ZBED6 binding motifs correlate with endogenous retroviruses and Syncytin genes

Retroviruses have infiltrated vertebrate germlines for millions of years as inherited endogenous retroviruses (ERVs). Mammalian genomes host large numbers of ERVs and transposable elements (TEs), including retrotransposons and DNA transposons, that contribute to genomic innovation and evolution as coopted genes and regulators of diverse functions. To explore features distinguishing coopted ERVs and TEs from other integrations, we focus on the potential role of ZBED6 and repeated ERV domestication as repurposed Syncytin genes. The placental mammal-specific ZBED6 is a DNA transposon-derived transcription regulator and we demonstrate that its binding motifs are associated with distinct Syncytins and that ZBED6 binding motifs are 2- to 3-fold more frequent in ERVs than in flanking DNA. Our observations suggest that ZBED6 could contribute an extended regulatory role of genomic expression, utilizing ERVs as platforms for genomic innovation and evolution.

Characterization of an eutherian gene cluster generated after transposon domestication identifies Bex3 as relevant for advanced neurological functions

BACKGROUND

One of the most unusual sources of phylogenetically restricted genes is the molecular domestication of transposable elements into a host genome as functional genes. Although these kinds of events are sometimes at the core of key macroevolutionary changes, their origin and organismal function are generally poorly understood.

RESULTS

Here, we identify several previously unreported transposable element domestication events in the human and mouse genomes. Among them, we find a remarkable molecular domestication that gave rise to a multigenic family in placental mammals, the Bex/Tceal gene cluster. These genes, which act as hub proteins within diverse signaling pathways, have been associated with neurological features of human patients carrying genomic microdeletions in chromosome X. The Bex/Tceal genes display neural-enriched patterns and are differentially expressed in human neurological disorders, such as autism and schizophrenia. Two different murine alleles of the cluster member Bex3 display morphological and physiopathological brain modifications, such as reduced interneuron number and hippocampal electrophysiological imbalance, alterations that translate into distinct behavioral phenotypes.

CONCLUSIONS

We provide an in-depth understanding of the emergence of a gene cluster that originated by transposon domestication and gene duplication at the origin of placental mammals, an evolutionary process that transformed a non-functional transposon sequence into novel components of the eutherian genome. These genes were integrated into existing signaling pathways involved in the development, maintenance, and function of the CNS in eutherians. At least one of its members, Bex3, is relevant for higher brain functions in placental mammals and may be involved in human neurological disorders.

ZBED1/DREF: A transcription factor that regulates cell proliferation

Maintenance of genomic diversity is critically dependent on gene regulation at the transcriptional level. This occurs via the interaction of regulatory DNA sequence motifs with DNA-binding transcription factors. The zinc finger, BED-type (ZBED) gene family contains major DNA-binding motifs present in human transcriptional factors. It encodes proteins that present markedly diverse regulatory functions. ZBED1 has similar structural and functional properties to its Drosophila homolog DNA replication-related element-binding factor (DREF) and plays a critical role in the regulation of transcription. ZBED1 regulates the expression of several genes associated with cell proliferation, including cell cycle regulation, chromatin remodeling and protein metabolism, and some genes associated with apoptosis and differentiation. In the present review, the origin, structure and functional role of ZBED1 were comprehensively assessed. In addition, the similarities and differences between ZBED1 and its Drosophila homolog DREF were highlighted, and future research directions, particularly in the area of clinical cancer, were discussed.

ZBED2 is an antagonist of interferon regulatory factor 1 and modifies cell identity in pancreatic cancer

Lineage plasticity is a prominent feature of pancreatic ductal adenocarcinoma (PDA) cells, which can occur via deregulation of lineage-specifying transcription factors. Here, we show that the zinc finger protein ZBED2 is aberrantly expressed in PDA and alters tumor cell identity in this disease. Unexpectedly, our epigenomic experiments reveal that ZBED2 is a sequence-specific transcriptional repressor of IFN-stimulated genes, which occurs through antagonism of IFN regulatory factor 1 (IRF1)-mediated transcriptional activation at cooccupied promoter elements. Consequently, ZBED2 attenuates the transcriptional output and growth arrest phenotypes downstream of IFN signaling in multiple PDA cell line models. We also found that ZBED2 is preferentially expressed in the squamous molecular subtype of human PDA, in association with inferior patient survival outcomes. Consistent with this observation, we show that ZBED2 can repress the pancreatic progenitor transcriptional program, enhance motility, and promote invasion in PDA cells. Collectively, our findings suggest that high ZBED2 expression is acquired during PDA progression to suppress the IFN response pathway and to promote lineage plasticity in this disease.

Plant NLRs with Integrated Domains: Unity Makes Strength

Integration of unusual domains in plant immune receptors is a widespread mechanism of NLR diversification enabling specific pathogen detection.

Candidate genes and pathways downstream of PAX8 involved in ovarian high-grade serous carcinoma

Understanding the biology and molecular pathogenesis of ovarian epithelial cancer (EOC) is key to developing improved diagnostic and prognostic indicators and effective therapies. Although research has traditionally focused on the hypothesis that high-grade serous carcinoma (HGSC) arises from the ovarian surface epithelium (OSE), recent studies suggest that additional sites of origin exist and a substantial proportion of cases may arise from precursor lesions located in the Fallopian tubal epithelium (FTE). In FTE cells, the transcription factor PAX8 is a marker of the secretory cell lineage and its expression is retained in 96% of EOC. We have recently reported that PAX8 is involved in the tumorigenic phenotype of ovarian cancer cells. In this study, to uncover genes and pathways downstream of PAX8 involved in ovarian carcinoma we have determined the molecular profiles of ovarian cancer cells and in parallel of Fallopian tube epithelial cells by means of a silencing approach followed by an RNA-seq analysis. Interestingly, we highlighted the involvement of pathways like WNT signaling, epithelial-mesenchymal transition, p53 and apoptosis. We believe that our analysis has led to the identification of candidate genes and pathways regulated by PAX8 that could be additional targets for the therapy of ovarian carcinoma.

Transcriptional activity of PIF and Pong-like Class II transposable elements in Triticeae

Background

Transposable elements are major contributors to genome size and variability, accounting for approximately 70–80% of the maize, barley, and wheat genomes. PIF and Pong-like elements belong to two closely-related element families within the PIF/Harbinger superfamily of Class II (DNA) transposons. Both elements contain two open reading frames; one encodes a transposase (ORF2) that catalyzes transposition of the functional elements and their related non-autonomous elements, while the function of the second is still debated. In this work, we surveyed for PIF- and Pong-related transcriptional activity in 13 diploid Triticeae species, all of which have been previously shown to harbor extensive within-genome diversity of both groups of elements.

Results

The results revealed that PIF elements have considerable transcriptional activity in Triticeae, suggesting that they can escape the initial levels of plant cell control and are regulated at the post-transcriptional level. Phylogenetic analysis of 156 PIF cDNA transposase fragments along with 240 genomic partial transposase sequences showed that most, if not all, PIF clades are transcriptionally competent, and that multiple transposases coexisting within a single genome have the potential to act simultaneously. In contrast, we did not detect any transcriptional activity of Pong elements in any sample.

Conclusions

The lack of Pong element transcription shows that even closely related transposon families can exhibit wide variation in their transposase transcriptional activity within the same genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1028-6) contains supplementary material, which is available to authorized users.

Influence of Irradiated Peripheral Blood Mononuclear Cells on Both Ex Vivo Proliferation of Human Natural Killer Cells and Change in Cellular Property

Clinical studies with adoptive immunotherapy using allogeneic natural killer (NK) cells showed feasibility, but also limitation regarding the transfused absolute cell numbers. First promising results with peripheral blood mononuclear cells (PBMCs) as feeder cells to improve the final cell number need further optimization and investigation of the unknown controlling mechanism in the cross-talk to NK cells. We investigated the influence of irradiated autologous PBMCs to boost NK cell proliferation in the presence of OKT3 and IL-2. Our findings demonstrate a requirement for receptor-ligand interactions between feeders and NK cells to produce soluble factors that can sustain NK cell proliferation. Thus, both physical contact between feeder and NK cells, and soluble factors produced in consequence, are required to fully enhance NK cell ex vivo proliferation. This occurred with an indispensable role of the cross-talk between T cells, monocytes, and NK cells, while B cells had no further influence in supporting NK cell proliferation under these co-culture conditions. Moreover, gene expression analysis of highly proliferating and non-proliferating NK cells revealed important phenotypic changes on 5-day cultured NK cells. Actively proliferating NK cells have reduced Siglec-7 and -9 expression compared with non-proliferating and resting NK cells (day 0), independently of the presence of feeder cells. Interestingly, proliferating NK cells cultured with feeder cells contained increased frequencies of cells expressing RANKL, B7-H3, and HLA class II molecules, particularly HLA-DR, compared with resting NK cells or expanded with IL-2 only. A subset of HLA-DR expressing NK cells, co-expressing RANKL, and B7-H3 corresponded to the most proliferative population under the established co-culture conditions. Our results highlight the importance of the crosstalk between T cells, monocytes, and NK cells in autologous feeder cell-based ex vivo NK cell expansion protocols, and reveal the appearance of a highly proliferative subpopulation of NK cells (HLA-DR⁺RANKL⁺B7-H3⁺) with promising characteristics to extend the therapeutic potential of NK cells.

The LncRNA ZBED3-AS1 induces chondrogenesis of human synovial fluid mesenchymal stem cells

Human synovial fluid-derived mesenchymal stem cells (SFMSCs) have great potential for cartilage induction and are promising for cell-based strategies for articular cartilage repair. Many long non-coding RNAs (lncRNAs) regulate chondrogenesis of MSCs. We hypothesized that the divergent lncRNA ZBED3-AS1, which binds locally to chromatin, could promote the expression of zbed3, a novel Axin-interacting protein that activates Wnt/β-catenin signaling, involved in chondrogenesis. However, the function of ZBED3-AS1 in SFMSCs is unclear. In this study, the expression, biological function, and roles of ZBED3-AS1 in SFMSC chondrogenesis were examined by multilineage differentiation, flow cytometry, and gain-of-function studies. We found that ZBED3-AS1 promotes chondrogenesis. Furthermore, ZBED3-AS1 could directly increase zbed3 expression. Finally, the wnt-inhibitor DKK1 could reverse the stimulatory effect of ZBED3-AS1 on chondrogenesis. These findings demonstrate the role of a new lncRNA, ZBED3-AS1, in SFMSC chondrogenesis and may improve osteoarthritis treatment.

Effector Mimics and Integrated Decoys, the Never-Ending Arms Race between Rice and Xanthomonas oryzae

Plants are constantly challenged by a wide range of pathogens and have therefore evolved an array of mechanisms to defend against them. In response to these defense systems, pathogens have evolved strategies to avoid recognition and suppress plant defenses (Brown and Tellier, 2011). Three recent reports dealing with the resistance of rice to Xanthomonas oryzae have added a new twist to our understanding of this fascinating co-evolutionary arms race (Ji et al., 2016; Read et al., 2016; Triplett et al., 2016). They show that pathogens also develop sophisticated effector mimics to trick recognition.

TE-domestication and horizontal transfer in a putative Nef-AP1mu mimic of HLA-A cytoplasmic domain re-trafficking

Genes of the major histocompatibility complex (MHC; also called HLA in human) are polymorphic elements in the genomes of sharks to humans. Class-I and class-II MHC loci appear responsible for much of the genetic linkage to myriad disease states via the capacity to bind short (~8-15 a.a.) peptides of a given pathogen's proteome, or in some cases, the altered proteomes of cancerous cells, and even (in autoimmunity) certain nominal 'self' peptides (Janeway, 2004).(1) Unfortunately, little is known about how the canonical structure of the MHC-I/-II peptide-presenting gene evolved, particularly since beyond ~500 Mya (sharks) no paralogs exist.(2,3) We previously reported that HLA-A isotype alleles with the α1-helix, R65 motif, are wide-spread in phylogeny, but that the α 2-helix, H151R motif, has apparently segregated out of most species. Surprisingly, an uncharacterized orf in T. syrichta (Loc-103275158) encoded R151, but within a truncated A-23 like gene containing 5'- and 3'- footprints of the transposon (TE), tigger-1; the extant tarsier A-23 allele is totally missing exon-3 and part-of exon-4; together, suggesting TE-mediated inactivation of an intact/ancestral A-23 allele (Murray, 2015a).(4) The unique Loc-103275158 orf encodes a putative 15-exon transcript with no apparent paralogs throughout phylogeny. However, an HLA-A11 like gene in M. leucophaeus with a shortened C-terminal domain, and an HLA-A like orf in C. atys with two linked α1/α2/α3 domains, both contain a second transmembrane segment, which is conserved in Loc-103275158. Thus, we could model the putative protein with its Nef-like tail domain docked to its MHC-I like α3 domain (i.e., on the same side of a membrane). This modeled tertiary structure is strikingly similar to the solved structure of the Nef:MHC-I CD:AP1mu transporter (Jia, 2012).(5) Nef:AP1mu binds the CD of MHC-I in trafficking MHC-I away from the trans-golgi and into the endocytic pathway in HIV-1 infected cells. The CD loop of the Loc-103275158 provisional protein conserved the nominal MHC-I CD tyrosine phosphorylation site, and it has an N-terminal SH3 domain that we docked in one conformation to its internal Nef-like domain. Here, we suggest that phosphorylation of the protein's CD-loop signals an exchange between the internal Nef-like domain and a lentiviral-Nef for binding the N-terminal SH3 domain - freeing the Nef-like domain to bind MHC-I CD. Since the 5'-tigger sequence encodes part of the pseudo α1/α2 MHC-I domain, and the 3'-tigger part of the Nef-like domain, we speculate that transposition proceeded phylogenetically disparate horizontal transfers, involving adjacent 5'- and 3'- parasitic footprints, which we also found in the Loc-103275158 orf.

hAT Transposable Elements

hAT transposons are ancient in their origin and they are widespread across eukaryote kingdoms. They can be present in large numbers in many genomes. However, only a few active forms of these elements have so far been discovered indicating that, like all transposable elements, there is selective pressure to inactivate them. Nonetheless, there have been sufficient numbers of active hAT elements and their transposases characterized that permit an analysis of their structure and function. This review analyzes these and provides a comparison with the several domesticated hAT genes discovered in eukaryote genomes. Active hAT transposons have also been developed as genetic tools and understanding how these may be optimally utilized in new hosts will depend, in part, on understanding the basis of their function in genomes.

Multidisciplinary investigation links backward-speech trait and working memory through genetic mutation

Case studies of unusual traits can provide unique snapshots of the effects of modified systems. In this study, we report on an individual from a Serbian family with the ability to rapidly, accurately and voluntarily speak backwards. We consider psychological, neural and genetic correlates of this trait to identify specific relevant neural mechanisms and new molecular pathways for working memory and speech-related tasks. EEG data suggest that the effect of word reversal precedes semantic integration of visually presented backward-words, and that event-related potentials above the frontal lobe are affected by both word reversal and the maintenance of backward-words in working memory. fMRI revealed that the left fusiform gyrus may facilitate the production of backward-speech. Exome sequencing identified three novel coding variants of potential significance in the RIC3, RIPK1 and ZBED5 genes. Taken together, our data suggest that, in this individual, the ability to speak backwards is afforded by an extraordinary working memory capacity. We hypothesise that this is served by cholinergic projections from the basal forebrain to the frontal cortex and supported by visual semantic loops within the left fusiform gyrus and that these neural processes may be mediated by a genetic mutation in RIC3; a chaperone for nicotinic acetylcholine receptors.

Transcriptional modulator ZBED6 affects cell cycle and growth of human colorectal cancer cells

The transcription factor ZBED6 (zinc finger, BED-type containing 6) is a repressor of IGF2 whose action impacts development, cell proliferation, and growth in placental mammals. In human colorectal cancers, IGF2 overexpression is mutually exclusive with somatic mutations in PI3K signaling components, providing genetic evidence for a role in the PI3K pathway. To understand the role of ZBED6 in tumorigenesis, we engineered and validated somatic cell ZBED6 knock-outs in the human colorectal cancer cell lines RKO and HCT116. Ablation of ZBED6 affected the cell cycle and led to increased growth rate in RKO cells but reduced growth in HCT116 cells. This striking difference was reflected in the transcriptome analyses, which revealed enrichment of cell-cycle-related processes among differentially expressed genes in both cell lines, but the direction of change often differed between the cell lines. ChIP sequencing analyses displayed enrichment of ZBED6 binding at genes up-regulated in ZBED6-knockout clones, consistent with the view that ZBED6 modulates gene expression primarily by repressing transcription. Ten differentially expressed genes were identified as putative direct gene targets, and their down-regulation by ZBED6 was validated experimentally. Eight of these genes were linked to the Wnt, Hippo, TGF-β, EGF receptor, or PI3K pathways, all involved in colorectal cancer development. The results of this study show that the effect of ZBED6 on tumor development depends on the genetic background and the transcriptional state of its target genes.

Discovery of novel genes derived from transposable elements using integrative genomic analysis

Complex eukaryotes contain millions of transposable elements (TEs), comprising large fractions of their nuclear genomes. TEs consist of structural, regulatory, and coding sequences that are ordinarily associated with transposition, but that occasionally confer on the organism a selective advantage and may thereby become exapted. Exapted transposable element genes (ETEs) are known to play critical roles in diverse systems, from vertebrate adaptive immunity to plant development. Yet despite their evident importance, most ETEs have been identified fortuitously and few systematic searches have been conducted, suggesting that additional ETEs may await discovery. To explore this possibility, we develop a comprehensive systematic approach to searching for ETEs. We use TE-specific conserved domains to identify with high precision genes derived from TEs and screen them for signatures of exaptation based on their similarities to reference sets of known ETEs, conventional (non-TE) genes, and TE genes across diverse genetic attributes including repetitiveness, conservation of genomic location and sequence, and levels of expression and repressive small RNAs. Applying this approach in the model plant Arabidopsis thaliana, we discover a surprisingly large number of novel high confidence ETEs. Intriguingly, unlike known plant ETEs, several of the novel ETE families form tandemly arrayed gene clusters, whereas others are relatively young. Our results not only identify novel TE-derived genes that may have practical applications but also challenge the notion that TE exaptation is merely a relic of ancient life, instead suggesting that it may continue to fundamentally drive evolution.

CGGBP1--an indispensable protein with ubiquitous cytoprotective functions

The human genome contains multiple stretches of CGG trinucleotide repeats, which act as transcription- and translation-regulatory elements but at the same time form secondary structures that impede replication and give rise to sites of chromosome fragility. Proteins binding to such DNA elements may be involved in divergent cellular processes such as transcription, DNA damage, and epigenetic state of the chromatin. We review here the work done on CGG repeats and associated proteins with special focus on a factor called CGGBP1. CGGBP1 presents with an interesting example of factors that do not have any single dedicated function, but participate indispensably in multiple processes. Both experimental results and data from cancer genome sequencing have revealed that any alteration in CGGBP1 that compromises its function is not tolerated by normal or cancer cells alike. Based upon a large amount of published data, information from databases, and unpublished results, we decipher in this review how CGGBP1 is a classic example of the 'one factor, divergent functions' paradigm of cytoprotection. By taking cues from the studies on CGGBP1, more such factors can be discovered for a better understanding of the evolution of mechanisms of cellular survival.

Maize and millet transcription factors annotated using comparative genomic and transcriptomic data

Background

Transcription factors (TFs) contain DNA-binding domains (DBDs) and regulate gene expression by binding to specific DNA sequences. In addition, there are proteins, called transcription coregulators (TCs), which lack DBDs but can alter gene expression through interaction with TFs or RNA Polymerase II. Therefore, it is interesting to identify and classify the TFs and TCs in a genome. In this study, maize (Zea mays) and foxtail millet (Setaria italica), two important species for the study of C4 photosynthesis and kranz anatomy, were selected.

Result

We conducted a comprehensive genome-wide annotation of TFs and TCs in maize B73 and in two strains of foxtail millet, Zhang gu and Yugu1, and classified them into families. To gain additional support for our predictions, we searched for their homologous genes in Arabidopsis or rice and studied their gene expression level using RNA-seq and microarray data. We identified many new TF and TC families in these two species, and described some evolutionary and functional aspects of the 9 new maize TF families. Moreover, we detected many pseudogenes and transposable elements in current databases. In addition, we examined tissue expression preferences of TF and TC families and identified tissue/condition-specific TFs and TCs in maize and millet. Finally, we identified potential C4-related TF and TC genes in maize and millet.

Conclusions

Our results significantly expand current TF and TC annotations in maize and millet. We provided supporting evidence for our annotation from genomic and gene expression data and identified TF and TC genes with tissue preference in expression. Our study may facilitate the study of regulation of gene expression, tissue morphogenesis, and C4 photosynthesis in maize and millet. The data we generated in this study are available at http://sites.google.com/site/jjlmmtf.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-818) contains supplementary material, which is available to authorized users.

Transcription factor ZBED6 mediates IGF2 gene expression by regulating promoter activity and DNA methylation in myoblasts

Zinc finger, BED-type containing 6 (ZBED6) is an important transcription factor in placental mammals, affecting development, cell proliferation and growth. In this study, we found that the expression of the ZBED6 and IGF2 were upregulated during C2C12 differentiation. The IGF2 expression levels were negatively associated with the methylation status in beef cattle (P < 0.05). A luciferase assay for the IGF2 intron 3 and P3 promoter showed that the mutant-type 439 A-SNP-pGL3 in driving reporter gene transcription is significantly higher than that of the wild-type 439 G-SNP-pGL3 construct (P < 0.05). An over-expression assay revealed that ZBED6 regulate IGF2 expression and promote myoblast differentiation. Furthermore, knockdown of ZBED6 led to IGF2 expression change in vitro. Taken together, these results suggest that ZBED6 inhibits IGF2 activity and expression via a G to A transition disrupts the interaction. Thus, we propose that ZBED6 plays a critical role in myogenic differentiation.

Visual ModuleOrganizer: a graphical interface for the detection and comparative analysis of repeat DNA modules

BACKGROUND

DNA repeats, such as transposable elements, minisatellites and palindromic sequences, are abundant in sequences and have been shown to have significant and functional roles in the evolution of the host genomes. In a previous study, we introduced the concept of a repeat DNA module, a flexible motif present in at least two occurences in the sequences. This concept was embedded into ModuleOrganizer, a tool allowing the detection of repeat modules in a set of sequences. However, its implementation remains difficult for larger sequences.

RESULTS

Here we present Visual ModuleOrganizer, a Java graphical interface that enables a new and optimized version of the ModuleOrganizer tool. To implement this version, it was recoded in C++ with compressed suffix tree data structures. This leads to less memory usage (at least 120-fold decrease in average) and decreases by at least four the computation time during the module detection process in large sequences. Visual ModuleOrganizer interface allows users to easily choose ModuleOrganizer parameters and to graphically display the results. Moreover, Visual ModuleOrganizer dynamically handles graphical results through four main parameters: gene annotations, overlapping modules with known annotations, location of the module in a minimal number of sequences, and the minimal length of the modules. As a case study, the analysis of FoldBack4 sequences clearly demonstrated that our tools can be extended to comparative and evolutionary analyses of any repeat sequence elements in a set of genomic sequences. With the increasing number of sequences available in public databases, it is now possible to perform comparative analyses of repeated DNA modules in a graphic and friendly manner within a reasonable time period.

AVAILABILITY

Visual ModuleOrganizer interface and the new version of the ModuleOrganizer tool are freely available at: http://lcb.cnrs-mrs.fr/spip.php?rubrique313.

A 5'-regulatory region and two coding region polymorphisms modulate promoter activity and gene expression of the growth suppressor gene ZBED6 in cattle

Zinc finger, BED-type containing 6 (ZBED6) is an important transcription factor in placental mammals, affecting development, cell proliferation and growth. Polymorphisms in its promoter and coding regions are likely to impact ZBED6 transcription and growth traits. In this study, rapid amplification of 5’ cDNA ends (5'-RACE) analysis revealed two transcription start sites (TSS) for the bovine ZBED6 starting within exon 1 of the ZC3H11A gene (TSS-1) and upstream of the translation start codon of the ZBED6 gene (TSS-2). There was one SNP in the promoter and two missense mutations in the coding region of the bovine ZBED6 by sequencing of the pooled DNA samples (Pool-Seq, n = 100). The promoter and coding region are the key regions for gene function; polymorphisms in these regions can alter gene expression. Quantitative real-time PCR (qPCR) analysis showed that ZBED6 has a broad tissue distribution in cattle and is highly expressed in skeletal muscle. Eleven promoter-detection vectors were constructed, which enabled the cloning of putative promoter sequences and analysis of ZBED6 transcriptional activity by luciferase reporter gene assays. The core region of the basal promoter of bovine ZBED6 is located within region -866 to -556. The activity of WT-826G-pGL3 in driving reporter gene transcription is significantly higher than that of the M-826A-pGL3 construct (P < 0.01). Analysis of gene expression patterns in homozygous full-sibling Chinese Qinchuan cattle showed that the mutant-type Hap-AGG exhibited a lower mRNA level than the wild-type Hap-GCA (P < 0.05) in longissimus dorsi muscle (LDM). Moreover, ZBED6 mRNA expression was low in C2C12 cells overexpressing the mutant-type ZBED6 (pcDNA3.1⁺-Hap-GG) (P < 0.01). Our results suggest that the polymorphisms in the promoter and coding regions may modulate the promoter activity and gene expression of bovine ZBED6 in the skeletal muscles of these cattle breeds.

Domesticated transposable element gene products in human cancer

The adaptation of transposable elements inserted within the genome to serve novel functions in a host cell, a process known as molecular domestication, is a widespread phenomenon in nature. Around fifty protein-coding genes in humans have arisen through this mechanism. Functional characterization of these domesticated genes has revealed involvement in a multitude of diverse cellular processes. Some of these functions are related to cellular activities and pathways known to be involved in cancer development. In this mini-review we discuss such roles of domesticated genes that may be aberrantly regulated in human cancer, as well as studies that have identified disrupted expression in tumors. We also describe studies that have provided definitive experimental evidence for transposable element-derived gene products in promoting tumorigenesis.