Comparative cytogenetics of microsatellite distribution in two tetra fishes Astyanax bimaculatus (Linnaeus, 1758) and Psalidodon scabripinnis (Jenyns, 1842)

Background

The main cytogenetic studies of the Characidae family comprise the genera Astyanax and Psalidodon involving the use of repetitive DNA probes. However, for the microsatellite classes, studies are still scarce and the function of these sequences in the genome of these individuals is still not understood. Thus, we aimed to analyze and compare the distribution of microsatellite sequences in the species Astyanax bimaculatus and Psalidodon scabripinnis.

Methods

We collected biopsies from the fins of A. bimaculatus and P. scabripinnis to perform cell culture, followed by chromosome extraction, and mapped the distribution of 14 microsatellites by FISH in both species.

Results and Discussion

The diploid number observed for both species was 2n = 50, with an acrocentric B microchromosome in A. bimaculatus and a metacentric B chromosome in P. scabripinnis. Regarding FISH, 11 probes hybridized in the karyotype of A. bimaculatus mainly in centromeric regions, and 13 probes hybridized in P. scabripinnis, mainly in telomeric regions, in addition to a large accumulation of microsatellite hybridization on its B chromosome.

Conclusion

Comparative FISH mapping of 14 microsatellite motifs revealed different patterns of distribution both in autosomes and supernumerary chromosomes of A. bimaculatus and P. scabripinnis, suggesting independent evolutionary processes in each of these species, representing excellent data on chromosome rearrangements and cytotaxonomy.

Draft genome of Brasenia schreberi, a worldwide distributed and endangered aquatic plant

OBJECTIVES

Brasenia is a monotypic genus in the family of Cabombaceae. The only species, B. schreberi, is a macrophyte distributed worldwide. Because it requires good water quality, it is endangered in China and other countries due to the deterioration of aquatic habitats. The young leaves and stems of B. schreberi are covered by thick mucilage, which has high medical value. As an allelopathic aquatic plant, it can also be used in the management of aquatic weeds. Here, we present its assembled and annotated genome to help shed light on medial and allelopathic substrates and facilitate their conservation.

DATA DESCRIPTION

Genomic DNA and RNA extracted from B. schreberi leaf tissues were used for whole genome and RNA sequencing using a Nanopore and/or MGI sequencer. The assembly was 1,055,148,839 bp in length, with 92 contigs and an N50 of 22,379,495 bp. The repetitive elements in the assembly were 555,442,205 bp. A completeness assessment of the assembly with BUSCO and compleasm indicated 88.4 and 90.9% completeness in the Eudicots database and 95.4 and 96.6% completeness in the Embryphyta database. Gene annotation revealed 67,747 genes that coded for 73,344 proteins.

Fast satellite DNA evolution in Nothobranchius annual killifishes

Satellite DNA (satDNA) is a rapidly evolving class of tandem repeats, with some monomers being involved in centromere organization and function. To identify repeats associated with (peri)centromeric regions, we investigated satDNA across Southern and Coastal clades of African annual killifishes of the genus Nothobranchius. Molecular cytogenetic and bioinformatic analyses revealed that two previously identified satellites, designated here as NkadSat01-77 and NfurSat01-348, are associated with (peri)centromeres only in one lineage of the Southern clade. NfurSat01-348 was, however, additionally detected outside centromeres in three members of the Coastal clade. We also identified a novel satDNA, NrubSat01-48, associated with (peri)centromeres in N. foerschi, N. guentheri, and N. rubripinnis. Our findings revealed fast turnover of satDNA associated with (peri)centromeres and different trends in their evolution in two clades of the genus Nothobranchius.

Assembly, annotation and analysis of the chloroplast genome of the Algarrobo tree Neltuma pallida (subfamily: Caesalpinioideae)

BACKGROUND

Neltuma pallida is a tree that grows in arid soils in northwestern Peru. As a predominant species of the Equatorial Dry Forest ecoregion, it holds significant economic and ecological value for both people and environment. Despite this, the species is severely threatened and there is a lack of genetic and genomic research, hindering the proposal of evidence-based conservation strategies.

RESULTS

In this work, we conducted the assembly, annotation, analysis and comparison of the chloroplast genome of a N. pallida specimen with those of related species. The assembled chloroplast genome has a length of 162,381 bp with a typical quadripartite structure (LSC-IRA-SSC-IRB). The calculated GC content was 35.97%. However, this is variable between regions, with a higher GC content observed in the IRs. A total of 132 genes were annotated, of which 19 were duplicates and 22 contained at least one intron in their sequence. A substantial number of repetitive sequences of different types were identified in the assembled genome, predominantly tandem repeats (> 300). In particular, 142 microsatellites (SSR) markers were identified. The phylogenetic reconstruction showed that N. pallida grouped with the other Neltuma species and with Prosopis cineraria. The analysis of sequence divergence between the chloroplast genome sequences of N. pallida, N. juliflora, P. farcta and Strombocarpa tamarugo revealed a high degree of similarity.

CONCLUSIONS

The N. pallida chloroplast genome was found to be similar to those of closely related species. With a size of 162,831 bp, it had the classical chloroplast quadripartite structure and GC content of 35.97%. Most of the 132 identified genes were protein-coding genes. Additionally, over 800 repetitive sequences were identified, including 142 SSR markers. In the phylogenetic analysis, N. pallida grouped with other Neltuma spp. and P. cineraria. Furthermore, N. pallida chloroplast was highly conserved when compared with genomes of closely related species. These findings can be of great potential for further diversity studies and genetic improvement of N. pallida.

A new type of tandem repeats in Myotis petax (Chiroptera, Vespertilionidae) mitochondrial control region

BACKGROUND

Tandem repeats in mitochondrial DNA control region are known to different animal taxa, including bat species of the family Vespertilionidae. The long R1-repeats in the bat ETAS-domain are often presented in a variable copy number and may exhibit both inter-individual and intra-individual sequence diversity. The function of repeats in the control region is still unclear, but it has been shown that repetitive sequences in some animal groups (shrews, cats and sheep) may include parts of ETAS1 and ETAS2 conservative blocks of mitochondrial DNA.

METHODS AND RESULTS

Analysis of the control region sequences for 31 Myotis petax specimens allowed the identification of the inter-individual variability and clarification of the composition of the R1-repeats. The copy number of the R1-repeats varies from 4 to 7 in individuals. The specimens examined do not exhibit a size heteroplasmy previously described for Myotis species. The unusual short 30 bp R1-repeats have been detected in M. petax for the first time. The ten specimens from Amur Region and Primorsky Territory have one or two copies of these additional repeats.

CONCLUSIONS

It was determined that the R1-repeats in M. petax control region consist of parts of the ETAS1 and ETAS2 blocks. The origin of the additional repeats seems to be related to the 51 bp deletion in the central part of the R1-repeat unit and subsequent duplication. Comparison of repetitive sequences in the control region of closely-related Myotis species identified the occurrence of incomplete repeats also resulting from the short deletions, but distinct from additional repeats of M. petax.

Direct detection of alpha satellite DNA with single-base resolution by using abasic Peptide Nucleic Acids and Fluorescent in situ Hybridization

The detection of repetitive sequences with single-base resolution is becoming increasingly important aiming to understand the biological implications of genomic variation in these sequences. However, there is a lack of techniques to experimentally validate sequencing data from repetitive sequences obtained by Next-Generation Sequencing methods, especially in the case of Single-Nucleotide Variations (SNVs). That is one of the reasons why repetitive sequences have been poorly studied and excluded from most genomic studies. Therefore, in addition to sequencing data, there is an urgent need for efficient validation methods of genomic variation in these sequences. Herein we report the development of chemFISH, an alternative method for the detection of SNVs in repetitive sequences. ChemFISH is an innovative method based on dynamic chemistry labelling and abasic Peptide Nucleic Acid (PNA) probes to detect in situ the α-satellite DNA, organized in tandem repeats, with single-base resolution in a direct and rapid reaction. With this approach, we detected by microscopy the α-satellite DNA in a variety of human cell lines, we quantified the detection showing a low coefficient of variation among samples (13.16%-25.33%) and we detected single-base specificity with high sensitivity (82.41%-88.82%). These results indicate that chemFISH can serve as a rapid method to validate previously detected SNVs in sequencing data, as well as to find novel SNVs in repetitive sequences. Furthermore, the versatile chemistry behind chemFISH can lead to develop novel molecular assays for the in situ detection of nucleic acids.

SINE-B1 Distribution and Chromosome Rearrangements in the South American Proechimys gr. goeldii (Echimyidae, Rodentia)

Proechimys species are remarkable for their extensive chromosome rearrangements, representing a good model to understand genome evolution. Herein, we cytogenetically analyzed 3 different cytotypes of Proechimys gr. goeldii to assess their evolutionary relationship. We also mapped the transposable element SINE-B1 on the chromosomes of P. gr. goeldii in order to investigate its distribution among individuals and evaluate its possible contribution to karyotype remodeling in this species. SINE-B1 showed a dispersed distribution along chromosome arms and was also detected at the pericentromeric regions of some chromosomes, including pair 1 and the sex chromosomes, which are involved in chromosome rearrangements. In addition, we describe a new cytotype for P. gr. goeldii, reinforcing the significant role of gross chromosomal rearrangements during the evolution of the genus. The results of FISH with SINE-B1 suggest that this issue should be more deeply investigated for a better understanding of its role in the mechanisms involved in the wide variety of Proechimys karyotypes.

New methodology for repetitive sequences identification in human X and Y chromosomes

Repetitive DNA sequences occupy the major proportion of DNA in the human genome and even in the other species' genomes. The importance of each repetitive DNA type depends on many factors: structural and functional roles, positions, lengths and numbers of these repetitions are clear examples. Conserving such DNA sequences or not in different locations in the chromosome remains a challenge for researchers in biology. Detecting their location despite their great variability and finding novel repetitive sequences remains a challenging task. To side-step this problem, we developed a new method based on signal and image processing tools. In fact, using this method we could find repetitive patterns in DNA images regardless of the repetition length. This new technique seems to be more efficient in detecting new repetitive sequences than bioinformatics tools. In fact, the classical tools present limited performances especially in case of mutations (insertion or deletion). However, modifying one or a few numbers of pixels in the image doesn't affect the global form of the repetitive pattern. As a consequence, we generated a new repetitive patterns database which contains tandem and dispersed repeated sequences. The highly repetitive sequences, we have identified in X and Y chromosomes, are shown to be located in other human chromosomes or in other genomes. The data we have generated is then taken as input to a Convolutional neural network classifier in order to classify them. The system we have constructed is efficient and gives an average of 94.4% as recognition score.

How diverse is heterochromatin in the Caesalpinia group? Cytogenomic characterization of Erythrostemon hughesii Gagnon & G.P. Lewis (Leguminosae: Caesalpinioideae)

Cytogenomic characterization of Erythrostemon hughesii reveals a heterogeneity of repeats in its subtelomeric heterochromatin. Comparative analyses with other Caesalpinia group species reveal a significant reduction in the abundance of Ty3-gypsy/Chromovirus Tekay retrotransposons during its evolution. In numerically stable karyotypes, repetitive DNA variability is one of the main causes of genome and chromosome variation and evolution. Species from the Caesalpinia group (Leguminosae) are karyotypically characterized by 2n = 24, with small chromosomes and highly variable CMA⁺ heterochromatin banding patterns that correlate with environmental variables. Erythrostemon hughesii differs from other species of the group examined to date for having subtelomeric CMA⁺ bands; this contrasts with most species in the group which have proximal bands. Here we analyse the repeatome of E. hughesii using genome skimming and chromosomal mapping approaches to characterize the identity of the most abundant repetitive elements and their physical location. The repetitive fraction of E. hughesii comprises 28.73% of the genome. The most abundant elements were retrotransposons (RT) with long terminal repeats (LTR-RT; 9.76%) and satellite DNAs (7.83%). Within the LTR-RTs, the most abundant lineages were: Ty1/copia-Ale (1%), Ty3/gypsy CRM (0.88%) and Ty3/gypsy Athila (0.75%). Using fluorescent in situ hybridization four satellite DNAs and several LTR-RT elements were shown to be present in most subtelomeric CMA⁺ bands. These results highlight how the repeatome in E. hughesii, a species from Oaxaca state in Mexico, is clearly distinct from Northeast Brazilian species of the Caesalpinia group, mainly due to its high diversity of repeats in its subtelomeric heterochromatic bands and low amount of LTR-RT Ty3/gypsy-Tekay elements. Comparative sequence analysis of Tekay elements from different species is congruent with a clade-specific origin of this LTR-RT after the divergence of the Caesalpinia group. We hypothesize that repeat-rich heterochromatin may play a role in leading to faster genomic divergence between individuals, increasing speciation and diversification.

Heterochromatin evolution in Arachis investigated through genome-wide analysis of repetitive DNA

The most conspicuous difference among chromosomes and genomes in Arachis species, the patterns of heterochromatin, was mainly modeled by differential amplification of different members of one superfamily of satellite DNAs. Divergence in repetitive DNA is a primary driving force for genome and chromosome evolution. Section Arachis is karyotypically diverse and has six different genomes. Arachis glandulifera (D genome) has the most asymmetric karyotype and the highest reproductive isolation compared to the well-known A and B genome species. These features make A. glandulifera an interesting model species for studying the main repetitive components that accompanied the genome and chromosome diversification in the section. Here, we performed a genome-wide analysis of repetitive sequences in A. glandulifera and investigated the chromosome distribution of the identified satellite DNA sequences (satDNAs). LTR retroelements, mainly the Ty3-gypsy families "Fidel/Feral" and "Pipoka/Pipa", were the most represented. Comparative analyses with the A and B genomes showed that many of the previously described transposable elements (TEs) were differently represented in the D genome, and that this variation accompanied changes in DNA content. In addition, four major satDNAs were characterized. Agla_CL8sat was the major component of pericentromeric heterochromatin, while Agla_CL39sat, Agla_CL69sat, and Agla_CL122sat were found in heterochromatic and/or euchromatic regions. Even though Agla_CL8sat belong to a different family than that of the major satDNA (ATR-2) found in the heterochromatin of the A, K, and F genomes, both satDNAs are members of the same superfamily. This finding suggests that closely related satDNAs of an ancestral library were differentially amplified leading to the major changes in the heterochromatin patterns that accompanied the karyotype and genome differentiation in Arachis.

Genotyping of Candida albicans isolated from animals using 25S ribosomal DNA and ALT repeats polymorphism in repetitive sequence

Background and Purpose:

Candida albicans is the most prevalent Candida species isolated from animals. Candidiasis can be systemic in animals or may affect a single organ, such as the mouth, urinary tract, and skin. The aim of the present study was to determine the genetic diversity of C. albicans isolated from different animals and investigate the presence of a relationship between host specificity and genetic typing of C. albicans.

Materials and Methods:

For the purpose of the study, DNA extraction was performed on 27 clinical isolates of C. albicans obtained from animals. Subsequently, they were subjected to 25S ribosomal DNA amplification and ALT repeats in repetitive sequences (RPSs). The minimum inhibitory concentrations of fluconazole, ketoconazole, clotrimazole, nystatin, amphotericin B, and caspofungin were determined using the microdilution method based on the Clinical and Laboratory Standards Institute M27-S4 standard.

Results:

Out of 27 C. albicans strains, 11, 6, 5, and 5 cases were recognized as genotypes A (40.8%), E (22.2%), B (18.5%), and C (18.5%), respectively, through amplification using AS-I, which revealed 17 different types of C. albicans. By combining the two typing methods, 27 C. albicans strains were finally divided into 22 genotypes.

Conclusion:

Different genotypes showed genetic diversity among the C. albicans strains isolated from animal sources. The results revealed no special genotype relationship according to the host, anatomical source of isolation, and antifungal susceptibility.

Progerin expression induces a significant downregulation of transcription from human repetitive sequences in iPSC-derived dopaminergic neurons

Repetitive DNA sequences represent about half of the human genome. They have a central role in human biology, especially neurobiology, but are notoriously difficult to study. The purpose of this study was to quantify the transcription from repetitive sequences in a progerin-expressing cellular model of neuronal aging. Progerin is a nuclear protein causative of the Hutchinson-Gilford progeria syndrome that is also incrementally expressed during the normal aging process. A dedicated pipeline of analysis allowed to quantify transcripts containing repetitive sequences from RNAseq datasets oblivious of their genomic localization, tolerating a sufficient degree of mutational noise, all with low computational requirements. The pipeline has been applied to a published panel of RNAseq datasets derived from a well-established and well-described cellular model of aging of dopaminergic neurons. Progerin expression strongly downregulated the transcription from all the classes of repetitive sequences: satellites, long and short interspersed nuclear elements, human endogenous retroviruses, and DNA transposon. The Alu element represented by far the principal source of transcript originating either from repetitive sequences or from canonical coding genes; it was expressed on average at 192,493.5 reads per kilobase million (RPKM) (SE = 21,081.3) in the control neurons and dropped to 43,760.1 RPKM (SE = 5315.0) in the progerin-expressing neurons, being significant downregulated (p = 0.0005). The results highlighted a global perturbation of transcripts derived from repetitive sequences in a cellular model of aging and provided a direct link between progerin expression and alteration of transcription from human repetitive elements.

The genome sequence of 'Kurodagosun', a major carrot variety in Japan and China, reveals insights into biological research and carrot breeding

Carrot (Daucus carota L.) is one of the most economically important root vegetables in the world, providing numerous nutrients for human health. China is the largest country of carrot production in the world, and 'Kurodagosun' has been a major carrot variety in China. Carrot material used in this study was the inbred line 'DC-27', which was derived by forced selfing from 'Kurodagosun'. To understand the genetic system and plant-specific genes of 'Kurodagosun', we report the draft genome sequence of carrot 'DC-27' assembled using a combination of Roche454 and HiSeq 2000 sequencing technologies to achieve 32-fold genome coverage. A total of 31,891 predicted genes were identified. These assembled sequences provide candidate genes involved in biological processes including stress response and carotenoid biosynthesis. Genomic sequences corresponding to 371.6 Mb was less than 473 Mb, which is the estimated genome size. The availability of a draft sequence of the 'DC-27' genome advances knowledge on the biological research and breeding of carrot, as well as other Apiaceae plants. The 'DC-27' genome sequence data also provide a new resource to explore the evolution of other higher plants.

Genomic Organization of Repetitive DNAs and Differentiation of an XX/XY Sex Chromosome System in the Amazonian Puffer Fish, Colomesus asellus (Tetraodontiformes)

The genus Colomesus is the sole representative of the family Tetraodontidae in the Amazon region. Here, Colomesus asellus was analyzed using conventional and molecular cytogenetic protocols. Its diploid chromosome number is 2n = 46 with 12 meta-, 10 submeta-, 16 subtelo-, and 8 acrocentric chromosomes and a fundamental number of FN = 84. An XX/XY sex chromosome system was identified. Mapping of 18S rDNA correlated with the nucleolus organizer regions (Ag-NORs) in the short arms of the 2 X chromosomes in females and in the Y chromosome in males. C-banding revealed heterochromatin in the centromeric regions of all chromosomes, except for pair 3. Prominent sex chromosome-specific heterochromatin amplification was observed, covering the short arms of the Y chromosome almost entirely. FISH with telomeric and tropomyosin (tpm1) sequences, respectively, revealed terminal signals in all chromosomes. The analysis of extended DNA fibers confirmed the colocalization and the interspersed pattern of the telomeric and tpm1 sequences. Thus, this study highlights the remarkable evolutionary dynamism presented by the Amazonian puffer fish regarding the differentiation of a heteromorphic XY sex chromosome system and a particular sex-specific amplification of rDNA sites. This is the first record of such an association in the Tetraodontidae family.

Genome size and identification of abundant repetitive sequences in Vallisneria spinulosa

Vallisneria spinulosa is a freshwater aquatic plant of ecological and economic importance. However, there is limited cytogenetic and genomics information on Vallisneria. In this study, we measured the nuclear DNA content of Vallisneria spinulosa by flow cytometry, performed a de novo assembly, and annotated repetitive sequences by using a combination of next-generation sequencing (NGS) and bioinformatics tools. The genome size of Vallisneria spinulosa is approximately 3,595 Mbp, in which nearly 60% of the genome consists of repetitive sequences. The majority of the repetitive sequences are LTR-retrotransposons comprising 43% of the genome. Although the amount of sequencing data used in this study was not sufficient for a whole-genome assembly, it could generate an overview of representative elements in the genome. These results will lay a new foundation for further studies on various species that belong to the Vallisneria genus.

Comparative cytogenetics in three Sciaenid species (Teleostei, Perciformes): evidence of interspecific chromosomal diversification

Background

Species belonging to the Sciaenidae family present a karyotype composed by 48 acrocentric chromosomes and are thus considered a striking example of chromosomal conservation. In this family, three species are extensively studied including Larimichthys crocea, Larimichthys polyactis and Nibea albiflora due to their importance in fishery and aquaculture in East Asia. Despite abundant data of population genetics available for some of them, cytogenetic information on these species is still scarce and obtained by conventional cytogenetic protocols. Therefore, a more detailed cytogenomic investigation was performed in these species to analyze their karyotype differentiation using conventional staining techniques and fluorescence in situ hybridization to map several repetitive DNAs.

Results

The three species showed a slight karyotype differentiation with 4sm + 2st + 42a in L. polyactis, 20st + 28a in L. crocea and 48a in N. albiflora. Additionally, the mapping of repetitive sequences further revealed a number of interspecific differences among them. Particularly, 18S and 5S rDNA sites showed syntenic arrangements in N. albiflora and non-syntenic arrangements in both Larimichthys species. The microsatellites (CA)₁₅ and (GA)₁₅ showed conspicuous terminal clusters in some chromosomes of all species. On the other hand, (CGG)₁₀ repeats, Rex6 elements and U2 snRNA displayed a scattered distribution on the chromosomes.

Conclusions

Although the three Sciaenid species examined displayed a general pattern of karyotypic conservatism, we explored chromosomal diversification among them. The diversificated karyotypic macrostructure is followed by intergeneric evolutionary diversification of the repetitive sequences. The data indicate some degree of intergeneric evolutionary diversification at chromosomal level, and suggest the evolutionary dynamics among Sciaenid species, higher than previously thought. The present cytogenetic data provide new insight into the chromosomal diversification in Sciaenidae, and contribute to inferring the chromosomal rearrangements and trends of karyotype evolution in this fish group.

MetaCRAST: reference-guided extraction of CRISPR spacers from unassembled metagenomes

Clustered regularly interspaced short palindromic repeat (CRISPR) systems are the adaptive immune systems of bacteria and archaea against viral infection. While CRISPRs have been exploited as a tool for genetic engineering, their spacer sequences can also provide valuable insights into microbial ecology by linking environmental viruses to their microbial hosts. Despite this importance, metagenomic CRISPR detection remains a major challenge. Here we present a reference-guided CRISPR spacer detection tool (Metagenomic CRISPR Reference-Aided Search Tool-MetaCRAST) that constrains searches based on user-specified direct repeats (DRs). These DRs could be expected from assembly or taxonomic profiles of metagenomes. We compared the performance of MetaCRAST to those of two existing metagenomic CRISPR detection tools-Crass and MinCED-using both real and simulated acid mine drainage (AMD) and enhanced biological phosphorus removal (EBPR) metagenomes. Our evaluation shows MetaCRAST improves CRISPR spacer detection in real metagenomes compared to the de novo CRISPR detection methods Crass and MinCED. Evaluation on simulated metagenomes show it performs better than de novo tools for Illumina metagenomes and comparably for 454 metagenomes. It also has comparable performance dependence on read length and community composition, run time, and accuracy to these tools. MetaCRAST is implemented in Perl, parallelizable through the Many Core Engine (MCE), and takes metagenomic sequence reads and direct repeat queries (FASTA or FASTQ) as input. It is freely available for download at https://github.com/molleraj/MetaCRAST.

In silico identification and characterization of stress and virulence associated repeats in Salmonella

So much genomic similarities yet causing different diseases, is like a paradox in Salmonella biology. Repeat is one of the probes that can explain such differences. Here, a comparative genomics approach is followed to identify and characterize repeats that might play role in adaptation and pathogenesis. Repeats are non-randomly distributed in the genomes except few typhoid causing strains. Perfect long repeats are rare compare to polymorphic ones and both are statistically consistent. Significant differences in repeat densities in stress related genes manifest its probable participation in survival and virulence. 573 and 1053 repeat loci have been identified which are exclusively associated with stress and virulent genes respectively. In Salmonella Typhi, an octameric VNTR locus is found in between acrD and yffB genes having more than 25 perfect copies across Salmonella Typhi but possesses only single copy in other serovars. This repeat can be used as a diagnostic probe for typhoid.

Role of specialized DNA polymerases in the limitation of replicative stress and DNA damage transmission

Replication stress is a strong and early driving force for genomic instability and tumor development. Beside replicative DNA polymerases, an emerging group of specialized DNA polymerases is involved in the technical assistance of the replication machinery in order to prevent replicative stress and its deleterious consequences. During S-phase, altered progression of the replication fork by endogenous or exogenous impediments induces replicative stress, causing cells to reach mitosis with genomic regions not fully duplicated. Recently, specific mechanisms to resolve replication intermediates during mitosis with the aim of limiting DNA damage transmission to daughter cells have been identified. In this review, we detail the two major actions of specialized DNA polymerases that limit DNA damage transmission: the prevention of replicative stress by non-B DNA replication and the recovery of stalled replication forks.

Whole genome sequencing of Chinese clearhead icefish, Protosalanx hyalocranius

Background: Chinese clearhead icefish, Protosalanx hyalocranius, is a representative icefish species with economic importance and special appearance. Due to its great economic value in China, the fish was introduced into Lake Dianchi and several other lakes from the Lake Taihu half a century ago. Similar to the Sinocyclocheilus cavefish, the clearhead icefish has certain cavefish-like traits, such as transparent body and nearly scaleless skin. Here, we provide the whole genome sequence of this surface-dwelling fish and generated a draft genome assembly, aiming at exploring molecular mechanisms for the biological interests. Findings: A total of 252.1 Gb of raw reads were sequenced. Subsequently, a novel draft genome assembly was generated, with the scaffold N50 reaching 1.163 Mb. The genome completeness was estimated to be 98.39 % by using the CEGMA evaluation. Finally, we annotated 19 884 protein-coding genes and observed that repeat sequences account for 24.43 % of the genome assembly. Conclusion: We report the first draft genome of the Chinese clearhead icefish. The genome assembly will provide a solid foundation for further molecular breeding and germplasm resource protection in Chinese clearhead icefish, as well as other icefishes. It is also a valuable genetic resource for revealing the molecular mechanisms for the cavefish-like characters.

Karyotypic evolution of the Medicago complex: sativa-caerulea-falcata inferred from comparative cytogenetic analysis

BACKGROUND

Polyploidy plays an important role in the adaptation and speciation of plants. The alteration of karyotype is a significant event during polyploidy formation. The Medicago sativa complex includes both diploid (2n = 2× = 16) and tetraploid (2n = 2× = 32) subspecies. The tetraploid M. ssp. sativa was regarded as having a simple autopolyploid origin from diploid ssp. caerulea, whereas the autopolyploid origin of tetraploid ssp. falcata from diploid form ssp. falcata is still in doubt. In this study, detailed comparative cytogenetic analysis between diploid to tetraploid species, as well as genomic affinity across different species in the M. sativa complex, were conducted based on comparative mapping of 11 repeated DNA sequences and two rDNA sequences by a fluorescence in situ hybridization (FISH) technique.

RESULTS

FISH patterns of the repeats in diploid subspecies caerulea were highly similar to those in tetraploid subspecies sativa. Distinctly different FISH patterns were first observed in diploid ssp. falcata, with only centromeric hybridizations using centromeric and multiple region repeats and a few subtelomeric hybridizations using subtelomeric repeats. Tetraploid subspecies falcata was unexpectedly found to possess a highly variable karyotype, which agreed with neither diploid ssp. falcata nor ssp. sativa. Reconstruction of chromosome-doubling process of diploid ssp. caerulea showed that chromosome changes have occurred during polyploidization process.

CONCLUSIONS

The comparative cytogenetic results provide reliable evidence that diploid subspecies caerulea is the direct progenitor of tetraploid subspecies sativa. And autotetraploid ssp. sativa has been suggested to undergo a partial diploidization by the progressive accumulation of chromosome structural rearrangements during evolution. However, the tetraploid subspecies falcata is far from a simple autopolyploid from diploid subspecies falcata although no obvious morphological change was observed between these two subspecies.

Combined sense-antisense Alu elements activate the EGFP reporter gene when stable transfection

Alu elements in the human genome are present in more than one million copies, accounting for 10% of the genome. However, the biological functions of most Alu repeats are unknown. In this present study, we detected the effects of Alu elements on EGFP gene expression using a plasmid system to find the roles of Alu elements in human genome. We inserted 5'-4TMI-Alus-CMV promoter-4TMI-Alus (or antisense Alus)-3' sequences into the pEGFP-C1 vector to construct expression vectors. We altered the copy number of Alus, the orientation of the Alus, and the presence of an enhancer (4TMI) in the inserted 5'-4TMI-Alus-CMV promoter-4TMI-Alus (or antisense Alus)-3' sequences. These expression vectors were stably transfected into HeLa cells, and EGFP reporter gene expression was determined. Our results showed that combined sense-antisense Alu elements activated the EGFP reporter gene in the presence of enhancers and stable transfection. The combined sense-antisense Alu vectors carrying four copies of Alus downstream of inserted CMV induced much stronger EGFP gene expression than two copies. Alus downstream of inserted CMV were replaced to AluJBs (having 76% homology with Alu) to construct expression vectors. We found that combined sense-antisense Alu (or antisense AluJB) vectors induced strong EGFP gene expression after stable transfection and heat shock. To further explore combined sense-antisense Alus activating EGFP gene expression, we constructed Tet-on system vectors, mini-C1-Alu-sense-sense and mini-C1-Alu-sense-antisense (EGFP gene was driven by mini-CMV). We found that combined sense-antisense Alus activated EGFP gene in the presence of reverse tetracycline repressor (rTetR) and doxycycline (Dox). Clone experiments showed that Mini-C1-Alu-sense-antisense vector had more positive cells than that of Mini-C1-Alu-sense-sense vector. The results in this paper proved that Alu repetitive sequences inhibited gene expression and combined sense-antisense Alus activated EGFP reporter gene when Alu transcribes, which suggests that Alus play roles in maintaining gene expression (silencing genes or activating genes) in human genome.

Genomic skimming for identification of medium/highly abundant transposable elements in Arundo donax and Arundo plinii

Transposable elements (TEs) are the most abundant genetic material for almost all eukaryotic genomes. Their effects on the host genomes range from an extensive size variation to the regulation of gene expression, altering gene function and creating new genes. Because of TEs pivotal contribute to the host genome structure and regulation, their identification and characterization provide a wealth of useful data for gaining an in-depth understanding of host genome functioning. The giant reed (Arundo donax) is a perennial rhizomatous C₃ grass, octadecaploid, with an estimated nuclear genome size of 2744 Mbp. It is a promising feedstock for second-generation biofuels and biomethane production. To identify and characterize the most repetitive TEs in the genomes of A. donax and its ancestral A. plinii species, we carried out low-coverage whole genome shotgun sequencing for both species. Using a de novo repeat identification approach, 33,041 and 28,237 non-redundant repetitive sequences were identified and characterized in A. donax and A. plinii genomes, representing 37.55 and 31.68% of each genome, respectively. Comparative phylogenetic analyses, including the major TE classes identified in A. donax and A. plinii, together with rice and maize TE paralogs, were carried out to understand the evolutionary relationship of the most abundant TE classes. Highly conserved copies of RIRE1-like Ty1-Copia elements were discovered in two Arundo spp. in which they represented nearly 3% of each genomic sequence. We identified and characterized the medium/highly repetitive TEs in two unexplored polyploid genomes, thus generating useful information for the study of the genomic structure, composition, and functioning of these two non-model species. We provided a valuable resource that could be exploited in any effort aimed at sequencing and assembling these two genomes.

Chromosome mapping of repetitive DNAs in sergeant major fishes (Abudefdufinae, Pomacentridae): a general view on the chromosomal conservatism of the genus

Species of the Abudefduf genus (sergeant-majors) are widely distributed in the Indian, Pacific and Atlantic oceans, with large schools inhabiting rocky coastal regions and coral reefs. This genus consists of twenty recognized species are of generalist habit, showing typical characteristics of colonizers. Some populations maintain gene flow between large oceanic areas, a condition that may influence their cytogenetic features. A number of species have been shown to be invaders and able to hybridize with local species. However, cytogenetic data in this genus are restricted to few species. In this way, the present study includes the chromosomal investigation, using conventional (Giemsa staining, Ag-NOR and C-banding) and molecular (in situ mapping of six different repetitive DNA classes) approaches in four Abudefduf species from different oceanic regions (A. bengalensis and A. sexfasciatus from the Indo-Pacific, A. vaigiensis from the Indian and A. saxatilis from the Atlantic oceans, respectively), to investigate the evolutionary events associated with the chromosomal diversification in this group. All species share a similar karyotype (2n = 48; NF = 52), except A. sexfasciatus (2n = 48; NF = 50), which possesses a characteristic pericentric inversion in the NOR-bearing chromosomal pair. Mapping of repetitive sequences suggests a chromosomal conservatism in this genus. The high karyotypic similarity between allopatric species of Abudefduf may be related to the success of natural viable hybrids among species with recent secondary contact.

The puzzling character of repetitive DNA in Phodopus genomes (Cricetidae, Rodentia)

Three novel repetitive DNA sequences are described, presenting a similar heterochromatic chromosomal location in two hamster species: Phodopus roborovskii and Phodopus sungorus (Cricetidae, Rodentia). Namely, two species-specific repetitive sequences (PROsat from P. roborovskii and PSUchr1sat from P. sungorus) surrounding a third one (PsatDNA), that is shared by both hamster genomes. Fiber-FISH analyses revealed that PROsat intermingles with PsatDNA in P. roborovskii and PSUchr1sat intermingles with PsatDNA in P. sungorus. A model explaining the evolution of this intricate chromosomal distribution is proposed, which can explain better the evolution of these very derivative genomes (in comparison to the ancestral Muroidea). The most plausible evolutionary scenario seems to be the expansion of a number of repeats into other's domain, most probably resulting in its intermingling, followed by the subsequent spread of these complex repeats from a single chromosomal location to other chromosomes. Evidences of an association between repetitive sequences and the chromosome evolution process were observed, namely for PROsat. Most probably, the evolutionary breakpoints that shaped PRO and PSU chromosomes (pericentric inversions and fusions) occurred within the boundaries of PROsat blocks in the ancestor. The repeats high diversity at the heterochromatic regions of Phodopus chromosomes, together with its complex organization, suggests that these species are important models for evolutionary studies, namely in the investigation of a possible relationship between repetitive sequences and the occurrence of chromosomal rearrangements and consequently, in genome evolution.

Sex chromosome evolution in snakes inferred from divergence patterns of two gametologous genes and chromosome distribution of sex chromosome-linked repetitive sequences

BACKGROUND

The discovery of differentially organized sex chromosome systems suggests that heteromorphic sex chromosomes evolved from a pair of homologous chromosomes. Whereas karyotypes are highly conserved in alethinophidian snakes, the degeneration status of the W chromosomes varies among species. The Z and W chromosomes are morphologically homomorphic in henophidian species, whereas in snakes belonging to caenophidian families the W chromosomes are highly degenerated. Snakes therefore are excellent animal models in which to study sex chromosome evolution. Herein, we investigated the differentiation processes for snake sex chromosomes using both coding and repetitive sequences. We analyzed phylogenetic relationships of CTNNB1 and WAC genes, localized to the centromeric and telomeric regions, respectively, of the long arms on snake sex chromosomes, and chromosome distribution of sex chromosome-linked repetitive sequences in several henophidian and caenophidian species.

RESULTS

Partial or full-length coding sequences of CTNNB1 and WAC were identified for Z homologs of henophidian species from Tropidophiidae, Boidae, Cylindrophiidae, Xenopeltidae, and Pythonidae, and for Z and W homologs of caenophidian species from Acrochordidae, Viperidae, Elapidae, and Colubridae. Female-specific sequences for the two genes were not found in the henophidian (boid and pythonid) species examined. Phylogenetic trees constructed using each gene showed that the Z and W homologs of the caenophidian species cluster separately. The repetitive sequence isolated from the W chromosome heterochromatin of the colubrid Elaphe quadrivirgata and a microsatellite motif (AGAT)8 were strongly hybridized with W chromosomes of the viperid and colubrid species examined.

CONCLUSION

Our phylogenetic analyses suggest that the cessation of recombination between the Z and W homologs of CTNNB1 and WAC predated the diversification of the caenophidian families. As the repetitive sequences on the W chromosomes were shared among viperid and colubrid species, heterochromatinization of the proto-W chromosome appears to have occurred before the splitting of these two groups. These results collectively suggest that differentiation of the proto-Z and proto-W chromosomes extended to wide regions on the sex chromosomes in the common ancestor of caenophidian families during a relatively short period.

Distinguishing highly similar gene isoforms with a clustering-based bioinformatics analysis of PacBio single-molecule long reads

Background

Gene isoforms are commonly found in both prokaryotes and eukaryotes. Since each isoform may perform a specific function in response to changing environmental conditions, studying the dynamics of gene isoforms is important in understanding biological processes and disease conditions. However, genome-wide identification of gene isoforms is technically challenging due to the high degree of sequence identity among isoforms. Traditional targeted sequencing approach, involving Sanger sequencing of plasmid-cloned PCR products, has low throughput and is very tedious and time-consuming. Next-generation sequencing technologies such as Illumina and 454 achieve high throughput but their short read lengths are a critical barrier to accurate assembly of highly similar gene isoforms, and may result in ambiguities and false joining during sequence assembly. More recently, the third generation sequencer represented by the PacBio platform offers sufficient throughput and long reads covering the full length of typical genes, thus providing a potential to reliably profile gene isoforms. However, the PacBio long reads are error-prone and cannot be effectively analyzed by traditional assembly programs.

Results

We present a clustering-based analysis pipeline integrated with PacBio sequencing data for profiling highly similar gene isoforms. This approach was first evaluated in comparison to de novo assembly of 454 reads using a benchmark admixture containing 10 known, cloned msg genes encoding the major surface glycoprotein of Pneumocystis jirovecii. All 10 msg isoforms were successfully reconstructed with the expected length (~1.5 kb) and correct sequence by the new approach, while 454 reads could not be correctly assembled using various assembly programs. When using an additional benchmark admixture containing 22 known P. jirovecii msg isoforms, this approach accurately reconstructed all but 4 these isoforms in their full-length (~3 kb); these 4 isoforms were present in low concentrations in the admixture. Finally, when applied to the original clinical sample from which the 22 known msg isoforms were cloned, this approach successfully identified not only all known isoforms accurately (~3 kb each) but also 48 novel isoforms.

Conclusions

PacBio sequencing integrated with the clustering-based analysis pipeline achieves high-throughput and high-resolution discrimination of highly similar sequences, and can serve as a new approach for genome-wide characterization of gene isoforms and other highly repetitive sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s13040-016-0090-8) contains supplementary material, which is available to authorized users.

A Gossypium BAC clone contains key repeat components distinguishing sub-genome of allotetraploidy cottons

Background

Dissecting genome organization is indispensable for further functional and applied studies. As genome sequences data shown, cotton genomes contain more than 60 % repetitive sequences, so study on repetitive sequences composition, structure, and distribution is the key step to dissect cotton genome.

Results

In this study, a bacterial artificial chromosome (BAC) clone enriched in repetitive sequences, was discovered initiatively by fluorescence in situ hybridization (FISH). FISHing with allotetraploidy cotton as target DNA, dispersed signals on most regions of all A sub-genome chromosomes, and only middle regions of all D sub-genome chromosomes were detected. Further FISHing with other cotton species bearing A or D genome as target DNA, specific signals were viewed. After BAC sequencing and bioinformational analysis, 129 repeat elements, size about 57,172 bp were found, accounting for more than 62 % of the BAC sequence (91,238 bp). Among them, a type of long terminal repeat-retrotransposon (LTR-RT), LTR/Gypsy was the key element causing the specific FISH results. Using the fragments of BAC matching with the identified Gypsy-like LTR as probes, the BAC-57I23-like FISH signals were reappeared. Running BLASTN, the fragments had good match with all chromosomes of G. arboreum (A₂) genome and A sub-genome of G. hirsutum (AD₁), and had relatively inferior match with all chromosomes of D sub-genome of AD₁, but had little match with the chromosomes of G. raimondii (D₅) genome, which was consistent with the FISH results.

Conclusion

A repeats-enriched cytogenetic marker to identify A and D sub-genomes of Gossypium was discovered by FISH. Combined sequences analysis with FISH verification, the assembly quality of repetitive sequences in the allotetraploidy cotton draft genome was assessed, and better chromosome belonging was verified. We also found the genomic distribution of the identified Gypsy-LTR-RT was similar to the distribution of heterochromatin. The expansion of this type of Gypsy-LTR-RT in heterochromatic regions may be one of the major reasons for the size gap between A and D genome. The findings showed here will help to understand the composition, structure, and evolution of cotton genome, and contribute to the further perfection of the draft genomes of cotton.

Visualization of Genomic Loci in Living Cells with a Fluorescent CRISPR/Cas9 System

The discovery that the RNA guided bacterial endonuclease Cas9 can be harnessed to target and manipulate user-defined genomic sequences has greatly influenced the field of genome engineering. Interestingly, a catalytically dead Cas9 (dCas9) can be employed as a targeted DNA-binding platform to alter gene expression. By fusing this dCas9 to eGFP, we and others could show that the CRISPR/Cas9 system can be further expanded to label and trace genomic loci in living cells. We demonstrated that by exchanging the sgRNA, dCas9-eGFP could be specifically directed to various heterochromatic sequences within the nucleus. Here, we provide a basic protocol for this versatile tool and describe how to verify new dCas9-eGFP targets.

Repetitive sequences in plant nuclear DNA: types, distribution, evolution and function

Repetitive DNA sequences are a major component of eukaryotic genomes and may account for up to 90% of the genome size. They can be divided into minisatellite, microsatellite and satellite sequences. Satellite DNA sequences are considered to be a fast-evolving component of eukaryotic genomes, comprising tandemly-arrayed, highly-repetitive and highly-conserved monomer sequences. The monomer unit of satellite DNA is 150-400 base pairs (bp) in length. Repetitive sequences may be species- or genus-specific, and may be centromeric or subtelomeric in nature. They exhibit cohesive and concerted evolution caused by molecular drive, leading to high sequence homogeneity. Repetitive sequences accumulate variations in sequence and copy number during evolution, hence they are important tools for taxonomic and phylogenetic studies, and are known as "tuning knobs" in the evolution. Therefore, knowledge of repetitive sequences assists our understanding of the organization, evolution and behavior of eukaryotic genomes. Repetitive sequences have cytoplasmic, cellular and developmental effects and play a role in chromosomal recombination. In the post-genomics era, with the introduction of next-generation sequencing technology, it is possible to evaluate complex genomes for analyzing repetitive sequences and deciphering the yet unknown functional potential of repetitive sequences.

A genetic program theory of aging using an RNA population model

Aging is a common characteristic of multicellular eukaryotes. Copious hypotheses have been proposed to explain the mechanisms of aging, but no single theory is generally acceptable. In this article, we refine the RNA population gene activating model (Lv et al., 2003) based on existing reports as well as on our own latest findings. We propose the RNA population model as a genetic theory of aging. The new model can also be applied to differentiation and tumorigenesis and could explain the biological significance of non-coding DNA, RNA, and repetitive sequence DNA. We provide evidence from the literature as well as from our own findings for the roles of repetitive sequences in gene activation. In addition, we predict several phenomena related to aging and differentiation based on this model.

Molecular characterization of Bacillus thuringiensis using rep-PCR

The genetic divergence of 65 strains of Bacillus thuringiensis (Bt) was determined using Rep-PCR. Based on the repetitive sequences the BOX primer was the most informative with 26 fragments, followed by ERIC (19) and REP (10), generating a total of 55 fragments. The dendogram shows that ten groups were formed when 45% was the average distance of the population: group 1 with 41,5% of the isolates, 33,8% of the isolates were distributed in other groups and 24,6% did not formed distinct group. 53,2% of the isolates from Embrapa are in the group 1, and 29,8% of the isolates are distributed in other groups. Bt strains from USDA and Institute Pasteur showed more variability.

RepEx: repeat extractor for biological sequences

Genomic sequences are far from being random but are made up of systematically ordered and information rich patterns. These repeated sequence patterns have been vastly utilized for their fundamental importance in understanding the genome function and organization. To this end, a comprehensive toolkit, RepEx, has been developed which extracts repeat (inverted, everted and mirror) patterns from the given genome sequence(s) without any constraints. The toolkit can also be used to fetch the inverted repeats present in the protein sequence(s). Further, it is capable of extracting exact and degenerate repeats with a user defined spacer intervals. It is remarkably more precise and sensitive when compared to the existing tools. An example with comprehensive case studies and a performance evaluation of the proposed toolkit has been presented to authenticate its efficiency and accuracy.