CenH3 distribution reveals extended centromeres in the model beetle Tribolium castaneum.

Centromeres are chromosomal domains essential for kinetochore assembly and correct chromosome segregation. Inconsistent in their underlying DNA sequences, centromeres are defined epigenetically by the presence of the centromere-specific histone H3 variant CenH3. Most of the analyzed eukaryotes have monocentric chromosomes in which CenH3 proteins deposit into a single, primary constriction visible at metaphase chromosomes. Contrary to monocentrics, evolutionary sporadic holocentric chromosomes lack a primary constriction and have kinetochore activity distributed along the entire chromosome length. In this work, we identified cCENH3 protein, the centromeric H3 histone of the coleopteran model beetle Tribolium castaneum. By ChIP-seq analysis we disclosed that cCENH3 chromatin assembles upon a repertoire of repetitive DNAs. cCENH3 in situ mapping revealed unusually elongated T. castaneum centromeres that comprise approximately 40% of the chromosome length. Being the longest insect regional centromeres evidenced so far, T. castaneum centromeres are characterized by metapolycentric structure composed of several individual cCENH3-containing domains. We suggest that the model beetle T. castaneum with its metapolycentromeres could represent an excellent model for further studies of non-canonical centromeres in insects.

Satellitome analyses in nematodes illuminate complex species history and show conserved features in satellite DNAs.

BACKGROUND: Satellite DNAs (satDNAs) are tandemly repeated non-coding DNA sequences that belong to the most abundant and the fastest evolving parts of the eukaryotic genome. A satellitome represents the collection of different satDNAs in a genome. Due to extreme diversity and methodological difficulties to characterize and compare satDNA collection in complex genomes, knowledge on their putative functional constraints and capacity to participate in genome evolution remains rather elusive. SatDNA transcripts have been detected in many species, however comparative studies of satDNA transcriptome between species are extremely rare. RESULTS: We conducted a genome-wide survey and comparative analyses of satellitomes among different closely related Meloidogyne spp. nematodes. The evolutionary trends of satDNAs suggest that each round of proposed polyploidization in the evolutionary history is concomitant with the addition of a new set of satDNAs in the satellitome of any particular Meloidogyne species. Successive incorporation of new sets of satDNAs in the genome along the process of polyploidization supports multiple hybridization events as the main factor responsible for the formation of these species. Through comparative analyses of 83 distinct satDNAs, we found a CENP-B box-like sequence motif conserved among 11 divergent satDNAs (similarity ranges from 36 to 74%). We also found satDNAs that harbor a splice leader (SL) sequence which, in spite of overall divergence, shows conservation across species in two putative functional regions, the 25-nt SL exon and the Sm binding site. Intra- and interspecific comparative expression analyses of the complete satDNA set in the analyzed Meloidogyne species revealed transcription profiles including a subset of 14 actively transcribed satDNAs. Among those, 9 show active transcription in every species where they are found in the genome and throughout developmental stages. CONCLUSIONS: Our results demonstrate the feasibility and power of comparative analysis of the non-coding repetitive genome for elucidation of the origin of species with a complex history. Although satDNAs generally evolve extremely quickly, the comparative analyses of 83 satDNAs detected in the analyzed Meloidogyne species revealed conserved sequence features in some satDNAs suggesting sequence evolution under selective pressure. SatDNAs that are actively transcribed in related genomes and throughout nematode development support the view that their expression is not stochastic.

The Centromere Histone Is Conserved and Associated with Tandem Repeats Sharing a Conserved 19-bp Box in the Holocentromere of Meloidogyne Nematodes.

Although centromeres have conserved function, centromere-specific histone H3 (CenH3) and centromeric DNA evolve rapidly. The centromere drive model explains this phenomenon as a consequence of the conflict between fast-evolving DNA and CenH3, suggesting asymmetry in female meiosis as a crucial factor. We characterized evolution of the CenH3 protein in three closely related, polyploid mitotic parthenogenetic species of the Meloidogyne incognita group, and in the distantly related meiotic parthenogen Meloidogyne hapla. We identified duplication of the CenH3 gene in a putative sexual ancestral Meloidogyne. We found that one CenH3 (αCenH3) remained conserved in all extant species, including in distant Meloidogyne hapla, whereas the other evolved rapidly and under positive selection into four different CenH3 variants. This pattern of CenH3 evolution in Meloidogyne species suggests the subspecialization of CenH3s in ancestral sexual species. Immunofluorescence performed on mitotic Meloidogyne incognita revealed a dominant role of αCenH3 on its centromere, whereas the other CenH3s have lost their function in mitosis. The observed αCenH3 chromosome distribution disclosed cluster-like centromeric organization. The ChIP-Seq analysis revealed that in M. incognita αCenH3-associated DNA dominantly comprises tandem repeats, composed of divergent monomers which share a completely conserved 19-bp long box. Conserved αCenH3-associated DNA is also confirmed in the related mitotic Meloidogyne incognita group species suggesting preservation of both centromere protein and DNA constituents. We hypothesize that the absence of centromere drive in mitosis might allow for CenH3 and its associated DNA to achieve an equilibrium in which they can persist for long periods of time.

Reference-Guided De Novo Genome Assembly of the Flour Beetle Tribolium freemani.

The flour beetle Tribolium freemani is a sibling species of the model organism and important pest Tribolium castaneum. The two species are so closely related that they can produce hybrid progeny, but the genetic basis of their differences has not been revealed. In this work, we sequenced the T. freemani genome by applying PacBio HiFi technology. Using the well-assembled T. castaneum genome as a reference, we assembled 262 Mb of the T. freemani genomic sequence and anchored it in 10 linkage groups corresponding to nine autosomes and sex chromosome X. The assembly showed 99.8% completeness of conserved insect genes, indicating a high-quality reference genome. Comparison with the T. castaneum assembly revealed that the main differences in genomic sequence between the two sibling species come from repetitive DNA, including interspersed and tandem repeats. In this work, we also provided the complete assembled mitochondrial genome of T. freemani. Although the genome assembly needs to be ameliorated in tandemly repeated regions, the first version of the T. freemani reference genome and the complete mitogenome presented here represent useful resources for comparative evolutionary studies of related species and for further basic and applied research on different biological aspects of economically important pests.

RhoA activation and nuclearization marks loss of chondrocyte phenotype in crosstalk with Wnt pathway.

De-differentiation comprises a major drawback for the use of autologous chondrocytes in cartilage repair. Here, we investigate the role of RhoA and canonical Wnt signaling in chondrocyte phenotype. Chondrocyte de-differentiation is accompanied by an upregulation and nuclear localization of RhoA. Effectors of canonical Wnt signaling including ß-catenin and YAP/TAZ are upregulated in de-differentiating chondrocytes in a Rho-dependent manner. Inhibition of Rho activation with C3 transferase inhibits nuclear localization of RhoA, induces expression of chondrogenic markers on 2D and enhances the chondrogenic effect of 3D culturing. Upregulation of chondrogenic markers by Rho inhibition is accompanied by loss of canonical Wnt signaling markers in 3D or on 2D whereas treatment of chondrocytes with Wnt-3a abrogates this effect. However, induction of canonical Wnt signaling inhibits chondrogenic markers on 2D but enhances chondrogenic re-differentiation on 2D with C3 transferase or in 3D. These data provide insights on the context-dependent role of RhoA and Wnt signaling in de-differentiation and on mechanisms to induce chondrogenic markers for therapeutic approaches.

The Genome Organization of 5S rRNA Genes in the Model Organism Tribolium castaneum and Its Sibling Species Tribolium freemani.

5S ribosomal DNAs (rDNAs) are arranged in tandem and are often under-represented in genome assemblies. In the present study, we performed a global and in-depth analysis of the 5S rDNAs in the model insect Tribolium castaneum and its closely related species Tribolium freemani. To accomplish this goal, we used our recently published genome assemblies based on Nanopore and PacBio long-read sequencing. Although these closely related species share the 5S rRNA gene sequence with high homology, they show a different organization of the 5S rDNA locus. Analysis of 5S rDNA arrays in T. castaneum revealed a typical tandemly repeated organization characterized by repeat units consisting of the 121 bp long 5S rRNA gene and the 71 bp long nontranscribed spacer (NTS). In contrast, T. freemani showed a much more complex organization of 5S rDNA arrays characterized by two patterns. The first is based on the association of 5S rRNA gene with arrays of a satellite DNA, representing the NTS sequence of the 5S rDNA genes in T. freemani. The second, more complex type is characterized by a somewhat less frequent occurrence of the 5S rRNA gene and its association with longer satellite DNA arrays that are regularly interrupted by Jockey-like retrotransposons. This organization, in which the ribosomal gene is associated with two completely different repetitive elements such as satellite DNAs and retrotransposons, suggests that the 5S rRNA gene, regardless of its crucial function in the genome, could be a subject of extremely dynamic genomic rearrangements.

Unzipped genome assemblies of polyploid root-knot nematodes reveal unusual and clade-specific telomeric repeats.

Using long-read sequencing, we assembled and unzipped the polyploid genomes of Meloidogyne incognita, M. javanica and M. arenaria, three of the most devastating plant-parasitic nematodes. We found the canonical nematode telomeric repeat to be missing in these and other Meloidogyne genomes. In addition, we find no evidence for the enzyme telomerase or for orthologs of C. elegans telomere-associated proteins, suggesting alternative lengthening of telomeres. Instead, analyzing our assembled genomes, we identify species-specific composite repeats enriched mostly at one extremity of contigs. These repeats are G-rich, oriented, and transcribed, similarly to canonical telomeric repeats. We confirm them as telomeric using fluorescent in situ hybridization. These repeats are mostly found at one single end of chromosomes in these species. The discovery of unusual and specific complex telomeric repeats opens a plethora of perspectives and highlights the evolutionary diversity of telomeres despite their central roles in senescence, aging, and chromosome integrity.

The Low-Copy-Number Satellite DNAs of the Model Beetle Tribolium castaneum.

The red flour beetle Tribolium castaneum is an important pest of stored agricultural products and the first beetle whose genome was sequenced. So far, one high-copy-number and ten moderate-copy-number satellite DNAs (satDNAs) have been described in the assembled part of its genome. In this work, we aimed to catalog the entire collection of T. castaneum satDNAs. We resequenced the genome using Illumina technology and predicted potential satDNAs via graph-based sequence clustering. In this way, we discovered 46 novel satDNAs that occupied a total of 2.1% of the genome and were, therefore, considered low-copy-number satellites. Their repeat units, preferentially 140-180 bp and 300-340 bp long, showed a high A + T composition ranging from 59.2 to 80.1%. In the current assembly, we annotated the majority of the low-copy-number satDNAs on one or a few chromosomes, discovering mainly transposable elements in their vicinity. The current assembly also revealed that many of the in silico predicted satDNAs were organized into short arrays not much longer than five consecutive repeats, and some of them also had numerous repeat units scattered throughout the genome. Although 20% of the unassembled genome sequence masked the genuine state, the predominance of scattered repeats for some low-copy satDNAs raises the question of whether these are essentially interspersed repeats that occur in tandem only sporadically, with the potential to be satDNA "seeds".

Isolation of High Molecular Weight DNA from the Model Beetle Tribolium for Nanopore Sequencing.

The long-read Nanopore sequencing has been recently applied for assembly of complex genomes and analysis of linear genome organization. The most critical factor for successful long-read sequencing is extraction of high molecular weight (HMW) DNA of sufficient purity and quantity. The challenges associated with input DNA quality are further amplified when working with extremely small insects with hard exoskeletons. Here, we optimized the isolation of HMW DNA from the model beetle Tribolium and tested for use in Nanopore sequencing. We succeeded in overcoming all the difficulties in HMW handling and library preparation that were encountered when using published protocols and commercial kits. Isolation of nuclei and subsequent purification of DNA on an anion-exchange chromatography column resulted in genomic HMW DNA that was efficiently relaxed, of optimal quality and in sufficient quantity for Nanopore MinION sequencing. DNA shearing increased average N50 read values up to 26 kb and allowed us to use a single flow cell in multiple library loads for a total output of more than 13 Gb. Although our focus was on T. castaneum and closely related species, we expect that this protocol, with appropriate modifications, could be extended to other insects, particularly beetles.

Sequence Composition Underlying Centromeric and Heterochromatic Genome Compartments of the Pacific Oyster Crassostrea gigas.

Segments of the genome enriched in repetitive sequences still present a challenge and are omitted in genome assemblies. For that reason, the exact composition of DNA sequences underlying the heterochromatic regions and the active centromeres are still unexplored for many organisms. The centromere is a crucial region of eukaryotic chromosomes responsible for the accurate segregation of genetic material. The typical landmark of centromere chromatin is the rapidly-evolving variant of the histone H3, CenH3, while DNA sequences packed in constitutive heterochromatin are associated with H3K9me3-modified histones. In the Pacific oyster Crassostrea gigas we identified its centromere histone variant, Cg-CenH3, that shows stage-specific distribution in gonadal cells. In order to investigate the DNA composition of genomic regions associated with the two specific chromatin types, we employed chromatin immunoprecipitation followed by high-throughput next-generation sequencing of the Cg-CenH3- and H3K9me3-associated sequences. CenH3-associated sequences were assigned to six groups of repetitive elements, while H3K9me3-associated-ones were assigned only to three. Those associated with CenH3 indicate the lack of uniformity in the chromosomal distribution of sequences building the centromeres, being also in the same time dispersed throughout the genome. The heterochromatin of C. gigas exhibited general paucity and limited chromosomal localization as predicted, with H3K9me3-associated sequences being predominantly constituted of DNA transposons.

Hypoxia regulates RhoA and Wnt/ß-catenin signaling in a context-dependent way to control re-differentiation of chondrocytes.

Cartilage tissue is avascular and hypoxic which regulates chondrocyte phenotype via stabilization of HIFs. Here, we investigated the role of hypoxia and HIFs in regulation of Rho and canonical Wnt signaling in chondrocytes. Our data demonstrates that hypoxia controls the expression of RhoA in chondrocytes in a context-dependent manner on the culturing conditions. Within a 3D microenvironment, hypoxia suppresses RhoA on which hypoxia-driven expression of chondrogenic markers depends. Conversely, hypoxia leads to upregulation of RhoA in chondrocytes on 2D with a failure in re-expression of chondrogenic markers. Similarly to RhoA, hypoxic regulation of Wnt/ß-catenin signaling depends on the microenvironment. Hypoxia downregulates ß-catenin within 3D hydrogels whereas it causes a potent increase on 2D. Hypoxia-induced suppression of canonical Wnt signaling in 3D contributes to the promotion of chondrogenic phenotype as induction of Wnt signaling abrogates the hypoxic re-differentiation of chondrocytes. Inhibiting Wnt/ß-catenin signaling via stabilization of Axin2 leads to a synergistic enhancement of hypoxia-induced expression of chondrogenic markers. The effects of hypoxia on Rho and Wnt/ß-catenin signaling are HIF-dependent as stabilizing HIFs under normoxia revealed similar effects on chondrocytes. The study reveals important insights on hypoxic signaling of chondrocytes and how hypoxia regulates cellular mechanisms depending on the cellular microenvironment.