Runs of homozygosity in European populations

Genomic data illuminates demography, genetic structure and selection of a popular dog breed

BACKGROUND

Genomic methods have proved to be important tools in the analysis of genetic diversity across the range of species and can be used to reveal processes underlying both short- and long-term evolutionary change. This study applied genomic methods to investigate population structure and inbreeding in a common UK dog breed, the Labrador Retriever.

RESULTS

We found substantial within-breed genetic differentiation, which was associated with the role of the dog (i.e. working, pet, show) and also with coat colour (i.e. black, yellow, brown). There was little evidence of geographical differentiation. Highly differentiated genomic regions contained genes and markers associated with skull shape, suggesting that at least some of the differentiation is related to human-imposed selection on this trait. We also found that the total length of homozygous segments (runs of homozygosity, ROHs) was highly correlated with inbreeding coefficient.

CONCLUSIONS

This study demonstrates that high-density genomic data can be used to quantify genetic diversity and to decipher demographic and selection processes. Analysis of genetically differentiated regions in the UK Labrador Retriever population suggests the possibility of human-imposed selection on craniofacial characteristics. The high correlation between estimates of inbreeding from genomic and pedigree data for this breed demonstrates that genomic approaches can be used to quantify inbreeding levels in dogs, which will be particularly useful where pedigree information is missing.

Exome sequencing in schizophrenic patients with high levels of homozygosity identifies novel and extremely rare mutations in the GABA/glutamatergic pathways

Inbreeding is a known risk factor for recessive Mendelian diseases and previous studies have suggested that it could also play a role in complex disorders, such as psychiatric diseases. Recent inbreeding results in the presence of long runs of homozygosity (ROHs) along the genome, which are also defined as autozygosity regions. Genetic variants in these regions have two alleles that are identical by descent, thus increasing the odds of bearing rare recessive deleterious mutations due to a homozygous state. A recent study showed a suggestive enrichment of long ROHs in schizophrenic patients, suggesting that recent inbreeding could play a role in the disease. To better understand the impact of autozygosity on schizophrenia risk, we selected, from a cohort of 180 Italian patients, seven subjects with extremely high numbers of large ROHs that were likely due to recent inbreeding and characterized the mutational landscape within their ROHs using Whole Exome Sequencing and, gene set enrichment analysis. We identified a significant overlap (17%; empirical p-value = 0.0171) between genes inside ROHs affected by low frequency functional homozygous variants (107 genes) and the group of most promising candidate genes mutated in schizophrenia. Moreover, in four patients, we identified novel and extremely rare damaging mutations in the genes involved in neurodevelopment (MEGF8) and in GABA/glutamatergic synaptic transmission (GAD1, FMN1, ANO2). These results provide insights into the contribution of rare recessive mutations and inbreeding as risk factors for schizophrenia. ROHs that are likely due to recent inbreeding harbor a combination of predisposing low-frequency variants and extremely rare variants that have a high impact on pivotal biological pathways implicated in the disease. In addition, this study confirms that focusing on patients with high levels of homozygosity could be a useful prioritization strategy for discovering new high-impact mutations in genetically complex disorders.

Population specific analysis of Yakut exomes

We studied the genetic diversity of the Yakut population using exome sequencing. We performed comparative analysis of the Yakut population and the populations that are included in the "1000 Genomes" project and we identified the alleles specific to the Yakut population. We showed, that the Yakuts population is a separate cluster between Europeans and East Asians.

Recurrent rhinovirus infections in a child with inherited MDA5 deficiency

MDA5 is a cytosolic sensor of double-stranded RNA (ds)RNA including viral byproducts and intermediates. We studied a child with life-threatening, recurrent respiratory tract infections, caused by viruses including human rhinovirus (HRV), influenza virus, and respiratory syncytial virus (RSV). We identified in her a homozygous missense mutation in IFIH1 that encodes MDA5. Mutant MDA5 was expressed but did not recognize the synthetic MDA5 agonist/(ds)RNA mimic polyinosinic-polycytidylic acid. When overexpressed, mutant MDA5 failed to drive luciferase activity from the IFNB1 promoter or promoters containing ISRE or NF-κB sequence motifs. In respiratory epithelial cells or fibroblasts, wild-type but not knockdown of MDA5 restricted HRV infection while increasing IFN-stimulated gene expression and IFN-β/λ. However, wild-type MDA5 did not restrict influenza virus or RSV replication. Moreover, nasal epithelial cells from the patient, or fibroblasts gene-edited to express mutant MDA5, showed increased replication of HRV but not influenza or RSV. Thus, human MDA5 deficiency is a novel inborn error of innate and/or intrinsic immunity that causes impaired (ds)RNA sensing, reduced IFN induction, and susceptibility to the common cold virus.

Enrichment of low-frequency functional variants revealed by whole-genome sequencing of multiple isolated European populations

The genetic features of isolated populations can boost power in complex-trait association studies, and an in-depth understanding of how their genetic variation has been shaped by their demographic history can help leverage these advantageous characteristics. Here, we perform a comprehensive investigation using 3,059 newly generated low-depth whole-genome sequences from eight European isolates and two matched general populations, together with published data from the 1000 Genomes Project and UK10K. Sequencing data give deeper and richer insights into population demography and genetic characteristics than genotype-chip data, distinguishing related populations more effectively and allowing their functional variants to be studied more fully. We demonstrate relaxation of purifying selection in the isolates, leading to enrichment of rare and low-frequency functional variants, using novel statistics, DVxy and SVxy. We also develop an isolation-index (Isx) that predicts the overall level of such key genetic characteristics and can thus help guide population choice in future complex-trait association studies.

Isolated populations often have special genetic compositions that can be leveraged for genetic association studies. Here, Xue and colleagues generate and analyse 3,059 low-depth whole-genome sequences from eight European isolated populations and two matched general populations.

"Like sugar in milk": reconstructing the genetic history of the Parsi population

BACKGROUND

The Parsis are one of the smallest religious communities in the world. To understand the population structure and demographic history of this group in detail, we analyzed Indian and Pakistani Parsi populations using high-resolution genetic variation data on autosomal and uniparental loci (Y-chromosomal and mitochondrial DNA). Additionally, we also assayed mitochondrial DNA polymorphisms among ancient Parsi DNA samples excavated from Sanjan, in present day Gujarat, the place of their original settlement in India.

RESULTS

Among present-day populations, the Parsis are genetically closest to Iranian and the Caucasus populations rather than their South Asian neighbors. They also share the highest number of haplotypes with present-day Iranians and we estimate that the admixture of the Parsis with Indian populations occurred ~1,200 years ago. Enriched homozygosity in the Parsi reflects their recent isolation and inbreeding. We also observed 48% South-Asian-specific mitochondrial lineages among the ancient samples, which might have resulted from the assimilation of local females during the initial settlement. Finally, we show that Parsis are genetically closer to Neolithic Iranians than to modern Iranians, who have witnessed a more recent wave of admixture from the Near East.

CONCLUSIONS

Our results are consistent with the historically-recorded migration of the Parsi populations to South Asia in the 7th century and in agreement with their assimilation into the Indian sub-continent's population and cultural milieu "like sugar in milk". Moreover, in a wider context our results support a major demographic transition in West Asia due to the Islamic conquest.

A Genetic Population Isolate in The Netherlands Showing Extensive Haplotype Sharing and Long Regions of Homozygosity

Genetic isolated populations have features that may facilitate genetic analyses and can be leveraged to improve power of mapping genes to complex traits. Our aim was to test the extent to which a population with a former history of geographic isolation and religious endogamy, and currently with one of the highest fertility rates in The Netherlands, shows signs of genetic isolation. For this purpose, genome-wide genotype data was collected of 72 unrelated individuals from this population as well as in a sample of 104 random control subjects from The Netherlands. Additional reference data from different populations and population isolates was available through HapMap and the Human Genome Diversity Project. We performed a number of analyses to compare the genetic structure between these populations: we calculated the pairwise genetic distance between populations, examined the extent of identical-by-descent (IBD) sharing and estimated the effective population size. Genetic analysis of this population showed consistent patterns of a population isolate at all levels tested. We confirmed that this population is most closely related to the Dutch control subjects, and detected high levels of IBD sharing and runs of homozygosity at equal or even higher levels than observed in previously described population isolates. The effective population size of this population was estimated to be several orders of magnitude smaller than that of the Dutch control sample. We conclude that the geographic isolation of this population combined with rapid population growth has resulted in a genetic isolate with great potential value for future genetic studies.

The distribution of runs of homozygosity and selection signatures in six commercial meat sheep breeds

Domestication and the subsequent selection of animals for either economic or morphological features can leave a variety of imprints on the genome of a population. Genomic regions subjected to high selective pressures often show reduced genetic diversity and frequent runs of homozygosity (ROH). Therefore, the objective of the present study was to use 42,182 autosomal SNPs to identify genomic regions in 3,191 sheep from six commercial breeds subjected to selection pressure and to quantify the genetic diversity within each breed using ROH. In addition, the historical effective population size of each breed was also estimated and, in conjunction with ROH, was used to elucidate the demographic history of the six breeds. ROH were common in the autosomes of animals in the present study, but the observed breed differences in patterns of ROH length and burden suggested differences in breed effective population size and recent management. ROH provided a sufficient predictor of the pedigree inbreeding coefficient, with an estimated correlation between both measures of 0.62. Genomic regions under putative selection were identified using two complementary algorithms; the fixation index and hapFLK. The identified regions under putative selection included candidate genes associated with skin pigmentation, body size and muscle formation; such characteristics are often sought after in modern-day breeding programs. These regions of selection frequently overlapped with high ROH regions both within and across breeds. Multiple yet uncharacterised genes also resided within putative regions of selection. This further substantiates the need for a more comprehensive annotation of the sheep genome as these uncharacterised genes may contribute to traits of interest in the animal sciences. Despite this, the regions identified as under putative selection in the current study provide an insight into the mechanisms leading to breed differentiation and genetic variation in meat production.

Genome-wide mapping and estimation of inbreeding depression of semen quality traits in a cattle population

Inbreeding depression is known to affect quantitative traits such as male fertility and sperm quality, but the genetic basis for these associations is poorly understood. Most studies have been limited to examining how pedigree- or marker-derived genome-wide autozygosity is associated with quantitative phenotypes. In this study, we analyzed possible associations of genetic features of inbreeding depression with percentage of live spermatozoa and total number of spermatozoa in 19,720 ejaculates obtained from 554 Austrian Fleckvieh bulls during routine artificial insemination programs. Genome-wide inbreeding depression was estimated and genomic regions contributing to inbreeding depression were mapped. Inbreeding depression did affect total number of spermatozoa, and such depression was predicted by pedigree-based inbreeding levels and genome-wide inbreeding levels based on runs of homozygosity (ROH). Genome-wide inbreeding depression did not seem to affect percentage of live spermatozoa. A model incorporating genetic effects of the bull, environmental factors, and additive genetic and ROH status effects of individual single-nucleotide polymorphisms revealed genomic regions significantly associated with ROH status for total number of spermatozoa (4 regions) or percentage of live spermatozoa (5 regions). All but one region contains genes related to spermatogenesis and sperm morphology. These genomic regions contain genes affecting sperm morphogenesis and efficacy. The results highlight that next-generation sequencing may help explain some of the genetic factors contributing to inbreeding depression of sperm quality traits in Fleckvieh bulls.

GeneImp: Fast Imputation to Large Reference Panels Using Genotype Likelihoods from Ultralow Coverage Sequencing

We address the task of genotype imputation to a dense reference panel given genotype likelihoods computed from ultralow coverage sequencing as inputs. In this setting, the data have a high-level of missingness or uncertainty, and are thus more amenable to a probabilistic representation. Most existing imputation algorithms are not well suited for this situation, as they rely on prephasing for computational efficiency, and, without definite genotype calls, the prephasing task becomes computationally expensive. We describe GeneImp, a program for genotype imputation that does not require prephasing and is computationally tractable for whole-genome imputation. GeneImp does not explicitly model recombination, instead it capitalizes on the existence of large reference panels-comprising thousands of reference haplotypes-and assumes that the reference haplotypes can adequately represent the target haplotypes over short regions unaltered. We validate GeneImp based on data from ultralow coverage sequencing (0.5×), and compare its performance to the most recent version of BEAGLE that can perform this task. We show that GeneImp achieves imputation quality very close to that of BEAGLE, using one to two orders of magnitude less time, without an increase in memory complexity. Therefore, GeneImp is the first practical choice for whole-genome imputation to a dense reference panel when prephasing cannot be applied, for instance, in datasets produced via ultralow coverage sequencing. A related future application for GeneImp is whole-genome imputation based on the off-target reads from deep whole-exome sequencing.

Genomic insights into the population structure and history of the Irish Travellers

The Irish Travellers are a population with a history of nomadism; consanguineous unions are common and they are socially isolated from the surrounding, ‘settled’ Irish people. Low-resolution genetic analysis suggests a common Irish origin between the settled and the Traveller populations. What is not known, however, is the extent of population structure within the Irish Travellers, the time of divergence from the general Irish population, or the extent of autozygosity. Using a sample of 50 Irish Travellers, 143 European Roma, 2232 settled Irish, 2039 British and 6255 European or world-wide individuals, we demonstrate evidence for population substructure within the Irish Traveller population, and estimate a time of divergence before the Great Famine of 1845–1852. We quantify the high levels of autozygosity, which are comparable to levels previously described in Orcadian 1^st/2^nd cousin offspring, and finally show the Irish Traveller population has no particular genetic links to the European Roma. The levels of autozygosity and distinct Irish origins have implications for disease mapping within Ireland, while the population structure and divergence inform on social history.

Overcoming the dichotomy between open and isolated populations using genomic data from a large European dataset

Human populations are often dichotomized into “isolated” and “open” categories using cultural and/or geographical barriers to gene flow as differential criteria. Although widespread, the use of these alternative categories could obscure further heterogeneity due to inter-population differences in effective size, growth rate, and timing or amount of gene flow. We compared intra and inter-population variation measures combining novel and literature data relative to 87,818 autosomal SNPs in 14 open populations and 10 geographic and/or linguistic European isolates. Patterns of intra-population diversity were found to vary considerably more among isolates, probably due to differential levels of drift and inbreeding. The relatively large effective size estimated for some population isolates challenges the generalized view that they originate from small founding groups. Principal component scores based on measures of intra-population variation of isolated and open populations were found to be distributed along a continuum, with an area of intersection between the two groups. Patterns of inter-population diversity were even closer, as we were able to detect some differences between population groups only for a few multidimensional scaling dimensions. Therefore, different lines of evidence suggest that dichotomizing human populations into open and isolated groups fails to capture the actual relations among their genomic features.

Genetic diversity, extent of linkage disequilibrium and persistence of gametic phase in Canadian pigs

Background

Knowledge on the levels of linkage disequilibrium (LD) across the genome, persistence of gametic phase between breed pairs, genetic diversity and population structure are important parameters for the successful implementation of genomic selection. Therefore, the objectives of this study were to investigate these parameters in order to assess the feasibility of a multi-herd and multi-breed training population for genomic selection in important purebred and crossbred pig populations in Canada. A total of 3,057 animals, representative of the national populations, were genotyped with the Illumina Porcine SNP60 BeadChip (62,163 markers).

Results

The overall LD (r²) between adjacent SNPs was 0.49, 0.38, 0.40 and 0.31 for Duroc, Landrace, Yorkshire and Crossbred (Landrace x Yorkshire) populations, respectively. The highest correlation of phase (r) across breeds was observed between Crossbred animals and either Landrace or Yorkshire breeds, in which r was approximately 0.80 at 1 Mbp of distance. Landrace and Yorkshire breeds presented r ≥ 0.80 in distances up to 0.1 Mbp, while Duroc breed showed r ≥ 0.80 for distances up to 0.03 Mbp with all other populations. The persistence of phase across herds were strong for all breeds, with r ≥ 0.80 up to 1.81 Mbp for Yorkshire, 1.20 Mbp for Duroc, and 0.70 Mbp for Landrace. The first two principal components clearly discriminate all the breeds. Similar levels of genetic diversity were observed among all breed groups. The current effective population size was equal to 75 for Duroc and 92 for both Landrace and Yorkshire.

Conclusions

An overview of population structure, LD decay, demographic history and inbreeding of important pig breeds in Canada was presented. The rate of LD decay for the three Canadian pig breeds indicates that genomic selection can be successfully implemented within breeds with the current 60 K SNP panel. The use of a multi-breed training population involving Landrace and Yorkshire to estimate the genomic breeding values of crossbred animals (Landrace × Yorkshire) should be further evaluated. The lower correlation of phase at short distances between Duroc and the other breeds indicates that a denser panel may be required for the use of a multi-breed training population including Duroc.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0473-y) contains supplementary material, which is available to authorized users.

Genetic structure in the Sherpa and neighboring Nepalese populations

Background

We set out to describe the fine-scale population structure across the Eastern region of Nepal. To date there is relatively little known about the genetic structure of the Sherpa residing in Nepal and their genetic relationship with the Nepalese. We assembled dense genotype data from a total of 1245 individuals representing Nepal and a variety of different populations resident across the greater Himalayan region including Tibet, China, India, Pakistan, Kazakhstan, Uzbekistan, Tajikistan and Kirghizstan. We performed analysis of principal components, admixture and homozygosity.

Results

We identified clear substructure across populations resident in the Himalayan arc, with genetic structure broadly mirroring geographical features of the region. Ethnic subgroups within Nepal show distinct genetic structure, on both admixture and principal component analysis. We detected differential proportions of ancestry from northern Himalayan populations across Nepalese subgroups, with the Nepalese Rai, Magar and Tamang carrying the greatest proportions of Tibetan ancestry.

Conclusions

We show that populations dwelling on the Himalayan plateau have had a clear impact on the Northern Indian gene pool. We illustrate how the Sherpa are a remarkably isolated population, with little gene flow from surrounding Nepalese populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3469-5) contains supplementary material, which is available to authorized users.

Inferring Individual Inbreeding and Demographic History from Segments of Identity by Descent in Ficedula Flycatcher Genome Sequences

Individual inbreeding and historical demography can be estimated by analyzing runs of homozygosity (ROH), which are indicative of chromosomal segments of identity by descent (IBD). Such analyses have so far been rare in natural populations due to limited genomic resources. We analyzed ROH in whole genome sequences from 287 Ficedula flycatchers representing four species, with the objectives of evaluating the causes of genome-wide variation in the abundance of ROH and inferring historical demography. ROH were clearly more abundant in genomic regions with low recombination rate. However, this pattern was substantially weaker when ROH were mapped using genetic rather than physical single nucleotide polymorphism (SNP) coordinates in the genome. Empirical results and simulations suggest that high ROH abundance in regions of low recombination was partly caused by increased power to detect the very long IBD segments typical of regions with a low recombination rate. Simulations also showed that hard selective sweeps (but not soft sweeps or background selection) likely contributed to variation in the abundance of ROH across the genome. Comparisons of the abundance of ROH among several study populations indicated that the Spanish pied flycatcher population had the smallest historical effective population size (N_e) for this species, and that a putatively recently founded island (Baltic) population had the smallest historical N_e among the collared flycatchers. Analysis of pairwise IBD in Baltic collared flycatchers indicated that this population was founded <60 generations ago. This study provides a rare genomic glimpse into demographic history and the mechanisms underlying the genome-wide distribution of ROH.

Genome-wide analysis in endangered populations: a case study in Barbaresca sheep

Analysis of genomic data is becoming increasingly common in the livestock industry and the findings have been an invaluable resource for effective management of breeding programs in small and endangered populations. In this paper, with the goal of highlighting the potential of genomic analysis for small and endangered populations, genome-wide levels of linkage disequilibrium, measured as the squared correlation coefficient of allele frequencies at a pair of loci, effective population size, runs of homozygosity (ROH) and genetic diversity parameters, were estimated in Barbaresca sheep using Illumina OvineSNP50K array data. Moreover, the breed's genetic structure and its relationship with other breeds were investigated. Levels of pairwise linkage disequilibrium decreased with increasing distance between single nucleotide polymorphisms. An average correlation coefficient <0.25 was found for markers located up to 50 kb apart. Therefore, these results support the need to use denser single nucleotide polymorphism panels for high power association mapping and genomic selection efficiency in future breeding programs. The estimate of past effective population size ranged from 747 animals 250 generations ago to 28 animals five generations ago, whereas the contemporary effective population size was 25 animals. A total of 637 ROH were identified, most of which were short (67%) and ranged from 1 to 10 Mb. The genetic analyses revealed that the Barbaresca breed tended to display lower variability than other Sicilian breeds. Recent inbreeding was evident, according to the ROH analysis. All the investigated parameters showed a comparatively narrow genetic base and indicated an endangered status for Barbaresca. Multidimensional scaling, model-based clustering, measurement of population differentiation, neighbor networks and haplotype sharing distinguished Barbaresca from other breeds, showed a low level of admixture with the other breeds considered in this study, and indicated clear genetic differences compared with other breeds. Attention should be given to the conservation of Barbaresca due to its critical conservation status. In this context, genomic information may have a crucial role in management of small and endangered populations.

Automated Integration of a UPLC Glycomic Profile

Ultra-performance liquid chromatography (UPLC) is the established technology for accurate analysis of IgG Fc N-glycosylation due to its superior sensitivity, resolution, speed, and its capability to provide branch-specific information of glycan species. Correct and cost-efficient preprocessing of chromatographic data is the major prerequisite for subsequent analyses ranging from inference of structural isomers to biomarker discovery and prediction of humoral immune response from characterized changes in glycosylation. The complexity of glycomic chromatograms poses a number of challenges for developing automated data annotation and quantitation algorithms, which frequently necessitated manual or semi-manual approaches to preprocessing, most notably to peak detection and integration. Such procedures are meticulous and time-consuming, and may be a source of confounding due to their dependence on human labelers. Although liquid chromatography is a mature field and a number of methods have been developed for automatic peak detection outside the area of glycomics analysis, we found that hardly any of them are suitable for automatic integration of UPLC glycomic profiles without substantial modifications. In this chapter, we illustrate practical challenges of automatic peak detection of UPLC glycomics chromatograms. We outline a robust, semi-supervised method ACE (Automatic Chromatogram Extraction) for automated alignment and detection of glycan peaks in chromatograms, developed by Pharmatics Limited (UK) in collaboration with Genos Limited (Croatia). Application of the tool requires minimal human interference, which results in a significant reduction in the time and cost of IgG glycomics signal integration using Waters Acquity UPLC instrument (Milford, MA, USA) in several human cohorts with blind technical replicas.

Genetics of consanguinity and inbreeding in health and disease

CONTEXT

Inbreeding increases the level of homozygotes for autosomal recessive disorders and is the major objective in clinical studies. The prevalence of consanguinity and the degree of inbreeding vary from one population to another depending on ethnicity, religion, culture and geography. Global epidemiological studies have revealed that consanguineous unions have been significantly associated with increased susceptibility to various forms of inherited diseases.

OBJECTIVE

The study aimed to determine the role of consanguinity in human health and to highlight the associated risks for various diseases or disorders.

METHODS

PubMed and Google Scholar search engines were used to explore the published literature on consanguinity and its associated risks using the key words "consanguinity", "prevalence", "inbreeding depression", "coefficient of inbreeding", "child health", "mortality", "human health", "homozygosity" and "complex diseases" in different combinations. The studies were screened for eligibility on the basis of their epidemiological relevance.

RESULTS

This comprehensive assessment highlights the deleterious consequences in populations with a higher prevalence of consanguinity among different countries worldwide.

CONCLUSIONS

To avoid the inbreeding load there is the need to improve socioeconomic and educational status and to increase public awareness of reproductive health and anticipated deleterious effects. Pre-marital and pre-conception counselling of consanguineous populations should be an integral part of health policy to train people and make people aware of its harmful consequences. Furthermore, runs of homozygosity (ROH) and whole-exome sequencing (WES) are useful tools in exploring new genomic signatures for the cause of inbreeding depression.

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH ) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED ). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome-wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.

Structural and genomic variation in preterm birth

BACKGROUND

Runs of homozygosity (ROH) are consecutive homozygous genotypes, which may result from population inbreeding or consanguineous marriages. ROH enhance the expression of recessive traits.

METHODS

We mapped ROH in a case control study of women delivering at term compared with women delivering at or before 34 wk gestation. Gene sets known to be important in risk of preterm birth were examined for their overlap with identified ROH segments.

RESULTS

While we found no evidence of increased burden of ROH or copy number variations in mothers delivering at or before 34 wk compared with term, we identified 424 genome-wide 50 kb segments with significant difference in abundance of overlapping ROH segments in cases vs. controls, P < 0.05. These regions overlap 199 known genes. We found preterm birth associated genes (CXCR4, MYLK, PAK1) and genes shown to have an evolutionary link to preterm (CXCR4, PPP3CB, C6orf57, DUSP13, and SLC25A45) with significant differences in abundance of overlapping ROH blocks in cases vs. controls, P < 0.001.

CONCLUSION

We conclude, while we found no significant burden of ROH, we did identify genomic regions with significantly greater abundance of ROH blocks in women delivering preterm that overlapped genes known to be involved in preterm birth.

Genomics advances the study of inbreeding depression in the wild

Inbreeding depression (reduced fitness of individuals with related parents) has long been a major focus of ecology, evolution, and conservation biology. Despite decades of research, we still have a limited understanding of the strength, underlying genetic mechanisms, and demographic consequences of inbreeding depression in the wild. Studying inbreeding depression in natural populations has been hampered by the inability to precisely measure individual inbreeding. Fortunately, the rapidly increasing availability of high-throughput sequencing data means it is now feasible to measure the inbreeding of any individual with high precision. Here, we review how genomic data are advancing our understanding of inbreeding depression in the wild. Recent results show that individual inbreeding and inbreeding depression can be measured more precisely with genomic data than via traditional pedigree analysis. Additionally, the availability of genomic data has made it possible to pinpoint loci with large effects contributing to inbreeding depression in wild populations, although this will continue to be a challenging task in many study systems due to low statistical power. Now that reliably measuring individual inbreeding is no longer a limitation, a major focus of future studies should be to more accurately quantify effects of inbreeding depression on population growth and viability.

Meiotic recombination shapes precision of pedigree- and marker-based estimates of inbreeding

The proportion of an individual's genome that is identical by descent (GWIBD) can be estimated from pedigrees (inbreeding coefficient 'Pedigree F') or molecular markers ('Marker F'), but both estimators come with error. Assuming unrelated pedigree founders, Pedigree F is the expected proportion of GWIBD given a specific inbreeding constellation. Meiotic recombination introduces variation around that expectation (Mendelian noise) and related pedigree founders systematically bias Pedigree F downward. Marker F is an estimate of the actual proportion of GWIBD but it suffers from the sampling error of markers plus the error that occurs when a marker is homozygous without reflecting common ancestry (identical by state). We here show via simulation of a zebra finch and a human linkage map that three aspects of meiotic recombination (independent assortment of chromosomes, number of crossovers and their distribution along chromosomes) contribute to variation in GWIBD and thus the precision of Pedigree and Marker F. In zebra finches, where the genome contains large blocks that are rarely broken up by recombination, the Mendelian noise was large (nearly twofold larger s.d. values compared with humans) and Pedigree F thus less precise than in humans, where crossovers are distributed more uniformly along chromosomes. Effects of meiotic recombination on Marker F were reversed, such that the same number of molecular markers yielded more precise estimates of GWIBD in zebra finches than in humans. As a consequence, in species inheriting large blocks that rarely recombine, even small numbers of microsatellite markers will often be more informative about inbreeding and fitness than large pedigrees.

Novel Graphical Analyses of Runs of Homozygosity among Species and Livestock Breeds

Runs of homozygosity (ROH), uninterrupted stretches of homozygous genotypes resulting from parents transmitting identical haplotypes to their offspring, have emerged as informative genome-wide estimates of autozygosity (inbreeding). We used genomic profiles based on 698 K single nucleotide polymorphisms (SNPs) from nine breeds of domestic cattle (Bos taurus) and the European bison (Bison bonasus) to investigate how ROH distributions can be compared within and among species. We focused on two length classes: 0.5-15 Mb to investigate ancient events and >15 Mb to address recent events (approximately three generations). For each length class, we chose a few chromosomes with a high number of ROH, calculated the percentage of times a SNP appeared in a ROH, and plotted the results. We selected areas with distinct patterns including regions where (1) all groups revealed an increase or decrease of ROH, (2) bison differed from cattle, (3) one cattle breed or groups of breeds differed (e.g., dairy versus meat cattle). Examination of these regions in the cattle genome showed genes potentially important for natural and human-induced selection, concerning, for example, meat and milk quality, metabolism, growth, and immune function. The comparative methodology presented here permits visual identification of regions of interest for selection, breeding programs, and conservation.

CONSANGUINITY BY RANDOM ISONYMY AND SOCIOECONOMIC DEVELOPMENT IN ARGENTINA: A POPULATION STUDY

In human populations various flexible, labile and interdependent structures (genetic, demographic, socioeconomic) co-exist, each of which can be organized in an hierarchical order corresponding to administrative entities. The relationship between consanguinity, as estimated by random isonymy (F ST), and socioeconomic conditions was analysed at different levels of political and administrative organization in Argentina. From the surnames of 22,666,139 voters from the 2001 electoral roll, F ST was estimated for 510 Argentinian departments. Using a principal component analysis, a Socio-Demographic and Economic Indicator (SDEI), summarizing the effect of 22 socioeconomic and demographic variables at the departmental level, was computed. The relationship between departmental F ST and SDEI values was analysed for the whole nation and within regions using multiple regression analysis. The F ST presented a clinal distribution with the highest values in the north and west of the country, while SDEI expressed the opposite behaviour. A negative and significant correlation was observed between F ST and SDEI, accounting for 46% of the variation in consanguinity in Argentina. The strongest correlations of F ST with SDEI were observed in the Central, Patagonia and Cuyo regions, i.e. those with the highest values of SDEI and lowest values of F ST.

No Reliable Association between Runs of Homozygosity and Schizophrenia in a Well-Powered Replication Study

It is well known that inbreeding increases the risk of recessive monogenic diseases, but it is less certain whether it contributes to the etiology of complex diseases such as schizophrenia. One way to estimate the effects of inbreeding is to examine the association between disease diagnosis and genome-wide autozygosity estimated using runs of homozygosity (ROH) in genome-wide single nucleotide polymorphism arrays. Using data for schizophrenia from the Psychiatric Genomics Consortium (n = 21,868), Keller et al. (2012) estimated that the odds of developing schizophrenia increased by approximately 17% for every additional percent of the genome that is autozygous (β = 16.1, CI(β) = [6.93, 25.7], Z = 3.44, p = 0.0006). Here we describe replication results from 22 independent schizophrenia case-control datasets from the Psychiatric Genomics Consortium (n = 39,830). Using the same ROH calling thresholds and procedures as Keller et al. (2012), we were unable to replicate the significant association between ROH burden and schizophrenia in the independent PGC phase II data, although the effect was in the predicted direction, and the combined (original + replication) dataset yielded an attenuated but significant relationship between Froh and schizophrenia (β = 4.86,CI(β) = [0.90,8.83],Z = 2.40,p = 0.02). Since Keller et al. (2012), several studies reported inconsistent association of ROH burden with complex traits, particularly in case-control data. These conflicting results might suggest that the effects of autozygosity are confounded by various factors, such as socioeconomic status, education, urbanicity, and religiosity, which may be associated with both real inbreeding and the outcome measures of interest.

Dual T cell- and B cell-intrinsic deficiency in humans with biallelic RLTPR mutations

In two complementary papers, Casanova, Malissen, and collaborators report the discovery of human RLTPR deficiency, the first primary immunodeficiency of the human CD28 pathway in T cells. Together, the two studies highlight the important and largely (but not completely) overlapping roles of RLTPR in T and B cells of humans and mice.

Combined immunodeficiency (CID) refers to inborn errors of human T cells that also affect B cells because of the T cell deficit or an additional B cell–intrinsic deficit. In this study, we report six patients from three unrelated families with biallelic loss-of-function mutations in RLTPR, the mouse orthologue of which is essential for CD28 signaling. The patients have cutaneous and pulmonary allergy, as well as a variety of bacterial and fungal infectious diseases, including invasive tuberculosis and mucocutaneous candidiasis. Proportions of circulating regulatory T cells and memory CD4⁺ T cells are reduced. Their CD4⁺ T cells do not respond to CD28 stimulation. Their CD4⁺ T cells exhibit a "Th2" cell bias ex vivo and when cultured in vitro, contrasting with the paucity of "Th1," "Th17," and T follicular helper cells. The patients also display few memory B cells and poor antibody responses. This B cell phenotype does not result solely from the T cell deficiency, as the patients’ B cells fail to activate NF-κB upon B cell receptor (BCR) stimulation. Human RLTPR deficiency is a CID affecting at least the CD28-responsive pathway in T cells and the BCR-responsive pathway in B cells.

East Eurasian ancestry in the middle of Europe: genetic footprints of Steppe nomads in the genomes of Belarusian Lipka Tatars

Medieval era encounters of nomadic groups of the Eurasian Steppe and largely sedentary East Europeans had a variety of demographic and cultural consequences. Amongst these outcomes was the emergence of the Lipka Tatars—a Slavic-speaking Sunni-Muslim minority residing in modern Belarus, Lithuania and Poland, whose ancestors arrived in these territories via several migration waves, mainly from the Golden Horde. Our results show that Belarusian Lipka Tatars share a substantial part of their gene pool with Europeans as indicated by their Y-chromosomal, mitochondrial and autosomal DNA variation. Nevertheless, Belarusian Lipkas still retain a strong genetic signal of their nomadic ancestry, witnessed by the presence of common Y-chromosomal and mitochondrial DNA variants as well as autosomal segments identical by descent between Lipkas and East Eurasians from temperate and northern regions. Hence, we document Lipka Tatars as a unique example of former Medieval migrants into Central Europe, who became sedentary, changed language to Slavic, yet preserved their faith and retained, both uni- and bi-parentally, a clear genetic echo of a complex population interplay throughout the Eurasian Steppe Belt, extending from Central Europe to northern China.

Genomic analysis reveals depression due to both individual and maternal inbreeding in a free-living mammal population

There is ample evidence for inbreeding depression manifested as a reduction in fitness or fitness-related traits in the focal individual. In many organisms, fitness is not only affected by genes carried by the individual, but also by genes carried by their parents, for example if receiving parental care. While maternal effects have been described in many systems, the extent to which inbreeding affects fitness directly through the focal individual, or indirectly through the inbreeding coefficients of its parents, has rarely been examined jointly. The Soay sheep study population is an excellent system in which to test for both effects, as lambs receive extended maternal care. Here, we tested for both maternal and individual inbreeding depression in three fitness-related traits (birthweight and weight and hindleg length at 4 months of age) and three fitness components (first-year survival, adult annual survival and annual breeding success), using either pedigree-derived inbreeding or genomic estimators calculated using ~37 000 SNP markers. We found evidence for inbreeding depression in 4-month hindleg and weight, first-year survival in males, and annual survival and breeding success in adults. Maternal inbreeding was found to depress both birthweight and 4-month weight. We detected more instances of significant inbreeding depression using genomic estimators than the pedigree, which is partly explained through the increased sample sizes available. In conclusion, our results highlight that cross-generational inbreeding effects warrant further exploration in species with parental care and that modern genomic tools can be used successfully instead of, or alongside, pedigrees in natural populations.

Genome-wide autozygosity is associated with lower general cognitive ability

Inbreeding depression refers to lower fitness among offspring of genetic relatives. This reduced fitness is caused by the inheritance of two identical chromosomal segments (autozygosity) across the genome, which may expose the effects of (partially) recessive deleterious mutations. Even among outbred populations, autozygosity can occur to varying degrees due to cryptic relatedness between parents. Using dense genome-wide single-nucleotide polymorphism (SNP) data, we examined the degree to which autozygosity associated with measured cognitive ability in an unselected sample of 4854 participants of European ancestry. We used runs of homozygosity-multiple homozygous SNPs in a row-to estimate autozygous tracts across the genome. We found that increased levels of autozygosity predicted lower general cognitive ability, and estimate a drop of 0.6 s.d. among the offspring of first cousins (P=0.003-0.02 depending on the model). This effect came predominantly from long and rare autozygous tracts, which theory predicts as more likely to be deleterious than short and common tracts. Association mapping of autozygous tracts did not reveal any specific regions that were predictive beyond chance after correcting for multiple testing genome wide. The observed effect size is consistent with studies of cognitive decline among offspring of known consanguineous relationships. These findings suggest a role for multiple recessive or partially recessive alleles in general cognitive ability, and that alleles decreasing general cognitive ability have been selected against over evolutionary time.

Whole-exome sequencing to analyze population structure, parental inbreeding, and familial linkage

Principal component analysis (PCA), homozygosity rate estimations, and linkage studies in humans are classically conducted through genome-wide single-nucleotide variant arrays (GWSA). We compared whole-exome sequencing (WES) and GWSA for this purpose. We analyzed 110 subjects originating from different regions of the world, including North Africa and the Middle East, which are poorly covered by public databases and have high consanguinity rates. We tested and applied a number of quality control (QC) filters. Compared with GWSA, we found that WES provided an accurate prediction of population substructure using variants with a minor allele frequency > 2% (correlation = 0.89 with the PCA coordinates obtained by GWSA). WES also yielded highly reliable estimates of homozygosity rates using runs of homozygosity with a 1,000-kb window (correlation = 0.94 with the estimates provided by GWSA). Finally, homozygosity mapping analyses in 15 families including a single offspring with high homozygosity rates showed that WES provided 51% less genome-wide linkage information than GWSA overall but 97% more information for the coding regions. At the genome-wide scale, 76.3% of linked regions were found by both GWSA and WES, 17.7% were found by GWSA only, and 6.0% were found by WES only. For coding regions, the corresponding percentages were 83.5%, 7.4%, and 9.1%, respectively. With appropriate QC filters, WES can be used for PCA and adjustment for population substructure, estimating homozygosity rates in individuals, and powerful linkage analyses, particularly in coding regions.

Farming, Foreign Holidays, and Vitamin D in Orkney

Orkney, north of mainland Scotland, has the world's highest prevalence of multiple sclerosis (MS); vitamin D deficiency, a marker of low UV exposure, is also common in Scotland. Strong associations have been identified between vitamin D deficiency and MS, and between UV exposure and MS independent of vitamin D, although causal relationships remain to be confirmed. We aimed to compare plasma 25-hydroxyvitamin D levels in Orkney and mainland Scotland, and establish the determinants of vitamin D status in Orkney. We compared mean vitamin D and prevalence of deficiency in cross-sectional study data from participants in the Orkney Complex Disease Study (ORCADES) and controls in the Scottish Colorectal Cancer Study (SOCCS). We used multivariable regression to identify factors associated with vitamin D levels in Orkney. Mean (standard deviation) vitamin D was significantly higher among ORCADES than SOCCS participants (35.3 (18.0) and 31.7 (21.2), respectively). Prevalence of severe vitamin D deficiency was lower in ORCADES than SOCCS participants (6.6% to 16.2% p = 1.1 x 10(-15)). Older age, farming occupations and foreign holidays were significantly associated with higher vitamin D in Orkney. Although mean vitamin D levels are higher in Orkney than mainland Scotland, this masks variation within the Orkney population which may influence MS risk.

Evidence of Inbreeding in Hodgkin Lymphoma

Genome-wide association studies (GWASs) have identified several, mainly co-dominantly acting, single-nucleotide polymorphisms (SNPs) associated with Hodgkin lymphoma (HL). We searched for recessively acting disease loci by performing an analysis of runs of homozygosity (ROH) based on windows of homozygous SNP-blocks and by calculating genomic inbreeding coefficients on a SNP-wise basis. We used data from a previous GWAS with 906 cases and 1217 controls from a population with a long history of no matings between relatives. Ten recurrent ROHs were identified among 25 055 ROHs across all individuals but their association with HL was not genome-wide significant. All recurrent ROHs showed significant evidence for natural selection. As a novel finding genomic inbreeding among cases was significantly higher than among controls (P = 2.11*10-14) even after correcting for covariates. Higher inbreeding among the cases was mainly based on a group of individuals with a higher average length of ROHs per person. This result suggests a correlation of higher levels of inbreeding with higher cancer incidence and might reflect the existence of recessive alleles causing HL. Genomic inbreeding may result in a higher expression of deleterious recessive genes within a population.

Whole-exome sequencing in an isolated population from the Dalmatian island of Vis

We have whole-exome sequenced 176 individuals from the isolated population of the island of Vis in Croatia in order to describe exonic variation architecture. We found 290 577 single nucleotide variants (SNVs), 65% of which are singletons, low frequency or rare variants. A total of 25 430 (9%) SNVs are novel, previously not catalogued in NHLBI GO Exome Sequencing Project, UK10K-Generation Scotland, 1000Genomes Project, ExAC or NCBI Reference Assembly dbSNP. The majority of these variants (76%) are singletons. Comparable to data obtained from UK10K-Generation Scotland that were sequenced and analysed using the same protocols, we detected an enrichment of potentially damaging variants (non-synonymous and loss-of-function) in the low frequency and common variant categories. On average 115 (range 93–140) genotypes with loss-of-function variants, 23 (15–34) of which were homozygous, were identified per person. The landscape of loss-of-function variants across an exome revealed that variants mainly accumulated in genes on the xenobiotic-related pathways, of which majority coded for enzymes. The frequency of loss-of-function variants was additionally increased in Vis runs of homozygosity regions where variants mainly affected signalling pathways. This work confirms the isolate status of Vis population by means of whole-exome sequence and reveals the pattern of loss-of-function mutations, which resembles the trails of adaptive evolution that were found in other species. By cataloguing the exomic variants and describing the allelic structure of the Vis population, this study will serve as a valuable resource for future genetic studies of human diseases, population genetics and evolution in this population.

The Association Between Glycosylation of Immunoglobulin G and Hypertension: A Multiple Ethnic Cross-Sectional Study

More than half of all known proteins, and almost all membrane and extra-cellular proteins have oligosaccharide structures or glycans attached to them. Defects in glycosylation pathways are directly involved in at least 30 severe human diseases.A multiple center cross-sectional study (China, Croatia, and Scotland) was carried out to investigate the possible association between hypertension and IgG glycosylation. A hydrophilic interaction chromatography of fluorescently labeled glycans was used to analyze N-glycans attached to IgG in plasma samples from a total of 4757 individuals of Chinese Han, Croatian, and Scottish ethnicity.Five glycans (IgG with digalactosylated glycans) significantly differed in participants with prehypertension or hypertension compared to those with normal blood pressure, while additional 17 glycan traits were only significantly differed in participants with hypertension compared to those of normal blood pressure. These glycans were also significant correlated with systolic blood pressure (SBP) or diastolic blood pressure (DBP).The present study demonstrated for the 1st time an association between hypertension and IgG glycome composition. These findings suggest that the individual variation in N-glycosylation of IgG contributes to pathogenesis of hypertension, presumably via its effect on pro- and/or anti-inflammatory pathways.

Genomic signatures of population decline in the malaria mosquito Anopheles gambiae

Background

Population genomic features such as nucleotide diversity and linkage disequilibrium are expected to be strongly shaped by changes in population size, and might therefore be useful for monitoring the success of a control campaign. In the Kilifi district of Kenya, there has been a marked decline in the abundance of the malaria vector Anopheles gambiae subsequent to the rollout of insecticide-treated bed nets.

Methods

To investigate whether this decline left a detectable population genomic signature, simulations were performed to compare the effect of population crashes on nucleotide diversity, Tajima’s D, and linkage disequilibrium (as measured by the population recombination parameter ρ). Linkage disequilibrium and ρ were estimated for An. gambiae from Kilifi, and compared them to values for Anopheles arabiensis and Anopheles merus at the same location, and for An. gambiae in a location 200 km from Kilifi.

Results

In the first simulations ρ changed more rapidly after a population crash than the other statistics, and therefore is a more sensitive indicator of recent population decline. In the empirical data, linkage disequilibrium extends 100–1000 times further, and ρ is 100–1000 times smaller, for the Kilifi population of An. gambiae than for any of the other populations. There were also significant runs of homozygosity in many of the individual An. gambiae mosquitoes from Kilifi.

Conclusions

These results support the hypothesis that the recent decline in An. gambiae was driven by the rollout of bed nets. Measuring population genomic parameters in a small sample of individuals before, during and after vector or pest control may be a valuable method of tracking the effectiveness of interventions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1214-9) contains supplementary material, which is available to authorized users.

Runs of homozygosity and inbreeding in thyroid cancer

BACKGROUND

Genome-wide association studies (GWASs) have identified several single-nucleotide polymorphisms (SNPs) influencing the risk of thyroid cancer (TC). Most cancer predisposition genes identified through GWASs function in a co-dominant manner, and studies have not found evidence for recessively functioning disease loci in TC. Our study examines whether homozygosity is associated with an increased risk of TC and searches for novel recessively acting disease loci.

METHODS

Data from a previously conducted GWAS were used for the estimation of the proportion of phenotypic variance explained by all common SNPs, the detection of runs of homozygosity (ROH) and the determination of inbreeding to unravel their influence on TC.

RESULTS

Inbreeding coefficients were significantly higher among cases than controls. Association on a SNP-by-SNP basis was controlled by using the false discovery rate at a level of q* < 0.05, with 34 SNPs representing true differences in homozygosity between cases and controls. The average size, the number and total length of ROHs per person were significantly higher in cases than in controls. A total of 16 recurrent ROHs of rather short length were identified although their association with TC risk was not significant at a genome-wide level. Several recurrent ROHs harbor genes associated with risk of TC. All of the ROHs showed significant evidence for natural selection (iHS, Fst, Fay and Wu's H).

CONCLUSIONS

Our results support the existence of recessive alleles in TC susceptibility. Although regions of homozygosity were rather small, it might be possible that variants within these ROHs affect TC risk and may function in a recessive manner.

Homozygous NOTCH3 null mutation and impaired NOTCH3 signaling in recessive early-onset arteriopathy and cavitating leukoencephalopathy

Notch signaling is essential for vascular physiology. Neomorphic heterozygous mutations in NOTCH3, one of the four human NOTCH receptors, cause cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL). Hypomorphic heterozygous alleles have been occasionally described in association with a spectrum of cerebrovascular phenotypes overlapping CADASIL, but their pathogenic potential is unclear. We describe a patient with childhood-onset arteriopathy, cavitating leukoencephalopathy with cerebral white matter abnormalities presented as diffuse cavitations, multiple lacunar infarctions and disseminated microbleeds. We identified a novel homozygous c.C2898A (p.C966*) null mutation in NOTCH3 abolishing NOTCH3 expression and causing NOTCH3 signaling impairment. NOTCH3 targets acting in the regulation of arterial tone (KCNA5) or expressed in the vasculature (CDH6) were downregulated. Patient's vessels were characterized by smooth muscle degeneration as in CADASIL, but without deposition of granular osmiophilic material (GOM), the CADASIL hallmark. The heterozygous parents displayed similar but less dramatic trends in decrease in the expression of NOTCH3 and its targets, as well as in vessel degeneration. This study suggests a functional link between NOTCH3 deficiency and pathogenesis of vascular leukoencephalopathies.

Patterns and frequencies of acquired and constitutional uniparental isodisomies in pediatric and adult B-cell precursor acute lymphoblastic leukemia

Single nucleotide polymorphism (SNP) arrays are increasingly being used in clinical routine for genetic analysis of pediatric B-cell precursor acute lymphoblastic leukemias (BCP ALL). Because constitutional DNA is not readily available as a control at the time of diagnosis, it is important to be able to distinguish between acquired and constitutional aberrations in a diagnostic setting. In the present study we focused on uniparental isodisomies (UPIDs). SNP array analyses of 143 pediatric and 38 adult B-cell precursor acute lymphoblastic leukemias and matched remission samples revealed acquired whole chromosome or segmental UPIDs (wUPIDs, sUPIDs) in 32 cases (18%), without any age- or gender-related frequency differences. Acquired sUPIDs were larger than the constitutional ones (mean 35.3 Mb vs. 10.7 Mb; P < 0.0001) and were more often terminally located in the chromosomes (69% vs. 4.5%; P < 0.0001). Chromosomes 3, 5, and 9 were most often involved in acquired wUPIDs, whilst recurrent acquired sUPIDs targeted 6p, 9p, 9q, and 14q. The majority (56%) of sUPID9p was associated with homozygous CDKN2A deletions. In pediatric ALL, all wUPIDs were found in high hyperdiploid (51-67 chromosomes) cases and an extended analysis, also including unmatched diagnostic samples, revealed a higher frequency of wUPID-positivity in higher modal number (56-67 chromosomes) than in lower modal number (51-55 chromosomes) high hyperdiploid cases (34% vs. 11%; P = 0.04), suggesting different underlying mechanisms of formation of these subtypes of high hyperdiploidy. © 2016 Wiley Periodicals, Inc.

Comparative Analysis of Genome Diversity in Bullmastiff Dogs

Management and preservation of genomic diversity in dog breeds is a major objective for maintaining health. The present study was undertaken to characterise genomic diversity in Bullmastiff dogs using both genealogical and molecular analysis. Genealogical analysis of diversity was conducted using a database consisting of 16,378 Bullmastiff pedigrees from year 1980 to 2013. Additionally, a total of 188 Bullmastiff dogs were genotyped using the 170,000 SNP Illumina CanineHD Beadchip. Genealogical parameters revealed a mean inbreeding coefficient of 0.047; 142 total founders (f); an effective number of founders (fe) of 79; an effective number of ancestors (fa) of 62; and an effective population size of the reference population of 41. Genetic diversity and the degree of genome-wide homogeneity within the breed were also investigated using molecular data. Multiple-locus heterozygosity (MLH) was equal to 0.206; runs of homozygosity (ROH) as proportion of the genome, averaged 16.44%; effective population size was 29.1, with an average inbreeding coefficient of 0.035, all estimated using SNP Data. Fine-scale population structure was analysed using NETVIEW, a population analysis pipeline. Visualisation of the high definition network captured relationships among individuals within and between subpopulations. Effects of unequal founder use, and ancestral inbreeding and selection, were evident. While current levels of Bullmastiff heterozygosity, inbreeding and homozygosity are not unusual, a relatively small effective population size indicates that a breeding strategy to reduce the inbreeding rate may be beneficial.

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F ROH) in three local dairy cattle breeds. F ROH values were compared with F estimated from the genomic relationship matrix (F GRM), based on the difference between observed v. expected number of homozygous genotypes (F HOM) and the genomic homozygosity of individual i (F MOL i ). The molecular coancestry coefficient (f MOL ij ) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F HOM and F ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F ROH and other estimated F coefficients were generally lower ranged from 0.45 (F MOL i -F ROH) in Cinisara to 0.17 (F GRM-F ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.

Analysis of Shared Haplotypes amongst Palauans Maps Loci for Psychotic Disorders to 4q28 and 5q23-q31

To localize genetic variation affecting risk for psychotic disorders in the population of Palau, we genotyped DNA samples from 203 Palauan individuals diagnosed with psychotic disorders, broadly defined, and 125 control subjects using a genome-wide single nucleotide polymorphism array. Palau has unique features advantageous for this study: due to its population history, Palauans are substantially interrelated; affected individuals often, but not always, cluster in families; and we have essentially complete ascertainment of affected individuals. To localize risk variants to genomic regions, we evaluated long-shared haplotypes, ≥10 Mb, identifying clusters of affected individuals who share such haplotypes. This extensive sharing, typically identical by descent, was significantly greater in cases than population controls, even after controlling for relatedness. Several regions of the genome exhibited substantial excess of shared haplotypes for affected individuals, including 3p21, 3p12, 4q28, and 5q23-q31. Two of these regions, 4q28 and 5q23-q31, showed significant linkage by traditional LOD score analysis and could harbor variants of more sizeable risk for psychosis or a multiplicity of risk variants. The pattern of haplotype sharing in 4q28 highlights <i>PCDH10</i>, encoding a cadherin-related neuronal receptor, as possibly involved in risk.

TRAIP promotes DNA damage response during genome replication and is mutated in primordial dwarfism

DNA lesions encountered by replicative polymerases threaten genome stability and cell cycle progression. Here we report the identification of mutations in TRAIP, encoding an E3 RING ubiquitin ligase, in patients with microcephalic primordial dwarfism. We establish that TRAIP relocalizes to sites of DNA damage, where it is required for optimal phosphorylation of H2AX and RPA2 during S-phase in response to ultraviolet (UV) irradiation, as well as fork progression through UV-induced DNA lesions. TRAIP is necessary for efficient cell cycle progression and mutations in TRAIP therefore limit cellular proliferation, providing a potential mechanism for microcephaly and dwarfism phenotypes. Human genetics thus identifies TRAIP as a component of the DNA damage response to replication-blocking DNA lesions.

Genetics and psychotic disorders: A fresh look at consanguinity

Consanguineous unions refer to marriages between related individuals who share a common ancestor. These unions are still commonplace in certain regions of the world such as the southern coast of the Mediterranean, throughout the Middle East and South-East Asia. According to available data, couples of second cousins or closer and their offspring currently represent 10.4% of the world's population, thus resulting in increased frequencies of autosomal recessive disorders. Furthermore, consanguinity may be implicated in the increased frequency of multifactorial pathologies such as mental disorders. The few existing epidemiological studies in consanguineous and/or geographically isolated populations confirm that there is a significant association between consanguinity and mental disorders and a higher risk of schizophrenia or bipolar disorders among offspring from consanguineous couples. There exists a strong and complex genetic component in the predisposition to psychotic disorders that has been confirmed in numerous studies. However, the genetic basis of these disorders remains poorly understood. GWAS studies (Genome Wide Association Studies) over the past 10 years have identified a few weak associations, thus refuting the "common diseases-common variants" hypothesis. A model implicating numerous rare variants has been supported by the recent discovery of CNVs (Copy Number Variants) and their statistically significant association with psychiatric disorders such as schizophrenia, bipolar disorders and autism. The study of consanguineous families may contribute to identifying rare variants in homogenous populations who have conserved certain alleles. Major developments in molecular biology techniques would facilitate these studies as well as contributing to identifying major genes. These results emphasize the need for genetic counseling in high-risk communities and the importance of implementing preventive actions and raising awareness concerning the risk of consanguineous unions.

Dynamics of Dark-Fly Genome Under Environmental Selections

Environmental adaptation is one of the most fundamental features of organisms. Modern genome science has identified some genes associated with adaptive traits of organisms, and has provided insights into environmental adaptation and evolution. However, how genes contribute to adaptive traits and how traits are selected under an environment in the course of evolution remain mostly unclear. To approach these issues, we utilize “Dark-fly”, a Drosophila melanogaster line maintained in constant dark conditions for more than 60 years. Our previous analysis identified 220,000 single nucleotide polymorphisms (SNPs) in the Dark-fly genome, but did not clarify which SNPs of Dark-fly are truly adaptive for living in the dark. We found here that Dark-fly dominated over the wild-type fly in a mixed population under dark conditions, and based on this domination we designed an experiment for genome reselection to identify adaptive genes of Dark-fly. For this experiment, large mixed populations of Dark-fly and the wild-type fly were maintained in light conditions or in dark conditions, and the frequencies of Dark-fly SNPs were compared between these populations across the whole genome. We thereby detected condition-dependent selections toward approximately 6% of the genome. In addition, we observed the time-course trajectory of SNP frequency in the mixed populations through generations 0, 22, and 49, which resulted in notable categorization of the selected SNPs into three types with different combinations of positive and negative selections. Our data provided a list of about 100 strong candidate genes associated with the adaptive traits of Dark-fly.

The Molecular Karyotype of 25 Clinical-Grade Human Embryonic Stem Cell Lines

The application of human embryonic stem cell (hESC) derivatives to regenerative medicine is now becoming a reality. Although the vast majority of hESC lines have been derived for research purposes only, about 50 lines have been established under Good Manufacturing Practice (GMP) conditions. Cell types differentiated from these designated lines may be used as a cell therapy to treat macular degeneration, Parkinson’s, Huntington’s, diabetes, osteoarthritis and other degenerative conditions. It is essential to know the genetic stability of the hESC lines before progressing to clinical trials. We evaluated the molecular karyotype of 25 clinical-grade hESC lines by whole-genome single nucleotide polymorphism (SNP) array analysis. A total of 15 unique copy number variations (CNVs) greater than 100 kb were detected, most of which were found to be naturally occurring in the human population and none were associated with culture adaptation. In addition, three copy-neutral loss of heterozygosity (CN-LOH) regions greater than 1 Mb were observed and all were relatively small and interstitial suggesting they did not arise in culture. The large number of available clinical-grade hESC lines with defined molecular karyotypes provides a substantial starting platform from which the development of pre-clinical and clinical trials in regenerative medicine can be realised.

Inbreeding and homozygosity in breast cancer survival

Genome-wide association studies (GWASs) help to understand the effects of single nucleotide polymorphisms (SNPs) on breast cancer (BC) progression and survival. We performed multiple analyses on data from a previously conducted GWAS for the influence of individual SNPs, runs of homozygosity (ROHs) and inbreeding on BC survival. (I.) The association of individual SNPs indicated no differences in the proportions of homozygous individuals among short-time survivors (STSs) and long-time survivors (LTSs). (II.) The analysis revealed differences among the populations for the number of ROHs per person and the total and average length of ROHs per person and among LTSs and STSs for the number of ROHs per person. (III.) Common ROHs at particular genomic positions were nominally more frequent among LTSs than in STSs. Common ROHs showed significant evidence for natural selection (iHS, Tajima's D, Fay-Wu's H). Most regions could be linked to genes related to BC progression or treatment. (IV.) Results were supported by a higher level of inbreeding among LTSs. Our results showed that an increased level of homozygosity may result in a preference of individuals during BC treatment. Although common ROHs were short, variants within ROHs might favor survival of BC and may function in a recessive manner.

Genome-wide association studies identify genetic loci for low von Willebrand factor levels

Low von Willebrand factor (VWF) levels are associated with bleeding symptoms and are a diagnostic criterion for von Willebrand disease, the most common inherited bleeding disorder. To date, it is unclear which genetic loci are associated with reduced VWF levels. Therefore, we conducted a meta-analysis of genome-wide association studies to identify genetic loci associated with low VWF levels. For this meta-analysis, we included 31 149 participants of European ancestry from 11 community-based studies. From all participants, VWF antigen (VWF:Ag) measurements and genome-wide single-nucleotide polymorphism (SNP) scans were available. Each study conducted analyses using logistic regression of SNPs on dichotomized VWF:Ag measures (lowest 5% for blood group O and non-O) with an additive genetic model adjusted for age and sex. An inverse-variance weighted meta-analysis was performed for VWF:Ag levels. A total of 97 SNPs exceeded the genome-wide significance threshold of 5 × 10(-8) and comprised five loci on four different chromosomes: 6q24 (smallest P-value 5.8 × 10(-10)), 9q34 (2.4 × 10(-64)), 12p13 (5.3 × 10(-22)), 12q23 (1.2 × 10(-8)) and 13q13 (2.6 × 10(-8)). All loci were within or close to genes, including STXBP5 (Syntaxin Binding Protein 5) (6q24), STAB5 (stabilin-5) (12q23), ABO (9q34), VWF (12p13) and UFM1 (ubiquitin-fold modifier 1) (13q13). Of these, UFM1 has not been previously associated with VWF:Ag levels. Four genes that were previously associated with VWF levels (VWF, ABO, STXBP5 and STAB2) were also associated with low VWF levels, and, in addition, we identified a new gene, UFM1, that is associated with low VWF levels. These findings point to novel mechanisms for the occurrence of low VWF levels.