Tracking footprints of artificial selection in the dog genome

Selection against domestication alleles in introduced rabbit populations

Humans have moved domestic animals around the globe for thousands of years. These have occasionally established feral populations in nature, often with devastating ecological consequences. To understand how natural selection shapes re-adaptation into the wild, we investigated one of the most successful colonizers in history, the European rabbit. By sequencing the genomes of 297 rabbits across three continents, we show that introduced populations exhibit a mixed wild-domestic ancestry. We show that alleles that increased in frequency during domestication were preferentially selected against in novel natural environments. Interestingly, causative mutations for common domestication traits sometimes segregate at considerable frequencies if associated with less drastic phenotypes (for example, coat colour dilution), whereas mutations that are probably strongly maladaptive in nature are absent. Whereas natural selection largely targeted different genomic regions in each introduced population, some of the strongest signals of parallelism overlap genes associated with neuronal or brain function. This limited parallelism is probably explained by extensive standing genetic variation resulting from domestication together with the complex mixed ancestry of introduced populations. Our findings shed light on the selective and molecular mechanisms that enable domestic animals to re-adapt to the wild and provide important insights for the mitigation and management of invasive populations.

Historic dog Furs Unravel the Origin and Artificial Selection of Modern Nordic Lapphund and Elkhound dog Breeds

The origins and extreme morphological evolution of the modern dog breeds are poorly studied because the founder populations are extinct. Here, we analyse eight 100 to 200 years old dog fur samples obtained from traditional North Swedish clothing, to explore the origin and artificial selection of the modern Nordic Lapphund and Elkhound dog breeds. Population genomic analysis confirmed the Lapphund and Elkhound breeds to originate from the local dog population, and showed a distinct decrease in genetic diversity in agreement with intense breeding. We identified eleven genes under positive selection during the breed development. In particular, the MSRB3 gene, associated with breed-related ear morphology, was selected in all Lapphund and Elkhound breeds, and functional assays showed that a SNP mutation in the 3'UTR region suppresses its expression through miRNA regulation. Our findings demonstrate analysis of near-modern dog artifacts as an effective tool for interpreting the origin and artificial selection of the modern dog breeds.

Genetic characteristics and selection signatures between Southern Chinese local and commercial chickens

The introduction of exotic breeds and the cultivation of new lines by breeding companies have posed challenges to native chickens in South China, including loss of breed characteristics, decreased genetic diversity, and declining purity. Understanding the population genetic structure and genetic diversity of native chickens in South China is crucial for further advancements in breeding efforts. In this study, we analyzed the population genetic structure and genetic diversity of 321 individuals from 10 different breeds in South China. By comparing commercial chickens with native ones, we identified selection signatures occurring between local chickens and commercial breeds. The analysis of population genetic structure revealed that the native chicken populations in South China exhibited a considerable level of genetic diversity. Moreover, the commercial lines of Xiaobai chicken and Huangma chicken displayed even higher levels of genetic diversity, which distinguished them from other native varieties at the clustering level. However, certain individuals within these commercial varieties showed a discernible genetic relationship with the native populations. Notably, both commercial varieties also retained a significant degree of genetic similarity to their respective native counterparts. In order to investigate the genomic changes occurring during the commercialization of native chickens, we employed 4 methods (Fst, ROD, XPCLR, and XPEHH) to identify potential candidate regions displaying selective signatures in Southern Chinese native chicken population. A total of 168 (identified by Fst and ROD) and 86 (identified by XPCLR and XPEHH) overlapping genes were discovered. Functional annotation analysis revealed that these genes may be associated with reproduction and growth (SAMSN1, HYLS1, ROBO3, FGF14, PRSS23), musculoskeletal development (DNER, MYBPC1, DGKB, ORC1, KLF10), disease resistance and environmental adaptability (PUS3, CRB2, CALD1, USP15, SGCD, LTBP1), as well as egg production (ADGRB3, ACSF3). Overall, native chickens in South China harbor numerous selective sweep regions compared to commercial chickens, enriching valuable genomic resources for future genetic research and breeding conservation.

Evolutionary genomics of socially polymorphic populations of Pogonomyrmex californicus

BACKGROUND

Social insects vary considerably in their social organization both between and within species. In the California harvester ant, Pogonomyrmex californicus (Buckley 1867), colonies are commonly founded and headed by a single queen (haplometrosis, primary monogyny). However, in some populations in California (USA), unrelated queens cooperate not only during founding (pleometrosis) but also throughout the life of the colony (primary polygyny). The genetic architecture and evolutionary dynamics of this complex social niche polymorphism (haplometrosis vs pleometrosis) have remained unknown.

RESULTS

We provide a first analysis of its genomic basis and evolutionary history using population genomics comparing individuals from a haplometrotic population to those from a pleometrotic population. We discovered a recently evolved (< 200 k years), 8-Mb non-recombining region segregating with the observed social niche polymorphism. This region shares several characteristics with supergenes underlying social polymorphisms in other socially polymorphic ant species. However, we also find remarkable differences from previously described social supergenes. Particularly, four additional genomic regions not in linkage with the supergene show signatures of a selective sweep in the pleometrotic population. Within these regions, we find for example genes crucial for epigenetic regulation via histone modification (chameau) and DNA methylation (Dnmt1).

CONCLUSIONS

Altogether, our results suggest that social morph in this species is a polygenic trait involving a potential young supergene. Further studies targeting haplo- and pleometrotic individuals from a single population are however required to conclusively resolve whether these genetic differences underlie the alternative social phenotypes or have emerged through genetic drift.

Ethical Concerns about Fashionable Dog Breeding

The historical relationship between humans and dogs has involved selective breeding for various purposes, such as hunting, guarding, and service roles. However, over time, there has been a shift in preferences from functionality to aesthetics, which has influenced the diverse sizes, shapes, and coats of dog breeds. This review looks at fashionable dog breeding and questions the ethics of prioritising looks over health and behaviour. It aims to alert potential owners, breeders, and regulators to the importance of considering a dog's overall well-being, not just its appearance, which has resulted in fad breeding, leading to genetic disorders, health issues, and a loss of biodiversity. Ethical concerns arise from breeding brachycephalic breeds with respiratory conditions, inbreeding causing inherited disorders, and overbreeding popular breeds while shelter dogs remain unadopted. Additionally, the impact of cosmetic surgeries on popular dog breeds, as well as the neglect of behavioural traits in favour of physical characteristics and strict breeding practices are also considered. The current breeding model can have a negative impact on the emotional and cognitive well-being of dogs, resulting in issues such as aggression, anxiety, and other behavioural problems that can significantly reduce their overall quality of life. Unregulated breeding practices and the demand for rare breeds can lead to illegal breeding, compromising animal welfare. Prospective owners, veterinarians, kennel clubs, and legislators all need to play a responsible role in protecting animals.

Network of Inorganic Nanocrystals Can Swell: Study of Swelling-Induced Surface Instability

A unique organic-inorganic hybrid network composed of inorganic nanocores (ranging from semiconductors to metallic ones) interconnected through organic molecules can be produced by crosslinking the organic ligands of colloidal inorganic nanocrystals in assemblies. This work reports that this network, which is conventionally considered an inorganic film, can swell when exposed to a solvent because of the interaction between the solvent and the organic linkage within the network. Intriguingly, this work discovers that drying the solvent of the swollen organic-inorganic hybrid network can significantly affect the morphology owing to the swelling-induced compress stress, which is widely observed in various organic network systems. This work studies the surface instability of crosslinked organic-inorganic hybrid networks swollen by various organic solvents, which led to buckling delamination. Specifically, this work investigates the effects of the i) solvent-network interaction, ii) crosslinking density of the network, and iii) thickness of the film on the delamination behavior of the crosslinked network.

Longevity of companion dog breeds: those at risk from early death

The companion dog is one of the most phenotypically diverse species. Variability between breeds extends not only to morphology and aspects of behaviour, but also to longevity. Despite this fact, little research has been devoted to assessing variation in life expectancy between breeds or evaluating the potential for phylogenetic characterisation of longevity. Using a dataset of 584,734 unique dogs located within the UK, including 284,734 deceased, we present variation in longevity estimates within the following: parental lineage (purebred = 1 breed, crossbred ≥ 2 breeds), breed (n = 155), body size (large, medium, small), sex (male, female) and cephalic index (brachycephalic, mesocephalic, dolichocephalic). Survival estimates were then partitioned amongst phylogenetic clades: providing evidence that canine evolutionary history (via domestication and associated artificial selection) is associated with breed lifespan. This information provides evidence to inform discussions regarding pedigree health, whilst helping current/prospective owners, breeders, policy makers, funding bodies and welfare organisations improve decision making regarding canine welfare.

Integrative functional genomic analyses identify genetic variants influencing skin pigmentation in Africans

Skin color is highly variable in Africans, yet little is known about the underlying molecular mechanism. Here we applied massively parallel reporter assays to screen 1,157 candidate variants influencing skin pigmentation in Africans and identified 165 single-nucleotide polymorphisms showing differential regulatory activities between alleles. We combine Hi-C, genome editing and melanin assays to identify regulatory elements for MFSD12, HMG20B, OCA2, MITF, LEF1, TRPS1, BLOC1S6 and CYB561A3 that impact melanin levels in vitro and modulate human skin color. We found that independent mutations in an OCA2 enhancer contribute to the evolution of human skin color diversity and detect signals of local adaptation at enhancers of MITF, LEF1 and TRPS1, which may contribute to the light skin color of Khoesan-speaking populations from Southern Africa. Additionally, we identified CYB561A3 as a novel pigmentation regulator that impacts genes involved in oxidative phosphorylation and melanogenesis. These results provide insights into the mechanisms underlying human skin color diversity and adaptive evolution.

Negative Selection on a SOD1 Mutation Limits Canine Degenerative Myelopathy While Avoiding Inbreeding

Several hundred disease-causing mutations are currently known in domestic dogs. Breeding management is therefore required to minimize their spread. Recently, genetic methods such as direct-to-consumer testing have gained popularity; however, their effects on dog populations are unclear. Here, we aimed to evaluate the influence of genetic testing on the frequency of mutations responsible for canine degenerative myelopathy and assess the changes in the genetic structure of a Pembroke Welsh corgi population from Japan. Genetic testing of 5,512 dogs for the causative mutation in superoxide dismutase 1 (SOD1) (c.118G>A (p.E40K)) uncovered a recent decrease in frequency, plummeting from 14.5% (95/657) in 2019 to 2.9% (24/820) in 2022. Weir and Cockerham population differentiation (FST) based on genome-wide single-nucleotide polymorphism (SNP) of 117 selected dogs detected the SNP with the highest FST located in the intron of SOD1 adjacent to the c.118G>A mutation, supporting a selection signature on SOD1. Further genome-wide SNP analyses revealed no obvious changes in inbreeding levels and genetic diversity between the 2019 and 2022 populations. Our study highlights that genetic testing can help inform improved mating choices in breeding programs to reduce the frequency of risk variants and avoid inbreeding. This combined strategy could decrease the genetic risk of canine degenerative myelopathy, a fatal disease, within only a few years.

Integrated analysis strategy of genome-wide functional gene mining reveals DKK2 gene underlying meat quality in Shaziling synthesized pigs

BACKGROUND

Shaziling pig is a well-known indigenous breed in China who has superior meat quality traits. However, the genetic mechanism and genomic evidence underlying meat quality characteristics of Shaziling pigs are still unclear. To explore and investigate the germplasm characteristics of Shaziling pigs, we totally analyzed 67 individual's whole genome sequencing data for the first time (20 Shaziling pigs [S], 20 Dabasha pigs [DBS], 11 Yorkshire pigs [Y], 10 Berkshire pigs [BKX], 5 Basha pigs [BS] and 1 Warthog).

RESULTS

A total of 2,538,577 SNPs with high quality were detected and 9 candidate genes which was specifically selected in S and shared in S to DBS were precisely mined and screened using an integrated analysis strategy of identity-by-descent (IBD) and selective sweep. Of them, dickkopf WNT signaling pathway inhibitor 2 (DKK2), the antagonist of Wnt signaling pathway, was the most promising candidate gene which was not only identified an association of palmitic acid and palmitoleic acid quantitative trait locus in PigQTLdb, but also specifically selected in S compared to other 48 Chinese local pigs of 12 populations and 39 foreign pigs of 4 populations. Subsequently, a mutation at 12,726-bp of DKK2 intron 1 (g.114874954 A > C) was identified associated with intramuscular fat content using method of PCR-RFLP in 21 different pig populations. We observed DKK2 specifically expressed in adipose tissues. Overexpression of DKK2 decreased the content of triglyceride, fatty acid synthase and expression of relevant genes of adipogenic and Wnt signaling pathway, while interference of DKK2 got contrary effect during adipogenesis differentiation of porcine preadipocytes and 3T3-L1 cells.

CONCLUSIONS

Our findings provide an analysis strategy for mining functional genes of important economic traits and provide fundamental data and molecular evidence for improving pig meat quality traits and molecular breeding.

Is it now time to iron out the wrinkles? Health of Shar Pei dogs under primary veterinary care in the UK

BACKGROUND

The Shar Pei is a common dog breed with a distinctive appearance caused by hyaluronosis that has been linked with several health conditions. Anonymised primary-care veterinary clinical records were explored to extract data on the demography, common disorders and mortality of Shar Pei in the UK in 2013.

RESULTS

The study population of 455,557 dogs included 1913 (0.42%) Shar Pei. The mean adult bodyweight was 22.26 kg. The most prevalent fine-level precision disorders were entropion (prevalence 17.88%, 95% CI: 16.16-19.59), otitis externa (16.36%, 95% CI: 14.70-18.02), ear disorders (6.69%, 95% CI: 5.57-7.81), aggression (5.23%, 95% CI: 4.23-6.22), and pyoderma (4.29%, 95% CI: 3.38-5.19). The most prevalent disorder groups were ophthalmologic (prevalence = 22.27%, 95% CI: 20.40-24.13), dermatological (21.01%, 95% CI: 19.19-22.84), aural (18.66%, 95% CI: 16.92-20.41), traumatic injury (7.53%, 95% CI: 6.35-8.71) and behavioural (7.21%, 95% CI: 6.05-8.37). The median longevity of 190 Shar Pei that died during the study period was 7.28 years (IQR 5.04-10.05, range 0.04-15.04). Of 184 deaths with a recorded method of death, 157 (85.33%) deaths involved euthanasia and 27 (14.67%) deaths were unassisted. Among 136 (71.58%) deaths with a recorded biomedical cause of death, the most common causes of death at group level precision were neoplasia (15.44%, 95% CI: 9.37-21.51), renal disorders (13.24%, 95% CI: 7.54-18.93), and behavioural disorders (11.03%, 95% CI: 5.76-16.29).

CONCLUSIONS

Almost one fifth of Shar Pei receive veterinary care each year for entropion, a condition linked strongly with the extreme conformation of thickened and folded skin and bristly hair that characterises the Shar Pei breed. Several other common disorders are also linked to hyaluronosis. Current UK legislation can help support efforts to avoid breeding or acquiring animals with extreme conformations and to promote adequate veterinary care for already-owned animals with extreme conformations.

Genetic Origin and Introgression Pattern of Pingliang Red Cattle Revealed Using Genome-Wide SNP Analyses

The Pingliang red cattle, an outstanding indigenous resource in China, possesses an exceptional breeding value attributed to its tender meat and superior marbling quality. Currently, research efforts have predominantly concentrated on exploring its maternal origin and conducting conventional phenotypic studies. However, there remains a lack of comprehensive understanding regarding its genetic basis. To address this gap, we conducted a thorough whole-genome analysis to investigate the population structure, phylogenetic relationships, and gene flows of this breed using genomic SNP chip data from 17 bovine breeds. The results demonstrate that Pingliang red cattle have evolved distinct genetic characteristics unique to this breed, clearly distinguishing it from other breeds. Based on the analysis of the population structure and phylogenetic tree, it can be classified as a hybrid lineage between Bos taurus and Bos indicus. Furthermore, Pingliang red cattle display a more prominent B. taurus pedigree in comparison with Jinnan, Qinchuan, Zaosheng, Nanyang, and Luxi cattle. Moreover, this study also revealed closer genetic proximity within the Chinese indigenous cattle breed, particularly Qinchuan cattle, which shares the longest identical by descent (IBD) fragment with Pingliang red cattle. Gene introgression analysis shows that Pingliang red cattle have undergone gene exchange with South Devon and Red Angus cattle from Europe. Admixture analysis revealed that the proportions of East Asian taurine and Chinese indicine in the ancestry of Pingliang red cattle are approximately 52.44% and 21.00%, respectively, while Eurasian taurine, European taurine, and Indian indicine account for approximately 17.55%, 7.27%, and 1.74%. Our findings unveil distinct genetic characteristics in Pingliang red cattle and attribute their origin to B. taurus and B. indicus ancestry, as well as contributions from Qinchuan cattle, South Devon, and Red Angus.

A purebred South American breed showing high effective population size and independent breed ancestry: The Chilean Terrier

The Chilean Terrier is a known breed in Chile that has not been genetically assessed despite its distinctive color patterns, agility, and hardiness across the diversity of climates encountered within the Chilean landscape. The population structure and its relatedness with other breeds, as well as the actual origin of the breed, remain unknown. We estimated several population parameters using samples from individuals representing the distribution of the Chilean Terrier across the country. By utilizing the Illumina HD canine genotyping array, we computed the effective population size (Ne ), individual inbreeding, and relatedness to evaluate the genetic diversity of the breed. The results show that linkage disequilibrium was relatively low and decayed rapidly; in fact, Ne was very high when compared to other breeds, and similar to other American indigenous breeds (such as the Chihuahua with values of Ne near 500). These results are in line with the low estimates of genomic inbreeding and relatedness and the relatively large number of effective chromosome segments (Me  = 2467) obtained using the properties of the genomic relationship matrix. Between population analysis (cross-population extended haplotype homozygosity, di ) with other breeds such as the Jack Russell Terrier, the Peruvian-Inca Orchid, and the Chihuahua suggested that candidate regions harboring FGF5, PAX3, and ASIP, probably explained some morphological traits, such as the distinctive color pattern characteristic of the breed. When considering Admixture estimates and phylogenetic analysis, together with other breeds of American and European origin, the Chilean Terrier does not have a recent European ancestry. Overall, the results suggest that the breed has evolved independently in Chile from other terrier breeds, from an unknown European terrier ancestor.

Genetic variation and domestication of horses revealed by 10 chromosome-level genomes and whole-genome resequencing

Understanding the genetic variations of the horse (Equus caballus) genome will improve breeding conservation and welfare. However, genetic variations in long segments, such as structural variants (SVs), remain understudied. We de novo assembled 10 chromosome-level three-dimensional horse genomes, each representing a distinct breed, and analysed horse SVs using a multi-assembly approach. Our findings suggest that SVs with the accumulation of mammalian-wide interspersed repeats related to long interspersed nuclear elements might be a horse-specific mechanism to modulate genome-wide gene regulatory networks. We found that olfactory receptors were commonly loss and accumulated deleterious mutations, but no purge of deleterious mutations occurred during horse domestication. We examined the potential effects of SVs on the spatial structure of chromatin via topologically associating domains (TADs). Breed-specific TADs were significantly enriched by breed-specific SVs. We identified 4199 unique breakpoint-resolved novel insertions across all chromosomes that account for 2.84 Mb sequences missing from the reference genome. Several novel insertions might have potential functional consequences, as 519 appeared to reside within 449 gene bodies. These genes are primarily involved in pathogen recognition, innate immune responses and drug metabolism. Moreover, 37 diverse horses were resequenced. Combining this with public data, we analysed 97 horses through a comparative population genomics approach to identify the genetic basis underlying breed characteristics using Thoroughbreds as a case study. We provide new scientific evidence for horse domestication, an understanding of the genetic mechanism underlying the phenotypic evolution of horses, and a comprehensive genetic variation resource for further genetic studies of horses.

Whole-genome sequencing reveals a complex African population demographic history and signatures of local adaptation

We conduct high coverage (>30×) whole-genome sequencing of 180 individuals from 12 indigenous African populations. We identify millions of unreported variants, many predicted to be functionally important. We observe that the ancestors of southern African San and central African rainforest hunter-gatherers (RHG) diverged from other populations >200 kya and maintained a large effective population size. We observe evidence for ancient population structure in Africa and for multiple introgression events from "ghost" populations with highly diverged genetic lineages. Although currently geographically isolated, we observe evidence for gene flow between eastern and southern Khoesan-speaking hunter-gatherer populations lasting until ∼12 kya. We identify signatures of local adaptation for traits related to skin color, immune response, height, and metabolic processes. We identify a positively selected variant in the lightly pigmented San that influences pigmentation in vitro by regulating the enhancer activity and gene expression of PDPK1.

The genetic and evolutionary basis of gene expression variation in East Africans

BACKGROUND

Mapping of quantitative trait loci (QTL) associated with molecular phenotypes is a powerful approach for identifying the genes and molecular mechanisms underlying human traits and diseases, though most studies have focused on individuals of European descent. While important progress has been made to study a greater diversity of human populations, many groups remain unstudied, particularly among indigenous populations within Africa. To better understand the genetics of gene regulation in East Africans, we perform expression and splicing QTL mapping in whole blood from a cohort of 162 diverse Africans from Ethiopia and Tanzania. We assess replication of these QTLs in cohorts of predominantly European ancestry and identify candidate genes under selection in human populations.

RESULTS

We find the gene regulatory architecture of African and non-African populations is broadly shared, though there is a considerable amount of variation at individual loci across populations. Comparing our analyses to an equivalently sized cohort of European Americans, we find that QTL mapping in Africans improves the detection of expression QTLs and fine-mapping of causal variation. Integrating our QTL scans with signatures of natural selection, we find several genes related to immunity and metabolism that are highly differentiated between Africans and non-Africans, as well as a gene associated with pigmentation.

CONCLUSION

Extending QTL mapping studies beyond European ancestry, particularly to diverse indigenous populations, is vital for a complete understanding of the genetic architecture of human traits and can reveal novel functional variation underlying human traits and disease.

Analysis of dog breed diversity using a composite selection index

During breed development, domestic dogs have undergone genetic bottlenecks and sustained selective pressures, as a result distinctive genomic diversity occurs to varying degrees within and between breed groups. This diversity can be identified using standard methods or combinations of these methods. This study explored the application of a combined selection index, composite selection signals (CSS), derived from multiple methods to an existing genotype dataset from three breed groups developed in distinct regions of Asia: Qinghai-Tibet plateau dogs (adapted to living at altitude), Xi dogs (with superior running ability) and Mountain hounds (used for hunting ability). The CSS analysis confirmed top ranked genomic regions on CFA10 and CFA21 in Qinghai-Tibet plateau dogs, CFA1 in Xi dogs and CFA5 in Mountain hounds. CSS analysis identified additional significant genomic regions in each group, defined by a total of 1,397, 1,475 and 1,675 significant SNPs in the Qinghai-Tibetan Plateau dogs, Xi dogs and Mountain hounds, respectively. Chitinase 3 Like 1 (CHI3L1) and Leucine Rich Repeat Containing G Protein-Coupled Receptor 6 (LGR6) genes were located in the top ranked region on CFA7 (0.02-1 Mb) in the Qinghai-Tibetan Plateau dogs. Both genes have been associated with hypoxia responses or altitude adaptation in humans. For the Xi dogs, the top ranked region on CFA25 contained the Transient Receptor Potential Cation Channel Subfamily C Member 4 (TRPC4) gene. This calcium channel is important for optimal muscle performance during exercise. The outstanding signals in the Mountain dogs were on CFA5 with 213 significant SNPs that spanned genes involved in cardiac development, sight and generation of biochemical energy. These findings support the use of the combined index approach for identifying novel regions of genome diversity in dogs. As with other methods, the results do not prove causal links between these regions and phenotypes, but they may assist in focusing future studies that seek to identify functional pathways that contribute to breed diversity.

Large-Scale Chromosomal Changes Lead to Genome-Level Expression Alterations, Environmental Adaptation, and Speciation in the Gayal (Bos frontalis)

Determining the functional consequences of karyotypic changes is invariably challenging because evolution tends to obscure many of its own footprints, such as accumulated mutations, recombination events, and demographic perturbations. Here, we describe the assembly of a chromosome-level reference genome of the gayal (Bos frontalis) thereby revealing the structure, at base-pair-level resolution, of a telo/acrocentric-to-telo/acrocentric Robertsonian translocation (2;28) (T/A-to-T/A rob[2;28]). The absence of any reduction in the recombination rate or genetic introgression within the fusion region of gayal served to challenge the long-standing view of a role for fusion-induced meiotic dysfunction in speciation. The disproportionate increase noted in the distant interactions across pro-chr2 and pro-chr28, and the change in open-chromatin accessibility following rob(2;28), may, however, have led to the various gene expression irregularities observed in the gayal. Indeed, we found that many muscle-related genes, located synthetically on pro-chr2 and pro-chr28, exhibited significant changes in expression. This, combined with genome-scale structural variants and expression alterations in genes involved in myofibril composition, may have driven the rapid sarcomere adaptation of gayal to its rugged mountain habitat. Our findings not only suggest that large-scale chromosomal changes can lead to alterations in genome-level expression, thereby promoting both adaptation and speciation, but also illuminate novel avenues for studying the relationship between karyotype evolution and speciation.

Selection signature analysis and genome-wide divergence of South African Merino breeds from their founders

Merino sheep are a breed of choice across the world, popularly kept for their wool and mutton value. They are often reared as a pure breed or used in crossbreeding and are a common component in synthetic breed development. This study evaluated genetic diversity, population structure, and breed divergence in 279 animals of Merino and Merino-based sheep breeds in South Africa using the Illumina Ovine SNP 50K BeadChip. The sheep breeds analysed included the three Merino-derived breeds of Dohne Merino (n = 50); Meatmaster (n = 47); and Afrino (n = 52) and five presumed ancestral populations of Merinos (Merino (n = 46); South African Merino (n = 10); and South African Mutton Merino (n = 8)); and the non-Merino founding breeds of Damara (n = 20); Ronderib Afrikaner (n = 17); and Nguni (n = 29). Highest genetic diversity values were observed in the Dohne Merino (DM), with H _o = 0.39 ± 0.01, followed by the Meatmaster and South African Merino (SAM), with H _o = 0.37 ± 0.03. The level of inbreeding ranged from 0.0 ± 0.02 (DM) to 0.27 ± 0.05 (Nguni). Analysis of molecular variance (AMOVA) showed high within-population variance (>80%) across all population categories. The first principal component (PC1) separated the Merino, South African Mutton Merino (SAMM), DM, and Afrino (AFR) from the Meatmaster, Damara, Nguni, and Ronderib Afrikaner (RDA). PC2 aligned each Merino-derived breed with its presumed ancestors and separated the SAMM from the Merino and SAM. The iHS analysis yielded selection sweeps across the AFR (12 sweeps), Meatmaster (four sweeps), and DM (29 sweeps). Hair/wool trait genes such as FGF12; metabolic genes of ICA1, NXPH1, and GPR171; and immune response genes of IL22, IL26, IFNAR1, and IL10RB were reported. Other genes include HMGA, which was observed as selection signatures in other populations; WNT5A, important in the development of the skeleton and mammary glands; ANTXR2, associated with adaptation to variation in climatic conditions; and BMP2, which has been reported as strongly selected in both fat-tailed and thin-tailed sheep. The DM vs. SAMM shared all six sweep regions on chromosomes 1, 10, and 11 with AFR vs. SAMM. Genes such as FGF12 on OAR 1:191.3-194.7 Mb and MAP2K4 on OAR 11:28.6-31.3 Mb were observed. The selection sweep on chromosome 10 region 28.6-30.3 Mb harbouring the RXFP2 for polledness was shared between the DM vs. Merino, the Meatmaster vs. Merino, and the Meatmaster vs. Nguni. The DM vs. Merino and the Meatmaster vs. Merino also shared an Rsb-based selection sweep on chromosome 1 region 268.5-269.9 Mb associated with the Calpain gene, CAPN7. The study demonstrated some genetic similarities between the Merino and Merino-derived breeds emanating from common founding populations and some divergence driven by breed-specific selection goals. Overall, information regarding the evolution of these composite breeds from their founding population will guide future breed improvement programs and management and conservation efforts.

Out of the Qinghai-Tibet plateau: Genomic biogeography of the alpine monospecific genus Megadenia (Biscutelleae, Brassicaceae)

Numerous high-elevation alpine plants of the Qinghai-Tibet Plateau (QTP) also have disjunct distribution in adjacent low-altitude mountains. The out-of-QTP versus into-the-QTP hypothesis of alpine plants provide strong evidence for the highly disputed assumption of the massive ice sheet developed in the central plateau during the Last Glacial Maximum (LGM). In this study, we sequenced the genomes of most known populations of Megadenia, a monospecific alpine genus of Brassicaceae distributed primarily in the QTP, though rarely found in adjacent low-elevation mountains of north China and Russia (NC-R). All sequenced samples clustered into four geographic genetic groups: one pair was in the QTP and another was in NC-R. The latter pair is nested within the former, and these findings support the out-of-QTP hypothesis. Dating the four genetic groups and niche distribution suggested that Megadenia migrated out of the QTP to adjacent regions during the LGM. The NC-R group showed a decrease in the effective population sizes. In addition, the genes with high genetic divergences in the QTP group were mainly involved in habitat adaptations during low-altitude colonization. These findings reject the hypothesis of development massive ice sheets, and support glacial survival of alpine plants within, as well as further migration out of, the QTP.

Genomic correlates for migratory direction in a free-ranging cervid

Animal migrations are some of the most ubiquitous and one of the most threatened ecological processes globally. A wide range of migratory behaviours occur in nature, and this behaviour is not uniform among and within species, where even individuals in the same population can exhibit differences. While the environment largely drives migratory behaviour, it is necessary to understand the genetic mechanisms influencing migration to elucidate the potential of migratory species to cope with novel conditions and adapt to environmental change. In this study, we identified genes associated with a migratory trait by undertaking pooled genome-wide scans on a natural population of migrating mule deer. We identified genomic regions associated with variation in migratory direction, including FITM1, a gene linked to the formation of lipids, and DPPA3, a gene linked to epigenetic modifications of the maternal line. Such a genetic basis for a migratory trait contributes to the adaptive potential of the species and might affect the flexibility of individuals to change their behaviour in the face of changes in their environment.

Domestic dog lineages reveal genetic drivers of behavioral diversification

Selective breeding of domestic dogs has generated diverse breeds often optimized for performing specialized tasks. Despite the heritability of breed-typical behavioral traits, identification of causal loci has proven challenging due to the complexity of canine population structure. We overcome longstanding difficulties in identifying genetic drivers of canine behavior by developing a framework for understanding relationships between breeds and the behaviors that define them, utilizing genetic data for over 4,000 domestic, semi-feral, and wild canids and behavioral survey data for over 46,000 dogs. We identify ten major canine genetic lineages and their behavioral correlates and show that breed diversification is predominantly driven by non-coding regulatory variation. We determine that lineage-associated genes converge in neurodevelopmental co-expression networks, identifying a sheepdog-associated enrichment for interrelated axon guidance functions. This work presents a scaffold for canine diversification that positions the domestic dog as an unparalleled system for revealing the genetic origins of behavioral diversity.

Pleiotropic functions of chordin gene causing drastic morphological changes in ornamental goldfish

Breeders and fanciers have established many peculiar morphological phenotypes in ornamental goldfish. Among them, the twin-tail and dorsal-finless phenotypes have particularly intrigued early and recent researchers, as equivalent morphologies are extremely rare in nature. These two mutated phenotypes appeared almost simultaneously within a short time frame and were fixed in several strains. However, little is known about how these two different mutations could have co-occurred during such a short time period. Here, we demonstrate that the chordin gene, a key factor in dorsal-ventral patterning, is responsible not only for the twin-tail phenotype but also for the dorsal-finless phenotype. Our F2 backcrossing and functional analyses revealed that the penetrance/expressivity of the dorsal-finless phenotype can be suppressed by the wild-type allele of chdS. Based on these findings, we propose that chdSwt may have masked the expression of the dorsal-finless phenotype, acting as a capacitor buffering gene to allow accumulation of genetic mutations. Once this gene lost its original function in the twin-tail goldfish lineages, the dorsal-finless phenotype could be highly expressed. Thus, this study experimentally demonstrates that the rapid genetic fixation of morphological mutations during a short domestication time period may be related to the robustness of embryonic developmental mechanisms.

A ∼ 4.1 kb deletion in IRX1 gene upstream is completely associated with rumplessness in Piao chicken

Piao chicken, a Chinese indigenous rumpless chicken breed, lacks pygostyle, caudal vertebra, uropygial gland and tail feathers. The rumplessness in Piao chicken presents an autosomal dominant inheritance pattern. However, the molecular genetic mechanisms underlying the rumplessness in Piao chicken remains unclear. In this study, whole-genome resequencing was performed for 146 individuals from 10 chicken breeds, including 9 tailed chicken breeds and Piao rumpless breed. Tailbone CT scan for Piao chickens and WL chickens, revealed that some Piao chicken tails were normal in number, and for a few Piao chickens tail length and tail bone numbers were between the rumpless and the normal tailed chickens. The results showed that the rumpless phenotype has not been completely fixed in Piao chicken breed. Using selection signature analysis and structural variation detection, we found a 4174 bp deletion located in the upstream region of IRX1 gene on chromosome 2 related to rumpless phenotype. Structural variation genotyping showed that the deletion was present in all 32 rumpless Piao chickens (del/del, wild/del) and absent from all 112 tailed chickens included in the dataset for the other 9 breeds and 2 tailed Piao chickens (wild/wild). In summary, all rumpless Piao chickens tested here carry this deletion mutation, to show a complete linkage association with rumplessness trait. We suggested that the 4174 bp deletion could be causative for rumpless phenotype in Piao chicken since this is the only mutation to show the complete linkage disequilibrium with rumplessness on whole genome level across all of 146 chickens from the 10 breeds. This study could facilitate a better understanding of the genetic characteristics of Piao chicken.

Evolutionary endocrinology and the problem of Darwin's tangled bank

Like Darwin's tangled bank of biodiversity, the endocrine mechanisms that give rise to phenotypic diversity also exhibit nearly endless forms. This tangled bank of mechanistic diversity can prove problematic as we seek general principles on the role of endocrine mechanisms in phenotypic evolution. A key unresolved question is therefore: to what degree are specific endocrine mechanisms re-used to bring about replicated phenotypic evolution? Related areas of inquiry are booming in molecular ecology, but behavioral traits are underrepresented in this literature. Here, I leverage the rich comparative tradition in evolutionary endocrinology to evaluate whether and how certain mechanisms may be repeated hotspots of behavioral evolutionary change. At one extreme, mechanisms may be parallel, such that evolution repeatedly uses the same gene or pathway to arrive at multiple independent (or, convergent) origins of a particular behavioral trait. At the other extreme, the building blocks of behavior may be unique, such that outwardly similar phenotypes are generated via lineage-specific mechanisms. This review synthesizes existing case studies, phylogenetic analyses, and experimental evolutionary research on mechanistic parallelism in animal behavior. These examples show that the endocrine building blocks of behavior have some elements of parallelism across replicated evolutionary events. However, support for parallelism is variable among studies, at least some of which relates to the level of complexity at which we consider sameness (i.e. pathway vs. gene level). Moving forward, we need continued experimentation and better testing of neutral models to understand whether, how - and critically, why - mechanism A is used in one lineage and mechanism B is used in another. We also need continued growth of large-scale comparative analyses, especially those that can evaluate which endocrine parameters are more or less likely to undergo parallel evolution alongside specific behavioral traits. These efforts will ultimately deepen understanding of how and why hormone-mediated behaviors are constructed the way that they are.

Signatures of Convergent Evolution and Natural Selection at the Alcohol Dehydrogenase Gene Region are Correlated with Agriculture in Ethnically Diverse Africans

The alcohol dehydrogenase (ADH) family of genes encodes enzymes that catalyze the metabolism of ethanol into acetaldehyde. Nucleotide variation in ADH genes can affect the catalytic properties of these enzymes and is associated with a variety of traits, including alcoholism and cancer. Some ADH variants, including the ADH1B*48His (rs1229984) mutation in the ADH1B gene, reduce the risk of alcoholism and are under positive selection in multiple human populations. The advent of Neolithic agriculture and associated increase in fermented foods and beverages is hypothesized to have been a selective force acting on such variants. However, this hypothesis has not been tested in populations outside of Asia. Here, we use genome-wide selection scans to show that the ADH gene region is enriched for variants showing strong signals of positive selection in multiple Afroasiatic-speaking, agriculturalist populations from Ethiopia, and that this signal is unique among sub-Saharan Africans. We also observe strong selection signals at putatively functional variants in nearby lipid metabolism genes, which may influence evolutionary dynamics at the ADH region. Finally, we show that haplotypes carrying these selected variants were introduced into Northeast Africa from a West-Eurasian source within the last ∼2,000 years and experienced positive selection following admixture. These selection signals are not evident in nearby, genetically similar populations that practice hunting/gathering or pastoralist subsistence lifestyles, supporting the hypothesis that the emergence of agriculture shapes patterns of selection at ADH genes. Together, these results enhance our understanding of how adaptations to diverse environments and diets have influenced the African genomic landscape.

Experimental test of the contributions of initial variation and new mutations to adaptive evolution in a novel environment

Experimental evolution is an approach that allows researchers to study organisms as they evolve in controlled environments. Despite the growing popularity of this approach, there are conceptual gaps among projects that use different experimental designs. One such gap concerns the contributions to adaptation of genetic variation present at the start of an experiment and that of new mutations that arise during an experiment. The primary source of genetic variation has historically depended largely on the study organisms. In the long-term evolution experiment (LTEE) using Escherichia coli, for example, each population started from a single haploid cell, and therefore, adaptation depended entirely on new mutations. Most other microbial evolution experiments have followed the same strategy. By contrast, evolution experiments using multicellular, sexually reproducing organisms typically start with preexisting variation that fuels the response to selection. New mutations may also come into play in later generations of these experiments, but it is generally difficult to quantify their contribution in these studies. Here, we performed an experiment using E. coli to compare the contributions of initial genetic variation and new mutations to adaptation in a new environment. Our experiment had four treatments that varied in their starting diversity, with 18 populations in each treatment. One treatment depended entirely on new mutations, while the other three began with mixtures of clones, whole-population samples, or mixtures of whole-population samples from the LTEE. We tracked a genetic marker associated with different founders in two treatments. These data revealed significant variation in fitness among the founders, and that variation impacted evolution in the early generations of our experiment. However, there were no differences in fitness among the treatments after 500 or 2,000 generations in the new environment, despite the variation in fitness among the founders. These results indicate that new mutations quickly dominated, and eventually they contributed more to adaptation than did the initial variation. Our study thus shows that preexisting genetic variation can have a strong impact on early evolution in a new environment, but new beneficial mutations may contribute more to later evolution and can even drive some initially beneficial variants to extinction.

Assembly of a pangenome for global cattle reveals missing sequences and novel structural variations, providing new insights into their diversity and evolutionary history

A cattle pangenome representation was created based on the genome sequences of 898 cattle representing 57 breeds. The pangenome identified 83 Mb of sequence not found in the cattle reference genome, representing 3.1% novel sequence compared with the 2.71-Gb reference. A catalog of structural variants developed from this cattle population identified 3.3 million deletions, 0.12 million inversions, and 0.18 million duplications. Estimates of breed ancestry and hybridization between cattle breeds using insertion/deletions as markers were similar to those produced by single nucleotide polymorphism-based analysis. Hundreds of deletions were observed to have stratification based on subspecies and breed. For example, an insertion of a Bov-tA1 repeat element was identified in the first intron of the APPL2 gene and correlated with cattle breed geographic distribution. This insertion falls within a segment overlapping predicted enhancer and promoter regions of the gene, and could affect important traits such as immune response, olfactory functions, cell proliferation, and glucose metabolism in muscle. The results indicate that pangenomes are a valuable resource for studying diversity and evolutionary history, and help to delineate how domestication, trait-based breeding, and adaptive introgression have shaped the cattle genome.

Recombination landscape divergence between populations is marked by larger low-recombining regions in domesticated rye

The genomic landscape of recombination plays an essential role in evolution. Patterns of recombination are highly variable along chromosomes, between sexes, individuals, populations, and species. In many eukaryotes, recombination rates are elevated in sub-telomeric regions and drastically reduced near centromeres, resulting in large low-recombining (LR) regions. The processes of recombination are influenced by genetic factors, such as different alleles of genes involved in meiosis and chromatin structure, as well as external environmental stimuli like temperature and overall stress. In this work, we focused on the genomic landscapes of recombination in a collection of 916 rye (Secale cereale) individuals. By analysing population structure among individuals of different domestication status and geographic origin, we detected high levels of admixture, reflecting the reproductive biology of a self-incompatible, wind-pollinating grass species. We then analysed patterns of recombination in overlapping subpopulations, which revealed substantial variation in the physical size of LR regions, with a tendency for larger LR regions in domesticated subpopulations. Genome-wide association scans (GWAS) for LR region size revealed a major quantitative-trait-locus (QTL) at which, among 18 annotated genes, an ortholog of histone H4 acetyltransferase ESA1 was located. Rye individuals belonging to domesticated subpopulations showed increased synaptonemal complex length, but no difference in crossover frequency, indicating that only the recombination landscape is different. Furthermore, the genomic region harbouring rye ScESA1 showed moderate patterns of selection in domesticated subpopulations, suggesting that larger LR regions were indirectly selected for during domestication to achieve more homogeneous populations for agricultural use.

Impact of natural selection on global patterns of genetic variation and association with clinical phenotypes at genes involved in SARS-CoV-2 infection

Human genomic diversity has been shaped by both ancient and ongoing challenges from viruses. The current coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a devastating impact on population health. However, genetic diversity and evolutionary forces impacting host genes related to SARS-CoV-2 infection are not well understood. We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (angiotensin converting enzyme 2 [ACE2], transmembrane protease serine 2 [TMPRSS2], dipeptidyl peptidase 4 [DPP4], and lymphocyte antigen 6 complex locus E [LY6E]). We analyzed data from 2,012 ethnically diverse Africans and 15,977 individuals of European and African ancestry with electronic health records and integrated with global data from the 1000 Genomes Project. At ACE2, we identified 41 nonsynonymous variants that were rare in most populations, several of which impact protein function. However, three nonsynonymous variants (rs138390800, rs147311723, and rs145437639) were common among central African hunter-gatherers from Cameroon (minor allele frequency 0.083 to 0.164) and are on haplotypes that exhibit signatures of positive selection. We identify signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage compared with the chimpanzee genome. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19. Our study provides insights into global variation at host genes related to SARS-CoV-2 infection, which have been shaped by natural selection in some populations, possibly due to prior viral infections.

Capturing Genetic Diversity and Selection Signatures of the Endangered Kosovar Balusha Sheep Breed

There is a growing concern about the loss of animal genetic resources. The aim of this study was to analyze the genetic diversity and potential peculiarity of the endangered Kosovar sheep breed Balusha. For this purpose, a dataset consisting of medium-density SNP chip genotypes (39,879 SNPs) from 45 Balusha sheep was generated and compared with SNP chip genotypes from 29 individuals of a second Kosovar breed, Bardhoka. Publicly available SNP genotypes from 39 individuals of the relatively closely located sheep breeds Istrian Pramenka and Ruda were additionally included in the analyses. Analysis of heterozygosity, allelic richness and effective population size was used to assess the genetic diversity. Inbreeding was evaluated using two different methods (F_IS, F_ROH). The standardized F_ST (d_i) and cross-population extended haplotype homozygosity (XPEHH) methods were used to detect signatures of selection. We observed the lowest heterozygosity (H_O = 0.351) and effective population size (Ne₅ = 25, Ne₅₀ = 228) for the Balusha breed. The mean allelic richness levels (1.780–1.876) across all analyzed breeds were similar and also comparable with those in worldwide breeds. F_ROH estimates (0.023–0.077) were highest for the Balusha population, although evidence of decreased inbreeding was observed in F_IS results for the Balusha breed. Two Gene Ontology (GO) TERMs were strongly enriched for Balusha, and involved genes belonging to the melanogenesis and T cell receptor signaling pathways, respectively. This could result from selection for the special coat color pattern of Balusha (black head) and resistance to certain infectious diseases. The analyzed diversity parameters highlight the urgency to preserve the local Kosovar Balusha sheep as it is clearly distinguished from other sheep of Southeastern Europe, has the lowest diversity level and may harbor valuable genetic variants, e.g., for resistance to infectious diseases.

Resequencing and Signatures of Selective Scans Point to Candidate Genetic Variants for Hair Length Traits in Long-Haired and Normal-Haired Tianzhu White Yak

Tianzhu white yak is a rare local yak breed with a pure white coat in China. In recent years, breeders have discovered long-haired individuals characterized by long hair on the forehead in the Tianzhu white yak, and the length and density of the hair on these two parts of the body are higher than that of the normal Tianzhu white yak. To elucidate the genetic mechanism of hair length in Tianzhu white yak, we re-sequence the whole genome of long-haired Tianzhu White yak (LTWY) (n = 10) and normal Tianzhu White yak (NTWY) (n = 10). Then, fixation index (F ST), θπ ratio, cross-population composite likelihood ratio (XP-CLR), integrated haplotype score (iHS), cross-population extended haplotype homozygosity (XP-EHH), and one composite method, the de-correlated composite of multiple signals (DCMS) were performed to discover the loci and genes related to long-haired traits. Based on five single methods, we found two hotspots of 0.2 and 1.1 MB in length on chromosome 6, annotating two (FGF5, CFAP299) and four genes (ATP8A1, SLC30A9, SHISA3, TMEM33), respectively. Function enrichment analysis of genes in two hotspots revealed Ras signaling pathway, MAPK signaling pathway, PI3K-Akt signaling pathway, and Rap1 signaling pathway were involved in the process of hair length differences. Besides, the DCMS method further found that four genes (ACOXL, PDPK1, MAGEL2, CDH1) were associated with hair follicle development. Henceforth, our work provides novel genetic insights into the mechanisms of hair growth in the LTWY.

Genomic diversity and signatures of selection in meat and fancy rabbit breeds based on high-density marker data

Background

Domestication of the rabbit (Oryctolagus cuniculus) has led to a multi-purpose species that includes many breeds and lines with a broad phenotypic diversity, mainly for external traits (e.g. coat colours and patterns, fur structure, and morphometric traits) that are valued by fancy rabbit breeders. As a consequence of this human-driven selection, distinct signatures are expected to be present in the rabbit genome, defined as signatures of selection or selective sweeps. Here, we investigated the genome of three Italian commercial meat rabbit breeds (Italian Silver, Italian Spotted and Italian White) and 12 fancy rabbit breeds (Belgian Hare, Burgundy Fawn, Champagne d’Argent, Checkered Giant, Coloured Dwarf, Dwarf Lop, Ermine, Giant Grey, Giant White, Rex, Rhinelander and Thuringian) by using high-density single nucleotide polymorphism data. Signatures of selection were identified based on the fixation index (F_ST) statistic with different approaches, including single-breed and group-based methods, the latter comparing breeds that are grouped based on external traits (different coat colours and body sizes) and types (i.e. meat vs. fancy breeds).

Results

We identified 309 genomic regions that contained signatures of selection and that included genes that are known to affect coat colour (ASIP, MC1R and TYR), coat structure (LIPH), and body size (LCORL/NCAPG, COL11A1 and HOXD) in rabbits and that characterize the investigated breeds. Their identification proves the suitability of the applied methodologies for capturing recent selection events. Other regions included novel candidate genes that might contribute to the phenotypic variation among the analyzed breeds, including genes for pigmentation-related traits (EDNRA, EDNRB, MITF and OCA2) and body size, with a strong candidate for dwarfism in rabbit (COL2A1).

Conclusions

We report a genome-wide view of genetic loci that underlie the main phenotypic differences in the analyzed rabbit breeds, which can be useful to understand the shift from the domestication process to the development of breeds in O. cuniculus. These results enhance our knowledge about the major genetic loci involved in rabbit external traits and add novel information to understand the complexity of the genetic architecture underlying body size in mammals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00696-9.

Forensic DNA phenotyping: Canis familiaris breed classification and skeletal phenotype prediction using functionally significant skeletal SNPs and indels

This review highlights a novel application of breed identification and prediction of skeletal traits in forensic investigations using canine DNA evidence. Currently, genotyping methods used for canine breed classification involve the application of highly polymorphic short tandem repeats in addition to larger commercially available SNP arrays. Both applications face technical challenges. An additional approach to breed identification could be through genotyping SNPs and indels that characterise the array of skeletal differences displayed across domestic dog populations. Research has shown that a small number of genetic variants of large effect drive differences in skeletal phenotypes among domestic dog breeds. This feature makes functionally significant canine skeletal variants a cost-effective target for forensic investigators to classify individuals according to their breed. Further analysis of these skeletal variants would enable the prediction of external appearance. To date, functionally significant genes with genetic variants associated with differences in size, bulk, skull shape, ear shape, limb length, digit type, and tail morphology have been uncovered. Recommendations of a cost-effective genotyping method that can be readily designed and applied by forensic investigators have been given. Further advances to improve the field of canine skeletal forensic DNA phenotyping include the refinement of phenotyping methods, further biological validation of the skeletal genetic variants and establishing a publicly available database for storage of allele frequencies of the skeletal genetic variants in the wider domestic dog population.

Genomic Traces of the Fruit Fly Anastrepha obliqua Associated with Its Polyphagous Nature

Anastrepha obliqua (Macquart) (Diptera: Tephritidae) is an important pest in the neotropical region. It is considered a polyphagous insect, meaning it infests plants of different taxonomic families and readily colonizes new host plants. The change to new hosts can lead to diversification and the formation of host races. Previous studies investigating the effect of host plants on population structure and selection in Anastrepha obliqua have focused on the use of data from the mitochondrial DNA sequence and microsatellite markers of nuclear DNA, and there are no analyses at the genomic level. To better understand this issue, we used a pooled restriction site-associated DNA sequencing (pooled RAD-seq) approach to assess genomic differentiation and population structure across sympatric populations of Anastrepha obliqua that infest three host plants-Spondias purpurea (red mombin), Mangifera indica (mango) of the family Anacardiaceae and Averrhoa carambola (carambola) of the family Oxalidaceae-in sympatric populations of the species Anastrepha obliqua of Inter-Andean Valley of the Cauca River in southwestern Colombia. Our results show genomic differentiation of populations from carambola compared to mango and red mombin populations, but the genetic structure was mainly established by geography rather than by the host plant. On the other hand, we identified 54 SNPs in 23 sequences significantly associated with the use of the host plant. Of these 23 sequences, we identified 17 candidate genes and nine protein families, of which four protein families are involved in the nutrition of these flies. Future studies should investigate the adaptive processes undergone by phytophagous insects in the Neotropics, using fruit flies as a model and state-of-the-art molecular tools.

A single-nucleotide mutation within the TBX3 enhancer increased body size in Chinese horses

Chinese ponies are endemic to the mountainous areas of southwestern China and were first reported in the archaeological record at the Royal Tomb of Zhongshan King, Mancheng, dated to approximately ∼2,100 YBP.¹ Previous work has started uncovering the genetic basis of size variation in western ponies and horses, revealing a limited number of loci, including HMGA2,²LCORL/NCAPG,³ZFAT, and LASP1.^4,5 Whether the same genetic pathways also drive the small body size of Chinese ponies, which show striking anatomical differences to Shetland ponies,⁶ remains unclear.^2,7 To test this, we combined whole-genome sequences of 187 horses across China. Statistical analyses revealed top association between genetic variation at the T-box transcription factor 3 (TBX3) and the body size. Fine-scale analysis across an extended population of 189 ponies and 574 horses narrowed down the association to one A/G SNP at an enhancer region upstream of the TBX3 (ECA8:20,644,555, p = 2.34e−39). Luciferase assays confirmed the single-nucleotide G mutation upregulating TBX3 expression, and enhancer-knockout mice exhibited shorter limbs than wild-type littermates (p < 0.01). Re-analysis of ancient DNA data showed that the G allele, which is most frequent in modern horses, first occurred some ∼2,300 years ago and rose in frequency since. This supports selection for larger size in Asia from approximately the beginning of the Chinese Empire. Overall, this study characterized the causal regulatory mutation underlying small body size in Chinese ponies and revealed size as one of the main selection targets of past Chinese breeders.

•

One single A/G SNP in TBX3 enhancer region drives size variation in Chinese horses

•

The frequency of the G variant correlates positively with size in 763 horses

•

Cellular and mice models confirm it affects TBX3 transcription and the limb length

•

The G variant first occurred ∼2,300 years ago and rose in frequency since

Ethics in canine reproduction - a survey among veterinarians who provide canine reproductive services

OBJECTIVE

Several dog breeds suffer from health and welfare problems due to unfavorable anatomies and high prevalence of inherited diseases. Veterinarians should care for animal health and welfare but they may also be involved in breeding management including the use of reproductive techniques.

MATERIAL AND METHODS

An online survey was conducted to find out how veterinarians who are involved in small animal reproduction regard ethical issues, which services they provide, how they discuss and handle ethical concerns and talk about them with breeders.

RESULTS

In total, 83 participants mainly from Europe completed the online survey. Opinions on which interventions are ethical or not were heterogeneous. Most respondents (79.2 %) found it ethical to perform artificial insemination (AI) in dogs which did not reproduce naturally before. However, 62.7 % of participants stated that surgical insemination is not ethical. Elective cesarean sections are considered not ethical by approximately half of the participants. A closer look at responses to the questions of whether AI is ethical under specific conse quences suggests that most veterinarians are making ethical and value judgments in terms of what is justifiable or not. However, more than 80 % feel at least sometimes compromised by ethical conflicts.

CLINICAL RELEVANCE

This research intended to foster discussion of this topic. International strategies should be developed to support solving conflicts and dilemmas related to ethical issues in dog breeding and assisted reproduction, as well as diminishing physical disabilities and heritable diseases.

Advancing Genetic Selection and Behavioral Genomics of Working Dogs Through Collaborative Science

The ancient partnership between people and dogs is struggling to meet modern day needs, with demand exceeding our capacity to safely breed high-performing and healthy dogs. New statistical genetic approaches and genomic technology have the potential to revolutionize dog breeding, by transitioning from problematic phenotypic selection to methods that can preserve genetic diversity while increasing the proportion of successful dogs. To fully utilize this technology will require ultra large datasets, with hundreds of thousands of dogs. Today, dog breeders struggle to apply even the tools available now, stymied by the need for sophisticated data storage infrastructure and expertise in statistical genetics. Here, we review recent advances in animal breeding, and how a new approach to dog breeding would address the needs of working dog breeders today while also providing them with a path to realizing the next generation of technology. We provide a step-by-step guide for dog breeders to start implementing estimated breeding value selection in their programs now, and we describe how genotyping and DNA sequencing data, as it becomes more widely available, can be integrated into this approach. Finally, we call for data sharing among dog breeding programs as a path to achieving a future that can benefit all dogs, and their human partners too.

Interface of Human/Wildlife Interactions: An Example of a Bold Coyote (Canis latrans) in Atlanta, GA, USA

There is arguably no other North American species that better illustrates the complexities of the human-wildlife interface than the coyote. In this study, a melanistic coyote in metropolitan Atlanta, Georgia was exhibiting unusually bold behaviors that included encounters with humans, domestic dogs, and attempts to enter homes. After tracking this coyote (nicknamed Carmine) across a highly urbanized landscape with participatory science, including at least 80 publicly reported sightings, he was captured and relocated to a wildlife sanctuary. Genome-wide analyses revealed 92.8% coyote ancestry, 1.7% gray wolf ancestry, and 5.5% domestic dog ancestry. The dog alleles in Carmine’s genome were estimated to have been acquired by his ancestors 14–29 years ago. Despite his bold behavior, Carmine did not carry any mutations known to shape hypersociability in canines. He did, however, carry a single copy of the dominant mutation responsible for his melanistic coat color. This detailed study of Carmine dispels common assumptions about the reticent coyote personality and the origins of behavior. His unusual bold behavior created a higher level of human-coyote interaction. He now serves as a public ambassador for human-wildlife coexistence, urging the global community to reconsider mythologies about wildlife and promote coexistence with them in landscapes significantly altered by human activity in our rapidly changing world.

Impact of natural selection on global patterns of genetic variation, and association with clinical phenotypes, at genes involved in SARS-CoV-2 infection

We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (ACE2, TMPRSS2, DPP4, and LY6E). We analyzed novel data from 2,012 ethnically diverse Africans and 15,997 individuals of European and African ancestry with electronic health records, and integrated with global data from the 1000GP. At ACE2, we identified 41 non-synonymous variants that were rare in most populations, several of which impact protein function. However, three non-synonymous variants were common among Central African hunter-gatherers from Cameroon and are on haplotypes that exhibit signatures of positive selection. We identify strong signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19.

Impact of natural selection on global patterns of genetic variation, and association with clinical phenotypes, at genes involved in SARS-CoV-2 infection

We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection ( ACE2, TMPRSS2, DPP4 , and LY6E ). We analyzed novel data from 2,012 ethnically diverse Africans and 15,997 individuals of European and African ancestry with electronic health records, and integrated with global data from the 1000GP. At ACE2 , we identified 41 non-synonymous variants that were rare in most populations, several of which impact protein function. However, three non-synonymous variants were common among Central African hunter-gatherers from Cameroon and are on haplotypes that exhibit signatures of positive selection. We identify strong signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2 , we identified 13 amino acid changes that are adaptive and specific to the human lineage. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19.

Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.

Toward Genome-Based Selection in Asian Seabass: What Can We Learn From Other Food Fishes and Farm Animals?

Due to the steadily increasing need for seafood and the plateauing output of fisheries, more fish need to be produced by aquaculture production. In parallel with the improvement of farming methods, elite food fish lines with superior traits for production must be generated by selection programs that utilize cutting-edge tools of genomics. The purpose of this review is to provide a historical overview and status report of a selection program performed on a catadromous predator, the Asian seabass (Lates calcarifer, Bloch 1790) that can change its sex during its lifetime. We describe the practices of wet lab, farm and lab in detail by focusing onto the foundations and achievements of the program. In addition to the approaches used for selection, our review also provides an inventory of genetic/genomic platforms and technologies developed to (i) provide current and future support for the selection process; and (ii) improve our understanding of the biology of the species. Approaches used for the improvement of terrestrial farm animals are used as examples and references, as those processes are far ahead of the ones used in aquaculture and thus they might help those working on fish to select the best possible options and avoid potential pitfalls.

Copy number variation underlies complex phenotypes in domestic dog breeds and other canids

Extreme phenotypic diversity, a history of artificial selection, and socioeconomic value make domestic dog breeds a compelling subject for genomic research. Copy number variation (CNV) is known to account for a significant part of inter-individual genomic diversity in other systems. However, a comprehensive genome-wide study of structural variation as it relates to breed-specific phenotypes is lacking. We have generated whole genome CNV maps for more than 300 canids. Our data set extends the canine structural variation landscape to more than 100 dog breeds, including novel variants that cannot be assessed using microarray technologies. We have taken advantage of this data set to perform the first CNV-based genome-wide association study (GWAS) in canids. We identify 96 loci that display copy number differences across breeds, which are statistically associated with a previously compiled set of breed-specific morphometrics and disease susceptibilities. Among these, we highlight the discovery of a long-range interaction involving a CNV near MED13L and TBX3, which could influence breed standard height. Integration of the CNVs with chromatin interactions, long noncoding RNA expression, and single nucleotide variation highlights a subset of specific loci and genes with potential functional relevance and the prospect to explain trait variation between dog breeds.

Introgression, admixture, and selection facilitate genetic adaptation to high-altitude environments in cattle

Domestication and subsequent selection of cattle to form breeds and biological types that can adapt to different environments partitioned ancestral genetic diversity into distinct modern lineages. Genome-wide selection particularly for adaptation to extreme environments left detectable signatures genome-wide. We used high-density genotype data for 42 cattle breeds and identified the influence of Bos grunniens and Bos javanicus on the formation of Chinese indicine breeds that led to their divergence from India-origin zebu. We also found evidence for introgression, admixture, and migration in most of the Chinese breeds. Selection signature analyses between high-altitude (≥1800 m) and low-altitude adapted breeds (<1500 m) revealed candidate genes (ACSS2, ALDOC, EPAS1, EGLN1, NUCB2) and pathways that are putatively involved in hypoxia adaptation. Immunohistochemical, real-time PCR and CRISPR/cas9 ACSS2-knockout analyses suggest that the up-regulation of ACSS2 expression in the liver promotes the metabolic adaptation of cells to hypoxia via the hypoxia-inducible factor pathway. High altitude adaptation involved the introgression of alleles from high-altitude adapted yaks into Chinese Bos taurus taurus prior to their formation into recognized breeds and followed by selection. In addition to selection, adaptation to high altitude environments has been facilitated by admixture and introgression with locally adapted cattle populations.

Biosynthesis and Cellular Functions of Tartaric Acid in Grapevines

Tartaric acid (TA) is an obscure end point to the catabolism of ascorbic acid (Asc). Here, it is proposed as a "specialized primary metabolite", originating from carbohydrate metabolism but with restricted distribution within the plant kingdom and lack of known function in primary metabolic pathways. Grapes fall into the list of high TA-accumulators, with biosynthesis occurring in both leaf and berry. Very little is known of the TA biosynthetic pathway enzymes in any plant species, although recently some progress has been made in this space. New technologies in grapevine research such as the development of global co-expression network analysis tools and genome-wide association studies, should enable more rapid progress. There is also a lack of information regarding roles for this organic acid in plant metabolism. Therefore this review aims to briefly summarize current knowledge about the key intermediates and enzymes of TA biosynthesis in grapes and the regulation of its precursor, ascorbate, followed by speculative discussion around the potential roles of TA based on current knowledge of Asc metabolism, TA biosynthetic enzymes and other aspects of fruit metabolism.

Patterns of allele frequency differences among domestic cat breeds assessed by a 63K SNP array

Cats are ubiquitous companion animals that have been keenly associated with humans for thousands of years and only recently have been intentionally bred for aesthetically appealing coat looks and body forms. The intense selection on single gene phenotypes and the various breeding histories of cat breeds have left different marks on the genomes. Using a previously published 63K Feline SNP array dataset of twenty-six cat breeds, this study utilized a genetic differentiation-based method (di) to empirically identify candidate regions under selection. Defined as three or more overlapping (500Kb) windows of high levels of population differentiation, we identified a total of 205 candidate regions under selection across cat breeds with an average of 6 candidate regions per breed and an average size of 1.5 Mb per candidate region. Using the combined size of candidate regions of each breed, we conservatively estimate that a minimum of ~ 0.1-0.7% of the autosomal genome is potentially under selection in cats. As positive controls and tests of our methodology, we explored the candidate regions of known breed-defining genes (e.g., FGF5 for longhaired breeds) and we were able to detect the genes within candidate regions, each in its corresponding breed. For breed specific exploration of candidate regions under selection, eleven representative candidate regions were found to encompass potential candidate genes for several phenotypes such as brachycephaly of Persian (DLX6, DLX5, DLX2), curled ears of American Curl (MCRIP2, PBX1), and body-form of Siamese and Oriental (ADGRD1), which encourages further molecular investigations. The current assessment of the candidate regions under selection is empiric and detailed analyses are needed to rigorously disentangle effects of demography and population structure from artificial selection.

Genome-Scale Phylogenetic and Population Genetic Studies Provide Insight Into Introgression and Adaptive Evolution of Takifugu Species in East Asia

As a typical marine adaptive radiation species, most Takifugu species are widely distributed in East Asian offshore, which have diversified morphological characteristics and different ecological habits. The phylogenetic relationship and population structure of the Takifugu species was complicated because of incomplete lineage sorting, widespread hybridization and introgression. Therefore, to systematically clarify the phylogenetic relationships of Takifugu genus, explore the introgression and natural hybridization between different Takifugu species, and detect the selective signatures in the adaptive evolution of diversified traits, whole-genome resequencing was used in 122 Takifugu samples from 10 species. Phylogenetic analysis showed solid sister-group relationships between Takifugu bimaculatus and Takifugu flavidus, Takifugu oblongus, and Takifugu niphobles, Takifugu rubripes, and Takifugu obscurus, Takifugu xanthoptreus, and Takifugu ocellatus. Further admixture analysis indicated the divergence of T. obscurus population and the bidirectional gene flow between T. bimaculatus and T. flavidus. Using species-specific homozygous genetic variance sites, we detected the asymmetric introgression between T. bimaculatus and T. flavidus at China East sea and southern Taiwan Strait. By genome-scale genetic diversity scanning, we detected two copies of syt1, zar1 and tgfbr1 related to the semilunar reproduction rhythm in T. niphobles, involved in memory formation, embryo maturation and female reproduction. Furthermore, we also found lots of T. niphobles specific mutations in CDS region of circadian rhythm related genes and endocrine hormone genes. For Takifugu species, our research provides reliable genetic resources and results for the phylogeny, introgression, hybridization and adaptive evolution, and could be used as a guide for the formulation of the protection and proliferation release policies.