Camelid genomes reveal evolution and adaptation to desert environments

Camelids: new players in the international animal production context

The Camelidae family comprises the Bactrian camel (Camelus bactrianus), the dromedary camel (Camelus dromedarius), and four species of South American camelids: llama (Lama glama), alpaca (Lama pacos) guanaco (Lama guanicoe), and vicuña (Vicugna vicugna). The main characteristic of these species is their ability to cope with either hard climatic conditions like those found in arid regions (Bactrian and dromedary camels) or high-altitude landscapes like those found in South America (South American camelids). Because of such interesting physiological and adaptive traits, the interest for these animals as livestock species has increased considerably over the last years. In general, the main animal products obtained from these animals are meat, milk, and hair fiber, although they are also used for races and work among other activities. In the near future, climate change will likely decrease agricultural areas for animal production worldwide, particularly in the tropics and subtropics where competition with crops for human consumption is a major problem already. In such conditions, extensive animal production could be limited in some extent to semi-arid rangelands, subjected to periodical draughts and erratic patterns of rainfall, severely affecting conventional livestock production, namely cattle and sheep. In the tropics and subtropics, camelids may become an important protein source for humans. This article aims to review some of the recent literature about the meat, milk, and hair fiber production in the six existing camelid species highlighting their benefits and drawbacks, overall contributing to the development of camelid production in the framework of food security.

Deiminated proteins in extracellular vesicles and serum of llama (Lama glama)-Novel insights into camelid immunity

Peptidylarginine deiminases (PADs) are phylogenetically conserved calcium-dependent enzymes which post-translationally convert arginine into citrulline in target proteins in an irreversible manner, causing functional and structural changes in target proteins. Protein deimination causes generation of neo-epitopes, affects gene regulation and also allows for protein moonlighting. Furthermore, PADs have been found to be a phylogenetically conserved regulator for extracellular vesicle (EVs) release. EVs are found in most body fluids and participate in cellular communication via transfer of cargo proteins and genetic material. In this study, post-translationally deiminated proteins in serum and serum-EVs are described for the first time in camelids, using the llama (Lama glama L. 1758) as a model animal. We report a poly-dispersed population of llama serum EVs, positive for phylogenetically conserved EV-specific markers and characterised by TEM. In serum, 103 deiminated proteins were overall identified, including key immune and metabolic mediators including complement components, immunoglobulin-based nanobodies, adiponectin and heat shock proteins. In serum, 60 deiminated proteins were identified that were not in EVs, and 25 deiminated proteins were found to be unique to EVs, with 43 shared deiminated protein hits between both serum and EVs. Deiminated histone H3, a marker of neutrophil extracellular trap formation, was also detected in llama serum. PAD homologues were identified in llama serum by Western blotting, via cross reaction with human PAD antibodies, and detected at an expected 70 kDa size. This is the first report of deiminated proteins in serum and EVs of a camelid species, highlighting a hitherto unrecognized post-translational modification in key immune and metabolic proteins in camelids, which may be translatable to and inform a range of human metabolic and inflammatory pathologies.

Divergent selection and drift shape the genomes of two avian sister species spanning a saline-freshwater ecotone

The role of species divergence due to ecologically based divergent selection-or ecological speciation-in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline-freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes-the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)-to evaluate the influence of selective and demographic processes in shaping genome-wide patterns of divergence. Genome-wide divergence between these two recently diverged sister species was notably high (genome-wide F ST = 0.32). Against a background of high genome-wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism-all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.

Thermotolerance and plasticity of camel somatic cells exposed to acute and chronic heat stress

The Arabian camel is the largest known mammal that can survive in severe hot climatic conditions. We provide the molecular explanation for the thermotolerance of camel granulosa somatic cells after exposure to 45 °C for 2 (acute heat shock) or 20 h (chronic heat shock). The common features of the cellular responses to acute heat stress were the increase of heat shock proteins and DNA repair enzymes expression. Actin polymerization and Rho signaling were critically activated as a cellular defense against heat shock. Cells exposed to chronic heat shock showed altered cell architecture with a decrease in total detected proteins, metabolic enzymes, and cytoskeletal protein expression. Treatment with transforming growth factor beta (TGFβ) pathway inhibitor SB-431542 suppressed the morphological alterations of cells exposed to chronic heat shock. Moreover, during the recovery stage at 38 °C for 24 h, proteomic changes were partially restored with an exponential increase in HSP70 expression, and the cells restored their normal cellular morphology on the 9th day of recovery. Full proteomics data are available via ProteomeXchange with identifier PXD012159. The strategies of cellular defense and tolerance to both thermal conditions reflect the flexible adaptability of camel somatic cells to conserve life under extremely hot conditions.

Genome-wide characterization and analysis of microsatellite sequences in camelid species

Microsatellites or simple sequence repeats (SSRs) are among the genetic markers most widely utilized in research. This includes applications in numerous fields such as genetic conservation, paternity testing, and molecular breeding. Though ordered draft genome assemblies of camels have been announced, including for the Arabian camel, systemic analysis of camel SSRs is still limited. The identification and development of informative and robust molecular SSR markers are essential for marker assisted breeding programs and paternity testing. Here we searched and compared perfect SSRs with 1–6 bp nucleotide motifs to characterize microsatellites for draft genome sequences of the Camelidae. We analyzed and compared the occurrence, relative abundance, relative density, and guanine-cytosine (GC) content in four taxonomically different camelid species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. A total of 546762, 544494, 547974, and 437815 SSRs were mined, respectively. Mononucleotide SSRs were the most frequent in the four genomes, followed in descending order by di-, tetra-, tri-, penta-, and hexanucleotide SSRs. GC content was highest in dinucleotide SSRs and lowest in mononucleotide SSRs. Our results provide further evidence that SSRs are more abundant in noncoding regions than in coding regions. Similar distributions of microsatellites were found in all four species, which indicates that the pattern of microsatellites is conserved in family Camelidae.

The Camel Adaptive Immune Receptors Repertoire as a Singular Example of Structural and Functional Genomics

The adaptive immune receptors repertoire is highly plastic, with its ability to produce antigen-binding molecules and select those with high affinity for their antigen. Species have developed diverse genetic and structural strategies to create their respective repertoires required for their survival in the different environments. Camelids, until now, considered as a case of evolutionary innovation because of their only heavy-chain antibodies, represent a new mammalian model particularly useful for understanding the role of diversity in the immune system function. Here, we review the structural and functional characteristics and the current status of the genomic organization of camel immunoglobulins (IG) or antibodies, α/ß and γ/δ T cell receptors (TR), and major histocompatibility complex (MHC). In camelid humoral response, in addition to the conventional antibodies, there are IG with “only-heavy-chain” (no light chain, and two identical heavy gamma chains lacking CH1 and with a VH domain designated as VHH). The unique features of these VHH offer advantages in biotechnology and for clinical applications. The TRG and TRD rearranged variable domains of Camelus dromedarius (Arabian camel) display somatic hypermutation (SHM), increasing the intrinsic structural stability in the γ/δ heterodimer and influencing the affinity maturation to a given antigen similar to immunoglobulin genes. The SHM increases the dromedary γ/δ repertoire diversity. In Camelus genus, the general structural organization of the TRB locus is similar to that of the other artiodactyl species, with a pool of TRBV genes positioned at the 5’ end of three in tandem D-J-C clusters, followed by a single TRBV gene with an inverted transcriptional orientation located at the 3’ end. At the difference of TRG and TRD, the diversity of the TRB variable domains is not shaped by SHM and depends from the classical combinatorial and junctional diversity. The MHC locus is located on chromosome 20 in Camelus dromedarius. Cytogenetic and comparative whole genome analyses revealed the order of the three major regions “Centromere-ClassII-ClassIII-ClassI”. Unexpectedly low extent of polymorphisms and haplotypes was observed in all Old World camels despite different geographic origins.

The Major Histocompatibility Complex of Old World Camels-A Synopsis

This study brings new information on major histocompatibility complex (MHC) class III sub-region genes in Old World camels and integrates current knowledge of the MHC region into a comprehensive overview for Old World camels. Out of the MHC class III genes characterized, TNFA and the LY6 gene family showed high levels of conservation, characteristic for MHC class III loci in general. For comparison, an MHC class II gene TAP1, not coding for antigen presenting molecules but functionally related to MHC antigen presenting functions was studied. TAP1 had many SNPs, even higher than the MHC class I and II genes encoding antigen presenting molecules. Based on this knowledge and using new camel genomic resources, we constructed an improved genomic map of the entire MHC region of Old World camels. The MHC class III sub-region shows a standard organization similar to that of pig or cattle. The overall genomic structure of the camel MHC is more similar to pig MHC than to cattle MHC. This conclusion is supported by differences in the organization of the MHC class II sub-region, absence of functional DY genes, different organization of MIC genes in the MHC class I sub-region, and generally closer evolutionary relationships of camel and porcine MHC gene sequences analyzed so far.

Thermotolerance of camel (Camelus dromedarius) somatic cells affected by the cell type and the dissociation method

Researchers dealing with heat stress experiments use different cell kinds and use trypsin that has been reported to affect the cellular proteins of cultured cells. Therefore, we compared the effects of acute and chronic exposures to high temperature (45 °C) on camel skin fibroblast and granulosa cells. Primary culture of fibroblasts and granulosa cells tolerated the acute heat shock for 2 h; however, granulosa cells cultured for long duration (20 h) showed thermotolerance when compared with the fibroblasts. Moreover, the effect of cell dispersion method (trypsin and mechanical dissociation) on the thermotolerance of sub-cultured cells was examined. Trypsin altered the morphology of fibroblasts and granulosa cells exposed to 45 °C for 4 h. Moreover, trypsin significantly reduced the fibroblast and granulosa cell migration in the wound healing assay. The current results demonstrate that cell passaging and cell type can affect the thermotolerance of the cells; it also revealed that trypsin could alter the cellular response to the heat shock. We raise the demand for another alternative method for cell dispersion in experiments dealing with cellular responses to the heat shock.

Different maternal lineages revealed by ancient mitochondrial genome of Camelus bactrianus from China

Domestic Bactrian camel (Camelus bactrianus) used to be one of the most important livestock species in Chinese history, as well as the major transport carrier on the ancient Silk Road. However, archeological studies on Chinese C. bactrianus are still limited, and molecular biology research on this species is mainly focused on modern specimens. In this study, we retrieved the complete mitochondrial genome from a C. bactrianus specimen, which was excavated from northwestern China and dated at 1290-1180 cal. years before present (yBP). Phylogenetic analyses using 18 mitochondrial genomes indicated that the C. bactrianus clade was divided into two maternal lineages. The majority of samples originating from Iran to Japan and Mongolia belong to subclade A1, while our sample together with two Mongolian individuals formed the much smaller subclade A2. Furthermore, the divergence time of these two maternal lineages was estimated as 165 Kya (95% credibility interval 117-222 Kya), this might indicate that several different evolutionary lineages were incorporated into the domestic gene pool during the initial domestication process. Bayesian skyline plot (BSP) analysis suggest a slow increase in female effective population size of C. bactrianus from 5000 years ago, which corresponds to the beginning of domestication of C. bactrianus. The present study also revealed that there were extensive exchanges of genetic information among C. bactrianus populations in regions along the Silk Road.

Polymorphisms in MC1R and ASIP Genes are Associated with Coat Color Variation in the Arabian Camel

Pigmentation in mammals is primarily determined by the distribution of eumelanin and pheomelanin, the ratio of which is mostly controlled by the activity of melanocortin 1 receptor (MC1R) and agouti signaling protein (ASIP) genes. Using 91 animals from 10 Arabian camel populations, that included the 4 predominant coat color phenotypes observed in the dromedary (light brown, dark brown, black, and white), we investigated the effects of the MC1R and ASIP sequence variants and identified candidate polymorphisms associated with coat color variation. In particular, we identified a single nucleotide polymorphism (SNP), found in the coding region of MC1R (901C/T), linked to the white coat color, whereas a 1-bp deletion (23delT/T) and a SNP (25G/A) in exon 2 of ASIP are associated with both black and dark-brown coat colors. Our results also indicate support that the light-brown coat color is likely the ancestral coat color for the dromedary. These sequence variations at the MC1R and ASIP genes represent the first documented evidence of candidate polymorphisms associated with Mendelian traits in the dromedary.

Genome Diversity and Signatures of Selection for Production and Performance Traits in Dromedary Camels

Dromedary camels (Camelus dromedarius) are single-humped animals found throughout the deserts of Africa, the Arabian Peninsula, and the southwest of Asia. This well-adapted species is mainly used for milk and meat production, although some specific types exhibit superior running performance and are used in racing competitions. However, neither performance nor production camels are bred under intensive genomic selection programs with specific aims to improve these traits. In this study, the full genome sequence data of six camels from the Arabian Peninsula and the genotyping-by-sequencing data of 44 camels (29 packing and 15 racing) from Sudan were analyzed to assess their genome diversities, relationships, and candidate signatures of positive selection. Genome ADMIXTURE and principle component analyses indicate clear geographic separation between the Sudanese and the Arabian Peninsula camels, but with no population-specific genetic distinction within populations. Camel samples from the Arabian Peninsula show higher mean heterozygosity (0.560 ± 0.003) than those from Sudan (0.347 ± 0.003). Analyses of signatures of selection, using pooled heterozygosity (Hp) approach, in the Sudanese camels revealed 176, 189, and 308 candidate regions under positive selection in the combined and packing and racing camel populations, respectively. These regions host genes that might be associated with adaptation to arid environment, dairy traits, energy homeostasis, and chondrogenesis. Eight regions show high genetic differentiation, based on Fst analysis, between the Sudanese packing and racing camel types. Genes associated with chondrogenesis, energy balance, and urinary system development were found within these regions. Our results advocate for further detailed investigation of the genome of the dromedary camel to identify and characterize genes and variants associated with their valuable phenotypic traits. The results of which may support the development of breeding programs to improve the production and performance traits of this unique domesticated species.

Old World camels in a modern world - a balancing act between conservation and genetic improvement

Old World camels have served humans in cross‐continental caravans, transporting people and goods, connecting different cultures and providing milk, meat, wool and draught since their domestication around 3000–6000 years ago. In a world of modern transport and fast connectivity, these beasts of burden seem to be out‐dated. However, a growing demand for sustainable milk and meat production, especially in countries affected by climate change and increasing desertification, brings dromedaries (Camelus dromedarius) and Bactrian camels (Camelus bactrianus) back onstage and into the focus of animal breeders and scientists. In this review on the molecular genetics of these economically important species we give an overview about the evolutionary history, domestication and dispersal of Old World camels, whereas highlighting the need for conservation of wild two‐humped camels (Camelus ferus) as an evolutionarily unique and highly endangered species. We provide cutting‐edge information on the current molecular resources and on‐going sequencing projects. We cannot emphasise enough the importance of balancing the need for improving camel production traits with maintaining the genetic diversity in two domestic species with specific physiological adaptation to a desert environment.

Pleurotus eryngii Genomes Reveal Evolution and Adaptation to the Gobi Desert Environment

Pleurotus eryngii (King Oyster) is one of the most highly prized edible mushrooms. Among the diverse varieties within P. eryngii, P. eryngii var. eryngii is the commonest one, with a worldwide distribution, while P. eryngii var. ferulae is only distributed in Europe and China, and is especially adapted to the Gobi Desert in Xinjiang Autonomous Region of China. However, little is known about the genome-wide pattern of evolution and adaptation to the divergent environments of P. eryngii. Here, we present the high-quality genome sequences of P. eryngii var. eryngii strain PEE81 originating from Europe and P. eryngii var. ferulae strain PEF12 originating from the Gobi Desert of China. The assembled genome sizes of PEE81 and PEF12 were 53.6 and 48.0 Mbp, respectively, which are larger than other reported genomes in the genus Pleurotus. We propose that the selective amplification of long terminal repeat (LTR) retrotransposons increases the genome size of the genus Pleurotus, and may play a key role in driving their rapid species diversification. Molecular clock analyses of five Pleurotus species, namely PEE81, PEF12, P. tuoliensis, P. ostreatus and P. cf. floridanus suggest that the divergence estimates of the genus Pleurotus over time scales ranged from ∼4 to ∼38 million years ago (Mya), and PEE81 and PEF12 diverged at ∼13 Mya. The whole genome resequencing of 33 geographically diverse strains of P. eryngii var. eryngii and var. ferulae was then performed and the genome variation among and within these two populations were investigated. Comparative analyses of these two populations detected several candidate genes related to stress responses and DNA repair that are putatively involved in adaptation to the Gobi Desert environment. These findings offer insights into genome evolution of the genus Pleurotus and provide valuable genomic resources for King Oyster mushroom breeding.

Casein Gene Cluster in Camelids: Comparative Genome Analysis and New Findings on Haplotype Variability and Physical Mapping

The structure of casein genes has been fully understood in llamas, whereas in other camelids, this information is still incomplete. In fact, structure and polymorphisms have been identified in three (CSN1S1, αs1-CN; CSN2, β-CN; CSN3, κ-CN) out of four casein genes, whereas controversial information is available for the CSN1S2 (αs2-CN) in terms of structure and genetic diversity. Data from the genome analysis, whose assembly is available for feral camel, Bactrian, dromedary, and alpaca, can contribute to a better knowledge. However, a majority of the scaffolds available in GenBank are still unplaced, and the comparative annotation is often inaccurate or lacking.Therefore, the aims of this study are 1) to perform a comparative genome analysis and synthesize the literature data on camelids casein cluster; 2) to analyze the casein variability in two dromedary populations (Sudanese and Nigerian) using polymorphisms at CSN1S1 (c.150G > T), CSN2 (g.2126A > G), and CSN3 (g.1029T > C); and 3) to physically map the casein cluster in alpaca. Exon structures, gene and intergenic distances, large insertion/deletion events, SNPs, and microsatellites were annotated. In all camelids, the CSN1S2 consists of 17 exons, confirming the structure of llama CSN1S2 gene. The comparative analysis of the complete casein cluster (∼190kb) shows 12,818 polymorphisms. The most polymorphic gene is the CSN1S1 (99 SNPs in Bactrian vs. 248 in dromedary vs. 626 in alpaca). The less polymorphic is the CSN3 in the Bactrian (22 SNPs) and alpaca (301 SNPs), whereas it is the CSN1S2 in dromedary (79 SNPs). In the two investigated dromedary populations, the allele frequencies for the three markers are slightly different: the allele C at CSN1S1 is very rare in Nigerian (0.054) and Sudanese dromedaries (0.094), whereas the frequency of the allele G at CSN2 is almost inverted. Haplotype analysis evidenced GAC as the most frequent (0.288) and TGC as the rarest (0.005). The analysis of R-banding metaphases hybridized with specific probes mapped the casein genes on chromosome 2q21 in alpaca. These data deepen the information on the structure of the casein cluster in camelids and add knowledge on the cytogenetic map and haplotype variability.

Genomic Analyses Reveal Evidence of Independent Evolution, Demographic History, and Extreme Environment Adaptation of Tibetan Plateau Agaricus bisporus

Agaricus bisporus distributed in the Tibetan Plateau of China has high-stress resistance that is valuable for breeding improvements. However, its evolutionary history, specialization, and adaptation to the extreme Tibetan Plateau environment are largely unknown. Here, we performed de novo genome sequencing of a representative Tibetan Plateau wild strain ABM and comparative genomic analysis with the reported European strain H97 and H39. The assembled ABM genome was 30.4 Mb in size, and comprised 8,562 protein-coding genes. The ABM genome shared highly conserved syntenic blocks and a few inversions with H97 and H39. The phylogenetic tree constructed by 1,276 single-copy orthologous genes in nine fungal species showed that the Tibetan Plateau and European A. bisporus diverged ∼5.5 million years ago. Population genomic analysis using genome resequencing of 29 strains revealed that the Tibetan Plateau population underwent significant differentiation from the European and American populations and evolved independently, and the global climate changes critically shaped the demographic history of the Tibetan Plateau population. Moreover, we identified key genes that are related to the cell wall and membrane system, and the development and defense systems regulated A. bisporus adapting to the harsh Tibetan Plateau environment. These findings highlight the value of genomic data in assessing the evolution and adaptation of mushrooms and will enhance future genetic improvements of A. bisporus.

Genome-Wide Identification of Microsatellites and Transposable Elements in the Dromedary Camel Genome Using Whole-Genome Sequencing Data

Transposable elements (TEs) along with simple sequence repeats (SSRs) are prevalent in eukaryotic genome, especially in mammals. Repetitive sequences form approximately one-third of the camelid genomes, so study on this part of genome can be helpful in providing deeper information from the genome and its evolutionary path. Here, in order to improve our understanding regarding the camel genome architecture, the whole genome of the two dromedaries (Yazdi and Trodi camels) was sequenced. Totally, 92- and 84.3-Gb sequence data were obtained and assembled to 137,772 and 149,997 contigs with a N50 length of 54,626 and 54,031 bp in Yazdi and Trodi camels, respectively. Results showed that 30.58% of Yazdi camel genome and 30.50% of Trodi camel genome were covered by TEs. Contrary to the observed results in the genomes of cattle, sheep, horse, and pig, no endogenous retrovirus-K (ERVK) elements were found in the camel genome. Distribution pattern of DNA transposons in the genomes of dromedary, Bactrian, and cattle was similar in contrast with LINE, SINE, and long terminal repeat (LTR) families. Elements like RTE-BovB belonging to LINEs family in cattle and sheep genomes are dramatically higher than genome of dromedary. However, LINE1 (L1) and LINE2 (L2) elements cover higher percentage of LINE family in dromedary genome compared to genome of cattle. Also, 540,133 and 539,409 microsatellites were identified from the assembled contigs of Yazdi and Trodi dromedary camels, respectively. In both samples, di-(393,196) and tri-(65,313) nucleotide repeats contributed to about 42.5% of the microsatellites. The findings of the present study revealed that non-repetitive content of mammalian genomes is approximately similar. Results showed that 9.1 Mb (0.47% of whole assembled genome) of Iranian dromedary's genome length is made up of SSRs. Annotation of repetitive content of Iranian dromedary camel genome revealed that 9,068 and 11,544 genes contain different types of TEs and SSRs, respectively. SSR markers identified in the present study can be used as a valuable resource for genetic diversity investigations and marker-assisted selection (MAS) in camel-breeding programs.

Evolution, Systematics, and Biogeography of the Triatominae, Vectors of Chagas Disease

In this chapter, we review and update current knowledge about the evolution, systematics, and biogeography of the Triatominae (Hemiptera: Reduviidae)-true bugs that feed primarily on vertebrate blood. In the Americas, triatomines are the vectors of Trypanosoma cruzi, the etiological agent of Chagas disease. Despite declining incidence and prevalence, Chagas disease is still a major public health concern in Latin America. Triatomines occur also in the Old World, where vector-borne T. cruzi transmission has not been recorded. Triatomines evolved from predatory reduviid bugs, most likely in the New World, and diversified extensively across the Americas (including the Caribbean) and in parts of Asia and Oceania. Here, we first discuss our current understanding of how, how many times, and when the blood-feeding habit might have evolved among the Reduviidae. Then we present a summary of recent advances in the systematics of this diverse group of insects, with an emphasis on the contribution of molecular tools to the clarification of taxonomic controversies. Finally, and in the light of both up-to-date phylogenetic hypotheses and a thorough review of distribution records, we propose a global synthesis of the biogeography of the Triatominae. Over 130 triatomine species contribute to maintaining T. cruzi transmission among mammals (sometimes including humans) in almost every terrestrial ecoregion of the Americas. This means that Chagas disease will never be eradicated and underscores the fact that effective disease prevention will perforce require stronger, long-term vector control-surveillance systems.

Natural Killer Cell Receptor Genes in Camels: Another Mammalian Model

Due to production of special homodimeric heavy chain antibodies, somatic hypermutation of their T-cell receptor genes and unusually low diversity of their major histocompatibility complex genes, camels represent an important model for immunogenetic studies. Here, we analyzed genes encoding selected natural killer cell receptors with a special focus on genes encoding receptors for major histocompatibility complex (MHC) class I ligands in the two domestic camel species, Camelus dromedarius and Camelus bactrianus. Based on the dromedary genome assembly CamDro2, we characterized the genetic contents, organization, and variability of two complex genomic regions, the leukocyte receptor complex and the natural killer complex, along with the natural cytotoxicity receptor genes NCR1, NCR2, and NCR3. The genomic organization of the natural killer complex region of camels differs from cattle, the phylogenetically most closely related species. With its minimal set of KLR genes, it resembles this complex in the domestic pig. Similarly, the leukocyte receptor complex of camels is strikingly different from its cattle counterpart. With KIR pseudogenes and few LILR genes, it seems to be simpler than in the pig. The syntenies and protein sequences of the NCR1, NCR2, and NCR3 genes in the dromedary suggest that they could be human orthologues. However, only NCR1 and NCR2 have a structure of functional genes, while NCR3 appears to be a pseudogene. High sequence similarities between the two camel species as well as with the alpaca Vicugna pacos were observed. The polymorphism in all genes analyzed seems to be generally low, similar to the rest of the camel genomes. This first report on natural killer cell receptor genes in camelids adds new data to our understanding of specificities of the camel immune system and its functions, extends our genetic knowledge of the innate immune variation in dromedaries and Bactrian camels, and contributes to studies of natural killer cell receptors evolution in mammals.

Chromosome-Level Alpaca Reference Genome VicPac3.1 Improves Genomic Insight Into the Biology of New World Camelids

The development of high-quality chromosomally assigned reference genomes constitutes a key feature for understanding genome architecture of a species and is critical for the discovery of the genetic blueprints of traits of biological significance. South American camelids serve people in extreme environments and are important fiber and companion animals worldwide. Despite this, the alpaca reference genome lags far behind those available for other domestic species. Here we produced a chromosome-level improved reference assembly for the alpaca genome using the DNA of the same female Huacaya alpaca as in previous assemblies. We generated 190X Illumina short-read, 8X Pacific Biosciences long-read and 60X Dovetail Chicago^® chromatin interaction scaffolding data for the assembly, used testis and skin RNAseq data for annotation, and cytogenetic map data for chromosomal assignments. The new assembly VicPac3.1 contains 90% of the alpaca genome in just 103 scaffolds and 76% of all scaffolds are mapped to the 36 pairs of the alpaca autosomes and the X chromosome. Preliminary annotation of the assembly predicted 22,462 coding genes and 29,337 isoforms. Comparative analysis of selected regions of the alpaca genome, such as the major histocompatibility complex (MHC), the region involved in the Minute Chromosome Syndrome (MCS) and candidate genes for high-altitude adaptations, reveal unique features of the alpaca genome. The alpaca reference genome VicPac3.1 presents a significant improvement in completeness, contiguity and accuracy over VicPac2 and is an important tool for the advancement of genomics research in all New World camelids.

Large-scale ruminant genome sequencing provides insights into their evolution and distinct traits

The ruminants are one of the most successful mammalian lineages, exhibiting morphological and habitat diversity and containing several key livestock species. To better understand their evolution, we generated and analyzed de novo assembled genomes of 44 ruminant species, representing all six Ruminantia families. We used these genomes to create a time-calibrated phylogeny to resolve topological controversies, overcoming the challenges of incomplete lineage sorting. Population dynamic analyses show that population declines commenced between 100,000 and 50,000 years ago, which is concomitant with expansion in human populations. We also reveal genes and regulatory elements that possibly contribute to the evolution of the digestive system, cranial appendages, immune system, metabolism, body size, cursorial locomotion, and dentition of the ruminants.

Cellular and Molecular Adaptation of Arabian Camel to Heat Stress

To cope with the extreme heat stress and drought of the desert, the Arabian camel (Camelus dromedarius) has developed exceptional physiological and biochemical particularities. Previous reports focused mainly on the physiological features of Arabian camel and neglected its cellular and molecular characteristics. Heat shock proteins are suggested to play a key role in the protein homeostasis and thermotolerance. Therefore, we aim by this review to elucidate the implication of camel HSPs in its physiological adaptation to heat stress and compare them with HSPs in related mammalian species. Correlation of these molecules to the adaptive mechanisms in camel is of special importance to expand our understanding of the overall camel physiology and homeostasis.

Beyond the Big Five: Investigating Myostatin Structure, Polymorphism and Expression in Camelus dromedarius

Myostatin, a negative regulator of skeletal muscle mass in animals, has been shown to play a role in determining muscular hypertrophy in several livestock species, and a high degree of polymorphism has been previously reported for this gene in humans and cattle. In this study, we provide a characterization of the myostatin gene in the dromedary (Camelus dromedarius) at the genomic, transcript and protein level. The gene was found to share high structural and sequence similarity with other mammals, notably Old World camelids. 3D modeling highlighted several non-conservative SNP variants compared to the bovine, as well as putative functional variants involved in the stability of the myostatin dimer. NGS data for nine dromedaries from various countries revealed 66 novel SNPs, all of them falling either upstream or downstream the coding region. The analysis also confirmed the presence of three previously described SNPs in intron 1, predicted here to alter both splicing and transcription factor binding sites (TFBS), thus possibly impacting myostatin processing and/or regulation. Several putative TFBS were identified in the myostatin upstream region, some of them belonging to the myogenic regulatory factor family. Patterns of SNP distribution across countries, as suggested by Bayesian clustering of the nine dromedaries using the 69 SNPs, pointed to weak geographic differentiation, in line with known recurrent gene flow at ancient trading centers along caravan routes. Myostatin expression was investigated in a set of 8 skeletal muscles, both at transcript and protein level, via Digital Droplet PCR and Western Blotting, respectively. No significant differences were observed at the transcript level, while, at the protein level, the only significant differences concerned the promyostatin dimer (75 kDa), in four pair-wise comparisons, all involving the tensor fasciae latae muscle. Beside the mentioned band at 75 kDa, additional bands were observed at around 40 and 25 kDa, corresponding to the promyostatin monomer and the active C-terminal myostatin dimer, respectively. Their weaker intensity suggests that the unprocessed myostatin dimers could act as important reservoirs of slowly available myostatin forms. Under this assumption, the sequential cleavage steps may contribute additional layers of control within an already complex regulatory framework.

Comparative Analysis of the TRB Locus in the Camelus Genus

T cells can be separated into two major subsets based on the heterodimer that forms their T cell receptors. αβ T cells have receptors consisting of α and β chains, while γδ T cells are composed of γ and δ chains. αβ T cells play an essential role within the adaptive immune responses against pathogens. The recent genomic characterization of the Camelus dromedarius T cell receptor β (TRB) locus has allowed us to infer the structure of this locus from the draft genome sequences of its wild and domestic Bactrian congeners, Camelus ferus and Camelus bactrianus. The general structural organization of the wild and domestic Bactrian TRB locus is similar to that of the dromedary, with a pool of TRBV genes positioned at the 5′ end of D-J-C clusters, followed by a single TRBV gene located at the 3′ end with an inverted transcriptional orientation. Despite the fragmented nature of the assemblies, comparative genomics reveals the existence of a perfect co-linearity between the three Old World camel TRB genomic sequences, which enables the transfer of information from one sequence to another and the filling of gaps in the genomic sequences. A virtual camelid TRB locus is hypothesized with the presence of 33 TRBV genes distributed in 26 subgroups. Likewise, in the artiodactyl species, three in-tandem D-J-C clusters, each composed of one TRBD gene, six or seven TRBJ genes, and one TRBC gene, are placed at the 3′ end of the locus. As reported in the ruminant species, a group of four functional TRY genes at the 5′ end and only one gene at the 3′ end, complete the camelid TRB locus. Although the gene content is similar, differences are observed in the TRBV functional repertoire, and genes that are functional in one species are pseudogenes in the other species. Hence, variations in the functional repertoire between dromedary, wild and domestic Bactrian camels, rather than differences in the gene content, may represent the molecular basis explaining the disparity in the TRB repertoire between the Camelus species. Finally, our data contribute to the knowledge about the evolutionary history of Old World camelids.

A First Y-Chromosomal Haplotype Network to Investigate Male-Driven Population Dynamics in Domestic and Wild Bactrian Camels

Polymorphic markers on the male-specific part of the Y chromosome (MSY) provide useful information for tracking male genealogies. While maternal lineages are well studied in Old World camelids using mitochondrial DNA, the lack of a Y-chromosomal reference sequence hampers the analysis of male-driven demographics. Recently, a shotgun assembly of the horse MSY was generated based on short read next generation sequencing data. The haplotype network resulting from single copy MSY variants using the assembly as a reference revealed sufficient resolution to trace individual male lines in this species. In a similar approach we generated a 3.8 Mbp sized assembly of the MSY of Camelus bactrianus. The camel MSY assembly was used as a reference for variant calling using short read data from eight Old World camelid individuals. Based on 596 single nucleotide variants we revealed a Y-phylogenetic network with seven haplotypes. Wild and domestic Bactrian camels were clearly separated into two different haplogroups with an estimated divergence time of 26,999 ± 2,268 years. Unexpectedly, one wild camel clustered into the domestic Bactrian camels' haplogroup. The observation of a domestic paternal lineage within the wild camel population is concerning in view of the importance to conserve the genetic integrity of these highly endangered species in their natural habitat.

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.

Evaluation of SNP Genotyping in Alpacas Using the Bovine HD Genotyping Beadchip

Alpacas are one of four South American Camelid species living in the highlands of the Andes. Production of alpaca fiber contributes to the economy of the region and the livelihood of many rural families. Fiber quantity and quality are important and in need of a modern breeding program based on genomic selection to accelerate genetic gain. To achieve this is necessary to discover enough molecular markers, single nucleotide polymorphisms (SNPs) in particular, to provide genome coverage and facilitate genome wide association studies to fiber production characteristics. The aim of this study was to discover alpaca SNPs by genotyping forty alpaca DNA samples using the BovineHD Genotyping Beadchip. Data analysis was performed with GenomeStudio (Illumina) software. Because different filters and thresholds are reported in the literature we investigated the effects of no-call threshold (≥0.05, ≥0.15, and ≥0.25) and call frequency (≥0.9 and =1.0) in identifying positive SNPs. Average GC Scores, calculated as the average of the 10% and 50% GenCall scores for each SNP (≥0.70) and the GenTrain score ≥ 0.25 parameters were applied to all comparisons. SNPs with minor allele frequency (MAF) ≥ 0.05 or ≥ 0.01 were retained. Since detection of SNPs is based on the stable binding of oligonucleotide probes to the target DNA immediately adjacent to the variant nucleotide, all positive SNP flanking sequences showing perfect alignments between the bovine and alpaca genomes for the first 21 or 26 nucleotides flanking the variant nucleotide at either side were selected. Only SNPs localized in one scaffold were assumed unique. Unique SNPs identified in both reference genomes were kept and mapped on the Vicugna_pacos 2.0.2 genome. The effects of the no-call threshold ≥ 0.25, call frequency = 1 and average GC ≥ 0.7 were meaningful and identified 6756 SNPs of which 400 were unique and polymorphic (MAF ≥ 0.01). Assignment to alpaca chromosomes was possible for 292 SNPs. Likewise, 209 SNPs were localized in 202 alpaca gene loci and 29 of these share the same loci with the dromedary. Interestingly, 69 of 400 alpaca SNPs have 100% similarity with dromedary.

Comparative FISH-Mapping of MC1R, ASIP, and TYRP1 in New and Old World Camelids and Association Analysis With Coat Color Phenotypes in the Dromedary (Camelus dromedarius)

Melanocortin 1 receptor (MC1R), the agouti signaling protein (ASIP), and tyrosinase related protein 1 (TYRP1) are among the major regulators of pigmentation in mammals. Recently, MC1R and ASIP sequence variants were associated with white and black/dark brown coat colors, respectively, in the dromedary. Here we confirmed this association by independent sequencing and mutation discovery of MC1R and ASIP coding regions and by TaqMan genotyping in 188 dromedaries from Saudi Arabia and United States, including 38 black, 53 white, and 97 beige/brown/red animals. We showed that heterozygosity for a missense mutation c.901C > T in MC1R is sufficient for the white coat color suggesting a possible dominant negative effect. Likewise, we confirmed that the majority of black dromedaries were homozygous for a frameshift mutation in ASIP exon 2, except for 4 animals, which were heterozygous. In search for additional mutations underlying the black color, we identified another frameshift mutation in ASIP exon 4 and 6 new variants in MC1R including a significantly associated SNP in 3'UTR. In pursuit of sequence variants that may modify dromedary wild-type color from dark-reddish brown to light beige, we identified 4 SNPs and one insertion in TYRP1 non-coding regions. However, none of these were associated with variations in wild-type colors. Finally, the three genes were cytogenetically mapped in New World (alpaca) and Old World (dromedary and Bactrian camel) camelids. The MC1R was assigned to chr21, ASIP to chr19 and TYRP1 to chr4 in all 3 species confirming extensive conservation of camelid karyotypes. Notably, while the locations of ASIP and TYRP1 were in agreement with human-camelid comparative map, mapping MC1R identified a new evolutionary conserved synteny segment between camelid chromosome 21 and HSA16. The findings contribute to coat color genomics and the development of molecular tests in camelids and toward the chromosome level reference assemblies of camelid genomes.

Key genes differential expressions and pathway involved in salt and water-deprivation stresses for renal cortex in camel

Background

Camels possess the characteristics of salt- and drought-resistances, due to the long-time adaption to the living environment in desert. The camel resistance research on transcriptome is rare and deficient, especially reabsorption in renal cortex. Non-coding RNAs are normally considered as the RNA molecules that are not translated into proteins, their current roles remain mostly in regulation of information flux from DNA to protein, further on normal life activities and diseases. In order to reveal the mysterious veil of the post-transcriptional regulation of ncRNAs in renal cortex for the first time as far as we know, we designed and carried out the experiment of salt stress and water-deprivation stress in camel.

Results

By means of RNA-seq in renal cortex of Alxa Bactrian Camel (Camelus bactrianus), we identified certain significantly differential RNAs, including 4 novel lncRNAs, 11 miRNAs and 13 mRNAs under salt stress, 0 lncRNAs, 18 miRNAs and 14 mRNAs under water-deprivation stress. By data analysis, the response pathway of post-transcriptional regulation concerning salt and water-deprivation stresses was put forward, involving preventing sodium from entering the cell, purifying of water and compensating neutral amino acids by miR-193b, miR-542-5p interaction with SLC6A19 mRNA.

Conclusion

Based on the resistance-related lncRNAs, miRNAs, and mRNAs, we proposed the post-transcriptional regulation pathway to explain how camels respond to salt and water-deprivation stresses in the ncRNAs regulation level of renal cortex for the first time, thus hoping to provide a theoretical basis for therapy of disease that is similar to high blood pressure in humans.

Electronic supplementary material

The online version of this article (10.1186/s12867-019-0129-8) contains supplementary material, which is available to authorized users.

De novo comparative transcriptome analysis of a rare cicada, with identification of candidate genes related to adaptation to a novel host plant and drier habitats

Background

Although the importance of host plant chemistry in plant–insect interactions is widely recognized, our understanding about the genetic basis underlying the relationship between changes in midgut proteins and adaptation of plant-feeding insects to novel host plants and habitats is very limited. To address this knowledge gap, the transcriptional profiles of midguts among three populations of the cicada Subpsaltria yangi Chen were compared. Among which, the Hancheng (HC) and Fengxiang (FX) populations occurring in the Loess Plateau feed on Ziziphus jujuba Mill. var. spinosa (Bunge) Hu ex H. F. Chow, while the population occurring in a much drier habitat in the Helan (HL) Mountains is locally specialized on a chemically divergent plant, Ephedra lepidosperma C. Y. Cheng.

Results

Based on comparative analysis, 1826 (HL vs HC) differentially expressed genes (DEGs) and 723 DEGs (HL vs FX) were identified between the populations utilizing different host plants, including 20, 36, 2, 5 and 2 genes related to digestion, detoxification, oxidation-reduction, stress response and water-deprivation response, respectively, and 35 genes presumably associated with osmoregulation. However, only 183 DEGs were identified between the HC and FX populations, including two genes related to detoxification, two genes related to stress response, and one gene presumably associated with osmoregulation. These results suggest that the weakest expression differences were between the populations utilizing the same host plant and occurring in the closest habitats, which may help explain the metabolic mechanism of adaptation in S. yangi populations to novel host plants and new niches.

Conclusions

The observed differences in gene expression among S. yangi populations are consistent with the hypothesis that the host plant shift and habitat adaptation in the HL population was facilitated by differential regulation of genes related to digestion, detoxification, oxidation-reduction, stress response, water-deprivation response and osmoregulation. The results may inform future studies on the molecular mechanisms underlying the relationship between changes in midgut proteins and adaptation of herbivorous insects to novel host plants and new niches.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5547-y) contains supplementary material, which is available to authorized users.

From Desert to Medicine: A Review of Camel Genomics and Therapeutic Products

Camels have an important role in the lives of human beings, especially in arid regions, due to their multipurpose role and unique ability to adapt to harsh conditions. In spite of its enormous economic, cultural, and biological importance, the camel genome has not been widely studied. The size of camel genome is roughly 2.38 GB, containing over 20,000 genes. The unusual genetic makeup of the camel is the main reason behind its ability to survive under extreme environmental conditions. The camel genome harbors several unique variations which are being investigated for the treatment of several disorders. Various natural products from camels have also been tested and prescribed as adjunct therapy to control the progression of ailments. Interestingly, the camel employs unique immunological and molecular mechanisms against pathogenic agents and pathological conditions. Here, we broadly review camel classification, distribution and breed as well as recent progress in the determination of the camel genome, its size, genetic distribution, response to various physiological conditions, immunogenetics and the medicinal potential of camel gene products.

Cdrom Archive: A Gateway to Study Camel Phenotypes

Camels are livestock that exhibit unique morphological, biochemical, and behavioral traits, which arose by natural and artificial selection. Investigating the molecular basis of camel traits has been limited by: (1) the absence of a comprehensive record of morphological trait variation (e.g., diseases) and the associated mode of inheritance, (2) the lack of extended pedigrees of specific trait(s), and (3) the long reproductive cycle of the camel, which makes the cost of establishing and maintaining a breeding colony (i.e., monitoring crosses) prohibitively high. Overcoming these challenges requires (1) detailed documentation of phenotypes/genetic diseases and their likely mode of inheritance (and collection of related DNA samples), (2) conducting association studies to identify phenotypes/genetic diseases causing genetic variants (instead of classical linkage analysis, which requires extended pedigrees), and (3) validating likely causative variants by screening a large number of camel samples from different populations. We attempt to address these issues by establishing a systematic way of collecting camel DNA samples, and associated phenotypic information, which we call the "Cdrom Archive." Here, we outline the process of building this archive to introduce it to other camel researchers (as an example). Additionally, we discuss the use of this archive to study the phenotypic traits of Arabian Peninsula camel breeds (the "Mezayen" camels). Using the Cdrom Archive, we report variable phenotypic traits related to the coat (color, length, and texture), ear and tail lengths, along with other morphological measurements.

A Near Chromosome Assembly of the Dromedary Camel Genome

The dromedary camel is an economically and socially important species of livestock in many parts of the world, being used for transport and the production of milk and meat. Much like cattle and horses, the camel may be found in industrial farming conditions as well as used in sporting. Camel racing is a multi-million dollar industry, with some specimens being valued at upward of 9.5 million USD. Despite its apparent value to humans, the dromedary camel is a neglected species in genomics. While cattle and other domesticated species have had much attention in terms of genome assembly, the camel has only been assembled to scaffold level, which does not give a clear indication of the order or chromosomal location of sequenced fragments. In this study, the Reference Assistant Chromosome Assembly (RACA) algorithm was implemented to use read-pair information of camel scaffolds, aligned with the cattle and human genomes in order to organize and orient these scaffolds in a near-chromosome level assembly. This method generated 72 large size fragments (N50 54.36 Mb). These predicted chromosome fragments (PCFs) were then compared with comparative maps of camel and cytogenetic map of alpaca chromosomes, allowing us to further upgrade the assembly. This dromedary camel assembly will be an invaluable tool to verify future camel assemblies generated with chromatin conformation or/and long read technologies. This study provides the first near-chromosome assembly of the dromedary camel, thus adding this economically important species to a growing pool of knowledge regarding the genome structure of domesticated livestock.

Quality and quantity of dromedary camel DNA sampled from whole-blood, saliva, and tail-hair

Camels are livestock with unique adaptations to hot-arid regions. To effectively study camel traits, a biobank of camel DNA specimens with associated biological information is needed. We examined whole-blood, saliva (buccal swabs), and tail-hair follicle samples to determine which is the best source for establishing a DNA biobank. We inspected five amounts of each of whole-blood, buccal swabs, and tail-hair follicles in nine camels, both qualitatively via gel electrophoresis and quantitatively using a NanoDrop spectrophotometer. We also tested the effects of long term-storage on the quality and quantity of DNA, and measured the rate of degradation, by analyzing three buccal swab samples and 30 tail-hair follicles over a period of nine months. Good quality DNA, in the form of visible large size DNA bands, was extracted from all three sources, for all five amounts. The five volumes of whole-blood samples (20–100μl) provided ~0.4–3.6 μg, the five quantities of buccal swabs (1–5) produced ~0.1–12 μg, while the five amounts of tail-hair follicles (10–50) resulted in ~0.7–25 μg. No differences in the rate of degradation of buccal swab and tail-hair follicle DNA were detected, but there was clearly greater deterioration in the quality of DNA extracted from buccal swabs when compared to tail-hair follicles. We recommend using tail-hair samples for camel DNA biobanking, because it resulted in both an adequate quality and quantity of DNA, along with its ease of collection, transportation, and storage. Compared to its success in studies of other domesticated animals, we anticipate that using ~50 tail-hair follicles will provide sufficient DNA for sequencing or SNP genotyping.

Monomeric Camelus dromedarius GSTM1 at low pH is structurally more thermostable than its native dimeric form

Glutathione S‒transferases (GSTs) are multifunctional enzymes that play an important role in detoxification, cellular signalling, and the stress response. Camelus dromedarius is well-adapted to survive in extreme desert climate and it has GSTs, for which limited information is available. This study investigated the structure-function and thermodynamic properties of a mu-class camel GST (CdGSTM1) at different pH. Recombinant CdGSTM1 (25.7 kDa) was expressed in E. coli and purified to homogeneity. Dimeric CdGSTM1 dissociated into stable but inactive monomeric subunits at low pH. Conformational and thermodynamic changes during the thermal unfolding pathway of dimeric and monomeric CdGSTM1 were characterised via a thermal shift assay and dynamic multimode spectroscopy (DMS). The thermal shift assay based on intrinsic tryptophan fluorescence revealed that CdGSTM1 underwent a two-state unfolding pathway at pH 1.0-10.0. Its Tm value varied with varying pH. Another orthogonal technique based on far-UV CD also exhibited two-state unfolding in the dimeric and monomeric states. Generally, proteins tend to lose structural integrity and stability at low pH; however, monomeric CdGSTM1 at pH 2.0 was thermally more stable and unfolded with lower van't Hoff enthalpy. The present findings provide essential information regarding the structural, functional, and thermodynamic properties of CdGSTM1 at pH 1.0-10.0.

Genetic variants analysis of three dromedary camels using whole genome sequencing data

Whole genome wide identification and annotation of genetic variations in camels is in its first steps. The aim of this study was the identification of genome wide variants, functional annotations of them and enrichment analysis of affected genes using whole genome sequencing data of three dromedary camels. The genomes of two Iranian female dromedary camels that mostly used to produce meat and milk were sequenced to 41.9-fold and 38.6-fold coverage. A total of 4,727,238 single-nucleotide polymorphisms (SNPs) and 692,908 indels (insertions and deletions) were found by mapping raw reads to the dromedary reference assembly (GenBank Accession: GCA_000767585.1). In-silico functional annotation of the discovered variants in under study samples revealed that most SNPs (2,305,738; 48.78%) and indels (339,756; 49.03%) were located in intergenic regions. A comparison of the identified SNPs with those of the African camel (BioProject Accession: PRJNA269274) indicated that they had 993,474 SNPs in common. We found 15,168 non-synonymous SNPs in the shared variants of the three camels that could affect gene function and protein structure. Obtained results revealed that there were 7085, 6271 and 4688 non-synonymous SNPs among the 3436, 3058 and 2882 genes in the specific gene sets of Yazd dromedary, Trod dromedary and African dromedary, respectively. The list of genes predicted to be affected by non-synonymous variants in different individuals was subjected to gene ontology (GO) enrichment analysis.

Molecular footprints of inshore aquatic adaptation in Indo-Pacific humpback dolphin (Sousa chinensis)

The Indo-Pacific humpback dolphin, Sousa chinensis, being a member of cetaceans, had fully adapted to inshore waters. As a threatened marine mammal, little molecular information available for understanding the genetic basis of ecological adaptation. We firstly sequenced and obtained the draft genome map of S. chinensis. Phylogenetic analysis in this study, based on the single copy orthologous genes of the draft genome, is consistent with traditional phylogenetic classification. The comparative genomic analysis indicated that S. chinensis had 494 species-specific gene families, which involved immune, DNA repair and sensory systems associated with the potential adaption mechanism. We also identified the expansion and positive selection genes in S. chinensis lineage to investigate the potential adaptation mechanism. Our study provided the potential insight into the molecular bases of ecological adaptation in Indo-Pacific humpback dolphin and will be also valuable for future understanding the ecological adaptation and evolution of cetaceans at the genomic level.

An association between differential expression and genetic divergence in the Patagonian olive mouse (Abrothrix olivacea)

Recent molecular studies have found striking differences between desert-adapted species and model mammals regarding water conservation. In particular, aquaporin 4, a classical gene involved in water regulation of model species, is absent or not expressed in the kidneys of desert-adapted species. To further understand the molecular response to water availability, we studied the Patagonian olive mouse Abrothrix olivacea, a species with an unusually broad ecological tolerance that exhibits a great urine concentration capability. The species is able to occupy both the arid Patagonian steppe and the Valdivian and Magellanic forests. We sampled 95 olive mouse specimens from four localities (two in the steppe and two in the forests) and analysed both phenotypic variables and transcriptomic data to investigate the response of this species to the contrasting environmental conditions. The relative size of the kidney and the ratio of urine to plasma concentrations were, as expected, negatively correlated with annual rainfall. Expression analyses uncovered nearly 3,000 genes that were differentially expressed between steppe and forest samples and indicated that this species resorts to the "classical" gene pathways for water regulation. Differential expression across biomes also involves genes that involved in immune and detoxification functions. Overall, genes that were differentially expressed showed a slight tendency to be more divergent and to display an excess of intermediate allele frequencies, relative to the remaining loci. Our results indicate that both differential expression in pathways involved in water conservation and geographical allelic variation are important in the occupation of contrasting habitats by the Patagonian olive mouse.

Immunohistochemical expression of apoptosis-related biomarkers in normal tissues of camel (Camelus dromedarius): A survey in a desert-dwelling mammalian model

Programmed cell death is a fundamental event that takes place during organ development and plays an important role in cellular homeostasis. Since various body organs of the camel are under high ecological and physiological stress during food and water deprivation, desiccation, and the long exposure to solar radiation in these desert nomads, we aimed to examine the immunohistochemical expression of apoptosis-related biomarkers in some of its normal body organs to illustrate a basic track for further pathological investigation. Regarding apoptosis, the present study has revealed that the higher expression of cleaved caspase-9 (CC9) [initiator of the intrinsic pathway] and CC3 (effector caspase), and the scanty expression of CC8 (initiator of the extrinsic pathway), highlight the role of the caspase-dependent, intrinsic apoptotic pathway particularly in the intestines and lymphoid organs. The apoptosis- inducing factor (AIF)-immunoexpression was completely missing in the cell nuclei of the examined tissues, indicating the absence of the caspase-independent pathway. The nuclear overexpression of the phospho-histone H2AX (γ H2AX) and the occasional expression of single-stranded DNA, particularly among the CNS neurons, suggest an efficient, protective DNA-repair mechanism in such cells. Thus, despite efficient anti-apoptotic mechanisms intrinsic apoptotic pathways exists in brain, intestine and lymph organs of adult desert camels.

Genetic characterization of the oxytocin-neurophysin I gene (OXT) and its regulatory regions analysis in domestic Old and New World camelids

Oxytocin is a neurohypophysial peptide linked to a wide range of biological functions, including milk ejection, temperament and reproduction. Aims of the present study were a) the characterization of the OXT (Oxytocin-neurophysin I) gene and its regulatory regions in Old and New world camelids; b) the investigation of the genetic diversity and the discovery of markers potentially affecting the gene regulation. On average, the gene extends over 814 bp, ranging between 825 bp in dromedary, 811 bp in Bactrian and 810 bp in llama and alpaca. Such difference in size is due to a duplication event of 21 bp in dromedary. The main regulatory elements, including the composite hormone response elements (CHREs), were identified in the promoter, whereas the presence of mature microRNAs binding sequences in the 3'UTR improves the knowledge on the factors putatively involved in the OXT gene regulation, although their specific biological effect needs to be still elucidated. The sequencing of genomic DNA allowed the identification of 17 intraspecific polymorphisms and 69 nucleotide differences among the four species. One of these (MF464535:g.622C>G) is responsible, in alpaca, for the loss of a consensus sequence for the transcription factor SP1. Furthermore, the same SNP falls within a CpG island and it creates a new methylation site, thus opening future possibilities of investigation to verify the influence of the novel allelic variant in the OXT gene regulation. A PCR-RFLP method was setup for the genotyping and the frequency of the allele C was 0.93 in a population of 71 alpacas. The obtained data clarify the structure of OXT gene in domestic camelids and add knowledge to the genetic variability of a genomic region, which has received little investigation so far. These findings open the opportunity for new investigations, including association studies with productive and reproductive traits.

Body mass-specific Na+-K+-ATPase activity in the medullary thick ascending limb: implications for species-dependent urine concentrating mechanisms

In general, the mammalian whole body mass-specific metabolic rate correlates positively with maximal urine concentration (U_max) irrespective of whether or not the species have adapted to arid or mesic habitat. Accordingly, we hypothesized that the thick ascending limb (TAL) of a rodent with markedly higher whole body mass-specific metabolism than rat exhibits a substantially higher TAL metabolic rate as estimated by Na⁺-K⁺-ATPase activity and Na⁺-K⁺-ATPase α1-gene and protein expression. The kangaroo rat inner stripe of the outer medulla exhibits significantly higher mean Na⁺-K⁺-ATPase activity (~70%) compared with two rat strains (Sprague-Dawley and Munich-Wistar), extending prior studies showing rat activity exceeds rabbit. Furthermore, higher expression of Na⁺-K⁺-ATPase α1-protein (~4- to 6-fold) and mRNA (~13-fold) and higher TAL mitochondrial volume density (~20%) occur in the kangaroo rat compared with both rat strains. Rat TAL Na⁺-K⁺-ATPase α1-protein expression is relatively unaffected by body hydration status or, shown previously, by dietary Na⁺, arguing against confounding effects from two unavoidably dissimilar diets: grain-based diet without water (kangaroo rat) or grain-based diet with water (rat). We conclude that higher TAL Na⁺-K⁺-ATPase activity contributes to relationships between whole body mass-specific metabolic rate and high U_max. More vigorous TAL Na⁺-K⁺-ATPase activity in kangaroo rat than rat may contribute to its steeper Na⁺ and urea axial concentration gradients, adding support to a revised model of the urine concentrating mechanism, which hypothesizes a leading role for vigorous active transport of NaCl, rather than countercurrent multiplication, in generating the outer medullary axial osmotic gradient.

Repeat associated mechanisms of genome evolution and function revealed by the Mus caroli and Mus pahari genomes

Understanding the mechanisms driving lineage-specific evolution in both primates and rodents has been hindered by the lack of sister clades with a similar phylogenetic structure having high-quality genome assemblies. Here, we have created chromosome-level assemblies of the Mus caroli and Mus pahari genomes. Together with the Mus musculus and Rattus norvegicus genomes, this set of rodent genomes is similar in divergence times to the Hominidae (human-chimpanzee-gorilla-orangutan). By comparing the evolutionary dynamics between the Muridae and Hominidae, we identified punctate events of chromosome reshuffling that shaped the ancestral karyotype of Mus musculus and Mus caroli between 3 and 6 million yr ago, but that are absent in the Hominidae. Hominidae show between four- and sevenfold lower rates of nucleotide change and feature turnover in both neutral and functional sequences, suggesting an underlying coherence to the Muridae acceleration. Our system of matched, high-quality genome assemblies revealed how specific classes of repeats can play lineage-specific roles in related species. Recent LINE activity has remodeled protein-coding loci to a greater extent across the Muridae than the Hominidae, with functional consequences at the species level such as reproductive isolation. Furthermore, we charted a Muridae-specific retrotransposon expansion at unprecedented resolution, revealing how a single nucleotide mutation transformed a specific SINE element into an active CTCF binding site carrier specifically in Mus caroli, which resulted in thousands of novel, species-specific CTCF binding sites. Our results show that the comparison of matched phylogenetic sets of genomes will be an increasingly powerful strategy for understanding mammalian biology.