Museum genomics approach to study the taxonomy and evolution of Woolly-necked storks using historic specimens.

The accessibility of genomic tools in evolutionary biology has allowed for a thorough exploration of various evolutionary processes associated with adaptation and speciation. However, genomic studies in natural systems present numerous challenges, reflecting the inherent complexities of studying organisms in their native habitats. The utilization of museum specimens for genomics research has received increased attention in recent times, facilitated by advancements in ancient DNA techniques. In this study, we have utilized a museum genomics approach to analyze historic specimens of Woolly-necked storks (Ciconia spp.) and examine their genetic composition and taxonomic status and explore the evolutionary and adaptive trajectories of populations over the years. The Woolly-necked storks are distributed in Asia and Africa with a taxonomic classification that has been a matter of ambiguity. Asian and African Woollynecks were recently recognized as different species based on their morphological differences; however, their genomic validation was lacking. In this study, we have used ∼70-year-old museum samples for whole-genome population-scale sequencing. Our study has revealed that Asian and African Woollynecks are genetically distinct, consistent with the current taxonomic classification based on morphological features. However, we also found a high genetic divergence between the Asian subspecies Ciconia episcopus neglecta and Ciconia episcopus episcopus, suggesting this classification requires a detailed examination to explore processes of ongoing speciation. Because taxonomic classification directly impacts conservation efforts, and there is evidence of declining populations of Asian Woollynecks in Southeast Asia, our results highlight that population-scale studies are urgent to determine the genetic, ecological, and phylogenetic diversity of these birds.

Genomic Variation, Population History, and Long-Term Genetic Adaptation to High Altitudes in Tibetan Partridge (Perdix hodgsoniae).

Species residing across elevational gradients display adaptations in response to environmental changes such as oxygen availability, ultraviolet radiation, and temperature. Here, we study genomic variation, gene expression, and long-term adaptation in Tibetan Partridge (Perdix hodgsoniae) populations residing across the elevational gradient of the Tibetan Plateau. We generated a high-quality draft genome and used it to carry out downstream population genomic and transcriptomic analysis. The P. hodgsoniae populations residing across various elevations were genetically distinct, and their phylogenetic clustering was consistent with their geographic distribution. We identified possible evidence of gene flow between populations residing in <3,000 and >4,200 m elevation that is consistent with known habitat expansion of high-altitude populations of P. hodgsoniae to a lower elevation. We identified a 60 kb haplotype encompassing the Estrogen Receptor 1 (ESR1) gene, showing strong genetic divergence between populations of P. hodgsoniae. We identified six single nucleotide polymorphisms within the ESR1 gene fixed for derived alleles in high-altitude populations that are strongly conserved across vertebrates. We also compared blood transcriptome profiles and identified differentially expressed genes (such as GAPDH, LDHA, and ALDOC) that correlated with differences in altitude among populations of P. hodgsoniae. These candidate genes from population genomics and transcriptomics analysis were enriched for neutrophil degranulation and glycolysis pathways, which are known to respond to hypoxia and hence may contribute to long-term adaptation to high altitudes in P. hodgsoniae. Our results highlight Tibetan Partridges as a useful model to study molecular mechanisms underlying long-term adaptation to high altitudes.

Bimodal evolution of the killer cell Ig-like receptor (KIR) family in New World primates.

The immunoglobulin-like receptor (KIR) gene family in New World primates (Platyrrhini) has been characterized only in the owl monkey (Aotus sp.). To gain a better understanding of the KIR system in Platyrrhini, we analyzed a KIR haplotype in Ateles geoffroyi, and sequenced KIR complementary DNAs (cDNAs) from other three Atelidae species, Ateles hybridus, Ateles belzebuth, and Lagothrix lagotricha. Atelidae expressed a variable set of activating and inhibitory KIRs that diversified independently from their Catarrhini counterparts. They had a unique mechanism to generate activating receptors from inhibitory ones, involving a single nucleotide deletion in exon 7 and a change in the donor splice site of intron 7. The A. geoffroyi haplotype contained at least six gene models including a pseudogene, two coding inhibitory receptors, and three coding activating receptors. The centromeric region was in a tail-to-tail orientation with respect to the telomeric region. The owl monkey KIR haplotype shared this organization, and in phylogenetic trees, the centromeric genes clustered together with those of A. geoffroyi, whereas their telomeric genes clustered independently. KIR cDNAs from the other Atelidae species conformed to this pattern. Signatures of positive selection were found in residues predicted to interact with the major histocompatibility complex. Such signatures, however, primarily explained variability between paralogous genes but not between alleles in a locus. Atelidae, therefore, has expanded the KIR family in a bimodal fashion, where an inverted centromeric region has remained relatively conserved and the telomeric region has diversified by a rapid process of gene duplication and divergence, likely favored by positive selection for ligand binding.

Evolutionary patterns of killer cell Ig-like receptor genes in Old World monkeys.

Killer cell Ig-like receptors (KIRs) modulate the cytotoxic effects of Natural Killer cells. KIR genes are encoded in the Leucocyte Receptor Complex and are characterized by their high haplotypic diversity and polymorphism. The KIR system has been studied in only three species of Old World monkeys, the rhesus macaque, the cynomolgus macaque, and the sabaeus monkey, displaying a complexity rivaling that of hominids (human and apes). Here we analyzed bacterial artificial chromosome draft sequences spanning the KIR haplotype of three other Old World monkeys, the vervet monkey (Chlorocebus aethiops), the olive baboon (Papio anubis) and the colobus monkey (Colobus guereza). A total of 25 KIR gene models were identified in these species, predicted to encode receptors with 1, 2, and 3 extracellular Ig domains, all of them with long cytoplasmic domains having two putative ITIMs, although three had a positively charged residue in the transmembrane domain. Sequence and phylogenetic analyses showed that most Old World monkeys shared five classes of KIR loci: i) KIR2DL5/3DL20 in the most centromeric region, followed by ii) the single Ig domain-encoding locus KIR1D, iii) the pseudogene KIR2DP, iv) the conserved KIR2DL4, and v) the highly diversified KIR3DL/H loci in the telomeric half of the cluster. An exception to this pattern was the KIR haplotype of the colobus monkey that lacked the KIR1D, KIR2DP, and KIR2DL4 loci of the central region of the cluster. Thus, Old World monkeys display a broad spectrum of KIR haplotype variation that has been generated upon an ancestral haplotype architecture by gene duplication, gene deletion, and non-homologous recombination.