Genomic Conservation Research Articles

Genomic Insights and Conservation Priorities for Kongshan Cattle: A Whole-Genome Resequencing Approach

Kongshan Cattle, indigenous to Sichuan Province and recognized as China’s 56th local cattle breed in 2024, exhibit unique adaptations including superior resistance to harsh conditions. Despite a declining population due to the influx of foreign breeds, there is a significant focus on preserving their genetic diversity through advanced genomic techniques. This study utilized whole-genome resequencing, a cost-effective and information-rich method, to perform a comprehensive genetic assessment of the Kongshan Cattle. High-quality resequencing data yielded an average of 17.5 billion clean bases per sample, with high proportions of Q20 and Q30 bases and a balanced GC content. SNP analysis revealed an average of 14 million SNPs per sample, with a notable transition-to-transversion ratio and a significant portion of heterozygosity. Further analysis of genomic and coding regions identified substantial insertions and deletions, particularly in coding sequences affecting gene functionality. A detailed examination of these genetic variations highlighted genes, including NEIL2 and PNKP, which are integral to stress resistance pathways, indicating potential adaptive advantages. This study not only underscores the genetic diversity of Kongshan Cattle but also contributes to broader efforts in germplasm conservation.

Comparative analysis of amphibian genomes: An emerging resource for basic and applied research.

Amphibians are the most threatened group of vertebrates and are in dire need of conservation intervention to ensure their continued survival. They exhibit unique features including a high diversity of reproductive strategies, permeable and specialized skin capable of producing toxins and antimicrobial compounds, multiple genetic mechanisms of sex determination and in some lineages, the ability to regenerate limbs and organs. Although genomic approaches would shed light on these unique traits and aid conservation, sequencing and assembly of amphibian genomes has lagged behind other taxa due to their comparatively large genome sizes. Fortunately, the development of long-read sequencing technologies and initiatives has led to a recent burst of new amphibian genome assemblies. Although growing, the field of amphibian genomics suffers from the lack of annotation resources, tools for working with challenging genomes and lack of high-quality assemblies in multiple clades of amphibians. Here, we analyse 51 publicly available amphibian genomes to evaluate their usefulness for functional genomics research. We report considerable variation in genome assembly quality and completeness and report some of the highest transposable element and repeat contents of any vertebrate. Additionally, we detected an association between transposable element content and climatic variables. Our analysis provides evidence of conserved genome synteny despite the long divergence times of this group, but we also highlight inconsistencies in chromosome naming and orientation across genome assemblies. We discuss sequencing gaps in the phylogeny and suggest key targets for future sequencing endeavours. Finally, we propose increased investment in amphibian genomics research to promote their conservation.

CNVoyant a machine learning framework for accurate and explainable copy number variant classification

The precise classification of copy number variants (CNVs) presents a significant challenge in genomic medicine, primarily due to the complex nature of CNVs and their diverse impact on rare genetic diseases (RGDs). This complexity is compounded by the limitations of existing methods in accurately distinguishing between benign, uncertain, and pathogenic CNVs. Addressing this gap, we introduce CNVoyant, a machine learning-based multi-class framework designed to enhance the clinical significance classification of CNVs. Trained on a comprehensive dataset of 52,176 ClinVar entries across pathogenic, uncertain, and benign classifications, CNVoyant incorporates a broad spectrum of genomic features, including genome position, disease-gene annotations, dosage sensitivity, and conservation scores. Models to predict the clinical significance of copy number gains and losses were trained independently. Final models were selected after testing 29 machine learning architectures and 10,000 hyperparameter combinations each for deletions and duplications via fivefold cross-validation. We validate the performance of CNVoyant by leveraging a comprehensive set of 21,574 CNVs from the DECIPHER database, a highly regarded resource known for its extensive catalog of chromosomal imbalances linked to clinical outcomes. Compared to alternative approaches, CNVoyant shows marked improvements in precision-recall and ROC AUC metrics for binary pathogenic classifications while going one step further, offering multi-classification of clinical significance and corresponding SHAP explainability plots. Additionally, when provided germline CNV calls from real-world RGD cases with diagnostic CNV(s), CNVoyant correctly classified all diagnostic CNVs as having pathogenic significance with high confidence. This large-scale validation demonstrates CNVoyant’s superior accuracy and underscores its potential to aid genomic researchers and clinical geneticists in interpreting the clinical implications of real CNVs.

From sequence to significance: A thorough investigation of the distinctive genome features uncovered in C. Werkmanii strain NIB003

Citrobacter werkmanii (C. werkmanii), an opportunistic urinary bacterium that causes diarrhea, is poorly understood. Our research focuses on genetic features that are crucial to disease development, such as pathogenic interactions, antibiotic resistance, virulence genes and genetic variation. Following its morphological, biochemical, and molecular identification, the whole genome of C. werkmanii strain NIB003 was sequenced in Bangladesh for the first time. Despite having around 80% whole genome conservation, the research shows that the Bangladeshi strain forms a separate phylogenetic cluster. This emphasises the genetic variability within C. werkmanii, resulting in particular modifications at the strain level and changes in its ability to cause disease. The results of the genetic diversity analysis indicate that the Bangladeshi sequenced genome is more diverse than the other strains due to the existence of unique features, such as the presence of t-RNA binding domain and N-6 adenine-specific DNA methylases.

Sequencing of Kaposi's Sarcoma Herpesvirus (KSHV) genomes from persons of diverse ethnicities and provenances with KSHV-associated diseases demonstrate multiple infections, novel polymorphisms, and low intra-host variance.

Recently published near full-length KSHV genomes from a Cameroon Kaposi sarcoma case-control study showed strong evidence of viral recombination and mixed infections, but no sequence variations associated with disease. Using the same methodology, an additional 102 KSHV genomes from 76 individuals with KSHV-associated diseases have been sequenced. Diagnoses comprise all KSHV-associated diseases (KAD): Kaposi sarcoma (KS), primary effusion lymphoma (PEL), KSHV-associated large cell lymphoma (KSHV-LCL), a type of multicentric Castleman disease (KSHV-MCD), and KSHV inflammatory cytokine syndrome (KICS). Participants originated from 22 different countries, providing the opportunity to obtain new near full-length sequences of a wide diversity of KSHV genomes. These include near full-length sequence of genomes with KSHV K1 subtypes A, B, C, and F as well as subtype E, for which no full sequence was previously available. High levels of recombination were observed. Fourteen individuals (18%) showed evidence of infection with multiple KSHV variants (from two to four unique genomes). Twenty-six comparisons of sequences, obtained from various sampling sites including PBMC, tissue biopsies, oral fluids, and effusions in the same participants, identified near complete genome conservation between different biological compartments. Polymorphisms were identified in coding and non-coding regions, including indels in the K3 and K15 genes and sequence inversions here reported for the first time. One such polymorphism in KSHV ORF46, specific to the KSHV K1 subtype E2, encoded a mutation in the leucine loop extension of the uracil DNA glycosylase that results in alteration of biochemical functions of this protein. This confirms that KSHV sequence variations can have functional consequences warranting further investigation. This study represents the largest and most diverse analysis of KSHV genome sequences to date among individuals with KAD and provides important new information on global KSHV genomics.

Genomic insights into ecological adaptation of oaks revealed by phylogenomic analysis of multiple species

Understanding the ecological adaptation of tree species can not only reveal the evolutionary potential but also benefit biodiversity conservation under global climate change. Quercus is a keystone genus in Northern Hemisphere forests, and its wide distribution in diverse ecosystems and long evolutionary history make it an ideal model for studying the genomic basis of ecological adaptations. Here we used a newly sequenced genome of Quercus gilva, an evergreen oak species from East Asia, with 18 published Fagales genomes to determine how Fagaceae genomes have evolved, identify genomic footprints of ecological adaptability in oaks in general, as well as between evergreen and deciduous oaks. We found that oak species exhibited a higher degree of genomic conservation and stability, as indicated by the absence of large-scale chromosomal structural variations or additional whole-genome duplication events. In addition, we identified expansion and tandem repetitions within gene families that contribute to plant physical and chemical defense (e.g., cuticle biosynthesis and oxidosqualene cyclase genes), which may represent the foundation for the ecological adaptation of oak species. Circadian rhythm and hormone-related genes may regulate the habits of evergreen and deciduous oaks. This study provides a comprehensive perspective on the ecological adaptations of tree species based on phylogenetic, genome evolutionary, and functional genomic analyses.

Chromosome-scale Elaeis guineensis and E. oleifera assemblies: comparative genomics of oil palm and other Arecaceae.

Elaeis guineensis and E. oleifera are the two species of oil palm. E. guineensis is the most widely cultivated commercial species, and introgression of desirable traits from E. oleifera is ongoing. We report an improved E. guineensis genome assembly with substantially increased continuity and completeness, as well as the first chromosome-scale E. oleifera genome assembly. Each assembly was obtained by integration of long-read sequencing, proximity ligation sequencing, optical mapping, and genetic mapping. High interspecific genome conservation is observed between the two species. The study provides the most extensive gene annotation to date, including 46,697 E. guineensis and 38,658 E. oleifera gene predictions. Analyses of repetitive element families further resolve the DNA repeat architecture of both genomes. Comparative genomic analyses identified experimentally validated small structural variants between the oil palm species and resolved the mechanism of chromosomal fusions responsible for the evolutionary descending dysploidy from 18 to 16 chromosomes.

Genomic conservation of Mongolian horses promoted by preservation of the intangible cultural heritage of Naadam in Mongolia

AbstractUncontrolled crossbreeding is a major challenge to the conservation of landrace horses in East Asia. Understanding the factors driving this trend is crucial for effective conservation efforts. Here, we investigate the genomic makeup of 40 Mongolian Naadam racehorses and 21 Asian landrace horse breeds through analyzing whole‐genome resequencing and Y chromosome data. Our results show that crossbreeding practices are linked to horse‐racing traditions. Regions characterized by strong horse‐racing traditions and a lack of crossbreeding regulations exhibit significant levels of exotic genetic introgression, as observed in populations from Inner Mongolia and Central Asia. However, in Mongolia, despite having strong horse‐racing traditions, the implementation of policies aimed at preserving traditional horse‐racing culture effectively reduces exotic introgression. These results suggest that horse‐racing traditions are the main driver of crossbreeding practices. Our research highlights that the preservation of traditional values in landrace horses, achieved through carefully managed horse‐racing activities, can lead to successful conservation outcomes.

An evolutionary view of the Fusarium core genome

Fusarium, a member of the Ascomycota fungi, encompasses several pathogenic species significant to plants and animals. Some phytopathogenic species have received special attention due to their negative economic impact on the agricultural industry around the world. Traditionally, identification and taxonomic analysis of Fusarium have relied on morphological and phenotypic features, including the fungal host, leading to taxonomic conflicts that have been solved using molecular systematic technologies. In this work, we applied a phylogenomic approach that allowed us to resolve the evolutionary history of the species complexes of the genus and present evidence that supports the F. ventricosum species complex as the most basal lineage of the genus. Additionally, we present evidence that proposes modifications to the previous hypothesis of the evolutionary history of the F. staphyleae, F. newnesense, F. nisikadoi, F. oxysporum, and F. fujikuroi species complexes. Evolutionary analysis showed that the genome GC content tends to be lower in more modern lineages, in both, the whole-genome and core-genome coding DNA sequences. In contrast, genome size gain and losses are present during the evolution of the genus. Interestingly, core genome duplication events positively correlate with genome size. Evolutionary and genome conservation analysis supports the F3 hypothesis of Fusarium as a more compact and conserved group in terms of genome conservation. By contrast, outside of the F3 hypothesis, the most basal clades only share 8.8% of its genomic sequences with the F3 clade.

Comparative Mitochondrial Genomic Analysis of Longhorn Beetles (Coleoptera: Chrysomeloidea) with Phylogenetic Implications

Longhorn beetles (Cerambycidae) play a vital role in global ecosystems. Some of them contribute to nutrient cycling and pollination, while others, pose a threat to forestry production. Despite their ecological importance, there has been a lack of comprehensive analyses on the mitochondrial genomes of Cerambycidae beetles. In this study, we have conducted mitochondrial genome sequencing and annotation for four Cerambycidae beetles: Monochamus sutor, Monochamus guerryi, Monochamus galloprovincialis, and Monochamus latefasciatus. Our analysis revealed a high degree of conservation in these mitochondrial genomes, with rare gene rearrangements observed across the Cerambycidae family. Additionally, a notable bias towards AT content was identified, with most genes using ATN as the start codon and TAA as the stop codon. Except for trnS1, all tRNA genes showed typical cloverleaf secondary structures. Phylogenetic analysis using IQ-TREE and Phylobayes consistently produced congruent topologies. At the gene level analyses, our results highlighted the remarkable conservation of the COX1 gene. Furthermore, at the species level, we observed strong adaptability in the Spondylidinae and Lepturinae subfamilies. We also offer our insights into contentious aspects of the phylogeny. Overall, our research contributes to a deeper understanding of the phylogeny and evolution of Cerambycidae, laying the groundwork for future population genetic investigations.

Genomic and virulence analysis of in vitro cultured Cryptosporidium parvum.

Recent advances in the in vitro cultivation of Cryptosporidium parvum using hollow fiber bioreactor technology (HFB) have permitted continuous growth of parasites that complete all life cycle stages. The method provides access to all stages of the parasite and provides a method for non-animal production of oocysts for use in clinical trials. Here we examined the effect of long-term (>20 months) in vitro culture on virulence-factors, genome conservation, and in vivo pathogenicity of the host by in vitro cultured parasites. We find low-level sequence variation that is consistent with that observed in calf-passaged parasites. Further using a calf model infection, oocysts obtained from the HFB caused diarrhea of the same volume, duration and oocyst shedding intensity as in vivo passaged parasites.

Flexible genomic island conservation across freshwater and marine Methylophilaceae.

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

Insight into the current genomic diversity, conservation status and population structure of Tunisian Barbarine sheep breed.

Local livestock breeds play a crucial role in global biodiversity, connecting natural and human-influenced environments and contributing significantly to ecosystem services. While commercial breeds dominate industrial systems, local livestock breeds in developing countries, like Barbarine sheep in Tunisia, are vital for food security and community maintenance. The Tunisian Barbarine sheep, known for its adaptability and distinctive fat-tailed morphology, faces challenges due to historical crossbreeding. In this study, the Illumina Ovine SNP50K BeadChip array was used to perform a genome-wide characterization of Tunisian Barbarine sheep to investigate its genetic diversity, the genome structure, and the relationship within the context of Mediterranean breeds. The results show moderate genetic diversity and low inbreeding. Runs of Homozygosity analysis find genomic regions linked to important traits, including fat tail characteristics. Genomic relationship analysis shows proximity to Algerian thin-tailed breeds, suggesting crossbreeding impacts. Admixture analysis reveals unique genetic patterns, emphasizing the Tunisian Barbarine's identity within the Mediterranean context and its closeness to African breeds. Current results represent a starting point for the creation of monitoring and conservation plans. In summary, despite genetic dilution due to crossbreeding, the identification of genomic regions offers crucial insights for conservation. The study confirms the importance of preserving unique genetic characteristics of local breeds, particularly in the face of ongoing crossbreeding practices and environmental challenges. These findings contribute valuable insights for the sustainable management of this unique genetic reservoir, supporting local economies and preserving sheep species biodiversity.

Specific Patterns in Correlations of Super-Short Tandem Repeats (SSTRs) with G+C Content, Genic and Intergenic Regions, and Retrotransposons on All Human Chromosomes.

The specific characteristics of k-mer words (2 ≤ k ≤ 11) regarding genomic distribution and evolutionary conservation were recently found. Among them are, in high abundance, words with a tandem repeat structure (repeat unit length of 1 bp to 3 bp). Furthermore, there seems to be a class of extremely short tandem repeats (≤12 bp), so far overlooked, that are non-random-distributed and, therefore, may play a crucial role in the functioning of the genome. In the following article, the positional distributions of these motifs we call super-short tandem repeats (SSTRs) were compared to other functional elements, like genes and retrotransposons. We found length- and sequence-dependent correlations between the local SSTR density and G+C content, and also between the density of SSTRs and genes, as well as correlations with retrotransposon density. In addition to many general interesting relations, we found that SINE Alu has a strong influence on the local SSTR density. Moreover, the observed connection of SSTR patterns to pseudogenes and -exons might imply a special role of SSTRs in gene expression. In summary, our findings support the idea of a special role and the functional relevance of SSTRs in the genome.

Clonality and Diversity in the Soft Rot Dickeya solani Phytopathogen.

Bacterial diversity analyses often suffer from a bias due to sampling only from a limited number of hosts or narrow geographic locations. This was the case for the phytopathogenic species Dickeya solani, whose members were mainly isolated from a few hosts-potato and ornamentals-and from the same geographical area-Europe and Israel, which are connected by seed trade. Most D. solani members were clonal with the notable exception of the potato isolate RNS05.1.2A and two related strains that are clearly distinct from other D. solani genomes. To investigate if D. solani genomic diversity might be broadened by analysis of strains isolated from other environments, we analysed new strains isolated from ornamentals and from river water as well as strain CFBP 5647 isolated from tomato in the Caribbean island Guadeloupe. While water strains were clonal to RNS05.1.2A, the Caribbean tomato strain formed a third clade. The genomes of the three clades are highly syntenic; they shared almost 3900 protein families, and clade-specific genes were mainly included in genomic islands of extrachromosomal origin. Our study thus revealed both broader D. solani diversity with the characterisation of a third clade isolated in Latin America and a very high genomic conservation between clade members.

Bioactive glycans in a microbiome-directed food for children with malnutrition

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1–4. Here we analyse biospecimens from a randomized, controlled trial of a microbiome-directed complementary food (MDCF-2) that produced superior rates of weight gain compared with a calorically more dense conventional ready-to-use supplementary food in 12–18-month-old Bangladeshi children with moderate acute malnutrition4. We reconstructed 1,000 bacterial genomes (metagenome-assembled genomes (MAGs)) from the faecal microbiomes of trial participants, identified 75 MAGs of which the abundances were positively associated with ponderal growth (change in weight-for-length Z score (WLZ)), characterized changes in MAG gene expression as a function of treatment type and WLZ response, and quantified carbohydrate structures in MDCF-2 and faeces. The results reveal that two Prevotella copri MAGs that are positively associated with WLZ are the principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilizing the component glycans of MDCF-2. The predicted specificities of carbohydrate-active enzymes expressed by their polysaccharide-utilization loci are correlated with (1) the in vitro growth of Bangladeshi P. copri strains, possessing varying degrees of polysaccharide-utilization loci and genomic conservation with these MAGs, in defined medium containing different purified glycans representative of those in MDCF-2, and (2) the levels of faecal carbohydrate structures in the trial participants. These associations suggest that identifying bioactive glycan structures in MDCFs metabolized by growth-associated bacterial taxa will help to guide recommendations about their use in children with acute malnutrition and enable the development of additional formulations.

Topological structures and syntenic conservation in sea anemone genomes

There is currently little information about the evolution of gene clusters, genome architectures and karyotypes in early branching animals. Slowly evolving anthozoan cnidarians can be particularly informative about the evolution of these genome features. Here we report chromosome-level genome assemblies of two related anthozoans, the sea anemones Nematostella vectensis and Scolanthus callimorphus. We find a robust set of 15 chromosomes with a clear one-to-one correspondence between the two species. Both genomes show chromosomal conservation, allowing us to reconstruct ancestral cnidarian and metazoan chromosomal blocks, consisting of at least 19 and 16 ancestral linkage groups, respectively. We show that, in contrast to Bilateria, the Hox and NK clusters of investigated cnidarians are largely disintegrated, despite the presence of staggered hox/gbx expression in Nematostella. This loss of microsynteny conservation may be facilitated by shorter distances between cis-regulatory sequences and their cognate transcriptional start sites. We find no clear evidence for topologically associated domains, suggesting fundamental differences in long-range gene regulation compared to vertebrates. These data suggest that large sets of ancestral metazoan genes have been retained in ancestral linkage groups of some extant lineages; yet, higher order gene regulation with associated 3D architecture may have evolved only after the cnidarian-bilaterian split.

The genomes of Hercules beetles reveal putative adaptive loci and distinct demographic histories in pristine North American forests.

Beetles, despite their remarkable biodiversity and a long history of research, remain lacking in reference genomes annotated with structural variations in loci of adaptive significance. We sequenced and assembled high-quality chromosome-level genomes of four Hercules beetles which exhibit divergence in male horn size and shape and body colouration. The four Hercules beetle genomes were assembled to 11 pseudo-chromosomes, where the three genomes assembled using Nanopore data (Dynastes grantii, D. hyllus and D. tityus) were mapped to the genome assembled using PacBio + Hi-C data (D. maya). We demonstrated a striking similarity in genome structure among the four species. This conservative genome structure may be attributed to our use of the D. maya assembly as the reference; however, it is worth noting that such a conservative genome structure is a recurring phenomenon among scarab beetles. We further identified homologues of nine and three candidate-gene families that may be associated with the evolution of horn structure and body colouration respectively. Structural variations in Scr and Ebony2 were detected and discussed for their putative impacts on generating morphological diversity in beetles. We also reconstructed the demographic histories of the four Hercules beetles using heterozygosity information from the diploid genomes. We found that the demographic histories of the beetles closely recapitulated historical changes in suitable forest habitats driven by climate shifts.

Identification of a morphogene required for tapered filament termini in filamentous cyanobacteria.

Although the photosynthetic cyanobacteria are monophyletic, they exhibit substantial morphological diversity across species, and even within an individual species due to phenotypic plasticity in response to life cycles and environmental signals. This is particularly prominent among the multicellular filamentous cyanobacteria. One example of this is the appearance of tapering at the filament termini. However, the morphogenes controlling this phenotype and the adaptive function of this morphology are not well defined. Here, using the model filamentous cyanobacterium Nostoc punctiforme ATCC29133 (PCC73102), we identify tftA, a morphogene required for the development of tapered filament termini. The tftA gene is specifically expressed in developing hormogonia, motile trichomes where the tapered filament morphology is observed, and encodes a protein containing putative amidase_3 and glucosaminidase domains, implying a function in peptidoglycan hydrolysis. Deletion of tftA abolished filament tapering inidcating that TftA plays a role in remodelling the cell wall to produce tapered filaments. Genomic conservation of tftA specifically in filamentous cyanobacteria indicates this is likely to be a conserved mechanism among these organisms. Finally, motility assays indicate that filaments with tapered termini migrate more efficiently through dense substratum, providing a plausible biological role for this morphology.

DNA Barcoding Based Identification of Rosa x Damascene and Prunus dulcis Herbs Using ITS2 Barcoding Gene Amplification.

The considerable risk of adulteration in the herbs has raised commercial interest in the identification of medicinal herbs globally. DNA barcoding is the primary techniques for identifying the herbs at genetic level. This technique's key benefit is that it can identify the material's purity. This study focuses on the accurate identification of species utilizing Polymerase chain reaction-based nuclear universal internal transcribed spacer region (a barcode region) amplification and sequencing in 2 medically significant plants (Rosa damascene and Prunus dulcis) procured from the local herbal market of Karachi in the year 2022. Results suggest that the ribosomal nuclear ITS2 region of the selected plant species shows the 100 % identity with the reference genome, therefore it has shown a good rate of identification at the species level. The findings of this concludes that ITS2 is the novel standard barcode that involve in species identification, genomic conservation, and secure utilization of the medically significant plant species.