Nucleic Acid Binding Research Articles

Production of stable and pure ZC3H11A – An extensively disordered RNA binding protein

Human ZC3H11A is an RNA-binding zinc finger protein involved in mRNA export and required for the efficient growth of human nuclear replicating viruses. Its biochemical properties are largely unknown so our goal has been to produce the protein in a pure and stable form suitable for its characterization. This has been challenging since the protein is large (810 amino acids) and with only the N-terminal zinc finger domain (amino acids 1–86) being well structured, the remainder is intrinsically disordered. Our production strategies have encompassed recombinant expression of full-length, truncated and mutated ZC3H11A variants with varying purification tags and fusion proteins in several expression systems, with or without co-expression of chaperones and putative interaction partners. A range of purification schemes have been explored. Initially, only truncated ZC3H11A encompassing the zinc finger domain could successfully be produced in a stable form. It required recombinant expression in insect cells since expression in E. coli gave a protein that aggregated. To reduce problematic nucleic acid contaminations, Cys8, located in one of the zinc fingers, was substituted by Ala and Ser. Interestingly, this did not affect nucleic acid binding, but the full-length protein was stabilised while the truncated version was insoluble. Ultimately, we discovered that when using alkaline buffers (pH 9) for purification, full-length ZC3H11A expressed in Sf9 insect cells was obtained in a stable and >90 % pure form, and as a mixture of monomers, dimers, tetramers and hexamers. Many of the challenges experienced are consistent with its predicted structure and unusual charge distribution.

Role of cold shock proteins B and D in Aeromonas salmonicida subsp. salmonicida physiology and virulence in lumpfish (Cyclopterus lumpus).

Cold shock proteins (Csp) are pivotal nucleic acid binding proteins known for their crucial roles in the physiology and virulence of various bacterial pathogens affecting plant, insect, and mammalian hosts. However, their significance in bacterial pathogens of teleost fish remains unexplored. Aeromonas salmonicida subsp. salmonicida (hereafter A. salmonicida) is a psychrotrophic pathogen and the causative agent of furunculosis in marine and freshwater fish. Four csp genes (cspB, cspD, cspA, and cspC) have been identified in the genome of A. salmonicida J223 (wild type). Here, we evaluated the role of DNA binding proteins, CspB and CspD, in A. salmonicida physiology and virulence in lumpfish (Cyclopterus lumpus). A. salmonicida ΔcspB, ΔcspD, and the double ΔcspBΔcspD mutants were constructed and characterized. A. salmonicida ΔcspB and ΔcspBΔcspD mutants showed a faster growth at 28°C, and reduced virulence in lumpfish. A. salmonicida ΔcspD showed a slower growth at 28°C, biofilm formation, lower survival in low temperatures and freezing conditions (-20°C, 0°C, and 4°C), deficient in lipopolysaccharide synthesis, and low virulence in lumpfish. Additionally, ΔcspBΔcspD mutants showed less survival in the presence of bile compared to the wild type. Transcriptome analysis revealed that 200, 37, and 921 genes were differentially expressed in ΔcspB, ΔcspD, and ΔcspBΔcspD, respectively. In ΔcspB and ΔcspBΔcspD virulence genes in the chromosome and virulence plasmid were downregulated. Our analysis indicates that CspB and CspD mostly act as a transcriptional activator, influencing cell division (e.g., treB), virulence factors (e.g., aexT), and ultimately virulence.

Nucleic acid-induced NADase activation of a short Sir2-associated prokaryotic Argonaute system

Inhibition of nucleic acid targets is mediated by Argonaute (Ago) proteins guided by RNA or DNA. Although the mechanisms underpinning the functions of eukaryotic and "long" prokaryotic Ago proteins (pAgos) are well understood, those for short pAgos remain enigmatic. Here, we determine two cryoelectron microscopy structures of short pAgos in association with the NADase-domain-containing protein Sir2-APAZ from Geobacter sulfurreducens (GsSir2/Ago): the guide RNA-target DNA-loaded GsSir2/Ago quaternary complex (2.58Å) and the dimer of the quaternary complex (2.93Å). These structures show that the nucleic acid binding causes profound conformational changes that result in disorder or partial dissociation of the Sir2 domain, suggesting that it adopts a NADase-active conformation. Subsequently, two RNA-/DNA-loaded GsSir2/Ago complexes form a dimer through their MID domains, further enhancing NADase activity through synergistic effects. The findings provide a structural basis for short-pAgo-mediated defense against invading nucleic acids.

Structure of the Hir histone chaperone complex

The evolutionarily conserved HIRA/Hir histone chaperone complex and ASF1a/Asf1 co-chaperone cooperate to deposit histone (H3/H4)2 tetramers on DNA for replication-independent chromatin assembly. The molecular architecture of the HIRA/Hir complex and its mode of histone deposition have remained unknown. Here, we report the cryo-EM structure of the S.cerevisiae Hir complex with Asf1/H3/H4 at 2.9-6.8Å resolution. We find that the Hir complex forms an arc-shaped dimer with a Hir1/Hir2/Hir3/Hpc2 stoichiometry of 2/4/2/4. The core of the complex containing two Hir1/Hir2/Hir2 trimers and N-terminal segments of Hir3 forms a central cavity containing two copies of Hpc2, with one engaged by Asf1/H3/H4, in a suitable position to accommodate a histone (H3/H4)2 tetramer, while the C-terminal segments of Hir3 harbor nucleic acid binding activity to wrap DNA around the Hpc2-assisted histone tetramer. The structure suggests a model for how the Hir/Asf1 complex promotes the formation of histone tetramers for their subsequent deposition onto DNA.

Unveiling potential inhibitors targeting the nucleocapsid protein of SARS-CoV-2: Structural insights into their binding sites

The Nucleocapsid (N) protein of SARS-CoV-2 plays a crucial role in viral replication and pathogenesis, making it an attractive target for developing antiviral therapeutics. In this study, we used differential scanning fluorimetry to establish a high-throughput screening method for identifying high-affinity ligands of N-terminal domain of the N protein (N-NTD). We screened an FDA-approved drug library of 1813 compounds and identified 102 compounds interacting with N-NTD. The screened compounds were further investigated for their ability to inhibit the nucleic-acid binding activity of the N protein using electrophoretic mobility-shift assays. We have identified three inhibitors, Ceftazidime, Sennoside A, and Tannic acid, that disrupt the N protein's interaction with RNA probe. Ceftazidime and Sennoside A exhibited nano-molar range binding affinities with N protein, determined through surface plasmon resonance. The binding sites of Ceftazidime and Sennoside A were investigated using [1H, 15N]-heteronuclear single quantum coherence (HSQC) NMR spectroscopy. Ceftazidime and Sennoside A bind to the putative RNA binding site of the N protein, thus providing insights into the inhibitory mechanism of these compounds. These findings will contribute to the development of novel antiviral agents targeting the N protein of SARS-CoV-2.

An α-ketoglutarate conformational switch controls iron accessibility, activation, and substrate selection of the human FTO protein

The fat mass and obesity-associated fatso (FTO) protein is a member of the Alkb family of dioxygenases and catalyzes oxidative demethylation of N6-methyladenosine (m6A), N1-methyladenosine (m1A), 3-methylthymine (m3T), and 3-methyluracil (m3U) in single-stranded nucleic acids. It is well established that the catalytic activity of FTO proceeds via two coupled reactions. The first reaction involves decarboxylation of alpha-ketoglutarate (αKG) and formation of an oxyferryl species. In the second reaction, the oxyferryl intermediate oxidizes the methylated nucleic acid to reestablish Fe(II) and the canonical base. However, it remains unclear how binding of the nucleic acid activates the αKG decarboxylation reaction and why FTO demethylates different methyl modifications at different rates. Here, we investigate the interaction of FTO with 5-mer DNA oligos incorporating the m6A, m1A, or m3T modifications using solution NMR, molecular dynamics (MD) simulations, and enzymatic assays. We show that binding of the nucleic acid to FTO activates a two-state conformational equilibrium in the αKG cosubstrate that modulates the O2 accessibility of the Fe(II) catalyst. Notably, the substrates that provide better stabilization to the αKG conformation in which Fe(II) is exposed to O2 are demethylated more efficiently by FTO. These results indicate that i) binding of the methylated nucleic acid is required to expose the catalytic metal to O2 and activate the αKG decarboxylation reaction, and ii) the measured turnover of the demethylation reaction (which is an ensemble average over the entire sample) depends on the ability of the methylated base to favor the Fe(II) state accessible to O2.

Weighted single-step genome-wide association study for direct and maternal genetic effects associated with birth and weaning weights in sheep

Body weight is an important economic trait for sheep meat production, and its genetic improvement is considered one of the main goals in the sheep breeding program. Identifying genomic regions that are associated with growth-related traits accelerates the process of animal breeding through marker-assisted selection, which leads to increased response to selection. In this study, we conducted a weighted single-step genome-wide association study (WssGWAS) to identify potential candidate genes for direct and maternal genetic effects associated with birth weight (BW) and weaning weight (WW) in Baluchi sheep. The data used in this research included 13,408 birth and 13,170 weaning records collected at Abbas-Abad Baluchi Sheep Breeding Station, Mashhad-Iran. Genotypic data of 94 lambs genotyped by Illumina 50K SNP BeadChip for 54,241 markers were used. The proportion of variance explained by genomic windows was calculated by summing the variance of SNPs within 1 megabase (Mb). The top 10 window genomic regions explaining the highest percentages of additive and maternal genetic variances were selected as candidate window genomic regions associated with body weights. Our findings showed that for BW, the top-ranked genomic regions (1 Mb windows) explained 4.30 and 4.92% of the direct additive and maternal genetic variances, respectively. The direct additive genetic variance explained by the genomic window regions varied from 0.31 on chromosome 1 to 0.59 on chromosome 8. The highest (0.84%) and lowest (0.32%) maternal genetic variances were explained by genomic windows on chromosome 10 and 17, respectively. For WW, the top 10 genomic regions explained 6.38 and 5.76% of the direct additive and maternal genetic variances, respectively. The highest and lowest contribution of direct additive genetic variances were 1.37% and 0.42%, respectively, both explained by genomic regions on chromosome 2. For maternal effects on WW, the highest (1.38%) and lowest (0.41%) genetic variances were explained by genomic windows on chromosome 2. Further investigation of these regions identified several possible candidate genes associated with body weight. Gene ontology analysis using the DAVID database identified several functional terms, such as translation repressor activity, nucleic acid binding, dehydroascorbic acid transporter activity, growth factor activity and SH2 domain binding.

Nucleic acid mediated activation of a short prokaryotic Argonaute immune system

A short prokaryotic Argonaute (pAgo) TIR-APAZ (SPARTA) defense system, activated by invading DNA to unleash its TIR domain for NAD(P)+ hydrolysis, was recently identified in bacteria. We report the crystal structure of SPARTA heterodimer in the absence of guide-RNA/target-ssDNA (2.66 Å) and a cryo-EM structure of the SPARTA oligomer (tetramer of heterodimers) bound to guide-RNA/target-ssDNA at nominal 3.15–3.35 Å resolution. The crystal structure provides a high-resolution view of SPARTA, revealing the APAZ domain as equivalent to the N, L1, and L2 regions of long pAgos and the MID domain containing a unique insertion (insert57). Cryo-EM structure reveals regions of the PIWI (loop10-9) and APAZ (helix αN) domains that reconfigure for nucleic-acid binding and decrypts regions/residues that reorganize to expose a positively charged pocket for higher-order assembly. The TIR domains amass in a parallel-strands arrangement for catalysis. We visualize SPARTA before and after RNA/ssDNA binding and uncover the basis of its active assembly leading to abortive infection.

SOFB is a comprehensive ensemble deep learning approach for elucidating and characterizing protein-nucleic-acid-binding residues

Proteins and nucleic-acids are essential components of living organisms that interact in critical cellular processes. Accurate prediction of nucleic acid-binding residues in proteins can contribute to a better understanding of protein function. However, the discrepancy between protein sequence information and obtained structural and functional data renders most current computational models ineffective. Therefore, it is vital to design computational models based on protein sequence information to identify nucleic acid binding sites in proteins. Here, we implement an ensemble deep learning model-based nucleic-acid-binding residues on proteins identification method, called SOFB, which characterizes protein sequences by learning the semantics of biological dynamics contexts, and then develop an ensemble deep learning-based sequence network to learn feature representation and classification by explicitly modeling dynamic semantic information. Among them, the language learning model, which is constructed from natural language to biological language, captures the underlying relationships of protein sequences, and the ensemble deep learning-based sequence network consisting of different convolutional layers together with Bi-LSTM refines various features for optimal performance. Meanwhile, to address the imbalanced issue, we adopt ensemble learning to train multiple models and then incorporate them. Our experimental results on several DNA/RNA nucleic-acid-binding residue datasets demonstrate that our proposed model outperforms other state-of-the-art methods. In addition, we conduct an interpretability analysis of the identified nucleic acid binding residue sequences based on the attention weights of the language learning model, revealing novel insights into the dynamic semantic information that supports the identified nucleic acid binding residues. SOFB is available at https://github.com/Encryptional/SOFB and https://figshare.com/articles/online_resource/SOFB_figshare_rar/25499452.

Transcriptional Response of Wolfberry to Infestation with the Endophytic Fusarium nematophilum Strain NQ8GII4.

Fusarium nematophilum NQ8GII4 is an endophytic fungus isolated from the root of healthy wolfberry (Lycium barbarum). Previous studies have reported that NQ8GII4 could dwell in wolfberry roots and enhance the defense responses in wolfberry against root rot, which is caused by F. oxysporum. To further elucidate the molecular mechanism of wolfberry disease resistance induced by NQ8GII4, in the present study, we adopted RNA sequencing analysis to profile the transcriptome of wolfberry response to NQ8GII4 infestation over a time course of 3 and 7 days postinoculation. Gene ontology enrichment analysis revealed that differentially expressed genes (DEGs) were enriched in biological regulation, response to stimulus, signaling, detoxification, immune system process, transporter activity, electron carrier activity, transcription factor activity, nucleic acid binding transcription factor, and antioxidant activity. Through Kyoto Encyclopedia of Genes and Genomes analysis, it was found that many of these DEGs were enriched in pathways related to plant-pathogen interactions, hormone signal transduction, and the phenylpropanoid biosynthesis pathway in wolfberry. This result suggested that innate immunity, phytohormone signaling, and numerous phenylpropanoid compounds comprise a complex defense network in wolfberry. Chloroplast 50S ribosomal proteins were consistently located at the core position of the response in wolfberry following infestation with NQ8GII4 analyzed by the protein-protein interaction network. This study elucidated the molecular mechanism underlying the interaction between NQ8GII4 and wolfberry, clarified the wolfberry immune response network to endophytic fungi infestation, identified candidate resistance genes in wolfberry, and provided a fundamental date for subsequent work.

Transcription factors ERα and Sox2 have differing multiphasic DNA- and RNA-binding mechanisms.

Many transcription factors (TFs) have been shown to bind RNA, leading to open questions regarding the mechanism(s) of this RNA binding and its role in regulating TF activities. Here, we use biophysical assays to interrogate the k on, k off, and K d for DNA and RNA binding of two model human TFs, ERα and Sox2. Unexpectedly, we found that both proteins exhibit multiphasic nucleic acid-binding kinetics. We propose that Sox2 RNA and DNA multiphasic binding kinetics can be explained by a conventional model for sequential Sox2 monomer association and dissociation. In contrast, ERα nucleic acid binding exhibited biphasic dissociation paired with novel triphasic association behavior, in which two apparent binding transitions are separated by a 10-20 min "lag" phase depending on protein concentration. We considered several conventional models for the observed kinetic behavior, none of which adequately explained all the ERα nucleic acid-binding data. Instead, simulations with a model incorporating sequential ERα monomer association, ERα nucleic acid complex isomerization, and product "feedback" on isomerization rate recapitulated the general kinetic trends for both ERα DNA and RNA binding. Collectively, our findings reveal that Sox2 and ERα bind RNA and DNA with previously unappreciated multiphasic binding kinetics, and that their reaction mechanisms differ with ERα binding nucleic acids via a novel reaction mechanism.

The nucleic acid binding protein SFPQ represses EBV lytic reactivation by promoting histone H1 expression.

Epstein-Barr virus (EBV) uses a biphasic lifecycle of latency and lytic reactivation to infect >95% of adults worldwide. Despite its central role in EBV persistence and oncogenesis, much remains unknown about how EBV latency is maintained. We used a human genome-wide CRISPR/Cas9 screen to identify that the nuclear protein SFPQ was critical for latency. SFPQ supported expression of linker histone H1, which stabilizes nucleosomes and regulates nuclear architecture, but has not been previously implicated in EBV gene regulation. H1 occupied latent EBV genomes, including the immediate early gene BZLF1 promoter. Upon reactivation, SFPQ was sequestered into sub-nuclear puncta, and EBV genomic H1 occupancy diminished. Enforced H1 expression blocked EBV reactivation upon SFPQ knockout, confirming it as necessary downstream of SFPQ. SFPQ knockout triggered reactivation of EBV in B and epithelial cells, as well as of Kaposi's sarcoma-associated herpesvirusin B cells, suggesting a conserved gamma-herpesvirus role. These findings highlight SFPQ as a major regulator of H1 expression and EBV latency.

Do Nucleic Acid–Binding Polymers (Nabps) May Be Reduce Proinflammatory Nucleic Acids Among Trauma-Associated Hemorrhagic Shock?

Do Nucleic Acid–Binding Polymers (Nabps) May Be Reduce Proinflammatory Nucleic Acids Among Trauma-Associated Hemorrhagic Shock?

Single-molecule imaging reveals a direct role of CTCF's zinc fingers in SA interaction and cluster-dependent RNA recruitment.

CTCF is a zinc finger protein associated with transcription regulation that also acts as a barrier factor for topologically associated domains (TADs) generated by cohesin via loop extrusion. These processes require different properties of CTCF-DNA interaction, and it is still unclear how CTCF's structural features may modulate its diverse roles. Here, we employ single-molecule imaging to study both full-length CTCF and truncation mutants. We show that CTCF enriches at CTCF binding sites (CBSs), displaying a longer lifetime than observed previously. We demonstrate that the zinc finger domains mediate CTCF clustering and that clustering enables RNA recruitment, possibly creating a scaffold for interaction with RNA-binding proteins like cohesin's subunit SA. We further reveal a direct recruitment and an increase of SA residence time by CTCF bound at CBSs, suggesting that CTCF-SA interactions are crucial for cohesin stability on chromatin at TAD borders. Furthermore, we establish a single-molecule T7 transcription assay and show that although a transcribing polymerase can remove CTCF from CBSs, transcription is impaired. Our study shows that context-dependent nucleic acid binding determines the multifaceted CTCF roles in genome organization and transcription regulation.

Studies on interactions of ruthenium(II) polypyridyl complexes [Ru(bpy)2(btip)]2+ and [Ru(dmb)2(btip)]2+ with double-stranded RNA poly(A)·poly(U) under molecular crowding conditions

Ruthenium(II) polypyridyl complexes have emerged as excellent nucleic acid binders based on the easy-to-construct three-dimensional structures and rich photophysical and photochemical properties. However, most studies are usually performed under dilute conditions, which might be disturbed by the crowding environment of cells. Herein, two ruthenium(II) polypyridyl complexes, [Ru(bpy)2(btip)]2+ (Ru1, bpy = 2,2′-bipyridine, btip = 2-benzo[b]thien-2-yl-1H-imidazo[4,5-f]-[1,10]phenanthroline) and [Ru(dmb)2(btip)]2+ (Ru2, dmb = 4,4′-dimethyl-2,2′-bipyridine) have been synthesized and characterized. Interactions of the two complexes with double-stranded RNA poly(A)·poly(U) in dilute solution and molecular crowding environment have been studied by spectroscopic techniques and viscosity measurements, respectively. Spectral titration and viscosity experiments show that both Ru1 and Ru2 combine with poly(A)·poly(U) through an intercalative mode in both environments, and the binding affinity in crowding environment is significantly lower than that in dilute solution, but still comparatively high. Moreover, thermal denaturation experiments show that the abilities of Ru1 and Ru2 to poly(A)·poly(U) are significantly reduced under crowding condition, which are attributed to the reduced binding affinity of the complexes to poly(A)·poly(U) in crowding environment.

The association of protein-bound methionine sulfoxide with proteomic basis for aging in beech seeds

BackgroundEuropean beech (Fagus sylvatica L.) trees produce seeds irregularly; therefore, it is necessary to store beech seeds for forestation. Despite the acquisition of desiccation tolerance during development, beech seeds are classified as intermediate because they lose viability during long-term storage faster than typical orthodox seeds. In this study, beech seeds stored for short (3 years) or long (20 years) periods under optimal conditions and displaying 92 and 30% germination capacity, respectively, were compared.ResultsAged seeds displayed increased membrane damage, manifested as electrolyte leakage and lipid peroxidation levels. Analyses have been based on embryonic axes, which contained higher levels of reactive oxygen species (ROS) and higher levels of protein-bound methionine sulfoxide (MetO) in aged seeds. Using label-free quantitative proteomics, 3,949 proteins were identified, of which 2,442 were reliably quantified pointing to 24 more abundant proteins and 35 less abundant proteins in beech seeds under long-term storage conditions. Functional analyses based on gene ontology annotations revealed that nucleic acid binding activity (molecular function), ribosome organization or biogenesis and transmembrane transport (cellular processes), translational proteins (protein class) and membranous anatomical entities (cellular compartment) were affected in aged seeds. To verify whether MetO, the oxidative posttranslational modification of proteins that can be reversed via the action of methionine sulfoxide reductase (Msr) enzymes, is involved in the aging of beech seeds, we identified and quantified 226 MetO-containing proteins, among which 9 and 19 exhibited significantly up- and downregulated MetO levels, respectively, in beech seeds under long-term storage conditions. Several Msr isoforms were identified and recognized as MsrA1-like, MsrA4, MsrB5 and MsrB5-like in beech seeds. Only MsrA1-like displayed decreased abundance in aged seeds.ConclusionsWe demonstrated that the loss of membrane integrity reflected in the elevated abundance of membrane proteins had a higher impact on seed aging progress than the MetO/Msr system. Proteome analyses enabled us to propose protein Sec61 and glyceraldehyde-3-phosphate dehydrogenase as potential longevity modulators in beech seeds.

4D label-free proteomics analysis of oxygen-induced retinopathy with or without anti-VEGF treatment

Oxygen-induced retinopathy (OIR) animal model is widely used for retinopathy of prematurity (ROP) researches. The purpose of this study was to identify proteins and related pathways of OIR with or without anti-vascular endothelial growth factor (VEGF) treatment, for use as biomarkers in diagnosing and treating ROP. Nine samples were subjected to proteomic analysis. Retina specimens were collected from 3 OIR mice, 3 OIR mice with anti-VEGF treatment and 3 normal mice (control group). Liquid chromatography-tandem mass spectrometry analysis was performed using the 4D label-free technique. Statistically significant differentially expressed proteins, gene ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway representations, InterPro (IPR) and protein interactions were analyzed. In total, 4585 unique proteins were identified as differentially expressed proteins (DEPs). Enrichment analysis of the GO and KEGG indicated functional clusters related to peptide biosynthetic and metabolic process, cellular macromolecule biosynthetic process and nucleic acid binding in OIR group. For anti-VEGF treatment group, DEPs were clustered in DNA replication, PI3K/Akt signaling pathway and Jak/STAT signaling pathway. Proteomic profiling is useful for the exploration of molecular mechanisms of OIR and mechanisms of anti-VEGF treatment. These findings may be useful for identification of novel biomarkers for ROP pathogenesis and treatment.

Fine Mapping of qAL5.2 Controlling Anther Length in Oryza sativa.

Anther length is the critical floral trait determining hybrid rice seed production and is controlled by many quantitative trait loci (QTL). However, the cloning of genes specifically controlling anther size has yet to be reported. Here, we report the fine mapping of qAL5.2 for anther size using backcross inbred lines (BILs) in the genetic background of Oryza sativa indica Huazhan (HZ). Gene chip analysis on the BC4F2 and BC5F1 population identified effective loci on Chr1, Chr5, and Chr8 and two genomic regions on Chr5, named qAL5.1 and qAL5.2. qAL5.2 was identified in both populations with LOD values of 17.54 and 10.19, which explained 35.73% and 25.1% of the phenotypic variances, respectively. Ultimately qAL5.2 was localized to a 73 kb region between HK139 and HK140 on chromosome 5. And we constructed two near-isogenic lines (NILs) for RNA-seq analysis, named NIL-qAL5.2HZ and NIL-qAL5.2KLY, respectively. The result of the GO enrichment analysis revealed that differential genes were significantly enriched in the carbohydrate metabolic process, extracellular region, and nucleic acid binding transcription, and KEGG enrichment analysis revealed that alpha-linolenic acid metabolism was significantly enriched. Meanwhile, candidate genes of qAL5.2 were analyzed in RNA-seq, and it was found that ORF8 is differentially expressed between NIL-qAL5.2HZ and NIL-qAL5.2KLY. The fine mapping of qAL5.2 conferring anther length will promote the breed improvement of the restorer line and understanding of the mechanisms driving crop mating patterns.

Genome-wide view and characterization of natural antisense transcripts in Cannabis Sativa L.

Natural Antisense Transcripts (NATs) are a kind of complex regulatory RNAs that play crucial roles in gene expression and regulation. However, the NATs in Cannabis Sativa L., a widely economic and medicinal plant rich in cannabinoids remain unknown. In this study, we comprehensively predicted C. sativa NATs genome-wide using strand-specific RNA sequencing (ssRNA-Seq) data, and validated the expression profiles by strand-specific quantitative reverse transcription PCR (ssRT-qPCR). Consequently, a total of 307 NATs were predicted in C. sativa, including 104 cis- and 203 trans- NATs. Functional enrichment analysis demonstrated the potential involvement of the C. sativa NATs in DNA polymerase activity, RNA-DNA hybrid ribonuclease activity, and nucleic acid binding. Finally, 18 cis- and 376 trans- NAT-ST pairs were predicted to produce 621 cis- and 5,679 trans- small interfering RNA (nat-siRNAs), respectively. These nat-siRNAs were potentially involved in the biosynthesis of cannabinoids and cellulose. All these results will shed light on the regulation of NATs and nat-siRNAs in C. sativa.

Definition of the binding specificity of the T7 bacteriophage primase by analysis of a protein binding microarray using a thermodynamic model.

Protein binding microarrays (PBM), SELEX, RNAcompeteand chromatin-immunoprecipitation have been intensively used to determine the specificity of nucleic acid binding proteins. While the specificity of proteins with pronounced sequence specificity is straightforward, the determination of the sequence specificity of proteins of modest sequence specificity is more difficult. In this work, an explorative data analysis workflow for nucleic acid binding data was developed that can be used by scientists that want to analyse their binding data. The workflow is based on a regressor realized in scikit-learn, the major machine learning module for the scripting language Python. The regressor is built on a thermodynamic model of nucleic acid binding and describes the sequence specificity with base- and position-specific energies. The regressor was used to determine the binding specificity of the T7 primase. For this, we reanalysed the binding data of the T7 primase obtained with a custom PBM. The binding specificity of the T7 primase agrees with the priming specificity (5'-GTC) and the template (5'-GGGTC) for the preferentially synthesized tetraribonucleotide primer (5'-pppACCC) but is more relaxed. The dominant contribution of two positions in the motif can be explained by the involvement of the initiating and elongating nucleotides for template binding.