Elasmobranch genome sequencing reveals evolutionary trends of vertebrate karyotype organization.

Genomic studies of vertebrate chromosome evolution have long been hindered by the scarcity of chromosome-scale DNA sequences of some key taxa. One of those limiting taxa has been the elasmobranchs (sharks and rays), which harbor species often with numerous chromosomes and enlarged genomes. Here, we report the chromosome-scale genome assembly for the zebra shark Stegostoma tigrinum, an endangered species that has a relatively small genome among sharks (3.71 Gb), as well as for the whale shark Rhincodon typus Our analysis using a male-female comparison identified an X Chromosome, the first genomically characterized shark sex chromosome. The X Chromosome harbors the Hox C cluster whose intact linkage has not been shown for an elasmobranch fish. The sequenced shark genomes show a gradualism of chromosome length with remarkable length-dependent characteristics-shorter chromosomes tend to have higher GC content, gene density, synonymous substitution rate, and simple tandem repeat content as well as smaller gene length and lower interspersed repeat content. We challenge the traditional binary classification of karyotypes as with and without so-called microchromosomes. Even without microchromosomes, the length-dependent characteristics persist widely in nonmammalian vertebrates. Our investigation of elasmobranch karyotypes underpins their unique characteristics and provides clues for understanding how vertebrate karyotypes accommodate intragenomic heterogeneity to realize a complex readout. It also paves the way to dissecting more genomes with variable sizes to be sequenced at high quality.

Diversity and Evolution of Highly Repetitive DNA Sequences Constituting Chromosome Site-Specific Heterochromatin in Two Gerbillinae Species.

Constitutive heterochromatin, consisting of repetitive sequences, diverges very rapidly; therefore, its nucleotide sequences and chromosomal distributions are often largely different, even between closely related species. The chromosome C-banding patterns of two Gerbillinae species, Meriones unguiculatus and Gerbillus perpallidus, vary greatly, even though they belong to the same subfamily. To understand the evolution of C-positive heterochromatin in these species, we isolated highly repetitive sequences, determined their nucleotide sequences, and characterized them using chromosomal and filter hybridization. We obtained a centromeric repeat (MUN-HaeIII) and a chromosome 13-specific repeat (MUN-EcoRI) from M. unguiculatus. We also isolated a centromeric/pericentromeric repeat (GPE-MBD) and an interspersed-type repeat that was predominantly amplified in the X and Y chromosomes (GPE-EcoRI) from G. perpallidus. GPE-MBD was found to contain a 17-bp motif that is essential for binding to the centromere-associated protein CENP-B. This indicates that it may play a role in the formation of a specified structure and/or function of centromeres. The nucleotide sequences of the three sequence families, except GPE-EcoRI, were conserved only in Gerbillinae. GPE-EcoRI was derived from the long interspersed nuclear elements 1 retrotransposon and showed sequence homology throughout Muridae and Cricetidae species, indicating that the repeat sequence occurred at least in the common ancestor of Muridae and Cricetidae. Due to a lack of assembly data of highly repetitive sequences constituting heterochromatin in whole-genome sequences of vertebrate species published to date, the knowledge obtained in this study provides useful information for a deep understanding of the evolution of repetitive sequences in not only rodents but also in mammals.

Molecular Cytogenetic Characterization of C-Band-Positive Heterochromatin of the Greater Long-Tailed Hamster (Tscherskia triton, Cricetinae).

The greater long-tailed hamster (Tscherskia triton, Cricetinae) has a unique karyotype (2n = 28), containing 11 pairs of acrocentric chromosomes with large C-band-positive centromeric heterochromatin blocks. To understand the origin and evolutionary process of heterochromatin in this species, we isolated novel families of chromosome site-specific highly repetitive DNA sequences from TaqI-digested genomic DNA and then characterized them by chromosome in situ and filter hybridization. The TaqI-families of repetitive sequences were classified into 2 types by their genome organization and chromosomal distribution: the 110-bp repeated sequence organized in large tandem arrays (as satellite DNA), localized to centromeric C-positive heterochromatin of acrocentric autosomes (chromosomes 1-11) and submetacentric X chromosome, and the 405-bp repeated sequence that was composed of 30-32-bp internal repeats, distributed in the pericentromeric region on the short arms of X and Y chromosomes. The repetitive sequences did not cross-hybridize with genomic DNA of any genera of Cricetinae (Mesocricetus, Cricetulus, and Phodopus). These results suggest that the 110-bp and 405-bp repeats rapidly diverged in the lineage of T. triton, evolving in a concerted manner among autosomes and X chromosome and within X and Y chromosomes, respectively. The 110-bp centromeric repeat contained a 17-bp motif in which 9 bases are essential for binding with the centromere-associated protein CENP-B, suggesting the possibility that the 110-bp major satellite DNA carrying the 17-bp motif may have a role in the formation of specified structure and/or function of centromeres in T. triton.

Genome evolution in the allotetraploid frog Xenopus laevis.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

Overexpression of microRNAs from the Gtl2-Rian locus contributes to postnatal death in mice.

The Dlk1-Dio3 imprinted domain functions in embryonic development but the roles of noncoding RNAs expressed from this domain remain unclear. We addressed this question by generating transgenic (TG) mice harbouring a BAC carrying IG-DMR (intergenic-differentially methylated region), Gtl2-DMR, Gtl2, Rtl1/Rtl1as, and part of Rian. High postnatal lethality (>85%) of the BAC-TG pups was observed in the maternally transmitted individuals (MAT-TG), but not following paternal transmission (PAT-TG). The DNA methylation status of IG-DMR and Gtl2-DMR in the BAC-allele was paternally imprinted similar to the genomic allele. The mRNA-Seq and miRNA-Seq analysis revealed marked expression changes in the MAT-TG, with 1,500 upregulated and 2,131 downregulated genes. The long noncoding RNAs and 12 miRNAs containing the BAC locus were markedly enhanced in the MAT-TG. We identified the 24 target genes of the overexpressed miRNAs and confirmed the downregulation in the MAT-TG. Notably, overexpression of mir770, mir493, and mir665 from Gtl2 in the MAT-TG embryos led to decreased expression of the 3 target genes, Col5a1, Pcgf2, and Clip2. Our results suggest that decreased expression of the 3 target genes concomitant with overexpression of the miRNAs within Gtl2 may be involved in the postnatal death in the MAT-TG. Because this imprinted domain is well conserved between mice and humans, the results of genetic and molecular analysis in mice hold important implications for related human disorders such as Temple syndrome.

Sex chromosome differentiation and the W- and Z-specific loci in Xenopus laevis.

Genetic sex-determining systems in vertebrates include two basic types of heterogamety; XX (female)/XY (male) and ZZ (male)/ZW (female) types. The African clawed frog Xenopus laevis has a ZZ/ZW-type sex-determining system. In this species, we previously identified a W-specific sex (female)-determining gene dmw, and specified W and Z chromosomes, which could be morphologically indistinguishable (homomorphic). In addition to dmw, we most recently discovered two genes, named scanw and ccdc69w, and one gene, named capn5z in the W- and Z-specific regions, respectively. In this study, we revealed the detail structures of the W/Z-specific loci and genes. Sequence analysis indicated that there is almost no sequence similarity between 278kb W-specific and 83kb Z-specific sequences on chromosome 2Lq32-33, where both the transposable elements are abundant. Synteny and phylogenic analyses indicated that all the W/Z-specific genes might have emerged independently. Expression analysis demonstrated that scanw and ccdc69w or capn5z are expressed in early differentiating ZW gonads or testes, thereby suggesting possible roles in female or male development, respectively. Importantly, the sex-determining gene (SDG) dmw might have been generated after allotetraploidization, thereby indicating the construction of the new sex-determining system by dmw after species hybridization. Furthermore, by direct genotyping, we confirmed that diploid WW embryos developed into normal female frogs, which indicate that the Z-specific region is not essential for female development. Overall, these findings indicate that sex chromosome differentiation has started, although no heteromorphic sex chromosomes are evident yet, in X. laevis. Homologous recombination suppression might have promoted the accumulation of mutations and transposable elements, and enlarged the W/Z-specific regions, thereby resulting in differentiation of the W/Z chromosomes.

Genome organization of the vg1 and nodal3 gene clusters in the allotetraploid frog Xenopus laevis.

Extracellular factors belonging to the TGF-ß family play pivotal roles in the formation and patterning of germ layers during early Xenopus embryogenesis. Here, we show that the vg1 and nodal3 genes of Xenopus laevis are present in gene clusters on chromosomes XLA1L and XLA3L, respectively, and that both gene clusters have been completely lost from the syntenic S chromosome regions. The presence of gene clusters and chromosome-specific gene loss were confirmed by cDNA FISH analyses. Sequence and expression analyses revealed that paralogous genes in the vg1 and nodal3 clusters on the L chromosomes were also altered compared to their Xenopus tropicalis orthologs. X. laevis vg1 and nodal3 paralogs have potentially become pseudogenes or sub-functionalized genes and are expressed at different levels. As X. tropicalis has a single vg1 gene on chromosome XTR1, the ancestral vg1 gene in X. laevis appears to have been expanded on XLA1L. Of note, two reported vg1 genes, vg1(S20) and vg1(P20), reside in the cluster on XLA1L. The nodal3 gene cluster is also present on X. tropicalis chromosome XTR3, but phylogenetic analysis indicates that nodal3 genes in X. laevis and X. tropicalis were independently expanded and/or evolved in concert within each cluster by gene conversion. These findings provide insights into the function and molecular evolution of TGF-ß family genes in response to allotetraploidization.

Molecular cloning and characterization of satellite DNA sequences from constitutive heterochromatin of the habu snake (Protobothrops flavoviridis, Viperidae) and the Burmese python (Python bivittatus, Pythonidae).

Highly repetitive DNA sequences of the centromeric heterochromatin provide valuable molecular cytogenetic markers for the investigation of genomic compartmentalization in the macrochromosomes and microchromosomes of sauropsids. Here, the relationship between centromeric heterochromatin and karyotype evolution was examined using cloned repetitive DNA sequences from two snake species, the habu snake (Protobothrops flavoviridis, Crotalinae, Viperidae) and Burmese python (Python bivittatus, Pythonidae). Three satellite DNA (stDNA) families were isolated from the heterochromatin of these snakes: 168-bp PFL-MspI from P. flavoviridis and 196-bp PBI-DdeI and 174-bp PBI-MspI from P. bivittatus. The PFL-MspI and PBI-DdeI sequences were localized to the centromeric regions of most chromosomes in the respective species, suggesting that the two sequences were the major components of the centromeric heterochromatin in these organisms. The PBI-MspI sequence was localized to the pericentromeric region of four chromosome pairs. The PFL-MspI and the PBI-DdeI sequences were conserved only in the genome of closely related species, Gloydius blomhoffii (Crotalinae) and Python molurus, respectively, although their locations on the chromosomes were slightly different. In contrast, the PBI-MspI sequence was also in the genomes of P. molurus and Boa constrictor (Boidae), and additionally localized to the centromeric regions of eight chromosome pairs in B. constrictor, suggesting that this sequence originated in the genome of a common ancestor of Pythonidae and Boidae, approximately 86 million years ago. The three stDNA sequences showed no genomic compartmentalization between the macrochromosomes and microchromosomes, suggesting that homogenization of the centromeric and/or pericentromeric stDNA sequences occurred in the macrochromosomes and microchromosomes of these snakes.

Molecular cloning and characterization of Siamese crocodile (Crocodylus siamensis) copper, zinc superoxide dismutase (CSI-Cu,Zn-SOD) gene.

Superoxide dismutase (SOD, EC 1.15.1.1) is an antioxidant enzyme found in all living cells. It regulates oxidative stress by breaking down superoxide radicals to oxygen and hydrogen peroxide. A gene coding for Cu,Zn-SOD was cloned and characterized from Siamese crocodile (Crocodylus siamensis; CSI). The full-length expressed sequence tag (EST) of this Cu,Zn-SOD gene (designated as CSI-Cu,Zn-SOD) contained 462bp encoding a protein of 154 amino acids without signal peptides, indicated as intracellular CSI-Cu,Zn-SOD. This agreed with the results from the phylogenetic tree, which indicated that CSI-Cu,Zn-SOD belonged to the intracellular Cu,Zn-SOD. Chromosomal location determined that the CSI-Cu,Zn-SOD was localized to the proximal region of the Siamese crocodile chromosome 1p. Several highly conserved motifs, two conserved signature sequences (GFHVHEFGDNT and GNAGGRLACGVI), and conserved amino acid residues for binding copper and zinc (His(47), His(49), His(64), His(72), His(81), Asp(84), and His(120)) were also identified in CSI-Cu,Zn-SOD. Real-time PCR analysis showed that CSI-Cu,Zn-SOD mRNA was expressed in all the tissues examined (liver, pancreas, lung, kidney, heart, and whole blood), which suggests a constitutively expressed gene in these tissues. Expression of the gene in Escherichia coli cells followed by purification yielded a recombinant CSI-Cu,Zn-SOD, with Km and Vmax values of 6.075mM xanthine and 1.4×10(-3)mmolmin(-1)mg(-1), respectively. This Vmax value was 40 times lower than native Cu,Zn-SOD (56×10(-3)mmolmin(-1)mg(-1)), extracted from crocodile erythrocytes. This suggests that cofactors, protein folding properties, or post-translational modifications were lost during the protein purification process, leading to a reduction in the rate of enzyme activity in bacterial expression of CSI-Cu,Zn-SOD.

Identification of the linkage group of the Z sex chromosomes of the sand lizard (Lacerta agilis, Lacertidae) and elucidation of karyotype evolution in lacertid lizards.

The sand lizard (Lacerta agilis, Lacertidae) has a chromosome number of 2n = 38, with 17 pairs of acrocentric chromosomes, one pair of microchromosomes, a large acrocentric Z chromosome, and a micro-W chromosome. To investigate the process of karyotype evolution in L. agilis, we performed chromosome banding and fluorescent in situ hybridization for gene mapping and constructed a cytogenetic map with 86 functional genes. Chromosome banding revealed that the Z chromosome is the fifth largest chromosome. The cytogenetic map revealed homology of the L. agilis Z chromosome with chicken chromosomes 6 and 9. Comparison of the L. agilis cytogenetic map with those of four Toxicofera species with many microchromosomes (Elaphe quadrivirgata, Varanus salvator macromaculatus, Leiolepis reevesii rubritaeniata, and Anolis carolinensis) showed highly conserved linkage homology of L. agilis chromosomes (LAG) 1, 2, 3, 4, 5(Z), 7, 8, 9, and 10 with macrochromosomes and/or macrochromosome segments of the four Toxicofera species. Most of the genes located on the microchromosomes of Toxicofera were localized to LAG6, small acrocentric chromosomes (LAG11-18), and a microchromosome (LAG19) in L. agilis. These results suggest that the L. agilis karyotype resulted from frequent fusions of microchromosomes, which occurred in the ancestral karyotype of Toxicofera and led to the disappearance of microchromosomes and the appearance of many small macrochromosomes.

Distribution of the T2-MITE Family Transposons in the Xenopus (Silurana) tropicalis Genome.

The T2 family of miniature inverted-repeat transposable elements (T2-MITE) is a prevalent MITE family found in both Xenopus(Silurana) tropicalis and X. laevis. Some subfamilies, particularly T2-A1 and T2-C, may have originated prior to the diversification of the 2 Xenopus lineages and currently include active members in X. tropicalis, whereas another subfamily, T2-E, may have lost its transposition activity even earlier. The distribution of each T2-MITE subfamily in X. tropicalis was investigated and compared to evaluate the evolutionary dynamics of the T2-MITE subfamilies. The subfamilies showed differences in chromosomal distribution, uniformity of insertion density on scaffolds, ratios of upstream to downstream insertions with respect to genes, and their distance from genes. Among these, the T2-C subfamily was interesting because it was frequently inserted upstream and close to genes and because genes with close insertions of this subfamily showed high correlations in spatial expression patterns. This unique distribution and long-lived transposition activity may reflect a mutual relationship evolved between this subfamily and the host.

Extraordinary Diversity in the Origins of Sex Chromosomes in Anurans Inferred from Comparative Gene Mapping.

Sex determination in frogs (anurans) is genetic and includes both male and female heterogamety. However, the origins of the sex chromosomes and their differentiation processes are poorly known. To investigate diversity in the origins of anuran sex chromosomes, we compared the chromosomal locations of sex-linked genes in 4 species: the African clawed frog (Xenopus laevis), the Western clawed frog (Silurana/X. tropicalis), the Japanese bell-ring frog (Buergeria buergeri), and the Japanese wrinkled frog (Rana rugosa). Comparative mapping data revealed that the sex chromosomes of X. laevis, X. tropicalis and R. rugosa are different chromosome pairs; however, the sex chromosomes of X. tropicalis and B. buergeri are homologous, although this may represent distinct evolutionary origins. We also examined the status of sex chromosomal differentiation in B. buergeri, which possesses heteromorphic ZW sex chromosomes, using comparative genomic hybridization and chromosome painting with DNA probes from the microdissected W chromosome. At least 3 rearrangement events have occurred in the proto-W chromosome: deletion of the nucleolus organizer region and a paracentric inversion followed by amplification of non-W-specific repetitive sequences.

A New Nomenclature of Xenopus laevis Chromosomes Based on the Phylogenetic Relationship to Silurana/Xenopus tropicalis.

Xenopus laevis (XLA) is an allotetraploid species which appears to have undergone whole-genome duplication after the interspecific hybridization of 2 diploid species closely related to Silurana/Xenopus tropicalis (XTR). Previous cDNA fluorescence in situ hybridization (FISH) experiments have identified 9 sets of homoeologous chromosomes in X. laevis, in which 8 sets correspond to chromosomes 1-8 of X. tropicalis (XTR1-XTR8), and the last set corresponds to a fusion of XTR9 and XTR10. In addition, recent X. laevis genome sequencing and BAC-FISH experiments support this physiological relationship and show no gross chromosome translocation in the X. laevis karyotype. Therefore, for the benefit of both comparative cytogenetics and genome research, we here propose a new chromosome nomenclature for X. laevis based on the phylogenetic relationship and chromosome length, i.e. XLA1L, XLA1S, XLA2L, XLA2S, and so on, in which the numbering of XLA chromosomes corresponds to that in X. tropicalis and the postfixes 'L' and 'S' stand for 'long' and 'short' chromosomes in the homoeologous pairs, which can be distinguished cytologically by their relative size. The last chromosome set is named XLA9L and XLA9S, in which XLA9 corresponds to both XTR9 and XTR10, and hence, to emphasize the phylogenetic relationship to X. tropicalis, XLA9_10L and XLA9_10S are also used as synonyms.

Chromosome size-correlated and chromosome size-uncorrelated homogenization of centromeric repetitive sequences in New World quails.

Many families of centromeric repetitive DNA sequences isolated from Struthioniformes, Galliformes, Falconiformes, and Passeriformes are localized primarily to microchromosomes. However, it is unclear whether chromosome size-correlated homogenization is a common characteristic of centromeric repetitive sequences in Aves. New World and Old World quails have the typical avian karyotype comprising chromosomes of two distinct sizes, and C-positive heterochromatin is distributed in centromeric regions of most autosomes and the whole W chromosome. We isolated six types of centromeric repetitive sequences from three New World quail species (Colinus virginianus, CVI; Callipepla californica, CCA; and Callipepla squamata, CSQ; Odontophoridae) and one Old World quail species (Alectoris chukar, ACH; Phasianidae), and characterized the sequences by nucleotide sequencing, chromosome in situ hybridization, and filter hybridization. The 385-bp CVI-MspI, 591-bp CCA-BamHI, 582-bp CSQ-BamHI, and 366-bp ACH-Sau3AI fragments exhibited tandem arrays of the monomer unit, and the 224-bp CVI-HaeIII and 135-bp CCA-HaeIII fragments were composed of minisatellite-like and microsatellite-like repeats, respectively. ACH-Sau3AI was a homolog of the chicken nuclear membrane repeat sequence, whose homologs are common in Phasianidae. CVI-MspI, CCA-BamHI, and CSQ-BamHI showed high homology and were specific to the Odontophoridae. CVI-MspI was localized to microchromosomes, whereas CVI-HaeIII, CCA-BamHI, and CSQ-BamHI were mapped to almost all chromosomes. CCA-HaeIII was localized to five pairs of macrochromosomes and most microchromosomes. ACH-Sau3AI was distributed in three pairs of macrochromosomes and all microchromosomes. Centromeric repetitive sequences may be homogenized in chromosome size-correlated and -uncorrelated manners in New World quails, although there may be a mechanism that causes homogenization of centromeric repetitive sequences primarily between microchromosomes, which is commonly observed in phasianid birds.

No Interstitial Telomeres on Autosomes but Remarkable Amplification of Telomeric Repeats on the W Sex Chromosome in the Sand Lizard (Lacerta agilis).

Telomeres are repeat (TTAGGG) n sequences that form terminal ends of chromosomes and have several functions, such as protecting the coding DNA from erosion at mitosis. Due to chromosomal rearrangements through evolutionary history (e.g., inversions and fusions), telomeric sequences are also found between the centromere and the terminal ends (i.e., at interstitial telomeric sites, ITSs). ITS telomere sequences have been implicated in heritable disease caused by genomic instability of ITS polymorphic variants, both with respect to copy number and sequence. In the sand lizard (Lacerta agilis), we have shown that telomere length is predictive of lifetime fitness in females but not males. To assess whether this sex specific fitness effect could be traced to ITSs differences, we mapped (TTAGGG) n sequences using fluorescence in situ hybridization in fibroblast cells cultured from 4 specimens of known sex. No ITSs could be found on autosomes in either sex. However, females have heterogametic sex chromosomes in sand lizards (ZW, 2n = 38) and the female W chromosome showed degeneration and remarkable (TTAGGG) n amplification, which was absent in the Z chromosomes. This work warrants further research on sex chromosome content, in particular of the degenerate W chromosome, and links to female fitness in sand lizards.

Genomic structures of the kW1 loci on the Z and W chromosomes in ratite birds: structural changes at an early stage of W chromosome differentiation.

The W chromosome of ratite birds shows minimal morphological differentiation and retains homology of genetic linkage and gene order with a substantial stretch of the Z chromosome; however, the molecular structure in the differentiated region is still not well known. The kW1 sequence was isolated from the kiwi as a W-specific DNA marker for PCR-based molecular sexing of ratite birds. In ratite W chromosomes, this sequence commonly contains a ∼200-bp deletion. To characterize the very early event of avian sex chromosome differentiation, we performed molecular cytogenetic analyses of kW1 and its flanking sequences in paleognathous and neognathous birds and reptiles. Female-specific repeats were found in the kW1-flanking sequence of the cassowary (Casuarius casuarius), and the repeats have been amplified in the pericentromeric region of the W chromosomes of ratites, which may have resulted from the cessation of meiotic recombination between the Z and W chromosomes at an early stage of sex chromosome differentiation. The presence of the kW1 sequence in neognathous birds and a crocodilian species suggests that the kW1 sequence was present in the ancestral genome of Archosauria; however, it disappeared in other reptilian taxa and several lineages of neognathous birds.

Karyotype evolution in monitor lizards: cross-species chromosome mapping of cDNA reveals highly conserved synteny and gene order in the Toxicofera clade.

The water monitor lizard (Varanus salvator macromaculatus (VSA), Platynota) has a chromosome number of 2n = 40: its karyotype consists of 16 macrochromosomes and 24 microchromosomes. To delineate the process of karyotype evolution in V. salvator macromaculatus, we constructed a cytogenetic map with 86 functional genes and compared it with those of the butterfly lizard (Leiolepis reevesii rubritaeniata (LRE); 2n = 36) and Japanese four-striped rat snake (Elaphe quadrivirgata (EQU); 2n = 36), members of the Toxicofera clade. The syntenies and gene orders of macrochromosomes were highly conserved between these species except for several chromosomal rearrangements: eight pairs of VSA macrochromosomes and/or chromosome arms exhibited homology with six pairs of LRE macrochromosomes and eight pairs of EQU macrochromosomes. Furthermore, the genes mapped to microchromosomes of three species were all located on chicken microchromosomes or chromosome 4p. No reciprocal translocations were found in the species, and their karyotypic differences were caused by: low frequencies of interchromosomal rearrangements, such as tandem fusions, or centric fissions/fusions between macrochromosomes and between macro- and microchromosomes; and intrachromosomal rearrangements, such as paracentric inversions or centromere repositioning. The chromosomal rearrangements that occurred in macrochromosomes of the Varanus lineage were also identified through comparative cytogenetic mapping of V. salvator macromaculatus and V. exanthematicus. Morphologic differences in chromosomes 6-8 between the two species could have resulted from pericentric inversion or centromere repositioning.

Unleashing diversity through flexibility: The evolutionary journey of sex chromosomes in amphibians and reptiles.

Sex determination systems have greatly diversified between amphibians and reptiles, with such as the different sex chromosome compositions within a single species and transition between temperature-dependent sex determination (TSD) and genetic sex determination (GSD). In most sex chromosome studies on amphibians and reptiles, the whole-genome sequence of Xenopous tropicalis and chicken have been used as references to compare the chromosome homology of sex chromosomes among each of these taxonomic groups, respectively. In the present study, we reviewed existing reports on sex chromosomes, including karyotypes, in amphibians and reptiles. Furthermore, we compared the identified genetic linkages of sex chromosomes in amphibians and reptiles with the chicken genome as a reference, which is believed to resemble the ancestral tetrapod karyotype. Our findings revealed that sex chromosomes in amphibians are derived from genetic linkages homologous to various chicken chromosomes, even among several frogs within single families, such as Ranidae and Pipidae. In contrast, sex chromosomes in reptiles exhibit conserved genetic linkages with chicken chromosomes, not only across most species within a single family, but also within closely related families. The diversity of sex chromosomes in amphibians and reptiles may be attributed to the flexibility of their sex determination systems, including the ease of sex reversal in these animals.

Highly species-specific centromeric repetitive DNA sequences in lizards: molecular cytogenetic characterization of a novel family of satellite DNA sequences isolated from the water monitor lizard (Varanus salvator macromaculatus, Platynota).

Two novel repetitive DNA sequences, VSAREP1 and VSAREP2, were isolated from the water monitor lizard (Varanus salvator macromaculatus, Platynota) and characterized using molecular cytogenetics. The respective lengths and guanine-cytosine (GC) contents of the sequences were 190 bp and 57.5% for VSAREP1 and 185 bp and 59.7% for VSAREP2, and both elements were tandemly arrayed as satellite DNA in the genome. VSAREP1 and VSAREP2 were each located at the C-positive heterochromatin in the pericentromeric region of chromosome 2q, the centromeric region of chromosome 5, and 3 pairs of microchromosomes. This suggests that genomic compartmentalization between macro- and microchromosomes might not have occurred in the centromeric repetitive sequences of V. salvator macromaculatus. These 2 sequences did only hybridize to genomic DNA of V. salvator macromaculatus, but no signal was observed even for other squamate reptiles, including Varanus exanthematicus, which is a closely related species of V. salvator macromaculatus. These results suggest that these sequences were differentiated rapidly or were specifically amplified in the V. salvator macromaculatus genome.

Squalomix: shark and ray genome analysis consortium and its data sharing platform.

The taxon Elasmobranchii (sharks and rays) contains one of the long-established evolutionary lineages of vertebrates with a tantalizing collection of species occupying critical aquatic habitats. To overcome the current limitation in molecular resources, we launched the Squalomix Consortium in 2020 to promote a genome-wide array of molecular approaches, specifically targeting shark and ray species. Among the various bottlenecks in working with elasmobranchs are their elusiveness and low fecundity as well as the large and highly repetitive genomes. Their peculiar body fluid composition has also hindered the establishment of methods to perform routine cell culturing required for their karyotyping. In the Squalomix consortium, these obstacles are expected to be solved through a combination of in-house cytological techniques including karyotyping of cultured cells, chromatin preparation for Hi-C data acquisition, and high fidelity long-read sequencing. The resources and products obtained in this consortium, including genome and transcriptome sequences, a genome browser powered by JBrowse2 to visualize sequence alignments, and comprehensive matrices of gene expression profiles for selected species are accessible through https://github.com/Squalomix/info.