Haplotype Resolution Research Articles

Addressing Technical Pitfalls in Pursuit of Molecular Factors That Mediate Immunoglobulin Gene Regulation.

The expressed Ab repertoire is a critical determinant of immune-related phenotypes. Ab-encoding transcripts are distinct from other expressed genes because they are transcribed from somatically rearranged gene segments. Human Abs are composed of two identical H and L chain polypeptides derived from genes in IGH locus and one of two L chain loci. The combinatorial diversity that results from Ab gene rearrangement and the pairing of different H and L chains contributes to the immense diversity of the baseline Ab repertoire. During rearrangement, Ab gene selection is mediated by factors that influence chromatin architecture, promoter/enhancer activity, and V(D)J recombination. Interindividual variation in the composition of the Ab repertoire associates with germline variation in IGH, implicating polymorphism in Ab gene regulation. Determining how IGH variants directly mediate gene regulation will require integration of these variants with other functional genomic datasets. In this study, we argue that standard approaches using short reads have limited utility for characterizing regulatory regions in IGH at haplotype resolution. Using simulated and chromatin immunoprecipitation sequencing reads, we define features of IGH that limit use of short reads and a single reference genome, namely 1) the highly duplicated nature of the DNA sequence in IGH and 2) structural polymorphisms that are frequent in the population. We demonstrate that personalized diploid references enhance performance of short-read data for characterizing mappable portions of the locus, while also showing that long-read profiling tools will ultimately be needed to fully resolve functional impacts of IGH germline variation on expressed Ab repertoires.

Improved haplotype resolution of highly duplicated MHC genes in a long-read genome assembly using MiSeq amplicons.

Long-read sequencing offers a great improvement in the assembly of complex genomic regions, such as the major histocompatibility complex (MHC) region, which can contain both tandemly duplicated MHC genes (paralogs) and high repeat content. The MHC genes have expanded in passerine birds, resulting in numerous MHC paralogs, with relatively high sequence similarity, making the assembly of the MHC region challenging even with long-read sequencing. In addition, MHC genes show rather high sequence divergence between alleles, making diploid-aware assemblers incorrectly classify haplotypes from the same locus as sequences originating from different genomic regions. Consequently, the number of MHC paralogs can easily be over- or underestimated in long-read assemblies. We therefore set out to verify the MHC diversity in an original and a haplotype-purged long-read assembly of one great reed warbler Acrocephalus arundinaceus individual (the focal individual) by using Illumina MiSeq amplicon sequencing. Single exons, representing MHC class I (MHC-I) and class IIB (MHC-IIB) alleles, were sequenced in the focal individual and mapped to the annotated MHC alleles in the original long-read genome assembly. Eighty-four percent of the annotated MHC-I alleles in the original long-read genome assembly were detected using 55% of the amplicon alleles and likewise, 78% of the annotated MHC-IIB alleles were detected using 61% of the amplicon alleles, indicating an incomplete annotation of MHC genes. In the haploid genome assembly, each MHC-IIB gene should be represented by one allele. The parental origin of the MHC-IIB amplicon alleles in the focal individual was determined by sequencing MHC-IIB in its parents. Two of five larger scaffolds, containing 6-19 MHC-IIB paralogs, had a maternal and paternal origin, respectively, as well as a high nucleotide similarity, which suggests that these scaffolds had been incorrectly assigned as belonging to different loci in the genome rather than as alternate haplotypes of the same locus. Therefore, the number of MHC-IIB paralogs was overestimated in the haploid genome assembly. Based on our findings we propose amplicon sequencing as a suitable complement to long-read sequencing for independent validation of the number of paralogs in general and for haplotype inference in multigene families in particular.

DelInsCaller: An Efficient Algorithm for Identifying Delins and Estimating Haplotypes from Long Reads with High Level of Sequencing Errors.

Delins, as known as complex indel, is a combined genomic structural variation formed by deleting and inserting DNA fragments at a common genomic location. Recent studies emphasized the importance of delins in cancer diagnosis and treatment. Although the long reads from PacBio CLR sequencing significantly facilitate delins calling, the existing approaches still encounter computational challenges from the high level of sequencing errors, and often introduce errors in genotyping and phasing delins. In this paper, we propose an efficient algorithmic pipeline, named delInsCaller, to identify delins on haplotype resolution from the PacBio CLR sequencing data. delInsCaller design a fault-tolerant method by calculating a variation density score, which helps to locate the candidate mutational regions under a high-level of sequencing errors. It adopts a base association-based contig splicing method, which facilitates contig splicing in the presence of false-positive interference. We conducted a series of experiments on simulated datasets, and the results showed that delInsCaller outperformed several state-of-the-art approaches, e.g., SVseq3, across a wide range of parameter settings, such as read depth, sequencing error rates, etc. delInsCaller often obtained higher f-measures than other approaches; specifically, it was able to maintain advantages at ~15% sequencing errors. delInsCaller was able to significantly improve the N50 values with almost no loss of haplotype accuracy compared with the existing approach as well.

Forensic assessment on the application of a virtual pool of 30 Y-STRs

Ideally, a panel of Y short tandem repeats (Y-STRs) should include markers with gene diversity (GD) ≥ 0.70 to guarantee a highly accurate and precise individual identification in all forensic situations where Y loci analysis is necessary to support the findings obtained with autosomal STRs. Since not all markers in a multiplex set reach the optimal GD value, to achieve a greater discrimination capacity (DC), the number of Y-STRs has been gradually expanded also with the addition of rapidly mutating Y-STRs (RM Y-STRs). The present study has tried to assess, based on the genetic diversity of each Y-STR and the cumulative haplotype diversity provided by the single Y-STR panels tested, whether to increase in the number of Y-STRs is effective. For this purpose, three commercially available forensic kits - the PowerPlex® Y23 System, the Yfiler™ Plus PCR Amplification Kit, and the ForenSeq™ DNA Signature Prep Kit - that allow for the simultaneous analysis of 22–25 Y chromosome markers (including Y-STRs and RM Y-STRs) have been used for typing 115 unrelated male individuals from North-East Italy employing two capillary electrophoresis (Applied Biosystems® 3130 Genetic Analyzer and SeqStudio Genetic Analyzer for HID) and one massively parallel sequencing (MiSeq FGx™ Forensic Genomics System) systems. Forensic parameters, such as gene diversity, cumulative haplotype diversity (cHD), and discrimination capacity, were determined for the three Y-STRs kits and the virtual pool of 30 Y-STRs obtained by the overall combination of the Y-loci of the three kits. From the findings, it emerges that twelve markers in the virtual panel of 30 Y-STRs, are characterized by a value of GD ≥ 0.70, and of these, eleven included in the Yfiler™ Plus kit seem to be sufficient to ensure the achievement of haplotype resolution. In light of this, it does not appear necessary to expand the number of Y-STRs in the currently available kits.

A novel set of polymorphic chloroplast microsatellite markers for northern red oak (Q. rubra L.)

AbstractOaks are model species due to their importance in various ecosystems, worldwide distribution, high economic value and emerging genomic resources. As such, knowledge on population differentiation across their distribution range is of high importance for sustainable forest management. As in most angiosperms, chloroplasts are maternally (via seeds) inherited in oak species and acorns are dispersed over comparably short distances. Consequently, chloroplast markers reveal comparatively high differentiation between populations, making them highly viable for the analysis of historic migration patterns and the certification of reproductive material. Despite the existence of various white oak (section: Quercus) chloroplast markers, red oak (section: Lobatae) chloroplast markers remain limited. Northern red oak is one of the most important North American oak species and a widespread non-native plantation tree species in European forests. We took advantage of chloroplast genomes of Q. rubra L. and related oak chloroplast genomes to develop a set of robust and easy-to-use chloroplast microsatellite primers for northern red oak. Furthermore, we tested transferability of those novel red oak primer pairs to Q. robur L. and Q. petraea Matt. Liebl. The new set of fifteen polymorphic chloroplast microsatellite markers revealed three additional northern red oak haplotypes after screening 80 northern red oak individuals representing seven haplotypes, identified based on formerly available markers. Therefore, they provide a higher resolution of haplotypes as compared to currently available markers.

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.

Improving the regional Y-STR haplotype resolution utilizing haplogroup-determining Y-SNPs and the application of machine learning in Y-SNP haplogroup prediction in a forensic Y-STR database: A pilot study on male Chinese Yunnan Zhaoyang Han population

Improving the resolution of the current widely used Y-chromosomal short tandem repeat (Y-STR) dataset is of great importance for forensic investigators, and the current approach is limited, except for the addition of more Y-STR loci. In this research, a regional Y-DNA database was investigated to improve the Y-STR haplotype resolution utilizing a Y-SNP Pedigree Tagging System that includes 24 Y-chromosomal single nucleotide polymorphism (Y-SNP) loci. This pilot study was conducted in the Chinese Yunnan Zhaoyang Han population, and 3473 unrelated male individuals were enrolled. Based on data on the male haplogroups under different panels, the matched or near-matching (NM) Y-STR haplotype pairs from different haplogroups indicated the critical roles of haplogroups in improving the regional Y-STR haplotype resolution. A classic median-joining network analysis was performed using Y-STR or Y-STR/Y-SNP data to reconstruct population substructures, which revealed the ability of Y-SNPs to correct misclassifications from Y-STRs. Additionally, population substructures were reconstructed using multiple unsupervised or supervised dimensionality reduction methods, which indicated the potential of Y-STR haplotypes in predicting Y-SNP haplogroups. Haplogroup prediction models were built based on nine publicly accessible machine-learning (ML) approaches. The results showed that the best prediction accuracy score could reach 99.71% for major haplogroups and 98.54% for detailed haplogroups. Potential influences on prediction accuracy were assessed by adjusting the Y-STR locus numbers, selecting Y-STR loci with various mutabilities, and performing data processing. ML-based predictors generally presented a better prediction accuracy than two available predictors (Nevgen and EA-YPredictor). Three tree models were developed based on the Yfiler Plus panel with unprocessed input data, which showed their strong generalization ability in classifying various Chinese Han subgroups (validation dataset). In conclusion, this study revealed the significance and application prospects of Y-SNP haplogroups in improving regional Y-STR databases. Y-SNP haplogroups can be used to discriminate NM Y-STR haplotype pairs, and it is important for forensic Y-STR databases to develop haplogroup prediction tools to improve the accuracy of biogeographic ancestry inferences.

A multi-dimensional evaluation of the ‘NIST 1032’ sample set across four forensic Y-STR multiplexes

This manuscript reports Y-chromosomal short tandem repeat (Y-STR) haplotypes for 1032 male U.S. population samples across 30 Y-STR loci characterized by three capillary electrophoresis (CE) length-based kits (PowerPlex Y23 System, Yfiler Plus PCR Amplification Kit, and Investigator Argus Y-28 QS Kit) and one sequence-based kit (ForenSeq DNA Signature Prep Kit): DYF387S1, DYS19, DYS385 a/b, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS449, DYS456, DYS458, DYS460, DYS481, DYS505, DYS518, DYS522, DYS533, DYS549, DYS570, DYS576, DYS612, DYS627, DYS635, DYS643, and Y-GATA-H4. The length-based Y-STR haplotypes include six loci that are not reported in the sequence-based kit (DYS393, DYS449, DYS456, DYS458, DYS518, and DYS627), whereas three loci included in the sequence-based kit are not present in length-based kits (DYS505, DYS522, and DYS612). For the latter, a custom multiplex was used to generate CE length-based data, allowing 1032 samples to be evaluated for concordance across the 30 Y-STR loci included in these four commercial Y-STR typing kits. Discordances between typing methods were analyzed further to assess underlying causes such as primer binding site mutations and flanking region insertions/deletions. Allele-level frequency and statistical information is provided for sequenced loci, excluding the multi-copy loci DYF387S1 and DYS385 a/b, for which locus-specific haplotype-level frequencies are provided instead. The resulting data reveals the degree of information gained through sequencing: 88% of sequenced Y-STR loci contain additional sequence-based alleles compared to length-based data, with the DYS389II locus containing the most additional alleles (51) observed by sequencing. Despite these allelic increases, only minimal improvement was observed in haplotype resolution by sequence, with all four commercial kits providing a similar ability to differentiate length-based haplotypes in this sample set. Finally, a subset of 369 male samples were compared to their corresponding additionally sequenced father samples, revealing the sequence basis for the 50 length-based changes observed, and no additional sequence-based mutations. GenBank accession numbers are reported for each unique sequence, and associated records are available in the STRSeq Y-Chromosomal STR Loci National Center for Biotechnology Information (NCBI) BioProject, accession PRJNA380347. Haplotype data is updated in the Y-STR Haplotype Reference Database (YHRD) for the ‘NIST 1032’ data set to now achieve the level of maximal haplotype of YHRD. All supplementary files including revisions to previously published Y-STR data are available in the NIST Public Data Repository: U.S. population data for human identification markers, DOI 10.18434/t4/1500024.

JWES: a new pipeline for whole genome/exome sequence data processing, management, and gene-variant discovery, annotation, prediction, and genotyping.

Whole genome and exome sequencing (WGS/WES) are the most popular next‐generation sequencing (NGS) methodologies and are at present often used to detect rare and common genetic variants of clinical significance. We emphasize that automated sequence data processing, management, and visualization should be an indispensable component of modern WGS and WES data analysis for sequence assembly, variant detection (SNPs, SVs), imputation, and resolution of haplotypes. In this manuscript, we present a newly developed findable, accessible, interoperable, and reusable (FAIR) bioinformatics‐genomics pipeline Java based Whole Genome/Exome Sequence Data Processing Pipeline (JWES) for efficient variant discovery and interpretation, and big data modeling and visualization. JWES is a cross‐platform, user‐friendly, product line application, that entails three modules: (a) data processing, (b) storage, and (c) visualization. The data processing module performs a series of different tasks for variant calling, the data storage module efficiently manages high‐volume gene‐variant data, and the data visualization module supports variant data interpretation with Circos graphs. The performance of JWES was tested and validated in‐house with different experiments, using Microsoft Windows, macOS Big Sur, and UNIX operating systems. JWES is an open‐source and freely available pipeline, allowing scientists to take full advantage of all the computing resources available, without requiring much computer science knowledge. We have successfully applied JWES for processing, management, and gene‐variant discovery, annotation, prediction, and genotyping of WGS and WES data to analyze variable complex disorders. In summary, we report the performance of JWES with some reproducible case studies, using open access and in‐house generated, high‐quality datasets.

Investigating Different DNA Methylation Patterns at the Resolution of Methylation Haplotypes.

Different DNA methylation patterns presented on different tissues or cell types are considered as one of the main reasons accounting for the tissue-specific gene expressions. In recent years, many methods have been proposed to identify differentially methylated regions (DMRs) based on the mixture of methylation signals from homologous chromosomes. To investigate the possible influence of homologous chromosomes on methylation analysis, this paper proposed a method (MHap) to construct methylation haplotypes for homologous chromosomes in CpG dense regions. Through comparing the methylation consistency between homologous chromosomes in different cell types, it can be found that majority of paired methylation haplotypes derived from homologous chromosomes are consistent, while a lower methylation consistency was observed in the breast cancer sample. It also can be observed that the hypomethylation consistency of differentiated cells is higher than that of the corresponding undifferentiated stem cells. Furthermore, based on the methylation haplotypes constructed on homologous chromosomes, a method (MHap_DMR) is developed to identify DMRs between differentiated cells and the corresponding undifferentiated stem cells, or between the breast cancer sample and the normal breast sample. Through comparing the methylation haplotype modes of DMRs in two cell types, the DNA methylation changing directions of homologous chromosomes in cell differentiation and cancerization can be revealed. The code is available at: https://github.com/xqpeng/MHap_DMR.

Stillbirth rates and their association with swine leucocyte antigen class II haplotypes in Microminipigs.

ObjectiveMicrominipig (MMP) is a miniature pig with an extra small body size for experimental use. In the present study, the occurrence of stillbirths and their genetic association with swine leukocyte antigen (SLA) class II haplotypes were evaluated in a population of MMPs.MethodsThe occurrences of stillbirth and genetic association with SLA class II haplotypes using 483 stillborn and 2,246 live piglets, and their parents were compared among the three groups of newborn piglet litters; an all stillborn (AS) group consisting of only stillborn piglet litters, a partial stillborn (PS) group consisting of stillborn and live piglet litters, and an all alive (AA) group consisting of only live piglet litters.ResultsThe incidence of stillborn piglets was 483/2,729 (17.7%). Distributions of litter sizes, numbers of stillborn piglets in a litter, parities, and gestation periods were distinct among the three groups. The frequencies of low resolution haplotype (Lr)-0.7 or Lr-0.23 were higher in the AS group than in the PS or AA groups. In sires, the frequency of Lr-0.7 associated with the AS group was significantly higher in the AS group than with the AA group. In dams, the frequency of Lr-0.23 was significantly higher in the AS group than in the PS or AA groups, whereas the frequency of Lr-0.7 was not significantly different.ConclusionThe incidence of stillborn piglets in MMPs appears to be higher than those in other pig breeds. Several traits related with stillbirths such as the number of stillborn piglets and parities of the AS group were different from those of the PS and AA groups. Specific SLA class II haplotypes were associated significantly with a high incidence of stillbirths and could be used as genetic markers to adopt breeding strategies to lower the rate of stillbirth in MMPs.

Genetic polymorphisms of 44 Y chromosomal genetic markers in the Inner Mongolia Han population and its genetic relationship analysis with other reference populations

Y chromosomal genetic markers in the non-recombining region are commonly used for human evolution research, familial searching, and forensic male differentiation since they strictly follow paternal inheritance. Y chromosomal short tandem repeats (Y-STRs) possess extraordinarily advantages in forensic applications because of their high polymorphisms and special genetic pattern. Here, we assessed the genetic diversities of 41 Y-STRs and three Y chromosomal insertion/deletion (Y-InDels) loci in the Chinese Inner Mongolia Han population; besides, genetic differentiation analyses among the studied Han population and other previously reported populations were conducted based on 27 same Y-STRs. Totally, 425 alleles were observed in 324 Inner Mongolia Han individuals for these Y-markers. Gene diversities of these Y-markers distributed from 0.0306 to 0.9634. The haplotype diversity and discriminatory capacity of these Y-markers in the Inner Mongolia Han population were 0.9999 and 0.98457, respectively. Haplotype resolution comparisons of different Y-marker groups in the studied Han population revealed that higher haplotype resolution could be achieved for these 44 Y-markers. Population genetic analyses of the Inner Mongolia Han population and other reference populations demonstrated that the studied Han population had relatively closer genetic affinities with Northern Han Chinese populations than Southern Han and other minority groups. To sum up, these 44 Y-markers can be utilized as a valuable tool for male differentiation in the Inner Mongolia Han population.

Haplotyping by linked-read sequencing (HLRS) of the genetic disease carriers for preimplantation genetic testing without a proband or relatives

BackgroundIn order to mitigate the risk of allele dropout (ADO) and ensure the accuracy of preimplantation genetic testing for monogenic disease (PGT-M), it is necessary to construct parental haplotypes. Typically, haplotype resolution is obtained by genotyping multiple polymorphic markers in both parents and a proband or a relative. Sometimes, single sperm typing, or tests on the polar bodies may also be useful. Nevertheless, this process is time-consuming. At present, there was no simple linkage analysis strategy for patients without affected relatives.MethodTo solve this problem, we established a haplotyping by linked-read sequencing (HLRS) method without the requirement for additional relatives. First, the haplotype of the genetic disease carriers in the family was constructed by linked-read sequencing, and then the informative single nucleotide polymorphisms (SNPs) in upstream and downstream mutation region were selected to construct the embryo haplotype and to determine whether the embryo was carrying the mutation. Two families were selected to validate this method; one with alpha thalassemia and the other with NDP gene disorder.ResultsThe haplotyping by linked-read sequencing (HLRS) method was successfully applied to construct parental haplotypes without recruiting additional family members; the method was also validated for PGT-M. The mutation carriers in these families were sequenced by linked-read sequencing, and their haplotypes were successfully phased. Adjacent SNPs of the mutation gene were identified. The informative SNPs were chosen for linkage analyses to identify the carrier embryos. For the alpha thalassemia family, a normal blastocyst was transferred to the uterus and the accuracy of PGT-M was confirmed by amniocentesis at 16 weeks of gestation.ConclusionsOur results suggest that HLRS can be applied for PGT-M of monogenic disorders or de novo mutations where the mutations haplotype cannot be determined due to absence of affected relatives.

A haplotype-resolved, de novo genome assembly for the wood tiger moth (Arctia plantaginis) through trio binning.

ABSTRACTBackgroundDiploid genome assembly is typically impeded by heterozygosity because it introduces errors when haplotypes are collapsed into a consensus sequence. Trio binning offers an innovative solution that exploits heterozygosity for assembly. Short, parental reads are used to assign parental origin to long reads from their F1 offspring before assembly, enabling complete haplotype resolution. Trio binning could therefore provide an effective strategy for assembling highly heterozygous genomes, which are traditionally problematic, such as insect genomes. This includes the wood tiger moth (Arctia plantaginis), which is an evolutionary study system for warning colour polymorphism.FindingsWe produced a high-quality, haplotype-resolved assembly for Arctia plantaginis through trio binning. We sequenced a same-species family (F1 heterozygosity ∼1.9%) and used parental Illumina reads to bin 99.98% of offspring Pacific Biosciences reads by parental origin, before assembling each haplotype separately and scaffolding with 10X linked reads. Both assemblies are contiguous (mean scaffold N50: 8.2 Mb) and complete (mean BUSCO completeness: 97.3%), with annotations and 31 chromosomes identified through karyotyping. We used the assembly to analyse genome-wide population structure and relationships between 40 wild resequenced individuals from 5 populations across Europe, revealing the Georgian population as the most genetically differentiated with the lowest genetic diversity.ConclusionsWe present the first invertebrate genome to be assembled via trio binning. This assembly is one of the highest quality genomes available for Lepidoptera, supporting trio binning as a potent strategy for assembling heterozygous genomes. Using our assembly, we provide genomic insights into the geographic population structure of A. plantaginis.

Development of novel Quercus rubra chloroplast genome CAPS markers for haplotype identification

Abstract Our main objective was to generate cost-effective chloroplast (cp) DNA markers that are easy to apply and to score. In combination with already published cpSSR markers they should increase haplotype resolution in populations. To discover new cpDNA markers, we sequenced 87-97 % of the entire chloroplast genome (except the second inverted repeat) of 8 trees representing different regions of the Quercus rubra L. natural range with 4,030X-6,297X coverage and assembled the genome sequences using the publicly available chloroplast genome of Quercus rubra L. as a reference. In total, 118 single nucleotide polymorphisms (SNPs) and 107 insertions or deletions (indels) were detected, and 15 cleaved amplified polymorphic sequence (CAPS) markers were developed for Q. rubra. Using these new markers together with five chloroplast microsatellite or simple sequence repeat (cpSSR) markers, we identified 10 haplotypes in our diversity panel of 19 Q. rubra populations. Specifically, two haplotypes based only on the cpSSR markers could now be separated in five haplotypes. These markers are useful to assess haplotype diversity with high resolution and are also transferable to a closely related species, Quercus ellipsoidalis E. J. Hill.

A study of genetic analysis using novel rapidly mutating Y-STR multiplex for Qatari population

Differentiating male lineages using non-recombining Y-chromosomal genetic markers is highly informative for tracing human migration and for forensic studies. The aim of this study was to develop allelic frequency to evaluate the resolution power of Rapidly Mutating Y-STR (RMY-STR) on Qatar population in comparison to Worldwide population. In Qatari population the haplotype resolution of 100 % was achieved. The set of 13 RMY-STR markers showed remarkable haplotype resolution power in the Qatar population, high gene diversity and sufficient robustness. The novel multiplex PCR reaction wasused to type 250 samples from unrelated Qatari males. A total of 250 different haplotypes were det, with no shared haplotypes between individuals. The overall haplotype diversity (HD) was 100 %. It was found that the markers which contributed the most toward high HD were DYF399S1 and DYF403S1a/b. FST values were generated and AMOVA test performed between Qatar population and Middle Eastern/several worldwide population data sets were undertaken for FST values; geography was found to account considerably for the pattern of population sub structuring, although some discrepancies were observed in the relationship between geographical and genetic distance.

Haplotype data for the 12 RM Y-STR loci in a Syrian population

Researches with RM Y-STRs have shown that these loci provide substantially higher haplotype diversity and haplotype discrimination capacity in worldwide populations when compared with the YSTRs commonly used in genetic forensics. The aim of this study was to develop an allelic frequency database for the Syrian population living in Turkey in order to obtain population data of 12 RM Y-STRs. A total of 80 unrelated males from the Syrian population living in Turkey were typed with 12 RM Y-STRs loci: DYF387S1, DYF399S1, DYF404S1, DYS449, DYS518, DYS526a/b, DYS547, DYS570, DYS576, DYS612, DYS626 and DYS627. The highest GD was observed for the locus DYF399S1 (0.91), followed by loci DYS449 (0.86) and DYS518 (0.83). RM Y-STR haplotype diversity was found 1.00 in these samples. Based on the results of this study, the RM YSTR loci showed remarkable haplotype resolution power in the Syrian population, high genetic diversity and, therefore, demonstrating their usefulness in forensic identification cases.

A new era of long-read sequencing for cancer genomics

Cancer is a disease largely caused by genomic aberrations. Utilizing many rapidly emerging sequencing technologies, researchers have studied cancer genomes to understand the molecular statuses of cancer cells and to reveal their vulnerabilities, such as driver mutations or gene expression. Long-read technologies enable us to identify and characterize novel types of cancerous mutations, including complicated structural variants in haplotype resolution. In this review, we introduce three representative platforms for long-read sequencing and research trends of cancer genomics with long-read data. Further, we describe that aberrant transcriptome and epigenome statuses, namely, fusion transcripts, as well as aberrant transcript isoforms and the phase information of DNA methylation, are able to be elucidated by long-read sequencers. Long-read sequencing may shed light on novel types of aberrations in cancer genomics that are being missed by conventional short-read sequencing analyses.

Off the street phasing (OTSP): no hassle haplotype phasing for molecular PGD applications.

Pre-implantation genetic diagnosis (PGD) for molecular disorders requires the construction of parental haplotypes. Classically, haplotype resolution ("phasing") is obtained by genotyping multiple polymorphic markers in both parents and at least one additional relative. However, this process is time-consuming, and immediate family members are not always available. The recent availability of massive genomic data for many populations promises to eliminate the needs for developing family-specific assays and for recruiting additional family members. In this study, we aimed to validate population-assisted haplotype phasing for PGD. Targeted sequencing of CFTR gene variants and ~ 1700 flanking polymorphic SNPs (± 2Mb) was performed on 54 individuals from 12 PGD families of (a) Full Ashkenazi (FA; n = 16), (b) mixed Ashkenazi (MA; n = 23 individuals with at least one Ashkenazi and one non-Ashkenazi grandparents), or (c) non-Ashkenazi (NA; n = 15) descent. Heterozygous genotype calls in each individual were phased using various whole genome reference panels and appropriate computational models. All computationally derived haplotype predictions were benchmarked against trio-based phasing. Using the Ashkenazi reference panel, phasing of FA was highly accurate (99.4% ± 0.2% accuracy); phasing of MA was less accurate (95.4% ± 4.5% accuracy); and phasing of NA was predictably low (83.4% ± 6.6% accuracy). Strikingly, for founder mutation carriers, our haplotyping approach facilitated near perfect phasing accuracy (99.9% ± 0.1% and 98.2% ± 2.8% accuracy for W1282X and delF508 carriers, respectively). Our results demonstrate the feasibility of replacing classical haplotype phasing with population-based phasing with uncompromised accuracy.

P073The common uncommon extended HLA haplotype: Do two weak associations make it strong?

Aim Genetic recombination occurs at specific hotspots within the MHC, leading to widely dispersed block-like structures; the alpha block containing HLA-A, beta block (HLA-C and HLA-B) and delta block (HLA-DR and HLA-DQ). Linkage disequilibrium (LD) is much stronger within than across these blocks. The HLA-A locus (most telomeric) has weaker LD with extended HLA haplotypes. Additionally, a recombination hotspot exists between HLA-DP genes and other class II loci, which leads to an HLA-DPB1 match rate of only 10–15% in 10/10 HLA-matched donor/recipient pairs. This report describes an extended 9 locus HLA haplotype found in 2 HCT and 1 kidney transplant recipients and their potential donors (8, 2 and 1 donors, respectively). Methods 1. Case review of 2 unrelated HCT recipients, both European Caucasian (CAU) and their related and unrelated potential donors. 2. Family study of HCT recipient and a sibling to identify the haplotype by segregation. 3. Determine frequency of observed haplotype in different populations through Haplostats.org (6 locus at 2 field resolution could yield higher estimate than the 9 locus 3 field resolution haplotype described here). 4. Searching lab database for incidence of this haplotype. Results All identified cases shared the following haplotype HLA-A*31:01:02∼C*07:02:01∼B*07:02:01∼DRB5*01:01:01∼DRB1*15:01:01∼DQA1*01:02:01∼DQB1*06:02:01∼DPA1*01:03:01∼DPB1*04:01:01. Although most alleles in this haplotype are fairly common, with HLA-A*31:01 least common (∼5%), Haplostats indicated that the haplotype HLA-A*31:01∼C*07:02∼B*07:02∼DRB5*01:01∼DRB1*15:01∼DQB1*06:02 ranks 61898th & 170477th in CAU & African Americans (AFA), respectively and is not seen in other populations. Searching our laboratory database for subjects A31 positive, 382 subjects out of 8340 typed between 2015 & 2017 were identified. Of those, 40 were positive for the haplotype HLA-A*31:01∼C*07:02∼B*07:02∼DRB5*01:01∼DRB1*15:01∼DQB1*06:02 and of these 17 were typed for HLA-DPB1. All were positive for the entire extended 9 locus haplotype and were either CAU (n = 6) or of unknown ethnicity (n = 11). Conclusions Although the 6 locus haplotype portion of this extended haplotype seems fairly uncommon, the extended 9 locus subset of that haplotype seems to be conserved especially in Caucasians.