Tetranucleotide Sequence Research Articles

Cytosine‐Functionalized Macrocyclic Receptor for Nucleotides: Combining Three Recognition Motifs

AbstractNucleoside monophosphates perform various essential functions and serve as second messengers in living systems. Despite the success of nucleotide‐selective probes based on proteins, the creation of synthetic receptors for nucleoside monophosphates remains a challenging task. In this paper, we explore the design of synthetic receptors based on a combination of three recognition motifs. We show that each binding motif contributes to the formation of the host‐guest complexes. The core macrocycle provides recognition of the phosphate residue of the nucleotide, while naphthalimide and the attached cytosine residue participate in the coordination of the nucleobase via π‐π and hydrogen bonding interactions. The presence of cytosine in the receptor structure substantially increases the overall affinity of receptors for nucleotides. A considerable increase (100‐fold) in binding affinities is observed for the tetranucleotide DNA sequences 5′‐GGGG ‐3′ and 5′‐CCCC‐3′, as compared to the previously reported system without the attached cytosine.

Genetic variability ofChamaerops humilis(Arecaceae) throughout its native range highlights two species movement pathways from its area of origin

AbstractThe European fan palm (Chamaerops humilis, Arecaceae) is the only native palm in continental Europe, providing ecosystem services that are hard to obtain from other species. However, its populations are declining in some areas due to anthropogenic effects including climate change. Knowledge of genetic variability among natural populations is needed to establish conservation plans, to prevent genetic contamination of native stands by cultivated germplasm and to exploit it as an ornamental species. However, information on the genetic similarities among C. humilis populations is scarce. The aims of this work were to study genetic structure in C. humilis using a set of specifically designed genetic markers and to highlight genetic similarities and their relationships with geographical proximity. We sampled 301 specimens from 42 natural populations throughout the distribution area and analysed these with ten di-, tri- and tetra-nucleotide simple sequence repeats. Relationships between genetic similarities and geographical distances were analysed and populations grouped according to a genetic, geographical or national clustering. We found lower variability in populations from the eastern half of the distribution, and this lower variability was accompanied by a stronger relationship between genetic differences and spatial proximity. In addition, we found that C. humilis probably showed two patterns of spread and further differentiation: one from Morocco to southern continental Spain and then to Portugal and the Balearic Islands, and one from Morocco to Algeria, Tunisia, Sicily and continental Italy. Populations from Sardinia and France showed similarities to those from Spain and Tunisia, respectively, and may have arisen from multiple colonization events. Our results support the hypothesis that isolation on large islands may have increased diversification of the species even if all populations shared the same founder. These results have important implications for both the ecological management and the conservation of the species.

Complete genome sequence of a novel victorivirus infecting cicada flower (Cordyceps chanhua).

Cicada flower, Cordyceps chanhua, is a precious edible and medicinal mushroom with uses in both medicine and food in China. In this study, Cordyceps chanhua strain RCEF5995 was found to be coinfected by a previously characterized alternavirus, Cordyceps chanhua alternavirus 1 (CcAV1), and a novel victorivirus, tentatively named "Cordyceps chanhua victorivirus 1" (CcV1). Molecular characterization of CcV1 showed that its complete genome is 5,232 nucleotides long with a GC content of 57.5%. Sequence analysis indicated that CcV1 contains two overlapping open reading frames (ORFs), ORF1 and ORF2, encoding a putative coat protein (CP) of 742 amino acids (aa) and a putative RNA-dependent RNA polymerase (RdRp) of 836 aa, respectively. The termination codon of the CP ORF overlaps with the initiation codon of the RdRp ORF at the tetranucleotide sequence AUGA. Homolog searches, sequence comparisons, and phylogenetic analysis based on deduced amino acid sequences of RdRp indicated that CcV1 is a new member of the genus Victorivirus, family Totiviridae.

Characterization of a novel victorivirus from Nigrospora chinensis, a fungus isolated from tobacco.

Here, we characterized a new mycovirus from the fungus Nigrospora chinensis, which was named "Nigrospora chinensis victorivirus 1" (NcVV1). The NcVV1 genome is 5283 bp in length, containing two continuous open reading frames (ORFs), ORF1 and ORF2. ORF1 and ORF2 were predicted to encode a putative coat protein (CP) and an RNA-dependent RNA polymerase (RdRp), respectively. The stop codon of ORF1 overlaps with the start codon of ORF2 by the tetranucleotide sequence AUGA. Phylogenetic analysis based on amino acid sequences of RdRp and CP indicated that NcVV1 clustered with members of the genus Victorivirus in the family Totiviridae. To our knowledge, this was the first report of a mycovirus infecting N. chinensis.

Complete genome sequence of a novel victorivirus infecting Aspergillus niger.

Aspergillus niger is an important filamentous phytopathogenic fungus with a broad host range. A novel double-stranded (ds) RNA mycovirus, named Aspergillus niger victorivirus 1 (AnV1), isolated from A. niger strain baiyun3.23-4, was sequenced and analyzed. The AnV1 genome is 5317 nucleotides long with a GC content of 56%. AnV1 contains two open reading frames (ORF1 and 2), overlapping at a tetranucleotide sequence (AUGA). ORF1 encodes a putative capsid protein (CP) of 778 amino acids (aa), while ORF2 potentially encodes a putative RNA-dependent RNA polymerase (RdRp) of 826 aa. Phylogenetic analysis indicated that AnV1 is a new member of the genus Victorivirus in the family Totiviridae. As far as we know, this is the first report of the complete genome sequence of a victorivirus infecting A. niger.

CTAG vs. GATC: Structural Basis for Representational Differences in Reverse Palindromic DNA Tetranucleotide Sequences

Genomic compositional biases are broadly acknowledged. For example, the tetranucleotide CTAG is extraordinarily rare in almost all bacterial, archaeal, eukaryotic genomes and some phage genomes. Moreover, CTAG is remarkably underrepresented relative to its reverse palindromic sequence GATC. To determine the origin of compositional biases, we analyze structural features of palindromic DNA duplexes derived in all-atom molecular simulations. We answer this question by comparing the structural features of 12-mer DNA duplexes comprising a central tetramer sequence of CTAG, GATC, and CGCG. This study uses unbiased and biased molecular dynamics simulations to compare the torsional angles of dinucleotides in a DNA duplex. Specifically, we use all-atom molecular dynamics and well-tempered metadynamics simulations to produce trajectories of the torsion angles between dinucleotides in the sequence and calculate the free energy surfaces using torsional collective variables. The data show a correlation between the sequence and the degree of the duplex unwinding. More significantly, data show that the CTAG sequence has a more unwinding effect on the DNA duplex than the GATC sequence. This conclusion suggests that the CTAG sequence bases are more accessible from the minor groove surface than the GATC bases. Consequently, the CTAG sequence appears more exposed to chemical modifications which may explain its extreme underrepresentation in the bacterial genome.

Immunological Features with DNA Microsatellite Alterations in Patients with Colorectal Cancer.

Competent human DNA mismatch repair (MMR) corrects DNA polymerase mistakes made during cell replication to maintain complete DNA fidelity in daughter cells; faulty DNA MMR occurs in the setting of inflammation and neoplasia, creating base substitutions (e.g. point mutations) and frameshift mutations at DNA microsatellite sequences in progeny cells. Frameshift mutations at DNA microsatellite sequences are a detected biomarker termed microsatellite instability (MSI) for human disease, as this marker can prognosticate and determine therapeutic approaches for patients with cancer. There are two types of MSI: MSI-High (MSI-H), defined by frameshifts at mono- and di-nucleotide microsatellite sequences, and elevated microsatellite alterations at selected tetranucleotide repeats or EMAST, defined by frameshifts in di- and tetranucleotide microsatellite sequences but not mononucleotide sequences. Patients with colorectal cancers (CRCs) manifesting MSI-H demonstrate improved survival over patients without an MSI-H tumor, driven by the generation of immunogenic neoantigens caused by novel truncated proteins from genes whose sequences contain coding microsatellites; these patients' tumors contain hundreds of somatic mutations, and show responsiveness to treatment with immune checkpoint inhibitors. Patients with CRCs manifesting EMAST demonstrate poor survival over patients without an EMAST tumor, and may be driven by a more dominant defect in double strand break repair attributed to the MMR protein MSH3 over its frameshift correcting function; these patients' tumors often have a component of inflammation (and are also termed inflammation-associated microsatellite alterations) and show less somatic mutations and lack coding mononucleotide frameshift mutations that seem to generate the neoantigens seen in the majority of MSI-H tumors. Overall, both types of MSI are biomarkers that can prognosticate patients with CRC, can be tested for simultaneously in marker panels, and informs the approach to specific therapy including immunotherapy for their cancers.

A multi-modal coarse grained model of DNA flexibility mappable to the atomistic level.

We present a new coarse grained method for the simulation of duplex DNA. The algorithm uses a generalized multi-harmonic model that can represent any multi-normal distribution of helical parameters, thus avoiding caveats of current mesoscopic models for DNA simulation and representing a breakthrough in the field. The method has been parameterized from accurate parmbsc1 atomistic molecular dynamics simulations of all unique tetranucleotide sequences of DNA embedded in long duplexes and takes advantage of the correlation between helical states and backbone configurations to derive atomistic representations of DNA. The algorithm, which is implemented in a simple web interface and in a standalone package reproduces with high computational efficiency the structural landscape of long segments of DNA untreatable by atomistic molecular dynamics simulations.

Tetranucleotide Microsatellite Mutational Behavior Assessed in Real Time: Implications for Future Microsatellite Panels.

Background & AimsFifty percent of colorectal cancers show elevated microsatellite alterations at selected tetranucleotide repeats (EMAST) and are associated with inflammation, metastasis, and poor patient outcome. EMAST results from interleukin 6–induced nuclear-to-cytosolic displacement of the DNA mismatch repair protein Mutated S Homolog 3, allowing frameshifts of dinucleotide and tetranucleotide but not mononucleotide microsatellites. Unlike mononucleotide frameshifts that universally shorten in length, we previously observed expansion and contraction frameshifts at tetranucleotide sequences. Here, we developed cell models to assess tetranucleotide frameshifts in real time.MethodsWe constructed plasmids containing native (AAAG)18 and altered-length ([AAAG]15 and [AAAG]12) human D9S242 locus that placed enhanced green fluorescent protein +1 bp/-1 bp out-of-frame for protein translation and stably transfected into DNA mismatch repair–deficient cells for clonal selection. We used flow cytometry to detect enhanced green fluorescent protein–positive cells to measure mutational behavior.ResultsFrameshift mutation rates were 31.6 to 71.1 × 10-4 mutations/cell/generation and correlated with microsatellite length (r2 = 0.986, P = .0375). Longer repeats showed modestly higher deletion over insertion rates, with both equivalent for shorter repeats. Accumulation of more deletion frameshifts contributed to a distinct mutational bias for each length (overall: 77.8% deletions vs 22.2% insertions), likely owing to continual deletional mutation of insertions. Approximately 78.9% of observed frameshifts were 1 AAAG repeat, 16.1% were 2 repeats, and 5.1% were 3 or more repeats, consistent with a slipped strand mispairing mutation model.ConclusionsTetranucleotide frameshifts show a deletion bias and undergo more than 1 deletion event via intermediates, with insertions converted into deletions. Tetranucleotide markers added to traditional microsatellite instability panels will be able to determine both EMAST and classic microsatellite instability, but needs to be assessed by multiple markers to account for mutational behavior and intermediates.

Profile and development of microsatellite primers for Acanthogobius ommaturus based on high-throughput sequencing technology

Acanthogobius ommaturus, a fish species of the Family Gobiidae, is a marine commercial fish perched on the bottom of seawater. In this study, Illumina high-throughput sequencing technology was applied to obtain the candidate microsatellite markers of A. ommaturus. A total of 4 746 microsatellite-rich fragments were found, of which 4 542 microsatellites are with primer fragments, containing 971 dinucleotide sequences, 2 643 trinucleotide sequences, 569 tetranucleotide sequences, 406 pentanucleotide sequences, and 212 hexanucleotide sequences. Based on the results of high-throughput sequencing, a total of 141 pairs of the microsatellite primers were designed and screened. And then 24 polymorphic primers were finally obtained by polyacrylamide gel electrophoresis. In total, 271 alleles were detected in the 24 pairs of primers. The number of alleles for different primers ranged from 5 to 19. The average number of effective alleles (Na) was 11.292; the average observed heterozygosity (Ho) of the 24 pairs of primers was 0.665, the average expected heterozygosity (He) was 0.880, and the average polymorphic information content was 0.846. All sites were highly polymorphic (PIC>0.50).

Complete genome sequence of a novel victorivirus isolated from the sesame charcoal rot fungus Macrophomina phaseolina.

Macrophomina phaseolina is an important phytopathogenic fungus with a broad host range. Here, the complete genome sequence of a novel victorivirus, tentatively named Macrophomina phaseolina victorivirus 1 (MpV1), was identified from strain 2012-019 of M. phaseolina. The MpV1 genome is 5,128 nucleotides long with a predicted GC content of 62%. Sequence analysis indicated that two open reading frames (ORF 1 and 2) overlap at a tetranucleotide AUGA sequence. Proteins encoded by ORF1 and ORF2 showed significant sequence similarity to coat proteins and the RNA-dependent RNA polymerases, respectively, of members of the family Totiviridae. Analysis of the genomic structure of MpV1, homolog searches of the deduced amino acid sequences, and phylogenetic analysis indicated that MpV1 is a new member of the genus Victorivirus. As far as we know, this is the first report of the full-length nucleotide sequence of the genome of a novel victorivirus that infects M. phaseolina.

Adenovirus VA RNAI Blocks ASC Oligomerization and Inhibits NLRP3 Inflammasome Activation.

Virus infected immune cells can rapidly respond to the invader by activating the inflammasome and as a consequence release proinflammatory cytokines and eventually die by pyroptosis. In human adenovirus-5 (Ad5) infected THP-1 cells, inhibition of NLRP3 inflammasome activation was demonstrated by a decreased secretion of HMGB1 and matured forms of caspase-1and IL-1ß. An Ad5 mutant virus defective in expression of the non-coding VA RNAI failed to inhibit the NLRP3 inflammasome and in addition displayed formation of ASC specks and increased cell lysis. Importantly, in vitro synthesized VA RNAI was able to inhibit the NLRP3 inflammasome activity in THP-1 cells in the absence of an Ad5 infection, suggesting that VA RNAI binding to PKR and blocking its function is sufficient for inhibition of the NLRP3 inflammasome. Although the inhibition of NLRP3 inflammasome activation required the phylogenetically conserved base paired tetranucleotide sequence in the central stem of VA RNAI, we demonstrate that PKR binding to VA RNAI primarily protected the apical stem, but not the tetranucleotide sequence itself. VA RNAI did not influence the interaction between PKR and NLRP3. In contrast, we describe a novel interaction between PKR and ASC and further show that VA RNAI inhibited ASC phosphorylation and oligomerization. Collectively, our results indicate a novel role for Ad5 VA RNAI as an inhibitor of NLRP3 inflammasome activation by targeting the cellular pro-inflammatory protein PKR.

The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules.

We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine–Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.

Genome-wide mining of microsatellites in king cobra (Ophiophagus hannah) and cross-species development of tetranucleotide SSR markers in Chinese cobra (Naja atra).

The complete genome sequence provides the opportunity for genome-wide and coding region analysis of SSRs in the king cobra and for cross-species identification of microsatellite markers in the Chinese cobra. In the Ophiophagus hannah genome, tetranucleotide repeats (38.03%) were the most abundant category, followed by dinucleotides (23.03%), pentanucleotides (13.07%), mononucleotides (11.78%), trinucleotides (11.49%) and hexanucleotides (2.6%). Twenty predominant motifs in the O. hannah genome were (A)n (C)n, (AC)n, (AG)n, (AT)n, (AGG)n, (AAT)n, (AAG)n, (AAC)n, (ATG)n, (ATAG)n, (AAGG)n, (ATCT)n, (CCTT)n, (ATTT)n, (AAAT)n, (AATAG)n, (ATTCT)n, (ATATGT)n, (AGATAT)n. In total, 4344 SSRs were found in coding sequences (CDSs). Tetranucleotides (52.79%) were the most abundant microsatellite type in CDS, followed by trinucleotides (28.50%), dinucleotides (11.02%), pentanucleotides (4.42%), mononucleotides (1.77%), and hexanucleotides (1.50%). A total of 984 CDSs containing microsatellites were assigned 11152 Gene Ontology (GO) functional terms. Gene Ontology (GO) analysis demonstrated that cellular process, cell and binding were the most frequent GO terms in biological process, cellular component and molecular function, respectively. Thirty-two novel highly polymorphic (PIC > 0.5) SSR markers for Naja atra were developed from cross-species amplification based on the tetranucleotide microsatellite sequences in the king cobra genome. The number of alleles (NA) per locus had between 3 and 11 alleles with an average of 6.5, the polymorphism information content (PIC) value ranged from 0.521 to 0.858 (average = 0.707), the observed heterozygosity (Ho) of 32 microsatellite loci ranged from 0.292 to 0.875 (mean = 0.678), the expected heterozygosity (HE) ranged from 0.561 to 0.889 (average = 0.761), and 3 microsatellite loci exhibited statistically significant departure from Hardy-Weinberg equilibrium (HWE) after Bonferroni correction (p < 0.003).

A Victorivirus and Two Novel Mitoviruses Co-Infected the Plant Pathogen Nigrospora oryzae.

Three dsRNAs, in sizes of approximately 2.5–5 kbp, were detected in the plant pathogenic fungus Nigrospora oryzae strain CS-7.5-4. Genomic analysis showed that the 5.0 kb dsRNA was a victorivirus named as Nigrospora oryzae victorivirus 2 (NoRV2). The genome of NoRV2 was 5166 bp in length containing two overlapping open reading frames (ORFs), ORF1 and ORF2. ORF1 was deduced to encode a coat protein (CP) showing homology to the CPs of viruses belonging to the Totiviridae family. The stop codon of ORF1 and the start codon of ORF2 were overlapped by the tetranucleotide sequence AUGA. ORF2 was predicted to encode an RNA-dependent RNA polymerase (RdRp), which was highly similar to the RdRps of victoriviruses. Virus-like particle examination demonstrated that the genome of NoRV2 was solely encapsidated by viral particles with a diameter of approximately 35 nm. The other two dsRNAs that were less than 3.0 kb were predicted to be the genomes of two mitoviruses, named as Nigrospora oryzae mitovirus 1 (NoMV1) and Nigrospora oryzae mitovirus 2 (NoMV2). Both NoMV1 and NoMV2 were A-U rich and with lengths of 2865 and 2507 bp, respectively. Mitochondrial codon usage inferred that each of the two mitoviruses contains a major large ORF encoding a mitoviral RdRp. Horizontal transfer experiments showed that the NoMV1 and NoMV2 could be cotransmitted horizontally via hyphal contact to other virus-free N. oryzae strains and causes phenotypic change to the recipient, such as an increase in growth rate. This is the first report of mitoviruses in N. oryzae.

Abstract 3368: Insertion and deletion frameshift rates and mutational spectra of tetranucleotide microsatellites in DNA mismatch repair-deficient human cells

Abstract Background. Oxidative DNA damage and inflammatory cytokines cause loss of function of the MMR protein MSH3. MSH3-MSH2 (MutSβ) function is required for correction of postreplicative polymerase errors in the form of frameshifts at tetranucleotide microsatellite sequences to prevent elevated microsatellite alterations at selected tetranucleotides (EMAST). EMAST is increasingly being used as a biomarker for MSH3 deficiency in subtyping colorectal cancers due to its association with metastasis and poor patient outcome. Here, we sought to understand frameshift mutational behavior of tetranucleotide sequences in the absence of MMR. Methods. We generated EGFP-based mutation reporter model systems containing the D9S242 tetranucleotide microsatellite at its native length of (AAAG)18 and at modified lengths of (AAAG)15 and (AAAG)12 to represent other endogenous microsatellite lengths. Plasmid vectors were constructed for each length that placed EGFP +1 bp or -1 bp out-of-frame for protein translation, enabling us to measure deletion and insertion frameshifts. These constructs were stably-integrated into MMR-deficient HCT116 cells, and monoclonal cell lines containing the constructs were isolated to measure mutation rates and spectra. Deletion and insertion frameshift mutations restored the EGFP reading frame to become detectable by flow cytometry. Genomic DNA from mutated cells were isolated and sequenced for insertion/deletion frameshift mutational biases. Results. Combined insertion/deletion frameshift mutation rates at D9S242 locus ranged from 31.6 X 10-4 to 71.1 X 10-4 mutations/cell/generation and was strongly correlated with longer length of tetranucleotide microsatellites (r2 = 0.986, P = 0.0375). We observed slightly higher deletion mutation rates over insertion rates at longer microsatellites, but equivalent insertion and deletion rates at shorter microsatellites. However, accumulation of more deletion frameshifts over time contributed to a distinct mutational bias for deletions for each length, likely due to continual frameshift mutation at insertions. About 79% of all observed frameshifts (insertion or deletion) were one- repeat (e.g., four nucleotides), 16% two-repeat, and 5% three or more repeat mutations that were consistent with a slipped strand mispairing mutation model. Conclusions. Our results suggest that intact MMR recognizes and repairs shorter length tetranucleotide frameshifts more efficiently than larger ones, and preferentially corrects loops on the template strand (preventing deletions) over loops on the newly-synthesized strand (preventing insertions) in human cells. MMR deficiency thus manifests in human tissue with a deletion mutation bias at tetranucleotide sequences (EMAST). Longer tetranucleotide microsatellites have increased propensities for larger DNA slippage errors. Citation Format: Maide O. Raeker, Jovan Pierre-Charles, John M. Carethers. Insertion and deletion frameshift rates and mutational spectra of tetranucleotide microsatellites in DNA mismatch repair-deficient human cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3368.

The sequence specificity of UV-induced DNA damage in a systematically altered DNA sequence

The sequence specificity of UV-induced DNA damage was investigated in a specifically designed DNA plasmid using two procedures: end-labelling and linear amplification. Absorption of UV photons by DNA leads to dimerisation of pyrimidine bases and produces two major photoproducts, cyclobutane pyrimidine dimers (CPDs) and pyrimidine(6-4)pyrimidone photoproducts (6-4PPs). A previous study had determined that two hexanucleotide sequences, 5′-GCTC*AC and 5′-TATT*AA, were high intensity UV-induced DNA damage sites. The UV clone plasmid was constructed by systematically altering each nucleotide of these two hexanucleotide sequences. One of the main goals of this study was to determine the influence of single nucleotide alterations on the intensity of UV-induced DNA damage. The sequence 5′-GCTC*AC was designed to examine the sequence specificity of 6-4PPs and the highest intensity 6-4PP damage sites were found at 5′-GTTC*CC nucleotides. The sequence 5′-TATT*AA was devised to investigate the sequence specificity of CPDs and the highest intensity CPD damage sites were found at 5′-TTTT*CG nucleotides. It was proposed that the tetranucleotide DNA sequence, 5’-YTC*Y (where Y is T or C), was the consensus sequence for the highest intensity UV-induced 6-4PP adduct sites; while it was 5′-YTT*C for the highest intensity UV-induced CPD damage sites. These consensus tetranucleotides are composed entirely of consecutive pyrimidines and must have a DNA conformation that is highly productive for the absorption of UV photons.

Targeting IS608 transposon integration to highly specific sequences by structure-based transposon engineering.

Transposable elements are efficient DNA carriers and thus important tools for transgenesis and insertional mutagenesis. However, their poor target sequence specificity constitutes an important limitation for site-directed applications. The insertion sequence IS608 from Helicobacter pylori recognizes a specific tetranucleotide sequence by base pairing, and its target choice can be re-programmed by changes in the transposon DNA. Here, we present the crystal structure of the IS608 target capture complex in an active conformation, providing a complete picture of the molecular interactions between transposon and target DNA prior to integration. Based on this, we engineered IS608 variants to direct their integration specifically to various 12/17-nt long target sites by extending the base pair interaction network between the transposon and the target DNA. We demonstrate in vitro that the engineered transposons efficiently select their intended target sites. Our data further elucidate how the distinct secondary structure of the single-stranded transposon intermediate prevents extended target specificity in the wild-type transposon, allowing it to move between diverse genomic sites. Our strategy enables efficient targeting of unique DNA sequences with high specificity in an easily programmable manner, opening possibilities for the use of the IS608 system for site-specific gene insertions.

Genome-wide mining of perfect microsatellites and tetranucleotide orthologous microsatellites estimates in six primate species

Advancement in genome sequencing and in silico mining tools have provided new opportunities for comparative primate genomics of microsatellites. The SSRs (simple sequence repeats) numbers were not correlated with the genome size (Pearson, r=0.310, p=0.550), and were positively correlated with the total length of SSRs (Pearson, r=0.992, p=0.00). A total of 224,289 tetranucleotide orthologous microsatellites families and 367 single-copy orthologous SSRs loci were found in six primate species by homologous alignment. The inner mutation types of single-copy orthologous SSRs loci included the copy number variance, point mutation, and chromosomal translocation. The accumulated repeat times and average length of tetranucleotide orthologous microsatellites in Rhinopithecus roxellana, Papio anubis and Macaca mulatta were longer than Homo sapiens and Pan troglodytes, which showed the tetranucleotide orthologous SSRs loci had more repeat times and longer average length on the branches with earlier divergence time, one exception may be Microcebus murinus as a primitive monkey with a smallest morphology in Malagasy. Our conclusion indicated that single-copy tetranucleotide orthologous SSRs sequences accumulated individual mutation more slowly through time in H. sapiens and P. troglodytes than in R. roxellanae, P. anubis and M. mulatta. However, such divergence wouldn't arise uniformly in all branches of the primate tree. A comparison of genomic sequence assemblages would offer remarkable insights about comparisons and contrasts, and the evolutionary processes of the microsatellites involved in human and nonhuman primate species.

Conserved intergenic sequences revealed by CTAG-profiling in Salmonella: thermodynamic modeling for function prediction

Highly conserved short sequences help identify functional genomic regions and facilitate genomic annotation. We used Salmonella as the model to search the genome for evolutionarily conserved regions and focused on the tetranucleotide sequence CTAG for its potentially important functions. In Salmonella, CTAG is highly conserved across the lineages and large numbers of CTAG-containing short sequences fall in intergenic regions, strongly indicating their biological importance. Computer modeling demonstrated stable stem-loop structures in some of the CTAG-containing intergenic regions, and substitution of a nucleotide of the CTAG sequence would radically rearrange the free energy and disrupt the structure. The postulated degeneration of CTAG takes distinct patterns among Salmonella lineages and provides novel information about genomic divergence and evolution of these bacterial pathogens. Comparison of the vertically and horizontally transmitted genomic segments showed different CTAG distribution landscapes, with the genome amelioration process to remove CTAG taking place inward from both terminals of the horizontally acquired segment.