Genome-wide Heterozygosity Research Articles

Strongly deleterious mutations influence reproductive output and longevity in an endangered population

Inbreeding depression has been documented in various fitness traits in a wide range of species and taxa, however, the mutational basis is not yet well understood. We investigate how putatively deleterious variation influences fitness and is shaped by individual ancestry by re-sequencing complete genomes of 37 individuals in a natural arctic fox (Vulpes lagopus) population subjected to both inbreeding depression and genetic rescue. We find that individuals with high proportion of homozygous loss of function genotypes (LoFs), which are predicted to exert a strong effect on fitness, generally have lower lifetime reproductive success and live shorter lives compared with individuals with lower proportion of LoFs. We also find that juvenile survival is negatively associated with the proportion of homozygous missense genotypes and positively associated with genome wide heterozygosity. Our results demonstrate that homozygosity of strongly and moderately deleterious mutations can be an important cause of trait specific inbreeding depression in wild populations, and mark an important step towards making more informed decisions using applied conservation genetics.

Chromosomal-level reference genome assembly of muskox (Ovibos moschatus) from Banks Island in the Canadian Arctic, a resource for conservation genomics

The muskox (Ovibos moschatus), an integral component and iconic symbol of arctic biocultural diversity, is under threat by rapid environmental disruptions from climate change. We report a chromosomal-level haploid genome assembly of a muskox from Banks Island in the Canadian Arctic Archipelago. The assembly has a contig N50 of 44.7 Mbp, a scaffold N50 of 112.3 Mbp, a complete representation (100%) of the BUSCO v5.2.2 set of 9225 mammalian marker genes and is anchored to the 24 chromosomes of the muskox. Tabulation of heterozygous single nucleotide variants in our specimen revealed a very low level of genetic diversity, which is consistent with recent reports of the muskox having the lowest genome-wide heterozygosity among the ungulates. While muskox populations are currently showing no overt signs of inbreeding depression, environmental disruptions are expected to strain the genomic resilience of the species. One notable impact of rapid climate change in the Arctic is the spread of emerging infectious and parasitic diseases in the muskox, as exemplified by the range expansion of muskox lungworms, and the recent fatal outbreaks of Erysipelothrix rhusiopathiae, a pathogen normally associated with domestic swine and poultry. As a genomics resource for conservation management of the muskox against existing and emerging disease modalities, we annotated the genes of the major histocompatibility complex on chromosome 2 and performed an initial assessment of the genetic diversity of this complex. This resource is further supported by the annotation of the principal genes of the innate immunity system, genes that are rapidly evolving and under positive selection in the muskox, genes associated with environmental adaptations, and the genes associated with socioeconomic benefits for Arctic communities such as wool (qiviut) attributes. These annotations will benefit muskox management and conservation.

Genetic and telomeric variability: Insights from a tropical avian hybrid zone.

Telomere lengths and telomere dynamics can correlate with lifespan, behaviour and individual quality. Such relationships have spurred interest in understanding variation in telomere lengths and their dynamics within and between populations. Many studies have identified how environmental processes can influence telomere dynamics, but the role of genetic variation is much less well characterized. To provide a novel perspective on how telomeric variation relates to genetic variability, we longitudinally sampled individuals across a narrow hybrid zone (n = 127 samples), wherein two Manacus species characterized by contrasting genome-wide heterozygosity interbreed. We measured individual (n = 66) and population (n = 3) differences in genome-wide heterozygosity and, among hybrids, amount of genetic admixture using RADseq-generated SNPs. We tested for population differences in telomere lengths and telomere dynamics. We then examined how telomere lengths and telomere dynamics covaried with genome-wide heterozygosity within populations. Hybrid individuals exhibited longer telomeres, on average, than individuals sampled in the adjacent parental populations. No population differences in telomere dynamics were observed. Within the parental population characterized by relatively low heterozygosity, higher genome-wide heterozygosity was associated with shorter telomeres and higher rates of telomere shortening-a pattern that was less apparent in the other populations. All of these relationships were independent of sex, despite the contrasting life histories of male and female manakins. Our study highlights how population comparisons can reveal interrelationships between genetic variation and telomeres, and how naturally occurring hybridization and genome-wide heterozygosity can relate to telomere lengths and telomere dynamics.

The last days of Aporia crataegi (L.) in Britain: Evaluating genomic erosion in an extirpated butterfly.

Current rates of habitat degradation and climate change are causing unprecedented declines in global biodiversity. Studies on vertebrates highlight how conservation genomics can be effective in identifying and managing threatened populations, but it is unclear how vertebrate-derived metrics of genomic erosion translate to invertebrates, with their markedly different population sizes and life histories. The Black-veined White butterfly (Aporia crataegi) was extirpated from Britain in the 1920s. Here, we sequenced historical DNA from 17 specimens collected between 1854 and 1924 to reconstruct demography and compare levels of genomic erosion between extirpated British and extant European mainland populations. We contrast these results using modern samples of the Common Blue butterfly (Polyommatus icarus); a species with relatively stable demographic trends in Great Britain. We provide evidence for bottlenecks in both these species around the period of post-glacial colonization of the British Isles. Our results reveal different demographic histories and Ne for both species, consistent with their fates in Britain, likely driven by differences in life history, ecology and genome size. Despite a difference, by an order of magnitude, in historical effective population sizes (Ne), reduction in genome-wide heterozygosity in A. crataegi was comparable to that in P. icarus. Symptomatic of A. crataegi's disappearance were marked increases in runs-of-homozygosity (RoH), potentially indicative of recent inbreeding, and accumulation of putatively mildly and weakly deleterious variants. Our results provide a rare glimpse of genomic erosion in a regionally extinct insect and support the potential use of genomic erosion metrics in identifying invertebrate populations or species in decline.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity.

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

P-560 Characterization of ploidy level and parental genome inheritance abnormalities in preimplantation human embryos

Abstract Study question Can preimplantation genetic testing (PGT) platforms with integrated genotyping aid the characterization of genome-wide chromosomal abnormalities in human embryos and factors contributing to their occurrence? Summary answer Signatures detected by genotype-criteria-based models in PGT pinpoint origin of haploidy, triploidy and genome-wide uniparental disomy (gwUPD), highlighting recombination deficiencies potentially underpinning these phenomena. What is known already While there is a growing body of knowledge on the frequency and origin of individual chromosome abnormalities, very little is known about genome-wide chromosomal abnormalities in human development. Indeed, genome-wide abnormalities are incompatible with normal human development, making the preimplantation stage the ideal time to characterize them. Novel PGT technologies provide the research and clinical opportunity to characterize both known ploidy-level abnormalities (triploid/haploid) and identify theoretical genome-wide abnormalities undetectable at later stages of development. Study design, size, duration PGT genotyping data from 96660 embryos were used to assess incidence of ploidy-level abnormalities and correlation with potential predisposing clinical factors. Parental origin of putative haploidy and triploidy was assessed on 14 and 41 embryos with parental genotypes available (trios), respectively. For the 41 triploid trios, meiotic origin and recombination rates were also evaluated. A novel genotype-criteria based model to detect theoretical gwUPD and meiotic origin was developed on a semi-contemporaneous 74009 female embryo dataset. Participants/materials, setting, methods PGT was performed with targeted next-generation sequencing, including 5000 genome-wide SNPs. Informative SNPs from trios were used to identify parental origin of haploidy/triploidy. Informative SNPs within 5cM of the centromere were used to evaluate meiotic origin and extended beyond 5cM to identify maternal recombination observing shifts between single (SPH) and biparental (BPH) homologs. gwUPD was predicted when low genome-wide heterozygosity, increased homozygosity in the centromeric region, and decreased similarity to sibling embryos were observed. Main results and the role of chance Genome wide ploidy abnormalities were detected in 1.1% (n = 1.063/96.660; CI 95%: 1-1,2%) of embryos. Haploidy was mainly of paternal origin (n = 12/14), whereas triploidy was exclusively maternal in origin (n = 41/41) and prominently due to meiosis II errors (n = 30/41; 73.2%; 95% CI:57.1-85.8%). Advanced maternal age was positively associated with an increased risk of triploid conceptions (OR = 1.059 per year; p < 0.001). Investigation of maternal recombination rates detected by three consecutive SNPs shifting between SPH and BPH (with a calculated genotyping error of 0.16%) in triploid embryos revealed an outlier group of 7 embryos with statistically lower (P < 0.001) identified crossover events (0-2 events genome-wide) compared to the average 23.7 detectable events genome-wide (SD:4.3). Recombinant gwUPD was predicted in 41 out of 74009 human embryos (0.05%; 95% CI:0.04-0.07%). These embryos showed a high homozygosity in pericentromeric regions, as a proof of each pair of chromosomes being the same parental homolog resulting from a faulty meiosis II, but with detectable heterozygosity regions throughout the rest of the genome, indicating that crossovers occurred between the two parental homologs during prophase I. Multifocal biopsies from the same embryos confirmed the reproducibility of identification of the number of recombination events in triploid embryos and the gwUPD prediction. Limitations, reasons for caution This assay genotypes only a portion of the genome and underestimation of recombination events is possible. Our current detection of gwUPD is limited to recombinant uniparental isodisomy of MII origin only. Therefore, the exact prevalence of these abnormal gwUPD configurations requires further understanding and analysis. Wider implications of the findings Genotyping in routine PGT can identify severe genome-wide abnormalities with associated gestational complications that would otherwise go undiagnosed. Comprehensive characterization of these abnormalities is not possible at later stages of development. Therefore, PGT offers a unique opportunity to delve into their biological origin and identify novel biomarkers for embryo competence. Trial registration number not applicable

The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus).

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

A chromosome-level genome assembly for the dugong (Dugong dugon).

The dugong (Dugong dugon) is a marine mammal widely distributed throughout the Indo-Pacific and the Red Sea, with a Vulnerable conservation status, and little is known about many of the more peripheral populations, some of which are thought to be close to extinction. We present a de novo high-quality genome assembly for the dugong from an individual belonging to the well-monitored Moreton Bay population in Queensland, Australia. Our assembly uses long-read PacBio HiFi sequencing and Omni-C data following the Vertebrate Genome Project pipeline to reach chromosome-level contiguity (24 chromosome-level scaffolds; 3.16 Gbp) and high completeness (97.9% complete BUSCOs). We observed relatively high genome-wide heterozygosity, which likely reflects historical population abundance before the last interglacial period, approximately 125,000 yr ago. Demographic inference suggests that dugong populations began declining as sea levels fell after the last interglacial period, likely a result of population fragmentation and habitat loss due to the exposure of seagrass meadows. We find no evidence for ongoing recent inbreeding in this individual. However, runs of homozygosity indicate some past inbreeding. Our draft genome assembly will enable range-wide assessments of genetic diversity and adaptation, facilitate effective management of dugong populations, and allow comparative genomics analyses including with other sirenians, the oldest marine mammal lineage.

Genomic analysis supports Cape Lion population connectivity prior to colonial eradication and extinction.

Cape lions (Panthera leo melanochaitus) formerly ranged throughout the grassland plains of the "Cape Flats" in what is today known as the Western Cape Province, South Africa. Cape lions were likely eradicated because of overhunting and habitat loss after European colonization. European naturalists originally described Cape lions as "black-maned lions" and claimed that they were phenotypically distinct. However, other depictions and historical descriptions of lions from the Cape report mixed or light coloration and without black or extensively developed manes. These findings suggest that, rather than forming a distinct population, Cape lions may have had phenotypic and genotypic variation similar to other African lions. Here we investigate Cape lion genome characteristics, population dynamics, and genetic distinctiveness prior to their extinction. We generated genomic data from 2 historic Cape lions to compare to 118 existing high-coverage mitogenomes, and low-coverage nuclear genomes of 53 lions from 13 African countries. We show that, before their eradication, lions from the Cape Flats had diverse mitogenomes and nuclear genomes that clustered with lions from both southern and eastern Africa. Cape lions had high genome-wide heterozygosity and low inbreeding coefficients, indicating that populations in the Cape Flats went extinct so rapidly that genomic effects associated with long-term small population size and isolation were not detectable. Our findings do not support the characterization of Cape lions as phylogeographically distinct, as originally put forth by some European naturalists, and illustrates how alternative knowledge systems, for example, Indigenous perspectives, could potentially further inform interpretations of species histories.

Heterogeneous patterns of heterozygosity loss in isolated populations of the threatened eastern barred bandicoot (Perameles gunnii).

Identifying and analysing isolated populations is critical for conservation. Isolation can make populations vulnerable to local extinction due to increased genetic drift and inbreeding, both of which should leave imprints of decreased genome-wide heterozygosity. While decreases in heterozygosity among populations are frequently investigated, fewer studies have analysed how heterozygosity varies among individuals, including whether heterozygosity varies geographically along lines of discrete population structure or with continuous patterns analogous to isolation by distance. Here we explore geographical patterns of differentiation and individual heterozygosity in the threatened eastern barred bandicoot (Perameles gunnii) in Tasmania, Australia, using genomic data from 85 samples collected between 2008 and 2011. Our analyses identified two isolated demes undergoing significant genetic drift, and several areas of fine-scale differentiation across Tasmania. We observed discrete genetic structures across geographical barriers and continuous patterns of isolation by distance, with little evidence of recent or historical migration. Using a recently developed analytical pipeline for estimating autosomal heterozygosity, we found individual heterozygosities varied within demes by up to a factor of two, and demes with low-heterozygosity individuals also still contained those with high heterozygosity. Spatial interpolation of heterozygosity scores clarified these patterns and identified the isolated Tasman Peninsula as a location where low-heterozygosity individuals were more common than elsewhere. Our results provide novel insights into the relationship between isolation-driven genetic structure and local heterozygosity patterns. These may help improve translocation efforts, by identifying populations in need of assistance, and by providing an individualised metric for identifying source animals for translocation.

Population genomics of the muskox' resilience in the near absence of genetic variation.

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Pangolin Genomes Offer Key Insights and Resources for the World's Most Trafficked Wild Mammals.

Pangolins form a group of scaly mammals that are trafficked at record numbers for their meat and purported medicinal properties. Despite their conservation concern, knowledge of their evolution is limited by a paucity of genomic data. We aim to produce exhaustive genomic resources that include 3,238 orthologous genes and whole-genome polymorphisms to assess the evolution of all eight extant pangolin species. Robust orthologous gene-based phylogenies recovered the monophyly of the three genera and highlighted the existence of an undescribed species closely related to Southeast Asian pangolins. Signatures of middle Miocene admixture between an extinct, possibly European, lineage and the ancestor of Southeast Asian pangolins, provide new insights into the early evolutionary history of the group. Demographic trajectories and genome-wide heterozygosity estimates revealed contrasts between continental versus island populations and species lineages, suggesting that conservation planning should consider intraspecific patterns. With the expected loss of genomic diversity from recent, extensive trafficking not yet realized in pangolins, we recommend that populations be genetically surveyed to anticipate any deleterious impact of the illegal trade. Finally, we produce a complete set of genomic resources that will be integral for future conservation management and forensic endeavors for pangolins, including tracing their illegal trade. These comprise the completion of whole-genomes for pangolins through the hybrid assembly of the first reference genome for the giant pangolin (Smutsia gigantea) and new draft genomes (∼43x-77x) for four additional species, as well as a database of orthologous genes with over 3.4 million polymorphic sites.

Assessing the genetic composition of cotton-top tamarins (Saguinus oedipus) before sweeping anthropogenic impact.

During the last century, the critically endangered cotton-top tamarin (Saguinus oedipus) has been threatened by multiple anthropogenic factors that drastically affected their habitat and population size. As the genetic impact of these pressures is largely unknown, this study aimed to establish a genetic baseline with the use of temporal sampling to determine the genetic makeup before detrimental anthropogenic impact. Genomes were resequenced from a combination of historical museum samples and modern wild samples at low-medium coverage, to unravel how the cotton-top tamarin population structure and genomic diversity may have changed during this period. Our data suggest two populations can be differentiated, probably separated historically by the mountain ranges of the Paramillo Massif in Colombia. Although this population structure persists in the current populations, modern samples exhibit genomic signals consistent with recent inbreeding, such as long runs of homozygosity and a reduction in genome-wide heterozygosity especially in the greater northeast population. This loss is likely the consequence of the population reduction following the mass exportation of cotton-top tamarins for biomedical research in the 1960s, coupled with the habitat loss this species continues to experience. However, current populations have not experienced an increase in genetic load. We propose that the historical genetic baseline established in this study can be used to provide insight into alteration in the modern population influenced by a drastic reduction in population size as well as providing background information to be used for future conservation decision-making for the species.

Contrasting genomic consequences of anthropogenic reintroduction and natural recolonization in high-arctic wild reindeer.

Anthropogenic reintroduction can supplement natural recolonization in reestablishing a species' distribution and abundance. However, both reintroductions and recolonizations can give rise to founder effects that reduce genetic diversity and increase inbreeding, potentially causing the accumulation of genetic load and reduced fitness. Most current populations of the endemic high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) originate from recent reintroductions or recolonizations following regional extirpations due to past overharvesting. We investigated and compared the genomic consequences of these two paths to reestablishment using whole-genome shotgun sequencing of 100 Svalbard reindeer across their range. We found little admixture between reintroduced and natural populations. Two reintroduced populations, each founded by 12 individuals around four decades (i.e. 8 reindeer generations) ago, formed two distinct genetic clusters. Compared to the source population, these populations showed only small decreases in genome-wide heterozygosity and increases in inbreeding and lengths of runs of homozygosity. In contrast, the two naturally recolonized populations without admixture possessed much lower heterozygosity, higher inbreeding and longer runs of homozygosity, possibly caused by serial population founder effects and/or fewer or more genetically related founders than in the reintroduction events. Naturally recolonized populations can thus be more vulnerable to the accumulation of genetic load than reintroduced populations. This suggests that in some organisms even small-scale reintroduction programs based on genetically diverse source populations can be more effective than natural recolonization in establishing genetically diverse populations. These findings warrant particular attention in the conservation and management of populations and species threatened by habitat fragmentation and loss.

The impact of habitat loss and population fragmentation on genomic erosion

Habitat loss and population fragmentation pose severe threats to biodiversity and the survival of many species. Population isolation and the decline in effective population size lead to increased genetic drift and inbreeding. In turn, this reduces neutral diversity, and it also affects the genetic load of deleterious mutations. Here, we analyse the effect of such genomic erosion by designing a spatially explicit, individual based model in SLiM, simulating the effects of the recorded habitat loss in Mauritius over the past ~ 250 years. We show that the loss of neutral diversity (genome-wide heterozygosity) was barely noticeable during the first 100 years of habitat loss. Changes to the genetic load took even more time to register, and they only became apparent circa 200 years after the start of habitat decline. Although a considerable number of deleterious mutations were lost by drift, others increased in frequency. The masked load was thus converted into a realised load, which compromised individual fitness and population viability after much of the native habitat had been lost. Importantly, genomic erosion continued after the metapopulation had stabilised at low numbers. Our study shows that historic habitat loss can pose a sustained threat to populations also in future generations, even without further habitat loss. The UN’s Decade on Ecosystem Restoration needs to lead to transformative change to save species from future extinction, and this requires the urgent restoration of natural habitats.

Genome-wide population structure, homozygosity, and heterozygosity patterns of Nero Siciliano pig in the framework of Italian and cosmopolitan breeds.

Analysis of genomic data is becoming more and more common for the effective management of livestock breeding programmes, even in the case of local populations. In this work, the genome-wide data of Nero Siciliano pig breed were compared to that of wild boar, Italian local and cosmopolitan breeds to investigate its genetic structure, and runs of homozygosity (ROH) and heterozygosity patterns. The Nero Siciliano has been reported to have the highest rate of genetic diversity among the Italian breeds, and a genetic variability comparable to that of the cosmopolitan breeds. Analyses of genomic structure and relationships underlined its proximity to wild boar, and an internal substructure probably linked to different family lines. The breed showed a low value of inbreeding estimated from ROH, and the highest diversity index among the Italian breeds, even if lower than that of the cosmopolitans. Four ROH islands in three chromosomes (SSC8, SSC11, and SSC14) and one heterozygosity-rich region (SSC1) were identified in Nero Siciliano, highlighting genomic regions related to productive QTL. Across breeds, SSC8 and SSC14 were the chromosomes with most ROH islands, with Mora Romagnola and wild boar showing the highest level of autozygosity. Chromosomes SSC2, SSC6, SSC8 and SSC13 showed the majority of runs of heterozygosity regions, mainly found in the cosmopolitan pig breeds, which reported several genes associated with health-related QTL. The outlined results can help to better identify the genomic profile of this local breed in order to plan matings, maintain adequate internal diversity and exploit the production system.

Increased heterozygosity in low-pass sequencing data allows identification of blood chimeras in cattle.

In about 90% of multiple pregnancies in cattle, shared blood circulation between fetuses leads to genetic chimerism in peripheral blood and can reduce reproductive performance in heterosexual co-twins. However, the early detection of heterosexual chimeras requires specialized tests. Here, we used low-pass sequencing data with a median coverage of 0.64× generated from blood samples of 322 F1 crosses between beef and dairy cattle and identified 20 putative blood chimeras through increased levels of genome-wide heterozygosity. In contrast, for 77 samples with routine SNP microarray data generated from hair bulbs of the same F1s, we found no evidence of chimerism, simultaneously observing high levels of genotype discordance with sequencing data. Fifteen out of 18 reported twins showed signs of blood chimerism, in line with previous reports, whereas the presence of five alleged singletons with strong signs of chimerism suggests that the in-utero death rate of co-twins is at the upper limit of former estimates. Together, our results show that low-pass sequencing data allow reliable screening for blood chimeras. They further affirm that blood is not recommended as a source of DNA for the detection of germline variants.

The rough-toothed dolphin genome provides new insights into the genetic mechanism of its rough teeth.

The rough-toothed dolphin (Steno bredanensis) is characterized by having teeth covered in finely wrinkled vertical ridges, which is a general manifestation of amelogenesis imperfecta. The rough surfaces are hypothesized to be an evolutionary morphological trait of feeding adaptation to increase the dolphin's grip on prey. Here, we assembled a rough-toothed dolphin genome and performed the comparative genomic analysis to reveal the genetic basis of the special enamel. Results showed that genes related to enamel development or dental diseases have undergone diversified adaptive changes that may shape the special enamel morphology of this dolphin species, including positive selection (CLDN19, PRKCE, SSUH2, and WDR72), rapid evolution (LAMB3), or unique amino acid substitutions (AMTN, ENAM, MMP20, and KLK4). Meanwhile, the historical demography of rough-toothed dolphin indicated several distinct population fluctuations associated with climate change. The genome-wide heterozygosity of this dolphin is in the middle of all published data for cetaceans. Although the population is considerable, there may be population or subspecies differentiation, and with the global warming and the increasing disturbance of human activities, we should pay more attention to protection in the future. Together, our study brings new insights into the genetic mechanisms that may have driven the evolution of the special enamel morphology in rough-toothed dolphins and provides the first results of genetic heterozygosity and population historical dynamics of this species, which have important guiding implications for the conservation of this dolphin species.

Genome diversity of Leishmania aethiopica.

Leishmania aethiopica is a zoonotic Old World parasite transmitted by Phlebotomine sand flies and causing cutaneous leishmaniasis in Ethiopia and Kenya. Despite a range of clinical manifestations and a high prevalence of treatment failure, L. aethiopica is one of the most neglected species of the Leishmania genus in terms of scientific attention. Here, we explored the genome diversity of L. aethiopica by analyzing the genomes of twenty isolates from Ethiopia. Phylogenomic analyses identified two strains as interspecific hybrids involving L. aethiopica as one parent and L. donovani and L. tropica respectively as the other parent. High levels of genome-wide heterozygosity suggest that these two hybrids are equivalent to F1 progeny that propagated mitotically since the initial hybridization event. Analyses of allelic read depths further revealed that the L. aethiopica - L. tropica hybrid was diploid and the L. aethiopica - L. donovani hybrid was triploid, as has been described for other interspecific Leishmania hybrids. When focusing on L. aethiopica, we show that this species is genetically highly diverse and consists of both asexually evolving strains and groups of recombining parasites. A remarkable observation is that some L. aethiopica strains showed an extensive loss of heterozygosity across large regions of the nuclear genome, which likely arose from gene conversion/mitotic recombination. Hence, our prospection of L. aethiopica genomics revealed new insights into the genomic consequences of both meiotic and mitotic recombination in Leishmania.

Chromosome-length genome assemblies and cytogenomic analyses of pangolins reveal remarkable chromosome counts and plasticity.

We report the first chromosome-length genome assemblies for three species in the mammalian order Pholidota: the white-bellied, Chinese, and Sunda pangolins. Surprisingly, we observe extraordinary karyotypic plasticity within this order and, in female white-bellied pangolins, the largest number of chromosomes reported in a Laurasiatherian mammal: 2n = 114. We perform the first karyotype analysis of an African pangolin and report a Y-autosome fusion in white-bellied pangolins, resulting in 2n = 113 for males. We employ a novel strategy to confirm the fusion and identify the autosome involved by finding the pseudoautosomal region (PAR) in the female genome assembly and analyzing the 3D contact frequency between PAR sequences and the rest of the genome in male and female white-bellied pangolins. Analyses of genetic variability show that white-bellied pangolins have intermediate levels of genome-wide heterozygosity relative to Chinese and Sunda pangolins, consistent with two moderate declines of historical effective population size. Our results reveal a remarkable feature of pangolin genome biology and highlight the need for further studies of these unique and endangered mammals.