Haplotype-based inference of the distribution of fitness effects

Recent genome sequencing studies with large sample sizes in humans have discovered a vast quantity of low-frequency variants, providing an important source of information to analyze how selection is acting on human genetic variation. In order to estimate the strength of natural selection acting on low-frequency variants, we have developed a likelihood-based method that uses the lengths of pairwise identity-by-state between haplotypes carrying low-frequency variants. We show that in some non-equilibrium populations (such as those that have had recent population expansions) it is possible to distinguish between positive or negative selection acting on a set of variants. With our new framework, one can infer a fixed selection intensity acting on a set of variants at a particular frequency, or a distribution of selection coefficients for standing variants and new mutations. We show an application of our method to the UK10K phased haplotype dataset of individuals.

Population Genetics of the Invasive Red Fox, Vulpes vulpes, in South-Eastern Australia

The use of genetic information in conservation biology has become more widespread with genetic information more readily available for non-model organisms. It has also been recognized that genetic information from invasive species can inform their management and control. The red fox poses a significant threat to Australian native fauna and the agricultural industry. Despite this, there are few recently published studies investigating the population genetics of foxes in Australia. This study investigated the population genetics of 94 foxes across the Illawarra and Shoalhaven regions of New South Wales, Australia. Diversity Array sequencing technology was used to genotype a large number of single nucleotide polymorphisms (N = 33,375). Moderate genetic diversity and relatedness were observed across the foxes sampled. Low to moderate levels of inbreeding, high-levels of identity-by-state values, as well as high identity-by-descent values were also found. There was limited evidence for population genetic structure among the foxes across the landscape sampled, supporting the presence of a single population across the study area. This indicates that there may be no barriers hindering fox dispersal across the landscape.

A likelihood ratio approach for identifying three-quarter siblings in genetic databases

The detection of family relationships in genetic databases is of interest in various scientific disciplines such as genetic epidemiology, population and conservation genetics, forensic science, and genealogical research. Nowadays, screening genetic databases for related individuals forms an important aspect of standard quality control procedures. Relatedness research is usually based on an allele sharing analysis of identity by state (IBS) or identity by descent (IBD) alleles. Existing IBS/IBD methods mainly aim to identify first-degree relationships (parent–offspring or full siblings) and second degree (half-siblings, avuncular, or grandparent–grandchild) pairs. Little attention has been paid to the detection of in-between first and second-degree relationships such as three-quarter siblings (3/4S) who share fewer alleles than first-degree relationships but more alleles than second-degree relationships. With the progressively increasing sample sizes used in genetic research, it becomes more likely that such relationships are present in the database under study. In this paper, we extend existing likelihood ratio (LR) methodology to accurately infer the existence of 3/4S, distinguishing them from full siblings and second-degree relatives. We use bootstrap confidence intervals to express uncertainty in the LRs. Our proposal accounts for linkage disequilibrium (LD) by using marker pruning, and we validate our methodology with a pedigree-based simulation study accounting for both LD and recombination. An empirical genome-wide array data set from the GCAT Genomes for Life cohort project is used to illustrate the method.

Plasmodium falciparum genomic surveillance reveals spatial and temporal trends, association of genetic and physical distance, and household clustering

ABSTRACTMolecular epidemiology using genomic data can help identify relationships between malaria parasite population structure, malaria transmission intensity, and ultimately help generate actionable data to assess the effectiveness of malaria control strategies. Genomic data, coupled with geographic information systems data, can further identify clusters or hotspots of malaria transmission, parasite genetic and spatial connectivity, and parasite movement by human or mosquito mobility over time and space. In this study, we performed longitudinal genomic surveillance in a cohort of 70 participants over four years from different neighborhoods and households in Thiès, Senegal—a region of exceptionally low malaria transmission (entomological inoculation rate (EIR) less than 1). Genetic identity (identity by state) was established using a 24 single nucleotide polymorphism molecular barcode and a multivariable linear regression model was used to establish genetic and spatial relationships. Our results show clustering of genetically similar parasites within households and a decline in genetic similarity of parasites with increasing distance. One household showed extremely high diversity and warrants further investigation as to the source of these diverse genetic types. This study illustrates the utility of genomic data with traditional epidemiological approaches for surveillance and detection of trends and patterns in malaria transmission not only by neighborhood but also by household. This approach can be implemented regionally and countrywide to strengthen and support malaria control and elimination efforts.

Probabilistic estimation of identity by descent segment endpoints and detection of recent selection

Most methods for fast detection of identity by descent (IBD) segments report identity by state segments without any quantification of the uncertainty in the endpoints and lengths of the IBD segments. We present a method for determining the posterior probability distribution of IBD segment endpoints. Our approach accounts for genotype errors, recent mutations, and gene conversions which disrupt DNA sequence identity within IBD segments. We find that our method’s estimates of uncertainty are well calibrated for homogeneous samples. We quantify endpoint uncertainty for 7.7 billion IBD segments from 408,883 individuals of White British ancestry in the UK Biobank, and use these IBD segments to find regions showing evidence of recent natural selection. We show that many spurious selection signals are eliminated by the use of unbiased estimates of IBD segment endpoints and a pedigree-based genetic map. Nine of the top ten regions with the greatest evidence for recent selection in our scan have been identified as selected in previous analyses using different approaches. Our computationally efficient method for quantifying IBD segment endpoint uncertainty is implemented in the open source ibd-ends software package.

Inbreeding and Population Structure

Inbreeding and its consequences are the main subject of Chapter 3, beginning with the concepts of identity by descent versus identity by state, the inbreeding coefficient F, genotype frequencies with inbreeding, and calculation of the inbreeding coefficient from pedigrees. Inbreeding and heterosis are discussed along with the effects of inbreeding in humans and other organisms, regular systems of mating (selfing and partial selfing, sib mating), and the utility of recombinant inbred lines. The chapter emphasizes the intimate connection between inbreeding and hierarchical population structure as measured by the F-statistics.

Identity-by-state-based haplotyping expands the application of comprehensive preimplantation genetic testing

STUDY QUESTION Is it possible to haplotype parents using parental siblings to leverage preimplantation genetic testing (PGT) for monogenic diseases and aneuploidy (comprehensive PGT) by genome-wide haplotyping? SUMMARY ANSWER We imputed identity-by-state (IBS) sharing of parental siblings to phase parental genotypes. WHAT IS KNOWN ALREADY Genome-wide haplotyping of preimplantation embryos is being implemented as a generic approach for genetic diagnosis of inherited single-gene disorders. To enable the phasing of genotypes into haplotypes, genotyping the direct family members of the prospective parent carrying the mutation is required. Current approaches require genotypes of either (i) both or one of the parents of the affected prospective parent or (ii) an affected or an unaffected child of the couple. However, this approach cannot be used when parents or children are not attainable, prompting an investigation into alternative phasing options. STUDY DESIGN, SIZE, DURATION This is a retrospective validation study, which applied IBS-based phasing of parental haplotypes in 56 embryos derived from 12 PGT families. Genome-wide haplotypes and copy number profiles generated for each embryo using the new phasing approach were compared with the reference PGT method to evaluate the diagnostic concordance. PARTICIPANTS/MATERIALS, SETTING, METHODS This study included 12 couples with a known hereditary genetic disorder, participating in the comprehensive PGT program and with at least one parental sibling available (e.g. brother and/or sister). Genotyping data from both prospective parents and the parental sibling(s) were used to perform IBS-based phasing and to trace the disease-associated alleles. The outcome of the IBS-based PGT was compared with the results of the clinically implemented reference haplotyping-based PGT method. MAIN RESULTS AND THE ROLE OF CHANCE IBS-based haplotyping was performed for 12 PGT families. In accordance with the theoretical prediction of allele sharing between sibling pairs, 6 out of 12 (50%) couples or 23 out of 56 embryos could be phased using parental siblings. In families where phasing was possible, haplotype calling in the locus of interest was 100% concordant between the reference PGT method and IBS-based approach using parental siblings. LARGE SCALE DATA N/A LIMITATIONS, REASONS FOR CAUTION Phasing of parental haplotypes will only be possible when the disease locus lies in an informative region (categorized as IBS1). Phasing prospective parents using relatives with reduced genetic relatedness as a reference (e.g. siblings) decreases the size and the occurrence of informative IBS1 regions, necessary for haplotype calling. By including more than one extended family member, the chance of obtaining IBS1 coverage in the interrogated locus can be increased. A pre-PGT work-up can define whether the carrier couple could benefit from this approach. WIDER IMPLICATIONS OF THE FINDINGS Phasing by relatives extends the potential of comprehensive PGT, since it allows the inclusion of couples who do not have access to the standard phasing references, such as parents or offspring. STUDY FUNDING/COMPETING INTEREST(S) The study was funded by the KU Leuven grant (C14/18/092), Research Foundation Flanders (FWO; GA09311N), Horizon 2020 innovation programme (WIDENLIFE, 692065) and Agilent Technologies. J.R.V., T.V. and M.Z.E. are co-inventors of a patent ZL910050-PCT/EP2011/060211-WO/2011/157846 ‘Methods for haplotyping single-cells’ and ZL913096-PCT/EP2014/068315-WO/2015/028576 ‘Haplotyping and copy number typing using polymorphic variant allelic frequencies’ licensed to Agilent Technologies. The other authors have no conflict of interest to declare.