The most frequent short sequences in non-coding DNA

Exploring objective feature sets in constructing the evolution relationship of animal genome sequences

BACKGROUND

Exploring evolution regularities of genome sequences and constructing more objective species evolution relationships at the genomic level are high-profile topics. Based on the evolution mechanism of genome sequences proposed in our previous research, we found that only the 8-mers containing CG or TA dinucleotides correlate directly with the evolution of genome sequences, and the relative frequency rather than the actual frequency of these 8-mers is more suitable to characterize the evolution of genome sequences.

RESULT

Therefore, two types of feature sets were obtained, they are the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers. The evolution relationships of mammals and reptiles were constructed by the relative frequency set of CG1 + CG2 8-mers, and two types of evolution relationships of insects were constructed by the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers respectively. Through comparison and analysis, we found that evolution relationships are consistent with the known conclusions. According to the evolution mechanism, we considered that the evolution relationship constructed by CG1 + CG2 8-mers reflects the evolution state of genome sequences in current time, and the evolution relationship constructed by TA1 + TA2 8-mers reflects the evolution state in the early stage.

CONCLUSION

Our study provides objective feature sets in constructing evolution relationships at the genomic level.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

Sequence analysis of wheat subtelomeres reveals a high polymorphism among homoeologous chromosomes

Bread wheat, Triticum aestivum L., is one of the most important crops in the world. Understanding its genome organization (allohexaploid; AABBDD; 2n = 6x = 42) is essential for geneticists and plant breeders. Particularly, the knowledge of how homologous chromosomes (equivalent chromosomes from the same genome) specifically recognize each other to pair at the beginning of meiosis, the cellular process to generate gametes in sexually reproducing organisms, is fundamental for plant breeding and has a big influence on the fertility of wheat plants. Initial homologous chromosome interactions contribute to specific recognition and pairing between homologues at the onset of meiosis. Understanding the molecular basis of these critical processes can help to develop genetic tools in a breeding context to promote interspecific chromosome associations in hybrids or interspecific genetic crosses to facilitate the transfer of desirable agronomic traits from related species into a crop like wheat. The terminal regions of chromosomes, which include telomeres and subtelomeres, participate in chromosome recognition and pairing. We present a detailed molecular analysis of subtelomeres of wheat chromosome arms 1AS, 4AS, 7AS, 7BS and 7DS. Results showed a high polymorphism in the subtelomeric region among homoeologues (equivalent chromosomes from related genomes) for all the features analyzed, including genes, transposable elements, repeats, GC content, predicted CpG islands, recombination hotspots and targeted sequence motifs for relevant DNA-binding proteins. These polymorphisms might be the molecular basis for the specificity of homologous recognition and pairing in initial chromosome interactions at the beginning of meiosis in wheat.

Single-strand annealing mediates the conservative repair of double-strand DNA breaks in homologous recombination-defective germ cells of Caenorhabditis elegans

A missense mutation in C. elegans RAD-54, a homolog of RAD54 that operates in the homologous recombination (HR) pathway, was found to decrease ATPase activity in vitro. The hypomorphic mutation caused hypersensitivity of C. elegans germ cells to double-strand DNA breaks (DSBs). Although the formation of RAD-51 foci at DSBs was normal in both the mutant and knockdown worms, their subsequent dissipation was slow. The rad-54-deficient phenotypes were greatly aggravated when combined with an xpf-1 mutation, suggesting a conservative role of single-strand annealing (SSA) for DSB repair in HR-defective worms. The phenotypes of doubly-deficient rad-54;xpf-1 worms were partially suppressed by a mutation of lig-4, a nonhomologous end-joining (NHEJ) factor. In summary, RAD-54 is required for the dissociation of RAD-51 from DSB sites in C. elegans germ cells. Also, NHEJ and SSA exert negative and positive effects, respectively, on genome stability when HR is defective.

The Poitiers School of Mathematical and Theoretical Biology: Besson-Gavaudan-Schützenberger's Conjectures on Genetic Code and RNA Structures

The French school of theoretical biology has been mainly initiated in Poitiers during the sixties by scientists like J. Besson, G. Bouligand, P. Gavaudan, M. P. Schützenberger and R. Thom, launching many new research domains on the fractal dimension, the combinatorial properties of the genetic code and related amino-acids as well as on the genetic regulation of the biological processes. Presently, the biological science knows that RNA molecules are often involved in the regulation of complex genetic networks as effectors, e.g., activators (small RNAs as transcription factors), inhibitors (micro-RNAs) or hybrids (circular RNAs). Examples of such networks will be given showing that (1) there exist RNA "relics" that have played an important role during evolution and have survived in many genomes, whose probability distribution of their sub-sequences is quantified by the Shannon entropy, and (2) the robustness of the dynamics of the networks they regulate can be characterized by the Kolmogorov-Sinaï dynamic entropy and attractor entropy.

Citosine-Adenine-Repeat Microsatellite of 11β-hydroxysteroid dehydrogenase 2 Gene in Hypertensive Children

BACKGROUND

The impairment of 11β-hydroxysteroid dehydrogenase type 2 enzyme (11βHSD2) results in an inefficient conversion of cortisol to cortisone, which triggers hypertension. Cytosine-adenine repeat (CA repeat) microsatellite has been associated with low HSD11B2 gene expression.

AIM

To determine whether the CA-repeat length in intron 1 affect the serum cortisol to cortisone (F/E) ratio and/or blood pressure (BP) levels in pediatric subjects.

SUBJECTS AND METHODS

Eighty-one hypertensive (HT) and 117 normotensive (NT) subjects participated in this study. We measured BP levels, as well as the F and E and F/E ratio in morning sera and 12-hour urine samples. The length of CA repeats was determined through fragment analysis. We compared the allele distribution between the HT and NT groups, and the patients were dichotomized into groups with short alleles (S) (<21 CA repeats) or long alleles (L), and also in groups according genotype (allele combination: S/S and S/L + L/L).

RESULTS

We found no differences in the distribution of CA-repeat allelic length between the NT and HT groups (P = 0.7807), and there was no correlation between the CA-repeat allelic length and BP (P = 0.1151) levels or the serum F/E ratio (P = 0.6778). However, the serum F/E ratio was higher in the HT group than in the NT group (P = 0.0251). The serum F/E ratio was associated with systolic BP index independent of body mass index only in HT group.

CONCLUSIONS

The CA-repeat length did not influence BP levels or serum F/E ratios in pediatric subjects. However, the serum F/E ratio was associated with BP, suggesting a role of 11βHSD2 in mineralocorticoid hypertension.

High evolutionary turnover of satellite families in Caenorhabditis

Background

The high density of tandem repeat sequences (satellites) in nematode genomes and the availability of genome sequences from several species in the group offer a unique opportunity to better understand the evolutionary dynamics and the functional role of these sequences. We take advantage of the previously developed SATFIND program to study the satellites in four Caenorhabditis species and investigate these questions.

Methods

The identification and comparison of satellites is carried out in three steps. First we find all the satellites present in each species with the SATFIND program. Each satellite is defined by its length, number of repeats, and repeat sequence. Only satellites with at least ten repeats are considered. In the second step we build satellite families with a newly developed alignment program. Satellite families are defined by a consensus sequence and the number of satellites in the family. Finally we compare the consensus sequence of satellite families in different species.

Results

We give a catalog of individual satellites in each species. We have also identified satellite families with a related sequence and compare them in different species. We analyze the turnover of satellites: they increased in size through duplications of fragments of 100-300 bases. It appears that in many cases they have undergone an explosive expansion. In C. elegans we have identified a subset of large satellites that have strong affinity for the centromere protein CENP-A. We have also compared our results with those obtained from other species, including one nematode and three mammals.

Conclusions

Most satellite families found in Caenorhabditis are species-specific; in particular those with long repeats. A subset of these satellites may facilitate the formation of kinetochores in mitosis. Other satellite families in C. elegans are either related to Helitron transposons or to meiotic pairing centers.

Electronic supplementary material

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

Altered LINE-1 Methylation in Mothers of Children with Down Syndrome

Down syndrome (DS, also known as trisomy 21) most often results from chromosomal nondisjunction during oogenesis. Numerous studies sustain a causal link between global DNA hypomethylation and genetic instability. It has been suggested that DNA hypomethylation might affect the structure and dynamics of chromatin regions that are critical for chromosome stability and segregation, thus favouring chromosomal nondisjunction during meiosis. Maternal global DNA hypomethylation has not yet been analyzed as a potential risk factor for chromosome 21 nondisjunction. This study aimed to asses the risk for DS in association with maternal global DNA methylation and the impact of endogenous and exogenous factors that reportedly influence DNA methylation status. Global DNA methylation was analyzed in peripheral blood lymphocytes by quantifying LINE-1 methylation using the MethyLight method. Levels of global DNA methylation were significantly lower among mothers of children with maternally derived trisomy 21 than among control mothers (P = 0.000). The combination of MTHFR C677T genotype and diet significantly influenced global DNA methylation (R2 = 4.5%, P = 0.046). The lowest values of global DNA methylation were observed in mothers with MTHFR 677 CT+TT genotype and low dietary folate. Although our findings revealed an association between maternal global DNA hypomethylation and trisomy 21 of maternal origin, further progress and final conclusions regarding the role of global DNA methylation and the occurrence of trisomy 21 are facing major challenges.

Conformations of p53 response elements in solution deduced using site-directed spin labeling and Monte Carlo sampling

The tumor suppressor protein p53 regulates numerous signaling pathways by specifically recognizing diverse p53 response elements (REs). Understanding the mechanisms of p53-DNA interaction requires structural information on p53 REs. However, such information is limited as a 3D structure of any RE in the unbound form is not available yet. Here, site-directed spin labeling was used to probe the solution structures of REs involved in p53 regulation of the p21 and Bax genes. Multiple nanometer distances in the p21-RE and BAX-RE, measured using a nucleotide-independent nitroxide probe and double-electron-electron-resonance spectroscopy, were used to derive molecular models of unbound REs from pools of all-atom structures generated by Monte-Carlo simulations, thus enabling analyses to reveal sequence-dependent DNA shape features of unbound REs in solution. The data revealed distinct RE conformational changes on binding to the p53 core domain, and support the hypothesis that sequence-dependent properties encoded in REs are exploited by p53 to achieve the energetically most favorable mode of deformation, consequently enhancing binding specificity. This work reveals mechanisms of p53-DNA recognition, and establishes a new experimental/computational approach for studying DNA shape in solution that has far-reaching implications for studying protein-DNA interactions.

A satellite explosion in the genome of holocentric nematodes

Centromere sequences in the genome are associated with the formation of kinetochores, where spindle microtubules grow in mitosis. Centromere sequences usually have long tandem repeats (satellites). In holocentric nematodes it is not clear how kinetochores are formed during mitosis; they are distributed throughout the chromosomes. For this reason it appeared of interest to study the satellites in nematodes in order to determine if they offer any clue on how kinetochores are assembled in these species. We have studied the satellites in the genome of six nematode species. We found that the presence of satellites depends on whether the nematode chromosomes are holocentric or monocentric. It turns out that holocentric nematodes are unique because they have a large number of satellites scattered throughout their genome. Their number, length and composition are different in each species: they apparently have very little evolutionary conservation. In contrast, no scattered satellites are found in the monocentric nematode Trichinella spiralis. It appears that the absence/presence of scattered satellites in the genome distinguishes monocentric from holocentric nematodes. We conclude that the presence of satellites is related to the holocentric nature of the chromosomes of most nematodes. Satellites may stabilize a higher order structure of chromatin and facilitate the formation of kinetochores. We also present a new program, SATFIND, which is suited to find satellite sequences.

Evolutionary divergence of intrinsic and trans-regulated nucleosome positioning sequences reveals plastic rules for chromatin organization

The packaging of eukaryotic genomes into nuclesomes plays critical roles in chromatin organization and gene regulation. Studies in Saccharomyces cerevisiae indicate that nucleosome occupancy is partially encoded by intrinsic antinucleosomal DNA sequences, such as poly(A) sequences, as well as by binding sites for trans-acting factors that can evict nucleosomes, such as Reb1 and the Rsc3/30 complex. Here, we use genome-wide nucleosome occupancy maps in 13 Ascomycota fungi to discover large-scale evolutionary reprogramming of both intrinsic and trans determinants of chromatin structure. We find that poly(G)s act as intrinsic antinucleosomal sequences, comparable to the known function of poly(A)s, but that the abundance of poly(G)s has diverged greatly between species, obscuring their antinucleosomal effect in low-poly(G) species such as S. cerevisiae. We also develop a computational method that uses nucleosome occupancy maps for discovering trans-acting general regulatory factor (GRF) binding sites. Our approach reveals that the specific sequences bound by GRFs have diverged substantially across evolution, corresponding to a number of major evolutionary transitions in the repertoire of GRFs. We experimentally validate a proposed evolutionary transition from Cbf1 as a major GRF in pre-whole-genome duplication (WGD) yeasts to Reb1 in post-WGD yeasts. We further show that the mating type switch-activating protein Sap1 is a GRF in S. pombe, demonstrating the general applicability of our approach. Our results reveal that the underlying mechanisms that determine in vivo chromatin organization have diverged and that comparative genomics can help discover new determinants of chromatin organization.

The distribution of alternating AT sequences in eukaryotic genomes suggests a role in homologous chromosome recognition in meiosis

There are general features of chromosome dynamics, such as homologue recognition in early meiosis, which are expected to involve related sequence motifs in non-coding DNA, with a similar distribution in different species. A search for such motifs is presented here. It has been carried out with the CONREPP programme. It has been found that short alternating AT sequences (10-20 bases) have a similar distribution in most eukaryotic organisms, with some exceptions related to unique meiotic features. All other microsatellite and repeat sequences vary significantly in different organisms. It is concluded that the unique structural features and uniform distribution of alternating AT sequences indicate that they may facilitate homologous chromosome pairing in the early preleptotene stage of meiosis. They may also play a role in the compaction of DNA in mitotic chromosomes.

Novel simple sequence repeats (SSRs) detected by ND-FISH in heterochromatin of Drosophila melanogaster

Background

In recent years, substantial progress has been made in understanding the organization of sequences in heterochromatin regions containing single-copy genes and transposable elements. However, the sequence and organization of tandem repeat DNA sequences, which are by far the majority fraction of D. melanogaster heterochromatin, are little understood.

Results

This paper reports that the heterochromatin, as well as containing long tandem arrays of pentanucleotide satellites (AAGAG, AAGAC, AATAT, AATAC and AACAC), is also enriched in other simple sequence repeats (SSRs) such as A, AC, AG, AAG, ACT, GATA and GACA. Non-denaturing FISH (ND-FISH) showed these SSRs to localize to the chromocentre of polytene chromosomes, and was used to map them on mitotic chromosomes. Different distributions were detected ranging from single heterochromatic clusters to complex combinations on different chromosomes. ND-FISH performed on extended DNA fibres, along with Southern blotting, showed the complex organization of these heterochromatin sequences in long tracts, and revealed subclusters of SSRs (several kilobase in length) flanked by other DNA sequences. The chromosomal characterization of C, AAC, AGG, AAT, CCG, ACG, AGC, ATC and ACC provided further detailed information on the SSR content of D. melanogaster at the whole genome level.

Conclusion

These data clearly show the variation in the abundance of different SSR motifs and reveal their non-random distribution within and between chromosomes. The greater representation of certain SSRs in D. melanogaster heterochromatin suggests that its complexity may be greater than previously thought.

A global view of 54,001 single nucleotide polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and their transferability to water buffalo

The Illumina BovineSNP50 BeadChip features 54,001 informative single nucleotide polymorphisms (SNPs) that uniformly span the entire bovine genome. Among them, 52,255 SNPs have locations assigned in the current genome assembly (Btau_4.0), including 19,294 (37%) intragenic SNPs (i.e., located within genes) and 32,961 (63%) intergenic SNPs (i.e., located between genes). While the SNPs represented on the Illumina Bovine50K BeadChip are evenly distributed along each bovine chromosome, there are over 14,000 genes that have no SNPs placed on the current BeadChip. Kernel density estimation, a non-parametric method, was used in the present study to identify SNP-poor and SNP-rich regions on each bovine chromosome. With bandwidth = 0.05 Mb, we observed that most regions have SNP densities within 2 standard deviations of the chromosome SNP density mean. The SNP density on chromosome X was the most dynamic, with more than 30 SNP-rich regions and at least 20 regions with no SNPs. Genotyping ten water buffalo using the Illumina BovineSNP50 BeadChip revealed that 41,870 of the 54,001 SNPs are fully scored on all ten water buffalo, but 6,771 SNPs are partially scored on one to nine animals. Both fully scored and partially/no scored SNPs are clearly clustered with various sizes on each chromosome. However, among 43,687 bovine SNPs that were successfully genotyped on nine and ten water buffalo, only 1,159 were polymorphic in the species. These results indicate that the SNPs sites, but not the polymorphisms, are conserved between two species. Overall, our present study provides a solid foundation to further characterize the SNP evolutionary process, thus improving understanding of within- and between-species biodiversity, phylogenetics and adaption to environmental changes.

The many twists and turns of DNA: template, telomere, tool, and target

If any proof were needed of DNA's versatile roles and use, it is certainly provided by the numerous depositions of new three-dimensional (3D) structures to the coordinate databanks (PDB, NDB) over the last two years. Quadruplex motifs involving G-repeats, adducted sequences and oligo-2'-deoxynucleotides (ODNs) with bound ligands are particularly well represented. In addition, structures of chemically modified DNAs (CNAs) and artificial analogs are yielding insight into stability, pairing properties, and dynamics, including those of the native nucleic acids. Besides being of significance for establishing diagnostic tools and in the analysis of protein-DNA interactions, chemical modification in conjunction with investigations of the structural consequences may yield novel nucleic acid-based therapeutics. DNA's predictable and highly specific pairing behavior makes it the material of choice for constructing 3D-nanostructures of defined architecture. Recently the first examples of DNA nanoparticle and self-assembled 3D-crystals were reported. Although the structures discussed in this review are all based either on X-ray crystallography or solution NMR, small angle X-ray scattering (SAXS), and cryoEM are proving to be useful approaches for the characterization of nanoscale DNA architecture.