Complete Haplotype Research Articles

Strainy: phasing and assembly of strain haplotypes from long-read metagenome sequencing.

Bacterial species in microbial communities are often represented by mixtures of strains, distinguished by small variations in their genomes. Short-read approaches can be used to detect small-scale variation between strains but fail to phase these variants into contiguous haplotypes. Long-read metagenome assemblers can generate contiguous bacterial chromosomes but often suppress strain-level variation in favor of species-level consensus. Here we present Strainy, an algorithm for strain-level metagenome assembly and phasing from Nanopore and PacBio reads. Strainy takes a de novo metagenomic assembly as input and identifies strain variants, which are then phased and assembled into contiguous haplotypes. Using simulated and mock Nanopore and PacBio metagenome data, we show that Strainy assembles accurate and complete strain haplotypes, outperforming current Nanopore-based methods and comparable with PacBio-based algorithms in completeness and accuracy. We then use Strainy to assemble strain haplotypes of a complex environmental metagenome, revealing distinct strain distribution and mutational patterns in bacterial species.

A telomere-to-telomere gap-free reference genome of Chionanthus retusus provides insights into the molecular mechanism underlying petal shape changes

Abstract Chionanthus retusus, an arbor tree of the Oleaceae family, is an ecologically and economically valuable ornamental plant for its remarkable adaptability in landscaping. During C. retusus breeding, we observed diverse floral shapes; however, no available genome for C. retusus has hindered the widespread identification of genes related to flower morphology. Thus, a de novo telomere-to-telomere (T2T) gap-free genome was generated. The assembly, incorporating high-coverage and long-read sequencing data, successfully yielded two complete haplotypes (687 and 683 Mb). The genome encompasses 42 864 predicted protein-coding genes, with all 46 telomeres and 23 centromeres in one haplotype. Whole genome duplication analysis revealed that C. retusus underwent one fewer event of whole-genome duplication after differentiation compared to other species in the Oleaceae family. Furthermore, flower vein diversity was the main reason for the differences in floral shapes. Auxin-related genes were responsible for petal shape formation on genome-based transcriptome analysis. Specifically, the removal and retention of the first intron in CrAUX/IAA20 resulted in the production of two transcripts, and the differences in the expression levels of CrAUX/IAA20 resulted in the variations of flower veins. Compared to transcripts lacking the first intron, transcripts with intron retention caused more severe decreases in the number and length of flower veins in transgenic Arabidopsis thaliana. Our findings will deepen our understanding of flower morphology development and provide important theoretical support for the cultivation of Oleaceae.

DNA structural features and variability of complete MHC locus sequences.

The major histocompatibility (MHC) locus, also known as the Human Leukocyte Antigen (HLA) genes, is located on the short arm of chromosome 6, and contains three regions (Class I, Class II and Class III). This 5Mbp locus is one of the most variable regions of the human genome, yet it also encodes a set of highly conserved and important proteins related to immunological response. Genetic variations in this region are responsible for more diseases than in the entire rest of the human genome. However, information on local structural features of the DNA is largely ignored. With recent advances in long-read sequencing technology, it is now becoming possible to sequence the entire 5Mbp MHC locus, producing complete diploid haplotypes of the whole region. Here, we describe structural maps based on the complete sequences from six different homozygous HLA cell lines. We find long-range structural variability in the different sequences for DNA stacking energy, position preference and curvature, variation in repeats, as well as more local changes in regions forming open chromatin structures, likely to influence gene expression levels. These structural maps can be useful in visualizing large scale structural variation across HLA types, in particular when this can be complemented with epigenetic signals.

Population structure and haplotype network analyses of Hyalomma anatolicum based on the large subunit ribosomal RNA (16S rRNA) gene.

Hyalomma anatolicum, an Anatolian hard tick is a well-recognized vector involved in the transmission of various pathogens to animals and humans. The present study elucidated the population structure and haplotype network of H. anatolicum based on the mitochondrial large subunit ribosomal RNA (16S rRNA) gene sequence. The population structure and haplotype network analysis of 75 sequences archived in the GenBank, including the 15 sequences generated herein, yielded 24 haplotypes. Haplotype 1 (Hap_1) was the predominant haplotype consisting of 45 sequences from India, China, Pakistan, Turkey, Egypt, Iraq, and Tajikistan. The complete haplotype network exhibited a stellate conformation, highlighting a recent population expansion. The overall dataset, together with the sequences corresponding to India, China, and Pakistan, showed a high haplotype (0.638 ± 0.065, 0.671 ± 0.103, 0.753 ± 0.099, and 0.854 ± 0.061, respectively) and low nucleotide (0.00407 ± 0.00090, 0.00525 ± 0.00196, 0.00680 ± 0.00233, and 0.00453 ± 0.00056, respectively) diversity, further emphasized a recent population expansion. The neutrality indices including Tajima's D, Fu and Li's D, and Fu and Li's F for the complete dataset (- 2.661, - 6.008, and - 5.649, respectively) as well as for the sequences from India (- 2.223, - 3.414, and - 3.567, respectively) were negative, suggesting deviation from neutrality and a recent population expansion. The present study provided novel insights into the population structure and haplotype networks of H. anatolicum based on the mitochondrial 16S rRNA gene, and the different tests inferred a low genetic differentiation and suggested a recent population expansion of this economically important tick species.

A precise and cost-efficient whole-genome haplotyping method without probands: preimplantation genetic testing analysis

Research questionIs there a precise and efficient haplotyping method to expand the application on preimplantation genetic testing (PGT)? DesignIn this study, 8 cell line families and 18 clinical families including 99 embryos, were used to construct whole-genome haplotyping based on link-read sequencing and optimized analytical workflow with correction algorithm. Correction algorithm was based on differentiation of assembly errors and homologous recombination, in which the main feature of parental assembly error was that all embryos (embryo number >=2) had breakpoints at the same chromosome position. ResultsWith our method, parental assembly errors and homologous recombination were accurately distinguished and corrected, and using the link-reads (that ≥30 kilobases (KB) were more than 25%), complete genome-wide parental haplotypes were constructed, and the consistency of the typing results of each chromosome with conventional method requiring other family members was more than 95%. Besides, the length of N50 contigs was from 11.03 to 16.2 million bases (MB), which was far beyond to N50 contigs from long-read sequencing (148 KB∼863 KB). The complete haplotype analysis of all embryos could be exhibited by our method, and revealed 100% concordance with the available diagnostic results obtained by conventional method requiring other family members. ConclusionsOur method has high clinical value as a precise and cost-efficient whole-genome haplotyping method without probands as part of preimplantation genetic testing (PGT) and other genetic research, which could promote the application of PGT to decrease the birth of children with genetic diseases and the development of linkage-related genetic research.

Assembly and analysis of the genome of Notholithocarpus densiflorus.

Tanoak (Notholithocarpus densiflorus) is an evergreen tree in the Fagaceae family found in California and southern Oregon. Historically, tanoak acorns were an important food source for Native American tribes, and the bark was used extensively in the leather tanning process. Long considered a disjunct relictual element of the Asian stone oaks (Lithocarpus spp.), phylogenetic analysis has determined that the tanoak is an example of convergent evolution. Tanoaks are deeply divergent from oaks (Quercus) of the Pacific Northwest and comprise a new genus with a single species. These trees are highly susceptible to "sudden oak death" (SOD), a plant pathogen (Phytophthora ramorum) that has caused widespread deaths of tanoaks. In this study, we set out to assemble the genome and perform comparative studies among a number of individuals that demonstrated varying levels of susceptibility to SOD. First, we sequenced and de novo assembled a draft reference genome of N. densiflorus using cobarcoded library processing methods and an MGI DNBSEQ-G400 sequencer. To increase the contiguity of the final assembly, we also sequenced Oxford Nanopore long reads to 30× coverage. To our knowledge, the draft genome reported here is one of the more contiguous and complete genomes of a tree species published to date, with a contig N50 of ∼1.2 Mb, a scaffold N50 of ∼2.1 Mb, and a complete gene score of 95.5% through BUSCO analysis. In addition, we sequenced 11 genetically distinct individuals and mapped these onto the draft reference genome, enabling the discovery of almost 25 million single nucleotide polymorphisms and ∼4.4 million small insertions and deletions. Finally, using cobarcoded data, we were able to generate a complete haplotype coverage of all 11 genomes.

AsmMix: an efficient haplotype-resolved hybrid de novo genome assembling pipeline.

Accurate haplotyping facilitates distinguishing allele-specific expression, identifying cis-regulatory elements, and characterizing genomic variations, which enables more precise investigations into the relationship between genotype and phenotype. Recent advances in third-generation single-molecule long read and synthetic co-barcoded read sequencing techniques have harnessed long-range information to simplify the assembly graph and improve assembly genomic sequence. However, it remains methodologically challenging to reconstruct the complete haplotypes due to high sequencing error rates of long reads and limited capturing efficiency of co-barcoded reads. We here present a pipeline, AsmMix, for generating both contiguous and accurate diploid genomes. It first assembles co-barcoded reads to generate accurate haplotype-resolved assemblies that may contain many gaps, while the long-read assembly is contiguous but susceptible to errors. Then two assembly sets are integrated into haplotype-resolved assemblies with reduced misassembles. Through extensive evaluation on multiple synthetic datasets, AsmMix consistently demonstrates high precision and recall rates for haplotyping across diverse sequencing platforms, coverage depths, read lengths, and read accuracies, significantly outperforming other existing tools in the field. Furthermore, we validate the effectiveness of our pipeline using a human whole genome dataset (HG002), and produce highly contiguous, accurate, and haplotype-resolved assemblies. These assemblies are evaluated using the GIAB benchmarks, confirming the accuracy of variant calling. Our results demonstrate that AsmMix offers a straightforward yet highly efficient approach that effectively leverages both long reads and co-barcoded reads for haplotype-resolved assembly.

Whole mitochondrial genome sequencing provides new insights into the phylogeography of loggerhead turtles (Caretta caretta) in the Mediterranean Sea

Population structure and phylogeography of the loggerhead sea turtle (Caretta caretta) have so far been assessed mainly by mitochondrial DNA (mtDNA) single-gene sequencing studies. However, phylogenetic relationships amongst matrilines, genetic characterisation of rookeries and mixed-stock analyses have suffered from the limited resolution obtained by comparison of relatively short sequences such as from the mtDNA control region. Whole mitogenome sequencing can significantly improve population genetics, particularly in marine organisms showing female natal philopatry. Despite mitogenomics becoming increasingly common in biodiversity monitoring and conservation, only a few complete mitogenomes are available for C. caretta. In this study, we sequenced the complete mtDNA of 61 loggerhead turtles sampled between 2008 and 2021 along the Italian coastline and central Mediterranean Sea. We assigned complete mtDNA haplotypes to dead embryos and bycatch samples, and introduced a first nomenclature for loggerhead mitogenomes. Analysis of mtDNA diversity, Maximum Parsimony and Bayesian phylogenetic reconstruction allowed improved resolution of lineages with respect to studies reporting on partial mtDNA control region sequence comparisons, and we were able to further inform previous analyses on loggerhead ancestry based on control region haplogroups. Overall, whole mitogenome analysis has potential for considerable improvement of evolutionary history and phylogeographic investigations as well as mixed-stock surveys of loggerhead turtles.

Single-cell individual full-length mtDNA sequencing by iMiGseq uncovers unexpected heteroplasmy shifts in mtDNA editing.

The ontogeny and dynamics of mtDNA heteroplasmy remain unclear due to limitations of current mtDNA sequencing methods. We developed individual Mitochondrial Genome sequencing (iMiGseq) of full-length mtDNA for ultra-sensitive variant detection, complete haplotyping, and unbiased evaluation of heteroplasmy levels, all at the individual mtDNA molecule level. iMiGseq uncovered unappreciated levels of heteroplasmic variants in single cells well below the conventional NGS detection limit and provided accurate quantitation of heteroplasmy level. iMiGseq resolved the complete haplotype of individual mtDNA in single oocytes and revealed genetic linkage of de novo mutations. iMiGseq detected sequential acquisition of detrimental mutations, including large deletions, in defective mtDNA in NARP/Leigh syndrome patient-derived induced pluripotent stem cells. iMiGseq identified unintended heteroplasmy shifts in mitoTALEN editing, while showing no appreciable level of unintended mutations in DdCBE-mediated mtDNA base editing. Therefore, iMiGseq could not only help elucidate the mitochondrial etiology of diseases, but also evaluate the safety of various mtDNA editing strategies.

Two haplotype-resolved, gap-free genome assemblies for Actinidia latifolia and Actinidia chinensis shed light on the regulatory mechanisms of vitamin C and sucrose metabolism in kiwifruit

Kiwifruit is a recently domesticated horticultural fruit crop with substantial economic and nutritional value, especially because of the high content of vitamin C in its fruit. In this study, we de novo assembled two telomere-to-telomere kiwifruit genomes from Actinidia chinensis var. ‘Donghong’ (DH) and Actinidia latifolia ‘Kuoye’ (KY), with total lengths of 608 327 852 and 640 561 626 bp for 29 chromosomes, respectively. With a burst of structural variants involving inversion, translocations, and duplications within 8.39 million years, the metabolite content of DH and KY exhibited differences in saccharides, lignans, and vitamins. A regulatory ERF098 transcription factor family has expanded in KY and Actinidia eriantha, both of which have ultra-high vitamin C content. With each assembly phased into two complete haplotypes, we identified allelic variations between two sets of haplotypes, leading to protein sequence variations in 26 494 and 27 773 gene loci and allele-specific expression of 4687 and 12 238 homozygous gene pairs. Synchronized metabolome and transcriptome changes during DH fruit development revealed the same dynamic patterns in expression levels and metabolite contents; free fatty acids and flavonols accumulated in the early stages, but sugar substances and amino acids accumulated in the late stages. The AcSWEET9b gene that exhibits allelic dominance was further identified to positively correlate with high sucrose content in fruit. Compared with wild varieties and other Actinidia species, AcSWEET9b promoters were selected in red-flesh kiwifruits that have increased fruit sucrose content, providing a possible explanation on why red-flesh kiwifruits are sweeter. Collectively, these two gap-free kiwifruit genomes provide a valuable genetic resource for investigating domestication mechanisms and genome-based breeding of kiwifruit.

Insights into cryptic speciation of quillworts in China

Cryptic species are commonly misidentified because of high morphological similarities to other species. One group of plants that may harbor large numbers of cryptic species is the quillworts (Isoëtes spp.), an ancient aquatic plant lineage. Although over 350 species of Isoëtes have been reported globally, only ten species have been recorded in China. The aim of this study is to better understand Isoëtes species diversity in China. For this purpose, we systematically explored the phylogeny and evolution of Isoëtes using complete chloroplast genome (plastome) data, spore morphology, chromosome number, genetic structure, and haplotypes of almost all Chinese Isoëtes populations. We identified three ploidy levels of Isoëtes in China—diploid (2n = 22), tetraploid (2n = 44), and hexaploid (2n = 66). We also found four megaspore and microspore ornamentation types in diploids, six in tetraploids, and three in hexaploids. Phylogenetic analyses confirmed that I. hypsophila as the ancestral group of the genus and revealed that Isoëtes diploids, tetraploids, and hexaploids do not form monophyletic clades. Most individual species possess a single genetic structure; however, several samples have conflicting positions on the phylogenetic tree based on SNPs and the tree based on plastome data. All 36 samples shared 22 haplotypes. Divergence time analysis showed that I. hypsophila diverged in the early Eocene (∼48.05 Ma), and most other Isoëtes species diverged 3–20 Ma. Additionally, different species of Isoëtes were found to inhabit different water systems and environments along the Yangtze River. These findings provide new insights into the relationships among Isoëtes species in China, where highly similar morphologic populations may harbor many cryptic species.

Large-Scale Complete Sequencing and Haplotyping of 1-10kb DNA Molecules Using Short Massively Parallel Reads.

In this chapter, we describe a simple, low-cost method for making many copies of a single DNA molecule (1-10kb in length) as a concatemer on a long DNA strand. This can enable applications requiring high-quality contiguous sequence and haplotype data from long single DNA molecules at large scale.

HaploDMF: viral haplotype reconstruction from long reads via deep matrix factorization.

Lacking strict proofreading mechanisms, many RNA viruses can generate progeny with slightly changed genomes. Being able to characterize highly similar genomes (i.e. haplotypes) in one virus population helps study the viruses' evolution and their interactions with the host/other microbes. High-throughput sequencing data has become the major source for characterizing viral populations. However, the inherent limitation on read length by next-generation sequencing makes complete haplotype reconstruction difficult. In this work, we present a new tool named HaploDMF that can construct complete haplotypes using third-generation sequencing (TGS) data. HaploDMF utilizes a deep matrix factorization model with an adapted loss function to learn latent features from aligned reads automatically. The latent features are then used to cluster reads of the same haplotype. Unlike existing tools whose performance can be affected by the overlap size between reads, HaploDMF is able to achieve highly robust performance on data with different coverage, haplotype number and error rates. In particular, it can generate more complete haplotypes even when the sequencing coverage drops in the middle. We benchmark HaploDMF against the state-of-the-art tools on simulated and real sequencing TGS data on different viruses. The results show that HaploDMF competes favorably against all others. The source code and the documentation of HaploDMF are available at https://github.com/dhcai21/HaploDMF. Supplementary data are available at Bioinformatics online.

HAT: haplotype assembly tool using short and error-prone long reads.

Haplotypes are the set of alleles co-occurring on a single chromosome and inherited together to the next generation. Because a monoploid reference genome loses this co-occurrence information, it has limited use in associating phenotypes with allelic combinations of genotypes. Therefore, methods to reconstruct the complete haplotypes from DNA sequencing data are crucial. Recently, several attempts have been made at haplotype reconstructions, but significant limitations remain. High-quality continuous haplotypes cannot be created reliably, particularly when there are few differences between the homologous chromosomes. Here, we introduce HAT, a haplotype assembly tool that exploits short and long reads along with a reference genome to reconstruct haplotypes. HAT tries to take advantage of the accuracy of short reads and the length of the long reads to reconstruct haplotypes. We tested HAT on the aneuploid yeast strain Saccharomyces pastorianus CBS1483 and multiple simulated polyploid datasets of the same strain, showing that it outperforms existing tools. https://github.com/AbeelLab/hat/. Supplementary data are available at Bioinformatics online.

Diversity of NKG2C genotypes in a European population: Conserved and recombinant haplotypes in the coding, promoter, and 3'-untranslated regions.

NK cells monitor altered molecular patterns in tumors and infected cells through an ample array of receptors. Two families of evolutionarily distant receptors have converged to enable human NK cells to sense levels of HLA class I ligands, frequently abnormal in altered cells. Whilst different forms of polymorphism are a hallmark of killer-cell immunoglobulin-like receptors and their classic HLA-A, B, and C ligands, genetic diversity of killer-cell lectin-like receptors for the non-classical HLA-E (CD94/NKG2 heterodimers) is less conspicuous and has attracted less attention. A common pattern of diversification in both receptor families is evolution of pairs of inhibitory and activating homologs for a common ligand, the genes encoding activating receptors being more frequently affected by copy number variation (CNV). This is exemplified by the gene encoding the activating NKG2C subunit (KLRC2 or NKG2C), which marks an NK-cell subpopulation that differentiates or expands in response to cytomegalovirus. We have studied NKG2C diversity in 240 South European individuals, using polymerase chain reaction and sequencing methods to assess both gene CNV and single-nucleotide polymorphisms (SNPs) affecting its promoter, coding and 3'-untranslated (3'UT) regions. Sequence analysis revealed eight common SNPs-one in the promoter, two in the coding sequence, and five in the 3'UT region. These SNPs associate strongly with each other, forming three conserved extended haplotypes (frequencies: 0.456, 0.221, and 0.117). Homo- and heterozygous combination of these, together with complete gene deletion (0.175) and additional haplotypes with frequencies lower than 0.015, generate a diversity of NKG2C genotypes of potential immunological importance.

Experimental method for haplotype phasing across the entire length of chromosome 21 in trisomy 21 cells using a chromosome elimination technique

Modern sequencing technologies produce a single consensus sequence without distinguishing between homologous chromosomes. Haplotype phasing solves this limitation by identifying alleles on the maternal and paternal chromosomes. This information is critical for understanding gene expression models in genetic disease research. Furthermore, the haplotype phasing of three homologous chromosomes in trisomy cells is more complicated than that in disomy cells. In this study, we attempted the accurate and complete haplotype phasing of chromosome 21 in trisomy 21 cells. To separate homologs, we established three corrected disomy cell lines (ΔPaternal chromosome, ΔMaternal chromosome 1, and ΔMaternal chromosome 2) from trisomy 21 induced pluripotent stem cells by eliminating one chromosome 21 utilizing the Cre-loxP system. These cells were then whole-genome sequenced by a next-generation sequencer. By simply comparing the base information of the whole-genome sequence data at the same position between each corrected disomy cell line, we determined the base on the eliminated chromosome and performed phasing. We phased 51,596 single nucleotide polymorphisms (SNPs) on chromosome 21, randomly selected seven SNPs spanning the entire length of the chromosome, and confirmed that there was no contradiction by direct sequencing.

Chromosome-Level Haplotype Assembly for Equus asinu.

Background: Haplotype provides significant insights into understanding genomes at both individual and population levels. However, research on many non-model organisms is still based on independent genetic variations due to the lack of haplotype. Results: We conducted haplotype assembling for Equus asinu, a non-model organism that plays a vital role in human civilization. We described the hybrid single individual assembled haplotype of the Dezhou donkey based on the high-depth sequencing data from single-molecule real-time sequencing (×30), Illumina short-read sequencing (×211), and high-throughput chromosome conformation capture (×56). We assembled a near-complete haplotype for the high-depth sequenced Dezhou donkey individual and a phased cohort for the resequencing data of the donkey population. Conclusion: Here, we described the complete chromosome-scale haplotype of the Dezhou donkey with more than a 99.7% phase rate. We further phased a cohort of 156 donkeys to form a donkey haplotype dataset with more than 39 million genetic variations.

The haplotype-resolved chromosome pairs of a heterozygous diploid African cassava cultivar reveal novel pan-genome and allele-specific transcriptome features.

BackgroundCassava (Manihot esculenta) is an important clonally propagated food crop in tropical and subtropical regions worldwide. Genetic gain by molecular breeding has been limited, partially because cassava is a highly heterozygous crop with a repetitive and difficult-to-assemble genome.FindingsHere we demonstrate that Pacific Biosciences high-fidelity (HiFi) sequencing reads, in combination with the assembler hifiasm, produced genome assemblies at near complete haplotype resolution with higher continuity and accuracy compared to conventional long sequencing reads. We present 2 chromosome-scale haploid genomes phased with Hi-C technology for the diploid African cassava variety TME204. With consensus accuracy >QV46, contig N50 >18 Mb, BUSCO completeness of 99%, and 35k phased gene loci, it is the most accurate, continuous, complete, and haplotype-resolved cassava genome assembly so far. Ab initio gene prediction with RNA-seq data and Iso-Seq transcripts identified abundant novel gene loci, with enriched functionality related to chromatin organization, meristem development, and cell responses. During tissue development, differentially expressed transcripts of different haplotype origins were enriched for different functionality. In each tissue, 20–30% of transcripts showed allele-specific expression (ASE) differences. ASE bias was often tissue specific and inconsistent across different tissues. Direction-shifting was observed in <2% of the ASE transcripts. Despite high gene synteny, the HiFi genome assembly revealed extensive chromosome rearrangements and abundant intra-genomic and inter-genomic divergent sequences, with large structural variations mostly related to LTR retrotransposons. We use the reference-quality assemblies to build a cassava pan-genome and demonstrate its importance in representing the genetic diversity of cassava for downstream reference-guided omics analysis and breeding.ConclusionsThe phased and annotated chromosome pairs allow a systematic view of the heterozygous diploid genome organization in cassava with improved accuracy, completeness, and haplotype resolution. They will be a valuable resource for cassava breeding and research. Our study may also provide insights into developing cost-effective and efficient strategies for resolving complex genomes with high resolution, accuracy, and continuity.

Reconstructing viral haplotypes using long reads.

Most RNA viruses lack strict proofreading during replication. Coupled with a high replication rate, some RNA viruses can form a virus population containing a group of genetically related but different haplotypes. Characterizing the haplotype composition in a virus population is thus important to understand viruses' evolution. Many attempts have been made to reconstruct viral haplotypes using next-generation sequencing (NGS) reads. However, the short length of NGS reads cannot cover distant single-nucleotide variants, making it difficult to reconstruct complete or near-complete haplotypes. Given the fast developments of third-generation sequencing technologies, a new opportunity has arisen for reconstructing full-length haplotypes with long reads. In this work, we developed a new tool, RVHaplo to reconstruct haplotypes for known viruses from long reads. We tested it rigorously on both simulated and real viral sequencing data and compared it against other popular haplotype reconstruction tools. The results demonstrated that RVHaplo outperforms the state-of-the-art tools for viral haplotype reconstruction from long reads. Especially, RVHaplo can reconstruct the rare (1% abundance) haplotypes that other tools usually missed. The source code and the documentation of RVHaplo are available at https://github.com/dhcai21/RVHaplo. Supplementary data are available at Bioinformatics online.

Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR–HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR–HLA interactions among all described worldwide populations, and that 83–97% of their KIR–HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR–HLA coevolution and its impact on human health and survival.