Disentangling river and swamp buffalo genetic diversity: initial insights from the 1000 Buffalo Genomes Project.

More people in the world depend on water buffalo for their livelihoods than on any other domesticated animals, but its genetics is still not extensively explored. The 1000 Buffalo Genomes Project (1000BGP) provides genetic resources for global buffalo population study and tools to breed more sustainable and productive buffaloes. Here we report the most contiguous swamp buffalo genome assembly (PCC_UOA_SB_1v2) with substantial resolution of telomeric and centromeric repeats, ∼4-fold more contiguous than the existing reference river buffalo assembly and exceeding a recently published male swamp buffalo genome. This assembly was used along with the current reference to align 140 water buffalo short-read sequences and produce a public genetic resource with an average of ∼41 million single nucleotide polymorphisms per swamp and river buffalo genome. Comparison of the swamp and river buffalo sequences showed ∼1.5% genetic differences, and estimated divergence time occurred 3.1 million years ago (95% CI, 2.6-4.9). The open science model employed in the 1000BGP provides a key genomic resource and tools for a species with global economic relevance.

Signals of positive selection in genomes of palearctic Myotis-bats coexisting with a fungal pathogen.

Disease can act as a driving force in shaping genetic makeup across populations, even species, if the impacts influence a particularly sensitive part of their life cycles. White-nose disease is caused by a fungal pathogen infecting bats during hibernation. The mycosis has caused massive population declines of susceptible species in North America, particularly in the genus Myotis. However, Myotis bats appear to tolerate infection in Eurasia, where the fungal pathogen has co-evolved with its bat hosts for an extended period of time. Therefore, with susceptible and tolerant populations, the fungal disease provides a unique opportunity to tease apart factors contributing to tolerance at a genomic level to and gain an understanding of the evolution of non-harmful in host-parasite interactions. To investigate if the fungal disease has caused adaptation on a genomic level in Eurasian bat species, we adopted both whole-genome sequencing approaches and a literature search to compile a set of 300 genes from which to investigate signals of positive selection in genomes of 11 Eurasian bats at the codon-level. Our results indicate significant positive selection in 38 genes, many of which have a marked role in responses to infection. Our findings suggest that white-nose syndrome may have applied a significant selective pressure on Eurasian Myotis-bats in the past, which can contribute their survival in co-existence with the pathogen. Our findings provide an insight on the selective pressure pathogens afflict on their hosts using methodology that can be adapted to other host-pathogen study systems.

The draft genomes of Crassostrea gasar and Crassostrea rhizophorae: key resources for leveraging oyster cultivation in the Southwest Atlantic.

OBJECTIVES: The two oyster species studied hold considerable economic importance for artisanal harvest (Crassostrea rhizophorae) and aquaculture (Crassostrea gasar). Their draft genomes will play an important role in the application of genomic methods such as RNAseq, population-based genomic scans aiming at addressing expression responses to pollution stress, adaptation to salinity and temperature variation, and will also permit investigating the genetic bases and enable marker-assisted selection of economically important traits like shell and mantle coloration and resistance to temperature and disease. DATA DESCRIPTION: The draft assembly size of Crassostrea gasar is 506 Mbp, and of Crassostrea rhizophorae is 584 Mbp with scaffolds N50 of 11,3 Mbp and 4,9 Mbp, respectively. The general masked bases by RepeatMasker in both genomes were highly similar using different datasets. The masked bases varied from 9.41% in C. gasar to 10.05% in C. rhizophorae and 42.85% in C. gasar to 44.44% in C. rhizophorae using Dfam and RepeatModeler datasets, respectively. Functional annotation with eggNog resulted in 34,693 annotated proteins in C. rhizophorae and 26,328 in C. gasar. BUSCO analysis shows that almost 99% of genes (5,295) are complete in relation to the mollusk orthologous genes dataset (mollusca_odb10).

Genomes were 'scrambled' when worms left the sea.

Chromosomal chaos may have aided their moves to fresh water and land.

XGRm: A Web Server for Interpreting Mouse Summary-level Genomic Data.

We introduce XGR-model (or XGRm), a web server made accessible at http://www.xgrm.pro, with the aim of meeting the increasing demand for effectively interpreting summary-level genomic data in model organisms. Currently, it hosts two enrichment analysers and two subnetwork analysers to support enrichment and subnetwork analyses for user-input mouse genomic data, whether gene-centric or genomic region-centric. The enrichment analysers identify ontology term enrichments for input genes (GElyser) or for genes linked from input genomic regions (RElyser). The subnetwork analysers rely on our previously established network algorithm to identify gene subnetworks from input gene-centric summary data (GSlyser) or from input region-centric summary data (RSlyser), leveraging network information about either functional interactions or pathway-derived interactions. Collectively, XGRm offers an all-in-one solution for gaining systems biology insights into summary-level genomic data in mice, underpinned by our commitment to regular updates as well as natural extensions to other model organisms.

Chromosome-level genome and population genomics of the intermediate horseshoe bat ( Rhinolophus affinis) reveal the molecular basis of virus tolerance in Rhinolophus and echolocation call frequency variation.

Horseshoe bats (genus Rhinolophus, family Rhinolophidae) represent an important group within chiropteran phylogeny due to their distinctive traits, including constant high-frequency echolocation, rapid karyotype evolution, and unique immune system. Advances in evolutionary biology, supported by high-quality reference genomes and comprehensive whole-genome data, have significantly enhanced our understanding of species origins, speciation mechanisms, adaptive evolutionary processes, and phenotypic diversity. However, genomic research and understanding of the evolutionary patterns of Rhinolophus are severely constrained by limited data, with only a single published genome of R. ferrumequinum currently available. In this study, we constructed a high-quality chromosome-level reference genome for the intermediate horseshoe bat ( R. affinis). Comparative genomic analyses revealed potential genetic characteristics associated with virus tolerance in Rhinolophidae. Notably, we observed expansions in several immune-related gene families and identified various genes functionally associated with the SARS-CoV-2 signaling pathway, DNA repair, and apoptosis, which displayed signs of rapid evolution. In addition, we observed an expansion of the major histocompatibility complex class II (MHC-II) region and a higher copy number of the HLA- DQB2 gene in horseshoe bats compared to other chiropteran species. Based on whole-genome resequencing and population genomic analyses, we identified multiple candidate loci (e.g., GLI3) associated with variations in echolocation call frequency across R. affinis subspecies. This research not only expands our understanding of the genetic characteristics of the Rhinolophus genus but also establishes a valuable foundation for future research.

Genome-wide identification, characterization, and expression analysis of the transient receptor potential gene family in mandarin fish Siniperca chuatsi.

BACKGROUND: Temperature is a crucial environmental determinant for the vitality and development of teleost fish, yet the underlying mechanisms by which they sense temperature fluctuations remain largely unexplored. Transient receptor potential (TRP) proteins, renowned for their involvement in temperature sensing, have not been characterized in teleost fish, especially regarding their temperature-sensing capabilities. RESULTS: In this study, a genome-wide analysis was conducted, identifying a total of 28 TRP genes in the mandarin fish Siniperca chuatsi. These genes were categorized into the families of TRPA, TRPC, TRPP, TRPM, TRPML, and TRPV. Despite notable variations in conserved motifs across different subfamilies, TRP family members shared common structural features, including ankyrin repeats and the TRP domain. Tissue expression analysis showed that each of these TRP genes exhibited a unique expression pattern. Furthermore, examination of the tissue expression patterns of ten selected TRP genes following exposure to both high and low temperature stress indicated the expression of TRP genes were responsive to temperatures changes. Moreover, the expression profiles of TRP genes in response to mandarin fish virus infections showed significant upregulation for most genes after Siniperca chuatsi rhabdovirus, mandarin fish iridovirus and infectious spleen and kidney necrosis virus infection. CONCLUSIONS: This study characterized the TRP family genes in mandarin fish genome-wide, and explored their expression patterns in response to temperature stress and virus infections. Our work will enhance the overall understanding of fish TRP channels and their possible functions.

Genome binding properties of Zic transcription factors underlie their changing functions during neuronal maturation.

BACKGROUND: The Zic family of transcription factors (TFs) promote both proliferation and maturation of cerebellar granule neurons (CGNs), raising the question of how a single, constitutively expressed TF family can support distinct developmental processes. Here we use an integrative experimental and bioinformatic approach to discover the regulatory relationship between Zic TF binding and changing programs of gene transcription during postnatal CGN differentiation. RESULTS: We first established a bioinformatic pipeline to integrate Zic ChIP-seq data from the developing mouse cerebellum with other genomic datasets from the same tissue. In newborn CGNs, Zic TF binding predominates at active enhancers that are co-bound by developmentally regulated TFs including Atoh1, whereas in mature CGNs, Zic TF binding consolidates toward promoters where it co-localizes with activity-regulated TFs. We then performed CUT&RUN-seq in differentiating CGNs to define both the time course of developmental shifts in Zic TF binding and their relationship to gene expression. Mapping Zic TF binding sites to genes using chromatin looping, we identified the set of Zic target genes that have altered expression in RNA-seq from Zic1 or Zic2 knockdown CGNs. CONCLUSIONS: Our data show that Zic TFs are required for both induction and repression of distinct, developmentally regulated target genes through a mechanism that is largely independent of changes in Zic TF binding. We suggest that the differential collaboration of Zic TFs with other TF families underlies the shift in their biological functions across CGN development.

Annelid Comparative Genomics and the Evolution of Massive Lineage-Specific Genome Rearrangement in Bilaterians.

The organization of genomes into chromosomes is critical for processes such as genetic recombination, environmental adaptation, and speciation. All animals with bilateral symmetry inherited a genome structure from their last common ancestor that has been highly conserved in some taxa but seemingly unconstrained in others. However, the evolutionary forces driving these differences and the processes by which they emerge have remained largely uncharacterized. Here, we analyze genome organization across the phylum Annelida using 23 chromosome-level annelid genomes. We find that while many annelid lineages have maintained the conserved bilaterian genome structure, the Clitellata, a group containing leeches and earthworms, possesses completely scrambled genomes. We develop a rearrangement index to quantify the extent of genome structure evolution and show that, compared to the last common ancestor of bilaterians, leeches and earthworms have among the most highly rearranged genomes of any currently sampled species. We further show that bilaterian genomes can be classified into two distinct categories-high and low rearrangement-largely influenced by the presence or absence, respectively, of chromosome fission events. Our findings demonstrate that animal genome structure can be highly variable within a phylum and reveal that genome rearrangement can occur both in a gradual, stepwise fashion, or rapid, all-encompassing changes over short evolutionary timescales.

A chromosome-level genome assembly of the Korean minipig (Sus scrofa).

Recent advancements in sequencing and genome assembly technologies have led to rapid generation of high-quality genome assemblies for various species and breeds. Despite the importance as minipigs an animal model in biomedical research, the construction of high-quality genome assemblies of minipigs still lags behind other pig breeds. To address this problem, we constructed a high-quality chromosome-level genome assembly of the Korean minipig (KMP) utilizing multiple different types of sequencing reads and reference genomes. The KMP assembly included 19 chromosome-level sequences with a total length of 2.52 Gb and N50 of 137 Mb. Comparative analyses with the pig reference genome (Sscrofa11.1) demonstrated comparable contiguity and completeness of the KMP assembly. Additionally, genome annotation analyses identified 22,666 protein-coding genes and repetitive elements occupying 40.10% of the genome. The KMP assembly and genome annotation provide valuable resources that can contribute to various future research on minipig and other pig breeds.

Validated DNA isolation method ensuring successful long-read sequencing of cattle semen genome.

Obtaining high-quality DNA suitable for long-read sequencing can be difficult for many types of tissues and cells, and it is a key step in current genomic studies. The challenge is even greater when it comes to isolating genomic DNA from mammalian spermatozoa, as DNA is tightly packed into a cell with a robust membrane rich in disulfide bonds. Here we describe a method for isolating high molecular weight DNA from Bovine commercial semen straws. This protocol includes a cleaning step to remove diluents and preservatives used for the long-term storage of the semen, which may affect long read sequencing. It is based on a simple salting-out method and avoid the use of spin columns, strong mixing or intensive centrifugation, in order to limit DNA fragmentation. However, we have adapted this protocol to facilitate the disruption of cell membranes and disulfide bonds with strong chaotropic and reducing agents. The average size of the fragments produced was approximately 49 kb, ranging from 25 to 85 kb, according to the femto pulse profiles.This method was used to isolate DNA from semen straws, more than 80 of them were successfully sequenced using the Continuous Long-Read (CLR) sequencing mode on the PacBio SequelII platform to study genome diversity and notably to detect large structural variations within genomes.

Chromosome-level genome assembly and annotation of the Rhabdophis nuchalis (Hubei keelback).

Rhabdophis nuchalis, a snake widely distributed in China, possesses a unique trait: glands beneath the skin on its neck and back, known as nucho-dorsal glands. These features make it a valuable subject for studying genetic diversity and the evolution of complex traits. In this study, we obtained a high-quality chromosome-level reference genome of R. nuchalis using MGI short-read sequencing, PacBio Revio long-read sequencing, and Hi-C sequencing techniques. The final assembly comprised 1.92 Gb of the R. nuchalis genome, anchored to 20 chromosomes (including 9 macrochromosomes and 11 microchromosomes), with a contig N50 of 104.79 Mb, a scaffold N50 of 204.96 Mb, and a BUSCO completeness of 97.50%. Additionally, we annotated a total of 1.09 Gb of repetitive sequences (which constitute 56.51% of the entire genome) and identified 22,057 protein-coding genes. This high-quality reference genome of R. nuchalis furnishes essential genomic data for comprehending the genetic diversity and evolutionary history of the species, as well as for facilitating species conservation efforts and comparative genomics studies.

Maptcha: an efficient parallel workflow for hybrid genome scaffolding.

BACKGROUND: Genome assembly, which involves reconstructing a target genome, relies on scaffolding methods to organize and link partially assembled fragments. The rapid evolution of long read sequencing technologies toward more accurate long reads, coupled with the continued use of short read technologies, has created a unique need for hybrid assembly workflows. The construction of accurate genomic scaffolds in hybrid workflows is complicated due to scale, sequencing technology diversity (e.g., short vs. long reads, contigs or partial assemblies), and repetitive regions within a target genome. RESULTS: In this paper, we present a new parallel workflow for hybrid genome scaffolding that would allow combining pre-constructed partial assemblies with newly sequenced long reads toward an improved assembly. More specifically, the workflow, called Maptcha, is aimed at generating long scaffolds of a target genome, from two sets of input sequences-an already constructed partial assembly of contigs, and a set of newly sequenced long reads. Our scaffolding approach internally uses an alignment-free mapping step to build a ⟨ contig,contig ⟩ graph using long reads as linking information. Subsequently, this graph is used to generate scaffolds. We present and evaluate a graph-theoretic "wiring" heuristic to perform this scaffolding step. To enable efficient workload management in a parallel setting, we use a batching technique that partitions the scaffolding tasks so that the more expensive alignment-based assembly step at the end can be efficiently parallelized. This step also allows the use of any standalone assembler for generating the final scaffolds. CONCLUSIONS: Our experiments with Maptcha on a variety of input genomes, and comparison against two state-of-the-art hybrid scaffolders demonstrate that Maptcha is able to generate longer and more accurate scaffolds substantially faster. In almost all cases, the scaffolds produced by Maptcha are at least an order of magnitude longer (in some cases two orders) than the scaffolds produced by state-of-the-art tools. Maptcha runs significantly faster too, reducing time-to-solution from hours to minutes for most input cases. We also performed a coverage experiment by varying the sequencing coverage depth for long reads, which demonstrated the potential of Maptcha to generate significantly longer scaffolds in low coverage settings ( 1 × - 10 × ).

Genome-wide detection of runs of homozygosity and heterozygosity in Tunchang pigs.

Tunchang pigs, mainly distributed throughout Hainan Province of China, are well-known for their superior meat quality, crude feed tolerance, and adaptability to high temperatures and humidity. Runs of homozygosity (ROH) can provide valuable information about the inbreeding coefficient in individuals and selection signals that may reveal candidate genes associated with key functional traits. Runs of heterozygosity (ROHet) are commonly associated with balance selection, which can help us understand the adaptive evolutionary history of domestic animals. In this study, we investigated ROHs and ROHets in 88 Tunchang pigs. We also compared the estimates of inbreeding coefficients in individuals calculated based on four methods. In summary, we detected a total of 16 ROH islands in our study, and 100 genes were found within ROH regions. These genes were correlated with economically important traits such as reproduction (e.g., SERPIND1, HIRA), meat quality (e.g., PI4KA, TBX1), immunity (e.g., ESS2, RANBP1), adaption to heat stress (TXNRD2 and DGCR8), and crude food tolerance (TRPM6). Moreover, we discovered 18 ROHet islands harbouring genes associated with reproduction (e.g., ARHGEF12, BMPR2), immune system (e.g., BRD4, DNMT3B). These findings may help us design effective breeding and conservation strategies for this unique breed.

Whole-genome scan for selection signature associated with temperature adaptation in Iranian sheep breeds.

The present study aimed to identify the selection signature associated with temperature adaptation in Iranian sheep breeds raised in cold and hot environments. The Illumina HD ovine SNP600K BeadChip genomic arrays were utilized to analyze 114 animals from eight Iranian sheep breeds, namely Ghezel, Afshari, Shall, Sanjabi, Lori-Bakhtiari, Karakul, Kermani, and Balochi. All animals were classified into two groups: cold-weather breeds and hot-weather breeds, based on the environments to which they are adapted and the regions where they have been raised for many years. The unbiased FST (Theta) and hapFLK tests were used to identify the selection signatures. The results revealed five genomic regions on chromosomes 2, 10, 11, 13, and 14 using the FST test, and three genomic regions on chromosomes 10, 14, and 15 using the hapFLK test to be under selection in cold and hot groups. Further exploration of these genomic regions revealed that most of these regions overlapped with genes previously identified to affect cold and heat stress, nervous system function, cell division and gene expression, skin growth and development, embryo and skeletal development, adaptation to hypoxia conditions, and the immune system. These regions overlapped with QTLs that had previously been identified as being associated with various important economic traits, such as body weight, skin color, and horn characteristics. The gene ontology and gene network analyses revealed significant pathways and networks that distinguished Iranian cold and hot climates sheep breeds from each other. We identified positively selected genomic regions in Iranian sheep associated with pathways related to cell division, biological processes, cellular responses to calcium ions, metal ions and inorganic substances. This study represents the initial effort to identify selective sweeps linked to temperature adaptation in Iranian indigenous sheep breeds. It may provide valuable insights into the genomic regions involved in climate adaptation in sheep.

Chromosome-level genome assembly of the glass catfish ( Kryptopterus vitreolus) reveals molecular clues to its transparent phenotype.

Glass catfish ( Kryptopterus vitreolus) are notable in the aquarium trade for their highly transparent body pattern. This transparency is due to the loss of most reflective iridophores and light-absorbing melanophores in the main body, although certain black and silver pigments remain in the face and head. To date, however, the molecular mechanisms underlying this transparent phenotype remain largely unknown. To explore the genetic basis of this transparency, we constructed a chromosome-level haplotypic genome assembly for the glass catfish, encompassing 32 chromosomes and 23 344 protein-coding genes, using PacBio and Hi-C sequencing technologies and standard assembly and annotation pipelines. Analysis revealed a premature stop codon in the putative albinism-related tyrp1b gene, encoding tyrosinase-related protein 1, rendering it a nonfunctional pseudogene. Notably, a synteny comparison with over 30 other fish species identified the loss of the endothelin-3 ( edn3b) gene in the glass catfish genome. To investigate the role of edn3b, we generated edn3b -/- mutant zebrafish, which exhibited a remarkable reduction in black pigments in body surface stripes compared to wild-type zebrafish. These findings indicate that edn3b loss contributes to the transparent phenotype of the glass catfish. Our high-quality chromosome-scale genome assembly and identification of key genes provide important molecular insights into the transparent phenotype of glass catfish. These findings not only enhance our understanding of the molecular mechanisms underlying transparency in glass catfish, but also offer a valuable genetic resource for further research on pigmentation in various animal species.

Into the Blue: Exploring genetic mechanisms behind the evolution of baleen whales.

Marine ecosystems are ideal for studying evolutionary adaptations involved in lineage diversification due to few physical barriers and reduced opportunities for strict allopatry compared to terrestrial ecosystems. Cetaceans (whales, dolphins, and porpoises) are a diverse group of mammals that successfully adapted to various habitats within the aquatic environment around 50 million years ago. While the overall adaptive transition from terrestrial to fully aquatic species is relatively well understood, the radiation of modern whales is still unclear. Here high-quality genomes derived from previously published data were used to identify genomic regions that potentially underpinned the diversification of baleen whales (Balaenopteridae). A robust molecular phylogeny was reconstructed based on 10,159 single copy and complete genes for eight mysticetes, seven odontocetes and two cetacean outgroups. Analysis of positive selection across 3,150 genes revealed that balaenopterids have undergone numerous idiosyncratic and convergent genomic variations that may explain their diversification. Genes associated with aging, survival and homeostasis were enriched in all species. Additionally, positive selection on genes involved in the immune system were disclosed for the two largest species, blue and fin whales. Such genes can potentially be ascribed to their morphological evolution, allowing them to attain greater length and increased cell number. Further evidence is presented about gene regions that might have contributed to the extensive anatomical changes shown by cetaceans, including adaptation to distinct environments and diets. This study contributes to our understanding of the genomic basis of diversification in baleen whales and the molecular changes linked to their adaptive radiation, thereby enhancing our understanding of cetacean evolution.

The genome sequence of the endemic Mexican common mustached Bat, Pteronotus mexicanus. Miller, 1902 [Mormoopidae; Pteronotus].

We describe here the first characterization of the genome of the bat Pteronotus mexicanus, an endemic species of Mexico, as part of the Mexican Bat Genome Project which focuses on the characterization and assembly of the genomes of endemic bats in Mexico. The genome was assembled from a liver tissue sample of an adult male from Jalisco, Mexico provided by the Texas Tech University Museum tissue collection. The assembled genome size was 1.9 Gb. The assembly of the genome was fitted in a framework of 110,533 scaffolds and 1,659,535 contigs. The ecological importance of bats such as P. mexicanus, and their diverse ecological roles, underscores the value of having complete genomes in addressing information gaps and facing challenges regarding their function in ecosystems and their conservation.

CARE-ing for Indigenous nonhuman genomic data - rethinking our approach.

Earlier this year, an international group of scientists published a paper in ScienceAdvances on the draft genome of the little bush moa (Anomalopteryx didiformis), one of about nine species of extinct flightless birds endemic to Aotearoa (New Zealand). The genome was sequenced from the ancient DNA of a "poorly provenanced" fossil bone acquired by the Royal Ontario Museum. It held important clues about the moa's evolutionary history and aspects of its biology.