Lipopolysaccharides of Herbaspirillum species and their relevance for bacterium-host interactions.

The lipopolysaccharides of Herbaspirillum lusitanum P6-12T (HlP6-12T) and H. frisingense GSF30T (HfGSF30T) was isolated by phenol-water extraction from bacterial cells and was characterized using chemical analysis and SDS-PAGE. It was shown that these bacteria produce LPSs that differ in their physicochemical properties and macromolecular organization. In this paper, the lipid A structure of the HlP6-12T LPS, was characterized through chemical analyses and matrix-assisted laser desorption ionization (MALDI) mass spectrometry. To prove the effect of the size of micelles on their bioavailability, we examined the activity of both LPSs toward the morphology of wheat seedlings. Analysis of the HlP6-12T and HfGSF30T genomes showed no significant differences between the operons that encode proteins involved in the biosynthesis of the lipids A and core oligosaccharides. The difference may be due to the composition of the O-antigen operon. HfGSF30T has two copies of the rfb operon, with the main one divided into two fragments. In contrast, the HlP6-12T genome contains only a single rfb-containing operon, and the other O-antigen operons are not comparable at all. The integrity of O-antigen-related genes may also affect LPS variability of. Specifically, we have observed a hairpin structure in the middle of the O-antigen glycosyltransferase gene, which led to the division of the gene into two fragments, resulting in incorrect protein synthesis and potential abnormalities in O-antigen production.

Chromosome-length genome assembly and karyotype of the endangered black-footed ferret (Mustela nigripes).

The black-footed ferret (Mustela nigripes) narrowly avoided extinction to become an oft-cited example of the benefits of intensive management, research, and collaboration to save a species through ex situ conservation breeding and reintroduction into its former range. However, the species remains at risk due to possible inbreeding, disease susceptibility, and multiple fertility challenges. Here, we report the de novo genome assembly of a male black-footed ferret generated through a combination of linked-read sequencing, optical mapping, and Hi-C proximity ligation. In addition, we report the karyotype for this species, which was used to anchor and assign chromosome numbers to the chromosome-length scaffolds. The draft assembly was ~2.5 Gb in length, with 95.6% of it anchored to 19 chromosome-length scaffolds, corresponding to the 2n = 38 chromosomes revealed by the karyotype. The assembly has contig and scaffold N50 values of 148.8 kbp and 145.4 Mbp, respectively, and is up to 96% complete based on BUSCO analyses. Annotation of the assembly, including evidence from RNA-seq data, identified 21,406 protein-coding genes and a repeat content of 37.35%. Phylogenomic analyses indicated that the black-footed ferret diverged from the European polecat/domestic ferret lineage 1.6 million yr ago. This assembly will enable research on the conservation genomics of black-footed ferrets and thereby aid in the further restoration of this endangered species.

Chromosome-Length Assembly of the Baikal Seal (Pusa sibirica) Genome Reveals a Historically Large Population Prior to Isolation in Lake Baikal.

Pusa sibirica, the Baikal seal, is the only extant, exclusively freshwater, pinniped species. The pending issue is, how and when they reached their current habitat-the rift lake Baikal, more than three thousand kilometers away from the Arctic Ocean. To explore the demographic history and genetic diversity of this species, we generated a de novo chromosome-length assembly, and compared it with three closely related marine pinniped species. Multiple whole genome alignment of the four species compared with their karyotypes showed high conservation of chromosomal features, except for three large inversions on chromosome VI. We found the mean heterozygosity of the studied Baikal seal individuals was relatively low (0.61 SNPs/kbp), but comparable to other analyzed pinniped samples. Demographic reconstruction of seals revealed differing trajectories, yet remarkable variations in Ne occurred during approximately the same time periods. The Baikal seal showed a significantly more severe decline relative to other species. This could be due to the difference in environmental conditions encountered by the earlier populations of Baikal seals, as ice sheets changed during glacial-interglacial cycles. We connect this period to the time of migration to Lake Baikal, which occurred ~3-0.3 Mya, after which the population stabilized, indicating balanced habitat conditions.

The Telomeric Repeats of HHV-6A Do Not Determine the Chromosome into Which the Virus Is Integrated.

Human herpes virus 6A (HHV-6A) is able to integrate into the telomeric and subtelomeric regions of human chromosomes representing chromosomally integrated HHV-6A (ciHHV-6A). The integration starts from the right direct repeat (DRR) region. It has been shown experimentally that perfect telomeric repeats (pTMR) in the DRR region are required for the integration, while the absence of the imperfect telomeric repeats (impTMR) only slightly reduces the frequency of HHV-6 integration cases. The aim of this study was to determine whether telomeric repeats within DRR may define the chromosome into which the HHV-6A integrates. We analysed 66 HHV-6A genomes obtained from public databases. Insertion and deletion patterns of DRR regions were examined. We also compared TMR within the herpes virus DRR and human chromosome sequences retrieved from the Telomere-to-Telomere consortium. Our results show that telomeric repeats in DRR in circulating and ciHHV-6A have an affinity for all human chromosomes studied and thus do not define a chromosome for integration.

The mixed liver and kidney transcriptome dataset of Darevskia valentini rock lizard.

OBJECTIVES: This study is performed in the frame of a bigger study dedicated to genomics and transcriptomics of parthenogenesis in vertebrates. Among vertebrates, obligate parthenogenesis was first described in the lizards of the genus Darevskia. In this genus, all found parthenogenetic species originated via interspecific hybridization. It remains unknown which genetic or genomic factors play a key role in the generation of parthenogenetic organisms. Comparative genomic and transcriptomic analysis of parthenogens and their parental species may elucidate this problem. Darevskia valentini is a paternal species for four (of seven) parthenogens of this genus, which we promote as a particularly important species for the generation of parthenogenetic forms. DATA DESCRIPTION: Total cellular RNA was isolated from kidney and liver tissues using the standard Trizol Tissue RNA Extraction protocol. Sequencing of transcriptome libraries prepared by random fragmentation of cDNA samples was performed on an Illumina HiSeq2500. Obtained raw sequences contained 117,6 million reads with the GC content of 47%. After preprocessing, raw data was assembled by Trinity and produced 491,482 contigs.

The first complete mitochondrial genome of Diadema antillarum (Diadematoida, Diadematidae) / El primer genoma mitocondrial completo de Diadema antillarum (Diadematoida, Diadematidae)

The mitochondrial genome of the long-spined black sea urchin, Diadema antillarum, was sequenced using Illumina next-generation sequencing technology. The complete mitogenome is 15,708 bp in length, containing two rRNA, 22 tRNA and 13 protein-coding genes, plus a noncoding control region of 133 bp. The nucleotide composition is 18.37% G, 23.79% C, 26.84% A and 30.99% T. The A + T bias is 57.84%. Phylogenetic analysis based on 12 complete mitochondrial genomes of sea urchins, including four species of the family Diadematidae, supported familial monophyly; however, the two Diadema species, D. antillarum and D. setosum were not recovered as sister taxa.

El genoma mitocondrial del erizo de mar negro de espinas largas, Diadema antillarum, se secuenció utilizando la tecnología de secuenciación de nueva generación de Illumina. El mitogenoma completo tiene un tamaño de 15,708 pb, que contiene dos ARNr, 22 ARNt y 13 genes codificadores de proteínas, además de una región de control no codificante de 133 pb. La composición de nucleótidos es 18.37% G, 23.79% C, 26.84% A y 30.99% T. El sesgo A+T es del 57.84%. El análisis filogenético basado en 12 genomas mitocondriales completos de erizos de mar, incluyendo cuatro especies de la familia Diadematidae, apoya la monofilia familiar. Sin embargo,  las dos especies de Diadema en este estudio,  D. antillarum y D. setosum no fueron identificadas como taxones hermanos.

First Genome of Rock Lizard Darevskia valentini Involved in Formation of Several Parthenogenetic Species.

The extant reptiles are one of the most diverse clades among terrestrial vertebrates and one of a few groups with instances of parthenogenesis. Due to the hybrid origin of parthenogenetic species, reference genomes of the parental species as well as of the parthenogenetic progeny are indispensable to explore the genetic foundations of parthenogenetic reproduction. Here, we report on the first genome assembly of rock lizard Darevskia valentini, a paternal species for several parthenogenetic lineages. The novel genome was used in the reconstruction of the comprehensive phylogeny of Squamata inferred independently from 7369 trees of single-copy orthologs and a supermatrix of 378 conserved proteins. We also investigated Hox clusters, the loci that are often regarded as playing an important role in the speciation of animal groups with drastically diverse morphology. We demonstrated that Hox clusters of D. valentini are invaded with transposons and contain the HoxC1 gene that has been considered to be lost in the amniote ancestor. This study provides confirmation for previous works and releases new genomic data that will contribute to future discoveries on the mechanisms of parthenogenesis as well as support comparative studies among reptiles.

Genome organization of major tandem repeats and their specificity for heterochromatin of macro- and microchromosomes in Japanese quail.

Tandemly repeated DNAs form heterochromatic regions of chromosomes, including the vital centromeric chromatin. Despite the progress in new genomic technologies, tandem repeats remain poorly deciphered and need targeted analysis in the species of interest. The Japanese quail is one of the highest-producing poultry species as well as a model organism. Its genome differs by a noticeable accumulation of heterochromatin, which led to an increase by 1/7 compared to the chicken genome size. Prominent heterochromatin blocks occupy the short arms of acrocentric macrochromosomes and of microchromosomes. We have applied de novo repeat finder approach to unassembled raw reads of the Japanese quail genome. We identified the 20 most common tandem repeats with the abundance >1 Mb, which represent about 4.8% of the genome. We found that tandem repeat CjapSAT primarily contributes to the centromeric regions of the macrochromosomes CJA1-8. Cjap31B together with previously characterized BglII makes up centromere regions of microchromosomes and W chromosome. Other repeats populate heterochromatin of microchromosomal short arms in unequal proportions, as revealed by fluorescence in situ hybridization. The Cjap84A, Cjap408A, and CjapSAT repeat sequences show similarities to retrotransposon motifs. This suggests that retroelements may have played a crucial role in the distribution of repeats throughout the Japanese quail genome.

Identification of Centromere-Specific Repeats in the Zebra Finch Genome.

Tandem repetitive sequences represent a significant part of many genomes but remain poorly characterized due to various methodological difficulties. Here, we describe the tandem repeat composition in the genome of zebra finch, Taeniopygia guttata, a species that has long served as an animal model, primarily in neurobiology and comparative genomics. Using available genome sequencing raw read datasets, we bioinformatically reconstructed consensus sequences of several tandem repeats and proved that the most abundant ones, Tgut191A and Tgut716A, are centromere-associated in chromosomes. Each centromeric region can have a different number of copies of each repeat, with Tgut716A enrichment in almost all microchromosomes and sex chromosomes. Sequences similar to Tgut191A and Tgut716A found in other Estrildidae and Viduidae species can be considered as candidate centromeric sequences, but this requires further cytogenetic verification.

The first complete mitochondrial genome of Diadema antillarum (Diadematoida, Diadematidae).

The mitochondrial genome of the long-spined black sea urchin, Diadema antillarum, was sequenced using Illumina next-generation sequencing technology. The complete mitogenome is 15,708 bp in length, containing two rRNA, 22 tRNA and 13 protein-coding genes, plus a noncoding control region of 133 bp. The nucleotide composition is 18.37% G, 23.79% C, 26.84% A and 30.99% T. The A + T bias is 57.84%. Phylogenetic analysis based on 12 complete mitochondrial genomes of sea urchins, including four species of the family Diadematidae, supported familial monophyly; however, the two Diadema species, D. antillarum and D. setosum were not recovered as sister taxa.

De novo transcriptome assembly and annotation of parthenogenetic lizard Darevskia unisexualis and its parental ancestors Darevskia valentini and Darevskia raddei nairensis.

Darevskia rock lizards include 29 sexual and seven parthenogenetic species of hybrid origin distributed in the Caucasus. All seven parthenogenetic species of the genus Darevskia were formed as a result of interspecific hybridization of only four sexual species. It remains unknown what are the main advantages of interspecific hybridization along with switching on parthenogenetic reproduction in evolution of reptiles. Data on whole transcriptome sequencing of parthenogens and their parental ancestors can provide value impact in solving this problem. Here we have sequenced ovary tissue transcriptomes from unisexual parthenogenetic lizard D. unisexualis and its parental bisexual ancestors to facilitate the subsequent annotation and to obtain the collinear characteristics for comparison with other lizard species. Here we report generated RNAseq data from total mRNA of ovary tissues of D. unisexualis, D. valentini and D. raddei with 58932755, 51634041 and 62788216 reads. Obtained RNA reads were assembled by Trinity assembler and 95141, 62123, 61836 contigs were identified with N50 values of 2409, 2801 and 2827 respectively. For further analysis top Gene Ontology terms were annotated for all species and transcript number was calculated. The raw data were deposited in the NCBI SRA database (BioProject PRJNA773939). The assemblies are available in Mendeley Data and can be accessed via doi:10.17632/rtd8cx7zc3.1.

Draft de novo Genome Assembly of the Elusive Jaguarundi, Puma yagouaroundi.

The Puma lineage within the family Felidae consists of 3 species that last shared a common ancestor around 4.9 million years ago. Whole-genome sequences of 2 species from the lineage were previously reported: the cheetah (Acinonyx jubatus) and the mountain lion (Puma concolor). The present report describes a whole-genome assembly of the remaining species, the jaguarundi (Puma yagouaroundi). We sequenced the genome of a male jaguarundi with 10X Genomics linked reads and assembled the whole-genome sequence. The assembled genome contains a series of scaffolds that reach the length of chromosome arms and is similar in scaffold contiguity to the genome assemblies of cheetah and puma, with a contig N50 = 100.2 kbp and a scaffold N50 = 49.27 Mbp. We assessed the assembled sequence of the jaguarundi genome using BUSCO, aligned reads of the sequenced individual and another published female jaguarundi to the assembled genome, annotated protein-coding genes, repeats, genomic variants and their effects with respect to the protein-coding genes, and analyzed differences of the 2 jaguarundis from the reference mitochondrial genome. The jaguarundi genome assembly and its annotation were compared in quality, variants, and features to the previously reported genome assemblies of puma and cheetah. Computational analyzes used in the study were implemented in transparent and reproducible way to allow their further reuse and modification.

Draft genome of Bugula neritina, a colonial animal packing powerful symbionts and potential medicines.

Many animal phyla have no representatives within the catalog of whole metazoan genome sequences. This dataset fills in one gap in the genome knowledge of animal phyla with a draft genome of Bugula neritina (phylum Bryozoa). Interest in this species spans ecology and biomedical sciences because B. neritina is the natural source of bioactive compounds called bryostatins. Here we present a draft assembly of the B. neritina genome obtained from PacBio and Illumina HiSeq data, as well as genes and proteins predicted de novo and verified using transcriptome data, along with the functional annotation. These sequences will permit a better understanding of host-symbiont interactions at the genomic level, and also contribute additional phylogenomic markers to evaluate Lophophorate or Lophotrochozoa phylogenetic relationships. The effort also fits well with plans to ultimately sequence all orders of the Metazoa.

Structural and genetic characterization of the colitose-containing O-specific polysaccharide from the lipopolysaccharide of Herbaspirillum frisingense GSF30T.

The lipopolysaccharide (LPS) of Herbaspirillum frisingense GSF30T (HfGSF30), a non-pathogenic diazotrophic endobiont, was isolated by phenol-water extraction from bacterial cells and was characterized by chemical analyses and SDS PAGE. The O-specific polysaccharide (OPS, O-antigen), obtained by mild acid hydrolysis of the LPS, was examined by sugar and methylation analysis, along with 1H and 13C NMR spectroscopy, including 2D 1H,1H COSY, 1H,1H TOCSY, 1H,1H ROESY, 1H,13C HSQC, and 1H,13C HMBC experiments. The OPS was found to consist of branched tetrasaccharide repeating units of the following structure: [Formula: see text] This structure is unique among the known bacterial polysaccharide structures. Analysis of the HfGSF30 genome showed that it contained a set of sequentially arranged operons (presumably a cluster of genes) associated with the O-antigen. Amino acid sequence analysis using the BLAST program demonstrated the specificity of this putative cluster for Herbaspirillum spp. The genes responsible for the biosynthesis of the OPS of HfGSF30 were dispersed in the genome, constituting small operons. A putative O-antigen gene cluster of HfGSF30 was identified and found to be consistent with the OPS structure.

Whole Genome Sequencing and Re-sequencing of the Sable Antelope (Hippotragus niger): A Resource for Monitoring Diversity in ex Situ and in Situ Populations.

Genome-wide assessment of genetic diversity has the potential to increase the ability to understand admixture, inbreeding, kinship and erosion of genetic diversity affecting both captive (ex situ) and wild (in situ) populations of threatened species. The sable antelope (Hippotragus niger), native to the savannah woodlands of sub-Saharan Africa, is a species that is being managed ex situ in both public (zoo) and private (ranch) collections in the United States. Our objective was to develop whole genome sequence resources that will serve as a foundation for characterizing the genetic status of ex situ populations of sable antelope relative to populations in the wild. Here we report the draft genome assembly of a male sable antelope, a member of the subfamily Hippotraginae (Bovidae, Cetartiodactyla, Mammalia). The 2.596 Gb draft genome consists of 136,528 contigs with an N50 of 45.5 Kbp and 16,927 scaffolds with an N50 of 4.59 Mbp. De novo annotation identified 18,828 protein-coding genes and repetitive sequences encompassing 46.97% of the genome. The discovery of single nucleotide variants (SNVs) was assisted by the re-sequencing of seven additional captive and wild individuals, representing two different subspecies, leading to the identification of 1,987,710 bi-allelic SNVs. Assembly of the mitochondrial genomes revealed that each individual was defined by a unique haplotype and these data were used to infer the mitochondrial gene tree relative to other hippotragine species. The sable antelope genome constitutes a valuable resource for assessing genome-wide diversity and evolutionary potential, thereby facilitating long-term conservation of this charismatic species.

White shark genome reveals ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability.

The white shark (Carcharodon carcharias; Chondrichthyes, Elasmobranchii) is one of the most publicly recognized marine animals. Here we report the genome sequence of the white shark and comparative evolutionary genomic analyses to the chondrichthyans, whale shark (Elasmobranchii) and elephant shark (Holocephali), as well as various vertebrates. The 4.63-Gbp white shark genome contains 24,520 predicted genes, and has a repeat content of 58.5%. We provide evidence for a history of positive selection and gene-content enrichments regarding important genome stability-related genes and functional categories, particularly so for the two elasmobranchs. We hypothesize that the molecular adaptive emphasis on genome stability in white and whale sharks may reflect the combined selective pressure of large genome sizes, high repeat content, high long-interspersed element retrotransposon representation, large body size, and long lifespans, represented across these two species. Molecular adaptation for wound healing was also evident, with positive selection in key genes involved in the wound-healing process, as well as Gene Ontology enrichments in fundamental wound-healing pathways. Sharks, particularly apex predators such as the white shark, are believed to have an acute sense of smell. However, we found very few olfactory receptor genes, very few trace amine-associated receptors, and extremely low numbers of G protein-coupled receptors. We did however, identify 13 copies of vomeronasal type 2 (V2R) genes in white shark and 10 in whale shark; this, combined with the over 30 V2Rs reported previously for elephant shark, suggests this gene family may underlie the keen odorant reception of chondrichthyans.

B Chromosomes of the Asian Seabass (Lates calcarifer) Contribute to Genome Variations at the Level of Individuals and Populations.

The Asian seabass (Lates calcarifer) is a bony fish from the Latidae family, which is widely distributed in the tropical Indo-West Pacific region. The karyotype of the Asian seabass contains 24 pairs of A chromosomes and a variable number of AT- and GC-rich B chromosomes (Bchrs or Bs). Dot-like shaped and nucleolus-associated AT-rich Bs were microdissected and sequenced earlier. Here we analyzed DNA fragments from Bs to determine their repeat and gene contents using the Asian seabass genome as a reference. Fragments of 75 genes, including an 18S rRNA gene, were found in the Bs; repeats represented 2% of the Bchr assembly. The 18S rDNA of the standard genome and Bs were similar and enriched with fragments of transposable elements. A higher nuclei DNA content in the male gonad and somatic tissue, compared to the female gonad, was demonstrated by flow cytometry. This variation in DNA content could be associated with the intra-individual variation in the number of Bs. A comparison between the copy number variation among the B-related fragments from whole genome resequencing data of Asian seabass individuals identified similar profiles between those from the South-East Asian/Philippines and Indian region but not the Australian ones. Our results suggest that Bs might cause variations in the genome among the individuals and populations of Asian seabass. A personalized copy number approach for segmental duplication detection offers a suitable tool for population-level analysis across specimens with low coverage genome sequencing.

Innovative assembly strategy contributes to understanding the evolution and conservation genetics of the endangered Solenodon paradoxus from the island of Hispaniola.

Solenodons are insectivores that live in Hispaniola and Cuba. They form an isolated branch in the tree of placental mammals that are highly divergent from other eulipothyplan insectivores The history, unique biology, and adaptations of these enigmatic venomous species could be illuminated by the availability of genome data. However, a whole genome assembly for solenodons has not been previously performed, partially due to the difficulty in obtaining samples from the field. Island isolation and reduced numbers have likely resulted in high homozygosity within the Hispaniolan solenodon (Solenodon paradoxus). Thus, we tested the performance of several assembly strategies on the genome of this genetically impoverished species. The string graph-based assembly strategy seemed a better choice compared to the conventional de Bruijn graph approach due to the high levels of homozygosity, which is often a hallmark of endemic or endangered species. A consensus reference genome was assembled from sequences of 5 individuals from the southern subspecies (S. p. woodi). In addition, we obtained an additional sequence from 1 sample of the northern subspecies (S. p. paradoxus). The resulting genome assemblies were compared to each other and annotated for genes, with an emphasis on venom genes, repeats, variable microsatellite loci, and other genomic variants. Phylogenetic positioning and selection signatures were inferred based on 4,416 single-copy orthologs from 10 other mammals. We estimated that solenodons diverged from other extant mammals 73.6 million years ago. Patterns of single-nucleotide polymorphism variation allowed us to infer population demography, which supported a subspecies split within the Hispaniolan solenodon at least 300 thousand years ago.

New high copy tandem repeat in the content of the chicken W chromosome.

The content of repetitive DNA in avian genomes is considerably less than in other investigated vertebrates. The first descriptions of tandem repeats were based on the results of routine biochemical and molecular biological experiments. Both satellite DNA and interspersed repetitive elements were annotated using library-based approach and de novo repeat identification in assembled genome. The development of deep-sequencing methods provides datasets of high quality without preassembly allowing one to annotate repetitive elements from unassembled part of genomes. In this work, we search the chicken assembly and annotate high copy number tandem repeats from unassembled short raw reads. Tandem repeat (GGAAA)n has been identified and found to be the second after telomeric repeat (TTAGGG)n most abundant in the chicken genome. Furthermore, (GGAAA)n repeat forms expanded arrays on the both arms of the chicken W chromosome. Our results highlight the complexity of repetitive sequences and update data about organization of sex W chromosome in chicken.

Genome-wide signatures of complex introgression and adaptive evolution in the big cats.

The great cats of the genus Panthera comprise a recent radiation whose evolutionary history is poorly understood. Their rapid diversification poses challenges to resolving their phylogeny while offering opportunities to investigate the historical dynamics of adaptive divergence. We report the sequence, de novo assembly, and annotation of the jaguar (Panthera onca) genome, a novel genome sequence for the leopard (Panthera pardus), and comparative analyses encompassing all living Panthera species. Demographic reconstructions indicated that all of these species have experienced variable episodes of population decline during the Pleistocene, ultimately leading to small effective sizes in present-day genomes. We observed pervasive genealogical discordance across Panthera genomes, caused by both incomplete lineage sorting and complex patterns of historical interspecific hybridization. We identified multiple signatures of species-specific positive selection, affecting genes involved in craniofacial and limb development, protein metabolism, hypoxia, reproduction, pigmentation, and sensory perception. There was remarkable concordance in pathways enriched in genomic segments implicated in interspecies introgression and in positive selection, suggesting that these processes were connected. We tested this hypothesis by developing exome capture probes targeting ~19,000 Panthera genes and applying them to 30 wild-caught jaguars. We found at least two genes (DOCK3 and COL4A5, both related to optic nerve development) bearing significant signatures of interspecies introgression and within-species positive selection. These findings indicate that post-speciation admixture has contributed genetic material that facilitated the adaptive evolution of big cat lineages.