Palindromic Structure Research Articles

Sequencing and characterization of the macronuclear rDNA minichromosome of the protozoan Tetrahymena pyriformis

Tetrahymena ribosomal DNA (rDNA) is an ideal system for studying eukaryotic DNA replication and gene transcription. In this study, we developed a new method to isolate rDNA from Tetrahymena cells and used it to sequence and annotate the complete 19,670 bp macronuclear rDNA minichromosome of Tetrahymena pyriformis, a species that lacks the germ-line micronucleus and is unable to undergo sexual reproduction. The key features of T. pyriformis and Tetrahymena thermophila rDNA sequences were then compared. Our results showed (i) the short inverted repeats (M repeats) essential for formation of rDNA minichromosome palindromic structure during sexual reproduction in Tetrahymena are highly conserved in T. pyriformis; (ii) in contrast to T. thermophila, which has two tandem domains that coordinately regulate rDNA replication, T. pyriformis has only a single domain; (iii) the 35S pre-rRNA precursor has 80.25% similarity between the two species; and (iv) the G + C content is higher in the transcribed region than the non-transcribed region in both species, but the GC-skew is more stable in T. pyriformis. The new isolation method and annotated information for the T. pyriformis rDNA minichromosome will provide a useful resource for studying DNA replication and chromosome copy number control in Tetrahymena.

Sequencing artifacts derived from a library preparation method using enzymatic fragmentation.

DNA fragmentation is a fundamental step during library preparation in hybridization capture-based, short-read sequencing. Ultra-sonication has been used thus far to prepare DNA of an appropriate size, but this method is associated with a considerable loss of DNA sample. More recently, studies have employed library preparation methods that rely on enzymatic fragmentation with DNA endonucleases to minimize DNA loss, particularly in nano-quantity samples. Yet, despite their wide use, the effect of enzymatic fragmentation on the resultant sequences has not been carefully assessed. Here, we used pairwise comparisons of somatic variants of the same tumor DNA samples prepared using ultrasonic and enzymatic fragmentation methods. Our analysis revealed a substantially larger number of recurrent artifactual SNVs/indels in endonuclease-treated libraries as compared with those created through ultrasonication. These artifacts were marked by palindromic structure in the genomic context, positional bias in sequenced reads, and multi-nucleotide substitutions. Taking advantage of these distinctive features, we developed a filtering algorithm to distinguish genuine somatic mutations from artifactual noise with high specificity and sensitivity. Noise cancelling recovered the composition of the mutational signatures in the tumor samples. Thus, we provide an informatics algorithm as a solution to the sequencing errors produced as a consequence of endonuclease-mediated fragmentation, highlighted for the first time in this study.

Birth, expansion, and death of VCY-containing palindromes on the human Y chromosome

BackgroundLarge palindromes (inverted repeats) make up substantial proportions of mammalian sex chromosomes, often contain genes, and have high rates of structural variation arising via ectopic recombination. As a result, they underlie many genomic disorders. Maintenance of the palindromic structure by gene conversion between the arms has been documented, but over longer time periods, palindromes are remarkably labile. Mechanisms of origin and loss of palindromes have, however, received little attention.ResultsHere, we use fiber-FISH, 10x Genomics Linked-Read sequencing, and breakpoint PCR sequencing to characterize the structural variation of the P8 palindrome on the human Y chromosome, which contains two copies of the VCY (Variable Charge Y) gene. We find a deletion of almost an entire arm of the palindrome, leading to death of the palindrome, a size increase by recruitment of adjacent sequence, and other complex changes including the formation of an entire new palindrome nearby. Together, these changes are found in ~ 1% of men, and we can assign likely molecular mechanisms to these mutational events. As a result, healthy men can have 1–4 copies of VCY.ConclusionsGross changes, especially duplications, in palindrome structure can be relatively frequent and facilitate the evolution of sex chromosomes in humans, and potentially also in other mammalian species.

On Some Special Property of The Farey Sequence

In this paper, some special property of the Farey sequence is discussed. We prove in each term of the Farey sequence, the sum of elements in the denominator is two times of the sum of elements in the numerator. We also prove that the Farey sequence contains a palindrome structure.

Male mice with large inversions or deletions of X-chromosome palindrome arms are fertile and express their associated genes during post-meiosis

Large (>10 kb) palindromic sequences are enriched on mammalian sex chromosomes. In mice, these palindromes harbor gene families (≥2 gene copies) expressed exclusively in post-meiotic testicular germ cells, a time when most single-copy sex-linked genes are transcriptionally repressed. This observation led to the hypothesis that palindromic structures or having ≥2 gene copies enable post-meiotic gene expression. We tested these hypotheses by using CRISPR to precisely engineer large (10’s of kb) inversions and deletions of X-chromosome palindrome arms for two regions that carry the mouse 4930567H17Rik and Mageb5 palindrome gene families. We found that 4930567H17Rik and Mageb5 gene expression is unaffected in mice carrying palindrome arm inversions and halved in mice carrying palindrome arm deletions. We assessed whether palindrome-associated genes were sensitive to reduced expression in mice carrying palindrome arm deletions. Male mice carrying palindrome arm deletions are fertile and show no defects in post-meiotic spermatogenesis. Together, these findings suggest palindromic structures on the sex chromosomes are not necessary for their associated genes to evade post-meiotic transcriptional repression and that these genes are not sensitive to reduced expression levels. Large sex chromosome palindromes may be important for other reasons, such as promoting gene conversion between palindrome arms.

Variation in human chromosome 21 ribosomal RNA genes characterized by TAR cloning and long-read sequencing.

Despite the key role of the human ribosome in protein biosynthesis, little is known about the extent of sequence variation in ribosomal DNA (rDNA) or its pre-rRNA and rRNA products. We recovered ribosomal DNA segments from a single human chromosome 21 using transformation-associated recombination (TAR) cloning in yeast. Accurate long-read sequencing of 13 isolates covering ∼0.82 Mb of the chromosome 21 rDNA complement revealed substantial variation among tandem repeat rDNA copies, several palindromic structures and potential errors in the previous reference sequence. These clones revealed 101 variant positions in the 45S transcription unit and 235 in the intergenic spacer sequence. Approximately 60% of the 45S variants were confirmed in independent whole-genome or RNA-seq data, with 47 of these further observed in mature 18S/28S rRNA sequences. TAR cloning and long-read sequencing enabled the accurate reconstruction of multiple rDNA units and a new, high-quality 44 838 bp rDNA reference sequence, which we have annotated with variants detected from chromosome 21 of a single individual. The large number of variants observed reveal heterogeneity in human rDNA, opening up the possibility of corresponding variations in ribosome dynamics.

Intracellular generation of single-strand template increases the knock-in efficiency by combining CRISPR/Cas9 with AAV.

Targeted integration of transgenes facilitates functional genomic research and holds prospect for gene therapy. The established microhomology-mediated end-joining (MMEJ)-based strategy leads to the precise gene knock-in with easily constructed donor, yet the limited efficiency remains to be further improved. Here, we show that single-strand DNA (ssDNA) donor contributes to efficient increase of knock-in efficiency and establishes a method to achieve the intracellular linearization of long ssDNA donor. We identified that the CRISPR/Cas9 system is responsible for breaking double-strand DNA (dsDNA) of palindromic structure in inverted terminal repeats (ITRs) region of recombinant adeno-associated virus (AAV), leading to the inhibition of viral second-strand DNA synthesis. Combing Cas9 plasmids targeting genome and ITR with AAV donor delivery, the precise knock-in of gene cassette was achieved, with 13-14% of the donor insertion events being mediated by MMEJ in HEK 293T cells. This study describes a novel method to integrate large single-strand transgene cassettes into the genomes, increasing knock-in efficiency by 13.6-19.5-fold relative to conventional AAV-mediated method. It also provides a comprehensive solution to the challenges of complicated production and difficult delivery with large exogenous fragments.

Combating Drug-Resistant Fungi with Novel Imperfectly Amphipathic Palindromic Peptides.

Antimicrobial peptides are an important weapon against invading pathogens and are potential candidates as novel antibacterial agents, but their antifungal activities are not fully developed. In this study, a set of imperfectly amphipathic peptides was developed based on the imperfectly amphipathic palindromic structure R n(XRXXXRX)R n ( n = 1, 2; X represents L, I, F, or W), and the engineered peptides exhibited high antimicrobial activities against all fungi and bacteria tested (including fluconazole-resistant Candida albicans), with geometric mean (GM) MICs ranging from 2.2 to 6.62 μM. Of such peptides, 13 (I6) (RRIRIIIRIRR-NH2) that was Ile rich in its hydrophobic face had the highest antifungal activity (GMfungi = 1.64 μM) while showing low toxicity and high salt and serum tolerance. It also had dramatic LPS-neutralizing propensity and a potent membrane-disruptive mechanism against microbial cells. In summary, these findings were useful for short AMPs design to combat the growing threat of drug-resistant fungal and bacterial infections.

Defining the binding determinants of Shewanella oneidensis OxyR: Implications for the link between the contracted OxyR regulon and adaptation

It is well-established that OxyR functions as a transcriptional activator of the peroxide stress response in bacteria, primarily based on studies on Escherichia coli Recent investigations have revealed that OxyRs of some other bacteria can regulate gene expression through both repression and activation or repression only; however, the underlying mechanisms remain largely unknown. Here, we demonstrated in γ-proteobacteriumShewanella oneidensis regulation of OxyR on expression of major catalase gene katB in a dual-control manner through interaction with a single site in the promoter region. Under non-stress conditions, katB expression was repressed by reduced OxyR (OxyRred), whereas when oxidized, OxyR (OxyRoxi) outcompeted OxyRred for the site because of substantially enhanced affinity, resulting in a graded response to oxidative stress, from repression to derepression to activation. The OxyR-binding motif is characterized as a combination of the E. coli motif (tetranucleotides spaced by heptanucleotide) and palindromic structure. We provided evidence to suggest that the S. oneidensis OxyR regulon is significantly contracted compared with those reported, probably containing only five members that are exclusively involved in oxygen reactive species scavenging and iron sequestering. These characteristics probably reflect the adapting strategy of the bacteria that S. oneidensis represents to thrive in redox-stratified microaerobic and anaerobic environments.

A Universal Classification System of Skin Substitutes Inspired by Factorial Design

The complexity of the dermal layer of skin means that damage to this section can result in permanent impairment of function. Partial or total dermal loss is a feature of deep burns and chronic wounds such as pressure sores or diabetic ulcers. The issues posed by traditional skin grafts have led to substantial research being carried out in the fields of tissue engineering and biomaterials science to develop a vast array of alternative skin substitutes. Given the large number of different materials, manufacturing methods, and techniques for implementation described for artificial skin substitutes, many classification systems have been created to simplify their categorization. Some of these systems are oriented toward clinicians while others toward researchers. However, none address the needs of both groups and none are intuitive. The creation of an effective classification system would be particularly helpful in the regulation, distribution, organization, and selection of skin substitutes. The aim of this review is to examine existing methods of classification of skin substitutes, and to propose a new system that uses an algorithm that is inspired by factorial design. Our system allows multiple factors to be simultaneously investigated or in this case, described, since all skin substitutes possess multiple characteristics: (1) cellularity (acellular or cellular), (2) layering (single layer or bilayer), (3) replaced region (epidermis, dermis, or both), (4) materials used (natural, synthetic, or both), and (5) permanence (temporary or permanent). The factors and levels are combined into an algorithm where all the possible combinations are shown. The multifactorial and palindromic structure of our system should enable all users to quickly understand the makeup of a selected skin substitute, or search for a skin substitute depending on their specific requirements. We feel that our proposed classification can be used by clinicians and biomedical researchers alike, which should be an advantage given the multidisciplinary nature of the tissue engineering field and the science that underpins the development of skin substitutes. We also touch upon some of the state-of-the-art skin substitutes that are commercially available or under development to demonstrate how our new method of classification might work.

SbcC-SbcD and ExoI process convergent forks to complete chromosome replication

SbcC-SbcD are the bacterial orthologs of Mre11-Rad50, a nuclease complex essential for genome stability, normal development, and viability in mammals. In vitro, these enzymes degrade long DNA palindromic structures. When inactivated along with ExoI in Escherichia coli, or Sae2 in eukaryotes, palindromic amplifications arise and propagate in cells. However, long DNA palindromes are not normally found in bacterial or human genomes, leaving the cellular substrates and function of these enzymes unknown. Here, we show that during the completion of DNA replication, convergent replication forks form a palindrome-like structural intermediate that requires nucleolytic processing by SbcC-SbcD and ExoI before chromosome replication can be completed. Inactivation of these nucleases prevents completion from occurring, and under these conditions, cells maintain viability by shunting the reaction through an aberrant recombinational pathway that leads to amplifications and instability in this region. The results identify replication completion as an event critical to maintain genome integrity and cell viability, demonstrate SbcC-SbcD-ExoI acts before RecBCD and is required to initiate the completion reaction, and reveal how defects in completion result in genomic instability.

Plasticity of the Mre11-Rad50-Xrs2-Sae2 nuclease ensemble in the processing of DNA-bound obstacles.

The budding yeast Mre11-Rad50-Xrs2 (MRX) complex and Sae2 function together in DNA end resection during homologous recombination. Here we show that the Ku complex shields DNA ends from exonucleolytic digestion but facilitates endonucleolytic scission by MRX with a dependence on ATP and Sae2. The incision site is enlarged into a DNA gap via the exonuclease activity of MRX, which is stimulated by Sae2 without ATP being present. RPA renders a partially resected or palindromic DNA structure susceptible to MRX-Sae2, and internal protein blocks also trigger DNA cleavage. We present models for how MRX-Sae2 creates entry sites for the long-range resection machinery.

Genome-wide motif predictions of BCARR-box in the amino-acid repressed genes of Lactobacillus helveticus CM4

BackgroundA BCARR (branched-chain amino acid responsive repressor) identified in proteolytic gene expressions in Lactobacillus helveticus is considered to negatively control transcriptions by binding to operator sites at the promoter regions in the presence of BCAAs. However, the distributions and regulatory potential of the BCARR in all genes repressed by BCAAs in CM4 remains unclear.ResultsA genome-wide search for the BCARR-box was conducted to clarify the contribution of BCARR in the regulation of amino acid metabolism in L. helveticus CM4. Among all 2174 genes of CM4, 390 genes repressed by amino acids were selected for the search of the BCARR-box. The annotated 33 genes among the 67 predicted BCARR-boxes were mainly linked to amino acid metabolism. The BCARR-boxes were mainly located adjacent to the −35 sequence of the promoter; however, the repressive effects in different locations were similar. Notably, the consensus BCARR-box motif, 5′-A1A2A3A4A5W6N7N8N9W10T11T12W13T14T15–3′, observed in highly repressed genes, revealed more frequent A-T base pairing and a lower free energy than that in lowly repressed genes. A MEME analysis also supported the lower frequency of T at positions 12, 14, 13 and 15 in the BCARR-box sequence of the lowly repressed gene group. These results reveal that genes with a more stable palindromic structure might be preferable targets for BCARR binding and result in higher repressions in the target gene expressions.ConclusionsOur genome-wide search revealed the involvement of the proteolytic system, transporter system and some transcriptional regulator systems in BCARR-box regulation in L. helveticus CM4.

Purcell effect in complex photonic structures based on optical Fibonacci lattices

The modification of probability of spontaneous emission for a dipole emitter that placed in Fibonacci lattice-based palindrome structures was examined. Localization properties of eigenstates, that appears in photonic band gaps of palindrome structure was investigated. It was shown, that for a dipole placed in centre of palindrome structure an increase in the spontaneous emission rate will be observed and, in addition, calculated value of Purcell factor is significantly greater than for critical localized states of usual Fibonacci lattices, and close to microcavity eigenmode case.

The complete chloroplast DNA sequence of Trapa maximowiczii Korsh. (Trapaceae), and comparative analysis with other Myrtales species

Trapaceae is an aquatic plant family widely distributed in and around Asia, Europe, and Africa. The fruits, while rich in starch and protein, also have edible and medicinal value. The first complete chloroplast genome in Trapaceae (Trapa maximowiczii Korsh.) has been determined. The genome size is 155,577bp, which includes two inverted repeats (IR) of 24,388bp, one large single copy (LSC) of 88,528bp, and one small single copy (SSC) of 18,273bp. The whole GC content is 36.4%. There are 111 genes in this genome, and 17 are duplicated in the IR. Fourteen genes have one intron, and rps12 and clpP have two introns. The SSR and repeat analysis revealed 61 SSRs, 22 repeat units, and 16 palindromic structures in this chloroplast genome. In a comparative analysis with five species, the detailed boundaries were compared between IR and LSC/SSC that are supposed to be closely related to T. maximowiczii and Lagerstroemia fauriei; gene infA has been lost only in T. maximowiczii and L. fauriei. The mVISTA analysis revealed several variable regions in the Myrtales genome including loci ycf1, ycf2, accD and rpl32-trnL, which can be useful in phylogenetic studies. The phylogenetic position of T. maximowiczii in the whole genome, LSC, SSC, and IR regions has been confirmed using 8 taxa from Myrtales, which is in agreement with other reports that place T. maximowiczii in the same clade as that of Lythraceae. These data represent valuable genetic resources for future research on Trapa population genetics.

Identification and Expression of Equine MER-Derived miRNAs.

MicroRNAs (miRNAs) are single-stranded, small RNAs (21–23 nucleotides) that function in gene silencing and translational inhibition via the RNA interference mechanism. Most miRNAs originate from host genomic regions, such as intergenic regions, introns, exons, and transposable elements (TEs). Here, we focused on the palindromic structure of medium reiteration frequencies (MERs), which are similar to precursor miRNAs. Five MER consensus sequences (MER5A1, MER53, MER81, MER91C, and MER117) were matched with paralogous transcripts predicted to be precursor miRNAs in the horse genome (equCab2) and located in either intergenic regions or introns. The MER5A1, MER53, and MER91C sequences obtained from RepeatMasker were matched with the eca-miR-544b, eca-miR-1302, and eca-miR-652 precursor sequences derived from Ensembl transcript database, respectively. Each precursor form was anticipated to yield two mature forms, and we confirmed miRNA expression in six different tissues (cerebrum, cerebellum, lung, spleen, adrenal gland, and duodenum) of one thorough-bred horse. MER5A1-derived miRNAs generally showed significantly higher expression in the lung than in other tissues. MER91C-derived miRNA-5p also showed significantly higher expression in the duodenum than in other tissues (cerebellum, lung, spleen, and adrenal gland). The MER117-overlapped expressed sequence tag generated polycistronic miRNAs, which showed higher expression in the duodenum than other tissues. These data indicate that horse MER transposons encode miR-NAs that are expressed in several tissues and are thought to have biological functions.

Analysis of the genome sequence of the pathogenic Muscovy duck parvovirus strain YY reveals a 14-nucleotide-pair deletion in the inverted terminal repeats

Genomic information about Muscovy duck parvovirus is still limited. In this study, the genome of the pathogenic MDPV strain YY was sequenced. The full-length genome of YY is 5075 nucleotides (nt) long, 57nt shorter than that of strain FM. Sequence alignment indicates that the 5' and 3' inverted terminal repeats (ITR) of strain YY contain a 14-nucleotide-pair deletion in the stem of the palindromic hairpin structure in comparison to strain FM and FZ91-30. The deleted region contains one "E-box" site and one repeated motif with the sequence "TTCCGGT" or "ACCGGAA". Phylogenetic trees constructed based the protein coding genes concordantly showed that YY, together with nine other MDPV isolates from various places, clustered in a separate branch, distinct from the branch formed by goose parvovirus (GPV) strains. These results demonstrate that, despite the distinctive deletion, the YY strain still belongs to the classical MDPV group. Moreover, the deletion of ITR may contribute to the genome evolution of MDPV under immunization pressure.

Analysis of cis and trans Requirements for DNA Replication at the Right-End Hairpin of the Human Bocavirus 1 Genome.

Parvoviruses are single-stranded DNA viruses that use the palindromic structures at the ends of the viral genome for their replication. The mechanism of parvovirus replication has been studied mostly in the dependoparvovirus adeno-associated virus 2 (AAV2) and the protoparvovirus minute virus of mice (MVM). Here, we used human bocavirus 1 (HBoV1) to understand the replication mechanism of bocaparvovirus. HBoV1 is pathogenic to humans, causing acute respiratory tract infections, especially in young children under 2 years old. By using the duplex replicative form of the HBoV1 genome in human embryonic kidney 293 (HEK293) cells, we identified the HBoV1 minimal replication origin at the right-end hairpin (OriR). Mutagenesis analyses confirmed the putative NS1 binding and nicking sites within the OriR. Of note, unlike the large nonstructural protein (Rep78/68 or NS1) of other parvoviruses, HBoV1 NS1 did not specifically bind OriR in vitro, indicating that other viral and cellular components or the oligomerization of NS1 is required for NS1 binding to the OriR. In vivo studies demonstrated that residues responsible for NS1 binding and nicking are within the origin-binding domain. Further analysis identified that the small nonstructural protein NP1 is required for HBoV1 DNA replication at OriR. NP1 and other viral nonstructural proteins (NS1 to NS4) colocalized within the viral DNA replication centers in both OriR-transfected cells and virus-infected cells, highlighting a direct involvement of NP1 in viral DNA replication at OriR. Overall, our study revealed the characteristics of HBoV1 DNA replication at OriR, suggesting novel characteristics of autonomous parvovirus DNA replication. Human bocavirus 1 (HBoV1) causes acute respiratory tract infections in young children. The duplex HBoV1 genome replicates in HEK293 cells and produces progeny virions that are infectious in well-differentiated airway epithelial cells. A recombinant AAV2 vector pseudotyped with an HBoV1 capsid has been developed to efficiently deliver the cystic fibrosis transmembrane conductance regulator gene to human airway epithelia. Here, we identified both cis-acting elements and trans-acting proteins that are required for HBoV1 DNA replication at the right-end hairpin in HEK293 cells. We localized the minimal replication origin, which contains both NS1 nicking and binding sites, to a 46-nucleotide sequence in the right-end hairpin. The identification of these essential elements of HBoV1 DNA replication acting both in cis and in trans will provide guidance to develop antiviral strategies targeting viral DNA replication at the right-end hairpin and to design next-generation recombinant HBoV1 vectors, a promising tool for gene therapy of lung diseases.

Single primer-mediated circular polymerase chain reaction for hairpin DNA cloning and plasmid editing

We developed and validated a universal polymerase chain reaction (PCR) method, single primer circular (SPC)–PCR, using single primer to simultaneously insert and amplify a short hairpin sequence into a vector with a high success rate. In this method, the hairpin structure is divided into two parts and fused into a vector by PCR. Then, a single primer is used to cyclize the chimera into a mature short hairpin RNA (shRNA) expression vector. It is not biased by loop length or palindromic structures. Six hairpin DNAs with short 4-nucleotide loops were successfully cloned. Moreover, SPC–PCR was also applied to plasmid editing within 3 h with a success rate higher than 95%.

ISC, a Novel Group of Bacterial and Archaeal DNA Transposons That Encode Cas9 Homologs.

ABSTRACTBacterial genomes encode numerous homologs of Cas9, the effector protein of the type II CRISPR-Cas systems. The homology region includes the arginine-rich helix and the HNH nuclease domain that is inserted into the RuvC-like nuclease domain. These genes, however, are not linked to cas genes or CRISPR. Here, we show that Cas9 homologs represent a distinct group of nonautonomous transposons, which we denote ISC (insertion sequences Cas9-like). We identify many diverse families of full-length ISC transposons and demonstrate that their terminal sequences (particularly 3′ termini) are similar to those of IS605 superfamily transposons that are mobilized by the Y1 tyrosine transposase encoded by the TnpA gene and often also encode the TnpB protein containing the RuvC-like endonuclease domain. The terminal regions of the ISC and IS605 transposons contain palindromic structures that are likely recognized by the Y1 transposase. The transposons from these two groups are inserted either exactly in the middle or upstream of specific 4-bp target sites, without target site duplication. We also identify autonomous ISC transposons that encode TnpA-like Y1 transposases. Thus, the nonautonomous ISC transposons could be mobilized in trans either by Y1 transposases of other, autonomous ISC transposons or by Y1 transposases of the more abundant IS605 transposons. These findings imply an evolutionary scenario in which the ISC transposons evolved from IS605 family transposons, possibly via insertion of a mobile group II intron encoding the HNH domain, and Cas9 subsequently evolved via immobilization of an ISC transposon.IMPORTANCE Cas9 endonucleases, the effectors of type II CRISPR-Cas systems, represent the new generation of genome-engineering tools. Here, we describe in detail a novel family of transposable elements that encode the likely ancestors of Cas9 and outline the evolutionary scenario connecting different varieties of these transposons and Cas9.