Diversity of selected toll-like receptor genes in cheetahs (Acinonyx jubatus) and African leopards (Panthera pardus pardus).

BackgroundThe growing world population amplifies the anthropogenic impact on wildlife globally. With shrinking habitats, wild populations are being pushed to co-exist in close proximity to humans, leading to an increased threat of infectious disease. Therefore, understanding the immune system of a species is key to assess its resilience in a changing environment. The innate immunity system (IIS) is the body’s first line of defense against pathogens. High variability in IIS-genes, such as the toll-like receptor (TLR) genes, appears to be associated with resistance to infectious diseases. However, few studies have investigated diversity in TLR genes in non-model organisms and drawn conclusions for the conservation of vulnerable species. Large predators are threatened globally, and their populations increasingly have been declining over the last decades. Big cats, such as leopards (Panthera pardus) and cheetahs (Acinonyx jubatus) are no exception to this trend and are listed as ‘vulnerable’ by the International Union for Conservation of Nature (IUCN) including several subspecies, e.g., A. j. venaticus and P. p. melas, that already face extinction. To better understand vulnerability in terms of immune genetic diversity in the two sympatric occurring species, we compared selected TLR genes (TLR2, TLR4, TLR6 and TLR8) between modern African leopards (P. p. pardus) and Southern African cheetahs (A. j. jubatus). ResultsOur study supports the previously detected high genetic diversity in African leopards and confirms genetic impoverishment in Southern African cheetahs. Despite notable differences, both species share some haplotypic similarities in the investigated TLRs. Moreover, our historic cheetah samples from all five subspecies showed levels of genetic diversity comparable to modern African leopards. By including historic cheetahs and samples from all known subspecies, we put the observed IIS diversity into an evolutionary context.ConclusionThe genetic diversity in the investigated TLR genes in modern Southern African cheetahs and in historic cheetahs is low compared to African leopards. However, according to previous studies, the low immune genetic diversity might not yet affect the health of this cheetah subspecies. Compared to historic cheetah data and other subspecies, a more recent population decline might explain the observed genetic impoverishment of TLR genes in modern Southern African cheetahs.<br/

Nucleotide diversity of functionally different groups of immune response genes in Old World camels based on newly annotated and reference-guided assemblies

Background Immune-response (IR) genes have an important role in the defense against highly variable pathogens, and therefore, diversity in these genomic regions is essential for species' survival and adaptation. Although current genome assemblies from Old World camelids are very useful for investigating genome-wide diversity, demography and population structure, they have inconsistencies and gaps that limit analyses at local genomic scales. Improved and more accurate genome assemblies and annotations are needed to study complex genomic regions like adaptive and innate IR genes. Results In this work, we improved the genome assemblies of the three Old World camel species - domestic dromedary and Bactrian camel, and the two-humped wild camel - via different computational methods. The newly annotated dromedary genome assembly CamDro3 served as reference to scaffold the NCBI RefSeq genomes of domestic Bactrian and wild camels. These upgraded assemblies were then used to assess nucleotide diversity of IR genes within and between species, and to compare the diversity found in immune genes and the rest of the genes in the genome. We detected differences in the nucleotide diversity among the three Old World camelid species and between IR gene groups, i.e., innate versus adaptive. Among the three species, domestic Bactrian camels showed the highest mean nucleotide diversity. Among the functionally different IR gene groups, the highest mean nucleotide diversity was observed in the major histocompatibility complex. Conclusions The new camel genome assemblies were greatly improved in terms of contiguity and increased size with fewer scaffolds, which is of general value for the scientific community. This allowed us to perform in-depth studies on genetic diversity in immunity-related regions of the genome. Our results suggest that differences of diversity across classes of genes appear compatible with a combined role of population history and differential exposures to pathogens, and consequent different selective pressures.Peer reviewe

Improving Illumina assemblies with Hi-C and long reads : An example with the North African dromedary

Researchers have assembled thousands of eukaryotic genomes using Illumina reads, but traditional mate-pair libraries cannot span all repetitive elements, resulting in highly fragmented assemblies. However, both chromosome conformation capture techniques, such as Hi-C and Dovetail Genomics Chicago libraries and long-read sequencing, such as Pacific Biosciences and Oxford Nanopore, help span and resolve repetitive regions and therefore improve genome assemblies. One important livestock species of arid regions that does not have a high-quality contiguous reference genome is the dromedary (Camelus dromedarius). Draft genomes exist but are highly fragmented, and a high-quality reference genome is needed to understand adaptation to desert environments and artificial selection during domestication. Dromedaries are among the last livestock species to have been domesticated, and together with wild and domestic Bactrian camels, they are the only representatives of the Camelini tribe, which highlights their evolutionary significance. Here we describe our efforts to improve the North African dromedary genome. We used Chicago and Hi-C sequencing libraries from Dovetail Genomics to resolve the order of previously assembled contigs, producing almost chromosome-level scaffolds. Remaining gaps were filled with Pacific Biosciences long reads, and then scaffolds were comparatively mapped to chromosomes. Long reads added 99.32 Mbp to the total length of the new assembly. Dovetail Chicago and Hi-C libraries increased the longest scaffold over 12-fold, from 9.71 Mbp to 124.99 Mbp and the scaffold N50 over 50-fold, from 1.48 Mbp to 75.02 Mbp. We demonstrate that Illumina de novo assemblies can be substantially upgraded by combining chromosome conformation capture and long-read sequencing.Peer reviewe

The potential and shortcomings of mitochondrial DNA analysis for cheetah conservation management

There are only about 7,100 adolescent and adult cheetahs (Acinonyx jubatus) remaining in the wild. With the majority occurring outside protected areas, their numbers are rapidly declining. Evidence-based conservation measures are essential for the survival of this species. Genetic data is routinely used to inform conservation strategies, e.g., by establishing conservation units (CU). A commonly used marker in conservation genetics is mitochondrial DNA (mtDNA). Here, we investigated the cheetah’s phylogeography using a large-scale mtDNA data set to refine subspecies distributions and better assign individuals to CUs. Our dataset mostly consisted of historic samples to cover the cheetah’s whole range as the species has been extinct in most of its former distribution. While our genetic data largely agree with geography-based subspecies assignments, several geographic regions show conflicting mtDNA signals. Our analyses support previous findings that evolutionary forces such as incomplete lineage sorting or mitochondrial capture likely confound the mitochondrial phylogeography of this species, especially in East and, to some extent, in Northeast Africa. We caution that subspecies assignments solely based on mtDNA should be treated carefully and argue for an additional standardized nuclear single nucleotide polymorphism (SNP) marker set for subspecies identification and monitoring. However, the detection of the A. j. soemmeringii specific haplogroup by a newly designed Amplification-Refractory Mutation System (ARMS) can already provide support for conservation measures.info:eu-repo/semantics/publishedVersio

Genomic Analyses Show Extremely Perilous Conservation Status of African and Asiatic cheetahs (Acinonyx jubatus).

We live in a world characterized by biodiversity loss and global environmental change. The extinction of large carnivores can have ramifying effects on ecosystems like an uncontrolled increase in wild herbivores, which in turn can have knock-on impacts on vegetation regeneration and communities. Cheetahs (Acinonyx jubatus) serve important ecosystem functions as apex predators; yet, they are quickly heading towards an uncertain future. Threatened by habitat loss, human-wildlife conflict and illegal trafficking, there are only approximately 7100 individuals remaining in nature. We present the most comprehensive genome-wide analysis of cheetah phylogeography and conservation genomics to date, assembling samples from nearly the entire current and past species' range. We show that their phylogeography is more complex than previously thought, and that East African cheetahs (A. j. raineyi) are genetically distinct from Southern African individuals (A. j. jubatus), warranting their recognition as a distinct subspecies. We found strong genetic differentiation between all classically recognized subspecies, thus refuting earlier findings that cheetahs show only little differentiation. The strongest differentiation was observed between the Asiatic and all the African subspecies. We detected high inbreeding in the Critically Endangered Iranian (A. j. venaticus) and North-western (A. j. hecki) subspecies, and show that overall cheetahs, along with snow leopards, have the lowest genome-wide heterozygosity of all the big cats. This further emphasizes the cheetah's perilous conservation status. Our results provide novel and important information on cheetah phylogeography that can support evidence-based conservation policy decisions to help protect this species. This is especially relevant in light of ongoing and proposed translocations across subspecies boundaries, and the increasing threats of illegal trafficking

Innate and adaptive immune genes associated with MERS-CoV infection in dromedaries

The recent SARS-CoV-2 pandemic has refocused attention to the betacoronaviruses, only eight years after the emergence of another zoonotic betacoronavirus, the Middle East respiratory syndrome coronavirus (MERS-CoV). While the wild source of SARS-CoV-2 may be disputed, for MERS-CoV, dromedaries are considered as source of zoonotic human infections. Testing 100 immune-response genes in 121 dromedaries from United Arab Emirates (UAE) for potential association with present MERS-CoV infection, we identified candidate genes with important functions in the adaptive, MHC-class I (HLA-A-24-like) and II (HLA-DPB1-like), and innate immune response (PTPN4, MAGOHB), and in cilia coating the respiratory tract (DNAH7). Some of these genes previously have been associated with viral replication in SARS-CoV-1/-2 in humans, others have an important role in the movement of bronchial cilia. These results suggest similar host genetic pathways associated with these betacoronaviruses, although further work is required to better understand the MERS-CoV disease dynamics in both dromedaries and humans

A Deadly Cargo: Gene Repertoire of Cytotoxic Effector Proteins in the Camelidae

Cytotoxic T cells and natural killer cells can kill target cells based on their expression and release of perforin, granulysin, and granzymes. Genes encoding these molecules have been only poorly annotated in camelids. Based on bioinformatic analyses of genomic resources, sequences corresponding to perforin, granulysin, and granzymes were identified in genomes of camelids and related ungulate species, and annotation of the corresponding genes was performed. A phylogenetic tree was constructed to study evolutionary relationships between the species analyzed. Re-sequencing of all genes in a panel of 10 dromedaries and 10 domestic Bactrian camels allowed analyzing their individual genetic polymorphisms. The data showed that all extant Old World camelids possess functional genes for two pore-forming proteins (PRF1, GNLY) and six granzymes (GZMA, GZMB, GZMH, GZMK, GZMM, and GZMO). All these genes were represented as single copies in the genome except the GZMH gene exhibiting interspecific differences in the number of loci. High protein sequence similarities with other camelid and ungulate species were observed for GZMK and GZMM. The protein variability in dromedaries and Bactrian camels was rather low, except for GNLY and chymotrypsin-like granzymes (GZMB, GZMH)