Runs of homozygosity in European populations

Genomic inbreeding coefficients based on the distribution of the length of runs of homozygosity in a closed line of Iberian pigs

BACKGROUND

The increasing availability of DNA markers provides new metrics of inbreeding based on single nucleotide polymorphisms (SNPs), i.e. molecular inbreeding or the proportion of runs of homozygosity (ROH), as alternatives to traditional pedigree-based inbreeding coefficients. However, none of these metrics incorporate the length of ROH as an indicator of recent inbreeding. Novel inbreeding coefficients that incorporate length of ROH as a random variable with an associated density are investigated.

METHODS

New inbreeding metrics based on the distribution of the length of ROH are proposed: (1) the Kolmolgorov-Smirnov test, (2) a function of the quantiles of the cumulative distribution function of an individual versus the population, and (3) fitting of an exponential distribution to ROH lengths (mean, variance, and the probability of drawing at random a ROH larger than a given threshold). The new inbreeding and pedigree-based metrics were compared using 217 sows of an Iberian line that belong to three groups: C1 (conservation), C2 (conservation derived from C1), and S (selected and derived from C1), with complete pedigrees and genotyped for 35,023 SNPs.

RESULTS

Correlations between pedigree-based and the new genomic inbreeding coefficients ranged from 0.22 to 0.72 but most ranged from 0.60 to 0.70. The correlation between quantile chromosomal inbreeding coefficients (using molecular information of just one chromosome at the time) and chromosomal length was 0.84 (SE = 0.14), supporting the hypothesis that these coefficients incorporate information on ROH length as an indication of recent inbreeding. Kolmogorov-Smirnov and exponential chromosomal inbreeding coefficients were also correlated with chromosomal length (0.57). Chromosome 1 had the largest quantile ROH inbreeding coefficient (largest ROH sizes), whereas chromosome 10 had the lowest (shortest ROH sizes). Selection for lean growth increased ROH-based inbreeding coefficients for group S when compared to unselected groups C1 and C2. At the chromosomal level, this comparison showed that the level of autozygosity and the length of ROH for most of the autosomes increased in the selection line.

CONCLUSIONS

Quantile and exponential probability inbreeding coefficients using ROH length as a random variable provide additional information about recent inbreeding compared to existing inbreeding coefficients such as molecular, pedigree-based or total ROH content inbreeding coefficients.

Long contiguous stretches of homozygosity spanning shortly the imprinted loci are associated with intellectual disability, autism and/or epilepsy

BACKGROUND

Long contiguous stretches of homozygosity (LCSH) (regions/runs of homozygosity) are repeatedly detected by single-nucleotide polymorphism (SNP) chromosomal microarrays. Providing important clues regarding parental relatedness (consanguinity), uniparental disomy, chromosomal recombination or rearrangements, LCSH are rarely considered as a possible epigenetic cause of neurodevelopmental disorders. Additionally, despite being relevant to imprinting, LCSH at imprinted loci have not been truly addressed in terms of pathogenicity. In this study, we examined LCSH in children with unexplained intellectual disability, autism, congenital malformations and/or epilepsy focusing on chromosomal regions which harbor imprinted disease genes.

RESULTS

Out of 267 cases, 14 (5.2 %) were found to have LCSH at imprinted loci associated with a clinical outcome. There were 5 cases of LCSH at 15p11.2, 4 cases of LCSH at 7q31.2, 3 cases of LCSH at 11p15.5, and 2 cases of LCSH at 7q21.3. Apart from a case of LCSH at 7q31.33q32.3 (~4 Mb in size), all causative LCSH were 1-1.5 Mb in size. Clinically, these cases were characterized by a weak resemblance to corresponding imprinting diseases (i.e., Silver-Russell, Beckwith-Wiedemann, and Prader-Willi/Angelman syndromes), exhibiting distinctive intellectual disability, autistic behavior, developmental delay, seizures and/or facial dysmorphisms. Parental consanguinity was detected in 8 cases (3 %), and these cases did not exhibit LCSH at imprinted loci.

CONCLUSIONS

This study demonstrates that shorter LCSH at chromosomes 7q21.3, 7q31.2, 11p15.5, and 15p11.2 occur with a frequency of about 5 % in the children with intellectual disability, autism, congenital malformations and/or epilepsy. Consequently, this type of epigenetic mutations appears to be the most common one among children with neurodevelopmental diseases. Finally, since LCSH less than 2.5-10 Mb in size are generally ignored in diagnostic SNP microarray studies, one can conclude that an important epigenetic cause of intellectual disability, autism or epilepsy is actually overlooked.

Runs of homozygosity reveal signatures of positive selection for reproduction traits in breed and non-breed horses

Background

Modern horses represent heterogeneous populations specifically selected for appearance and performance. Genomic regions under high selective pressure show characteristic runs of homozygosity (ROH) which represent a low genetic diversity. This study aims at detecting the number and functional distribution of ROHs in different horse populations using next generation sequencing data.

Methods

Next generation sequencing was performed for two Sorraia, one Dülmen Horse, one Arabian, one Saxon-Thuringian Heavy Warmblood, one Thoroughbred and four Hanoverian. After quality control reads were mapped to the reference genome EquCab2.70. ROH detection was performed using PLINK, version 1.07 for a trimmed dataset with 11,325,777 SNPs and a mean read depth of 12. Stretches with homozygous genotypes of >40 kb as well as >400 kb were defined as ROHs. SNPs within consensus ROHs were tested for neutrality. Functional classification was done for genes annotated within ROHs using PANTHER gene list analysis and functional variants were tested for their distribution among breed or non-breed groups.

Results

ROH detection was performed using whole genome sequences of ten horses of six populations representing various breed types and non-breed horses. In total, an average number of 3492 ROHs were detected in windows of a minimum of 50 consecutive homozygous SNPs and an average number of 292 ROHs in windows of 500 consecutive homozygous SNPs. Functional analyses of private ROHs in each horse revealed a high frequency of genes affecting cellular, metabolic, developmental, immune system and reproduction processes. In non-breed horses, 198 ROHs in 50-SNP windows and seven ROHs in 500-SNP windows showed an enrichment of genes involved in reproduction, embryonic development, energy metabolism, muscle and cardiac development whereas all seven breed horses revealed only three common ROHs in 50-SNP windows harboring the fertility-related gene YES1. In the Hanoverian, a total of 18 private ROHs could be shown to be located in the region of genes potentially involved in neurologic control, signaling, glycogen balance and reproduction. Comparative analysis of homozygous stretches common in all ten horses displayed three ROHs which were all located in the region of KITLG, the ligand of KIT known to be involved in melanogenesis, haematopoiesis and gametogenesis.

Conclusions

The results of this study give a comprehensive insight into the frequency and number of ROHs in various horses and their potential influence on population diversity and selection pressures. Comparisons of breed and non-breed horses suggest a significant artificial as well as natural selection pressure on reproduction performance in all types of horse populations.

Electronic supplementary material

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

Estimation of Recent and Ancient Inbreeding in a Small Endogamous Tunisian Community Through Genomic Runs of Homozygosity

Runs of homozygosity (ROHs) are extended genomic regions of homozygous genotypes that record populations' mating patterns in the past. We performed microarray genotyping on 15 individuals from a small isolated Tunisian community. We estimated the individual and population genome-wide level of homozygosity from data on ROH above 0.5 Mb in length. We found a high average number of ROH per individual (48.2). The smallest ROH category (0.5-1.49 Mb) represents 0.93% of the whole genome, while medium-size (1.5-4.99 Mb) and long-size ROH (≥5 Mb) cover 1.18% and 0.95%, respectively. We found that genealogical individual inbreeding coefficients (Fped ) based on three- to four-generation pedigrees are not reliable indicators of the current proportion of genome-wide homozygosity inferred from ROH (FROH ) either for 0.5 or 1.5 Mb ROH length thresholds, while identity-by-descent sharing is a function of shared coancestry. This study emphasizes the effect of reproductive isolation and a prolonged practice of consanguinity that limits the genetic heterogeneity. It also provides evidence of both recent and ancient parental relatedness contribution to the current level of genome-wide homozygosity in the studied population. These findings may be useful for evaluation of long-term effects of inbreeding on human health and for future applications of ROHs in identifying recessive susceptibility genes.

Somatic mosaicism for copy-neutral loss of heterozygosity and DNA copy number variations in the human genome

BACKGROUND

Somatic mosaicism denotes the presence of genetically distinct populations of somatic cells in one individual who has developed from a single fertilised oocyte. Mosaicism may result from a mutation that occurs during postzygotic development and is propagated to only a subset of the adult cells. Our aim was to investigate both somatic mosaicism for copy-neutral loss of heterozygosity (cn-LOH) events and DNA copy number variations (CNVs) in fully differentiated tissues.

RESULTS

We studied panels of tissue samples (11-12 tissues per individual) from four autopsy subjects using high-resolution Illumina HumanOmniExpress-12 BeadChips to reveal the presence of possible intra-individual tissue-specific cn-LOH and CNV patterns. We detected five mosaic cn-LOH regions >5 Mb in some tissue samples in three out of four individuals. We also detected three CNVs that affected only a portion of the tissues studied in one out of four individuals. These three somatic CNVs range from 123 to 796 kb and are also found in the general population. An attempt was made to explain the succession of genomic events that led to the observed somatic genetic mosaicism under the assumption that the specific mosaic patterns of CNV and cn-LOH changes reflect their formation during the postzygotic embryonic development of germinal layers and organ systems.

CONCLUSIONS

Our results give further support to the idea that somatic mosaicism for CNVs, and also cn-LOHs, is a common phenomenon in phenotypically normal humans. Thus, the examination of only a single tissue might not provide enough information to diagnose potentially deleterious CNVs within an individual. During routine CNV and cn-LOH analysis, DNA derived from a buccal swab can be used in addition to blood DNA to get information about the CNV/cn-LOH content in tissues of both mesodermal and ectodermal origin. Currently, the real frequency and possible phenotypic consequences of both CNVs and cn-LOHs that display somatic mosaicism remain largely unknown. To answer these questions, future studies should involve larger cohorts of individuals and a range of tissues.

Genetic analysis of the isolated Faroe Islands reveals SORCS3 as a potential multiple sclerosis risk gene

Background:

In search of the missing heritability in multiple sclerosis (MS), additional approaches adding to the genetic discoveries of large genome-wide association studies are warranted.

Objective:

The objective of this research paper is to search for rare genetic MS risk variants in the genetically homogenous population of the isolated Faroe Islands.

Methods:

Twenty-nine Faroese MS cases and 28 controls were genotyped with the HumanOmniExpressExome-chip. The individuals make up 1596 pair-combinations in which we searched for identical-by-descent shared segments using the PLINK-program.

Results:

A segment spanning 63 SNPs with excess case-case-pair sharing was identified (0.00173 < p > 0.00212). A haplotype consisting of 42 of the 63 identified SNPs which spanned the entire the Sortilin-related vacuolar protein sorting 10 domain containing receptor 3 ( SORCS3) gene had a carrier frequency of 0.34 in cases but was not present in any controls ( p = 0.0008).

Conclusion:

This study revealed an oversharing in case-case-pairs of a segment spanning 63 SNPs and the entire SORCS3. While not previously associated with MS, SORCS3 appears to be important in neuronal plasticity through its binding of neurotrophin factors and involvement in glutamate homeostasis. Although additional work is needed to scrutinise the genetic effect of the SORCS3-covering haplotype, this study suggests that SORCS3 may also be important in MS pathogenesis.

The Diagnostic Utility of Single Long Contiguous Stretches of Homozygosity in Patients without Parental Consanguinity

We present data from our clinical department's experience with chromosomal microarray analysis (CMA) regarding the diagnostic utility of 1 or 2 long contiguous stretches of homozygosity (LCSHs) in an outbred population. The study group consisted of 2,110 consecutive patients from 2011 to 2014 for whom CMA was performed. The minimum cut-off size for defining a homozygous stretch was 5 Mb. To focus on cases with no parental consanguinity, we further studied only patients in whom the total length of homozygous stretches did not exceed 28 Mb or 1% of the autosomal genome length. We identified 6 chromosomal regions where homozygous stretches appeared in at least 3 patients and excluded these from further analysis. In 2 out of 120 patients with an isolated finding of 1 or 2 non-recurrent LCSHs, a plausible candidate gene associated with their phenotype was identified within the homozygous stretch. In both of these cases, a pathogenic mutation was detected, leading to diagnoses of pyruvate kinase deficiency and Marinesco-Sjögren syndrome. To clarify whether previously found homozygous stretches could be important for the interpretation of genome-wide sequencing data, we report 7 cases in which homozygous stretches not encompassing a clinically associated gene were first found on CMA, followed by the diagnostic whole-exome sequencing. The diagnostic utility of single LCSHs, unlikely to be caused by uniparental disomy, is discussed in detail.

Runs of Homozygosity: Association with Coronary Artery Disease and Gene Expression in Monocytes and Macrophages

Runs of homozygosity (ROHs) are recognized signature of recessive inheritance. Contributions of ROHs to the genetic architecture of coronary artery disease and regulation of gene expression in cells relevant to atherosclerosis are not known. Our combined analysis of 24,320 individuals from 11 populations of white European ethnicity showed an association between coronary artery disease and both the count and the size of ROHs. Individuals with coronary artery disease had approximately 0.63 (95% CI: 0.4-0.8) excess of ROHs when compared to coronary-artery-disease-free control subjects (p = 1.49 × 10(-9)). The average total length of ROHs was approximately 1,046.92 (95% CI: 634.4-1,459.5) kb greater in individuals with coronary artery disease than control subjects (p = 6.61 × 10(-7)). None of the identified individual ROHs was associated with coronary artery disease after correction for multiple testing. However, in aggregate burden analysis, ROHs favoring increased risk of coronary artery disease were much more common than those showing the opposite direction of association with coronary artery disease (p = 2.69 × 10(-33)). Individual ROHs showed significant associations with monocyte and macrophage expression of genes in their close proximity-subjects with several individual ROHs showed significant differences in the expression of 44 mRNAs in monocytes and 17 mRNAs in macrophages when compared to subjects without those ROHs. This study provides evidence for an excess of homozygosity in coronary artery disease in outbred populations and suggest the potential biological relevance of ROHs in cells of importance to the pathogenesis of atherosclerosis.

Multicohort analysis of the maternal age effect on recombination

Several studies have reported that the number of crossovers increases with maternal age in humans, but others have found the opposite. Resolving the true effect has implications for understanding the maternal age effect on aneuploidies. Here, we revisit this question in the largest sample to date using single nucleotide polymorphism (SNP)-chip data, comprising over 6,000 meioses from nine cohorts. We develop and fit a hierarchical model to allow for differences between cohorts and between mothers. We estimate that over 10 years, the expected number of maternal crossovers increases by 2.1% (95% credible interval (0.98%, 3.3%)). Our results are not consistent with the larger positive and negative effects previously reported in smaller cohorts. We see heterogeneity between cohorts that is likely due to chance effects in smaller samples, or possibly to confounders, emphasizing that care should be taken when interpreting results from any specific cohort about the effect of maternal age on recombination.

The use and abuse of genetic marker-based estimates of relatedness and inbreeding

Genetic marker-based estimators remain a popular tool for measuring relatedness (r xy ) and inbreeding (F) coefficients at both the population and individual level. The performance of these estimators fluctuates with the number and variability of markers available, and the relatedness composition and demographic history of a population. Several methods are available to evaluate the reliability of the estimates of r xy and F, some of which are implemented in the program COANCESTRY. I used the simulation module in COANCESTRY since assess the performance of marker-based estimators of r xy and F in a species with very low genetic diversity, New Zealand's little spotted kiwi (Apteryx owenii). I also conducted a review of published papers that have used COANCESTRY as its release to assess whether and how the reliability of the estimates of r xy and F produced by genetic markers are being measured and reported in published studies. My simulation results show that even when the correlation between true (simulated) and estimated r xy or F is relatively high (Pearson's r = 0.66-0.72 and 0.81-0.85, respectively) the imprecision of the estimates renders them highly unreliable on an individual basis. The literature review demonstrates that the majority of studies do not report the reliability of marker-based estimates of r xy and F. There is currently no standard practice for selecting the best estimator for a given data set or reporting an estimator's performance. This could lead to experimental results being interpreted out of context and render the robustness of conclusions based on measures of r xy and F debatable.

Directional dominance on stature and cognition in diverse human populations

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10(-300), 2.1 × 10(-6), 2.5 × 10(-10) and 1.8 × 10(-10), respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months' less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

Estimation of inbreeding using pedigree, 50k SNP chip genotypes and full sequence data in three cattle breeds

Background

Levels of inbreeding in cattle populations have increased in the past due to the use of a limited number of bulls for artificial insemination. High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. Various estimators based on different sources, e.g., pedigree or genomic data, have been used to estimate inbreeding coefficients in cattle populations. However, the comparative advantage of using full sequence data to assess inbreeding is unknown. We used pedigree and genomic data at different densities from 50k to full sequence variants to compare how different methods performed for the estimation of inbreeding levels in three different cattle breeds.

Results

Five different estimates for inbreeding were calculated and compared in this study: pedigree based inbreeding coefficient (F_PED); run of homozygosity (ROH)-based inbreeding coefficients (F_ROH); genomic relationship matrix (GRM)-based inbreeding coefficients (F_GRM); inbreeding coefficients based on excess of homozygosity (F_HOM) and correlation of uniting gametes (F_UNI). Estimates using ROH provided the direct estimated levels of autozygosity in the current populations and are free effects of allele frequencies and incomplete pedigrees which may increase in inaccuracy in estimation of inbreeding. The highest correlations were observed between F_ROH estimated from the full sequence variants and the F_ROH estimated from 50k SNP (single nucleotide polymorphism) genotypes. The estimator based on the correlation between uniting gametes (F_UNI) using full genome sequences was also strongly correlated with F_ROH detected from sequence data.

Conclusions

Estimates based on ROH directly reflected levels of homozygosity and were not influenced by allele frequencies, unlike the three other estimates evaluated (F_GRM, F_HOM and F_UNI), which depended on estimated allele frequencies. F_PED suffered from limited pedigree depth. Marker density affects ROH estimation. Detecting ROH based on 50k chip data was observed to give estimates similar to ROH from sequence data. In the absence of full sequence data ROH based on 50k can be used to access homozygosity levels in individuals. However, genotypes denser than 50k are required to accurately detect short ROH that are most likely identical by descent (IBD).

Determining the genome-wide kinship coefficient seems unhelpful in distinguishing consanguineous couples with a high versus low risk for adverse reproductive outcome

Background

Offspring of consanguineous couples are at increased risk of congenital disorders. The risk increases as parents are more closely related. Individuals that have the same degree of relatedness according to their pedigree, show variable genomic kinship coefficients. To investigate whether we can differentiate between couples with high- and low risk for offspring with congenital disorders, we have compared the genomic kinship coefficient of consanguineous parents with a child affected with an autosomal recessive disorder with that of consanguineous parents with only healthy children, corrected for the degree of pedigree relatedness.

Methods

151 consanguineous couples (73 cases and 78 controls) from 10 different ethnic backgrounds were genotyped on the Affymetrix platform and passed quality control checks. After pruning SNPs in linkage disequilibrium, 57,358 SNPs remained. Kinship coefficients were calculated using three different toolsets: PLINK, King and IBDelphi, yielding five different estimates (IBDelphi, PLINK (all), PLINK (by population), King robust (all) and King homo (by population)). We performed a one-sided Mann Whitney test to investigate whether the median relative difference regarding observed and expected kinship coefficients is bigger for cases than for controls. Furthermore, we fitted a mixed effects linear model to correct for a possible population effect.

Results

Although the estimated degrees of genomic relatedness with the different toolsets show substantial variability, correlation measures between the different estimators demonstrated moderate to strong correlations. Controls have higher point estimates for genomic kinship coefficients. The one-sided Mann Whitney test did not show any evidence for a higher median relative difference for cases compared to controls. Neither did the regression analysis exhibit a positive association between case–control status and genomic kinship coefficient.

Conclusions

In this case–control setting, in which we compared consanguineous couples corrected for degree of pedigree relatedness, a higher degree of genomic relatedness was not significantly associated with a higher likelihood of having an affected child. Further translational research should focus on which parts of the genome and which pathogenic mutations couples are sharing. Looking at relatedness coefficients by determining genome-wide SNPs does not seem to be an effective measure for prospective risk assessment in consanguineous parents.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-015-0191-0) contains supplementary material, which is available to authorized users.

Homozygous loss-of-function variants in European cosmopolitan and isolate populations

Homozygous loss of function (HLOF) variants provide a valuable window on gene function in humans, as well as an inventory of the human genes that are not essential for survival and reproduction. All humans carry at least a few HLOF variants, but the exact number of inactivated genes that can be tolerated is currently unknown—as are the phenotypic effects of losing function for most human genes. Here, we make use of 1432 whole exome sequences from five European populations to expand the catalogue of known human HLOF mutations; after stringent filtering of variants in our dataset, we identify a total of 173 HLOF mutations, 76 (44%) of which have not been observed previously. We find that population isolates are particularly well suited to surveys of novel HLOF genes because individuals in such populations carry extensive runs of homozygosity, which we show are enriched for novel, rare HLOF variants. Further, we make use of extensive phenotypic data to show that most HLOFs, ascertained in population-based samples, appear to have little detectable effect on the phenotype. On the contrary, we document several genes directly implicated in disease that seem to tolerate HLOF variants. Overall HLOF genes are enriched for olfactory receptor function and are expressed in testes more often than expected, consistent with reduced purifying selection and incipient pseudogenisation.

The Relationship between Runs of Homozygosity and Inbreeding in Jersey Cattle under Selection

Inbreeding is often an inevitable outcome of strong directional artificial selection but on average it reduces population fitness with increased frequency of recessive deleterious alleles. Runs of homozygosity (ROH) representing genomic autozygosity that occur from mating between selected and genomically related individuals may be able to reveal the regions affecting fitness. To examine the influence of genomic autozygosity on fitness, we used a genome-wide association test to evaluate potential negative correlations between ROH and daughter pregnancy rate (DPR) or somatic cell score (SCS) in US Jersey cattle. In addition, relationships between changes of local ROH and inbreeding coefficients (F) were assessed to locate genomic regions with increased inbreeding. Despite finding some decreases in fertility associated with incremental increases in F, most emerging local ROH were not significantly associated with DPR or SCS. Furthermore, the analyses of ROH could be approximated with the most frequent haplotype(s), including the associations of ROH and F or traits. The analysis of the most frequent haplotype revealed that associations of ROH and fertility could be accounted for by the additive genetic effect on the trait. Thus, we suggest that a change of autozygosity is more likely to demonstrate footprints of selected haplotypes for production rather than highlight the possible increased local autozygosity of a recessive detrimental allele resulting from the mating between closely related animals in Jersey cattle.

Factors influencing success of clinical genome sequencing across a broad spectrum of disorders

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.

Recent genomic heritage in Scotland

BACKGROUND

The Generation Scotland Scottish Family Health Study (GS:SFHS) includes 23,960 participants from across Scotland with records for many health-related traits and environmental covariates. Genotypes at ~700 K SNPs are currently available for 10,000 participants. The cohort was designed as a resource for genetic and health related research and the study of complex traits. In this study we developed a suite of analyses to disentangle the genomic differentiation within GS:SFHS individuals to describe and optimise the sample and methods for future analyses.

RESULTS

We combined the genotypic information of GS:SFHS with 1092 individuals from the 1000 Genomes project and estimated their genomic relationships. Then, we performed Principal Component Analyses of the resulting relationships to investigate the genomic origin of different groups. We characterised two groups of individuals: those with a few sparse rare markers in the genome, and those with several large rare haplotypes which might represent relatively recent exogenous ancestors. We identified some individuals with likely Italian ancestry and a group with some potential African/Asian ancestry. An analysis of homozygosity in the GS:SFHS sample revealed a very similar pattern to other European populations. We also identified an individual carrying a chromosome 1 uniparental disomy. We found evidence of local geographic stratification within the population having impact on the genomic structure.

CONCLUSIONS

These findings illuminate the history of the Scottish population and have implications for further analyses such as the study of the contributions of common and rare variants to trait heritabilities and the evaluation of genomic and phenotypic prediction of disease.

Genomic analysis for managing small and endangered populations: a case study in Tyrol Grey cattle

Analysis of genomic data is increasingly becoming part of the livestock industry. Therefore, the routine collection of genomic information would be an invaluable resource for effective management of breeding programs in small, endangered populations. The objective of the paper was to demonstrate how genomic data could be used to analyse (1) linkage disequlibrium (LD), LD decay and the effective population size (Ne_LD); (2) Inbreeding level and effective population size (Ne_ROH) based on runs of homozygosity (ROH); (3) Prediction of genomic breeding values (GEBV) using small within-breed and genomic information from other breeds. The Tyrol Grey population was used as an example, with the goal to highlight the potential of genomic analyses for small breeds. In addition to our own results we discuss additional use of genomics to assess relatedness, admixture proportions, and inheritance of harmful variants. The example data set consisted of 218 Tyrol Grey bull genotypes, which were all available AI bulls in the population. After standard quality control restrictions 34,581 SNPs remained for the analysis. A separate quality control was applied to determine ROH levels based on Illumina GenCall and Illumina GenTrain scores, resulting into 211 bulls and 33,604 SNPs. LD was computed as the squared correlation coefficient between SNPs within a 10 mega base pair (Mb) region. ROHs were derived based on regions covering at least 4, 8, and 16 Mb, suggesting that animals had common ancestors approximately 12, 6, and 3 generations ago, respectively. The corresponding mean inbreeding coefficients (F_ROH) were 4.0% for 4 Mb, 2.9% for 8 Mb and 1.6% for 16 Mb runs. With an average generation interval of 5.66 years, estimated Ne_ROH was 125 (Ne_{ROH>16 Mb}), 186 (Ne_{ROH>8 Mb}) and 370 (Ne_{ROH>4 Mb}) indicating strict avoidance of close inbreeding in the population. The LD was used as an alternative method to infer the population history and the Ne. The results show a continuous decrease in Ne_LD, to 780, 120, and 80 for 100, 10, and 5 generations ago, respectively. Genomic selection was developed for and is working well in large breeds. The same methodology was applied in Tyrol Grey cattle, using different reference populations. Contrary to the expectations, the accuracy of GEBVs with very small within breed reference populations were very high, between 0.13–0.91 and 0.12–0.63, when estimated breeding values and deregressed breeding values were used as pseudo-phenotypes, respectively. Subsequent analyses confirmed the high accuracies being a consequence of low reliabilities of pseudo-phenotypes in the validation set, thus being heavily influenced by parent averages. Multi-breed and across breed reference sets gave inconsistent and lower accuracies. Genomic information may have a crucial role in management of small breeds, even if its primary usage differs from that of large breeds. It allows to assess relatedness between individuals, trends in inbreeding and to take decisions accordingly. These decisions would be based on the real genome architecture, rather than conventional pedigree information, which can be missing or incomplete. We strongly suggest the routine genotyping of all individuals that belong to a small breed in order to facilitate the effective management of endangered livestock populations.

Genomic prediction of complex human traits: relatedness, trait architecture and predictive meta-models

We explore the prediction of individuals' phenotypes for complex traits using genomic data. We compare several widely used prediction models, including Ridge Regression, LASSO and Elastic Nets estimated from cohort data, and polygenic risk scores constructed using published summary statistics from genome-wide association meta-analyses (GWAMA). We evaluate the interplay between relatedness, trait architecture and optimal marker density, by predicting height, body mass index (BMI) and high-density lipoprotein level (HDL) in two data cohorts, originating from Croatia and Scotland. We empirically demonstrate that dense models are better when all genetic effects are small (height and BMI) and target individuals are related to the training samples, while sparse models predict better in unrelated individuals and when some effects have moderate size (HDL). For HDL sparse models achieved good across-cohort prediction, performing similarly to the GWAMA risk score and to models trained within the same cohort, which indicates that, for predicting traits with moderately sized effects, large sample sizes and familial structure become less important, though still potentially useful. Finally, we propose a novel ensemble of whole-genome predictors with GWAMA risk scores and demonstrate that the resulting meta-model achieves higher prediction accuracy than either model on its own. We conclude that although current genomic predictors are not accurate enough for diagnostic purposes, performance can be improved without requiring access to large-scale individual-level data. Our methodologically simple meta-model is a means of performing predictive meta-analysis for optimizing genomic predictions and can be easily extended to incorporate multiple population-level summary statistics or other domain knowledge.

Recent artificial selection in U.S. Jersey cattle impacts autozygosity levels of specific genomic regions

BACKGROUND

Genome signatures of artificial selection in U.S. Jersey cattle were identified by examining changes in haplotype homozygosity for a resource population of animals born between 1953 and 2007. Genetic merit of this population changed dramatically during this period for a number of traits, especially milk yield. The intense selection underlying these changes was achieved through extensive use of artificial insemination (AI), which also increased consanguinity of the population to a few superior Jersey bulls. As a result, allele frequencies are shifted for many contemporary animals, and in numerous cases to a homozygous state for specific genomic regions. The goal of this study was to identify those selection signatures that occurred after extensive use of AI since the 1960, using analyses of shared haplotype segments or Runs of Homozygosity. When combined with animal birth year information, signatures of selection associated with economically important traits were identified and compared to results from an extended haplotype homozygosity analysis.

RESULTS

Overall, our results reveal that more recent selection increased autozygosity across the entire genome, but some specific regions increased more than others. A genome-wide scan identified more than 15 regions with a substantial change in autozygosity. Haplotypes found to be associated with increased milk, fat and protein yield in U.S. Jersey cattle also consistently increased in frequency.

CONCLUSIONS

The analyses used in this study was able to detect directional selection over the last few decades when individual production records for Jersey animals were available.

Surnames and Y-chromosomal markers reveal low relationships in Southern Spain

A sample of 416 males from western and eastern Andalusia has been jointly analyzed for surnames and Y-chromosome haplogroups and haplotypes. The observed number of different surnames was 222 (353 when the second surname of the Spanish system of naming is considered). The great majority of recorded surnames have a Castilian-Leonese origin, while Catalan or Basque surnames have not been found. A few Arab-related surnames appear but none discernible of Sephardic-Jewish descent. Low correlation among surnames with different population frequencies and Y-chromosome markers, at different levels of genetic resolution, has been observed in Andalusia. This finding could be explained mainly by the very low rate of monophyletic surnames because of the historical process of surname ascription and the resulting high frequencies of the most common Spanish surnames. The introduction of surnames in Spain during the Middle Ages coincided with Reconquest of the territories under Islamic rule, and Muslims and Jews progressively adopted the present male line surname system. Sampled surnames and Y-chromosome lineages fit well a power-law distribution and observed isonymy is very close to that of the general population. Besides, our data and results show that the reliability of the isonymy method should be questioned because of the high rate of polyphyletic surnames, even in small geographic regions and autochthonous populations. Random isonymy would be consistently dependent of the most common surname frequencies in the population.

Effective population size, extended linkage disequilibrium and signatures of selection in the rare dog breed lundehund

The Lundehund is an old dog breed with remarkable anatomical features including polydactyly in all four limbs and extraordinary flexibility of the spine. We genotyped 28 Lundehund using the canine Illumina high density beadchip to estimate the effective population size (N_e) and inbreeding coefficients as well as to identify potential regions of positive selection. The decay of linkage disequilibrium was slow with r² = 0.95 in 50 kb distance. The last 7-200 generations ago, Ne was at 10-13. An increase of N_e was noted in the very recent generations with a peak value of 19 for N_e at generation 4. The FROH estimated for 50-, 65- and 358-SNP windows were 0.87, 087 and 0.81, respectively. The most likely estimates for F_ROH after removing identical-by-state segments due to linkage disequilibria were at 0.80-0.81. The extreme loss of heterozygosity has been accumulated through continued inbreeding over 200 generations within a probably closed population with a small effective population size. The mean inbreeding coefficient based on pedigree data for the last 11 generations (F_Ped = 0.10) was strongly biased downwards due to the unknown coancestry of the founders in this pedigree data. The long-range haplotype test identified regions with genes involved in processes of immunity, olfaction, woundhealing and neuronal development as potential targets of selection. The genes QSOX2, BMPR1B and PRRX2 as well as MYOM1 are candidates for selection on the Lundehund characteristics small body size, increased number of digits per paw and extraordinary mobility, respectively.

Evolutionary behavioral genetics

We describe the scientific enterprise at the intersection of evolutionary psychology and behavioral genetics-a field that could be termed Evolutionary Behavioral Genetics-and how modern genetic data is revolutionizing our ability to test questions in this field. We first explain how genetically informative data and designs can be used to investigate questions about the evolution of human behavior, and describe some of the findings arising from these approaches. Second, we explain how evolutionary theory can be applied to the investigation of behavioral genetic variation. We give examples of how new data and methods provide insight into the genetic architecture of behavioral variation and what this tells us about the evolutionary processes that acted on the underlying causal genetic variants.

Infectious disease. Life-threatening influenza and impaired interferon amplification in human IRF7 deficiency

Severe influenza disease strikes otherwise healthy children and remains unexplained. We report compound heterozygous null mutations in IRF7, which encodes the transcription factor interferon regulatory factor 7, in an otherwise healthy child who suffered life-threatening influenza during primary infection. In response to influenza virus, the patient's leukocytes and plasmacytoid dendritic cells produced very little type I and III interferons (IFNs). Moreover, the patient's dermal fibroblasts and induced pluripotent stem cell (iPSC)-derived pulmonary epithelial cells produced reduced amounts of type I IFN and displayed increased influenza virus replication. These findings suggest that IRF7-dependent amplification of type I and III IFNs is required for protection against primary infection by influenza virus in humans. They also show that severe influenza may result from single-gene inborn errors of immunity.

Genome-wide scan of gastrointestinal nematode resistance in closed Angus population selected for minimized influence of MHC

Genetic markers associated with parasite indicator traits are ideal targets for study of marker assisted selection aimed at controlling infections that reduce herd use of anthelminthics. For this study, we collected gastrointestinal (GI) nematode fecal egg count (FEC) data from post-weaning animals of an Angus resource population challenged to a 26 week natural exposure on pasture. In all, data from 487 animals was collected over a 16 year period between 1992 and 2007, most of which were selected for a specific DRB1 allele to reduce the influence of potential allelic variant effects of the MHC locus. A genome-wide association study (GWAS) based on BovineSNP50 genotypes revealed six genomic regions located on bovine Chromosomes 3, 5, 8, 15 and 27; which were significantly associated (-log₁₀ p=4.3) with Box-Cox transformed mean FEC (BC-MFEC). DAVID analysis of the genes within the significant genomic regions suggested a correlation between our results and annotation for genes involved in inflammatory response to infection. Furthermore, ROH and selection signature analyses provided strong evidence that the genomic regions associated BC-MFEC have not been affected by local autozygosity or recent experimental selection. These findings provide useful information for parasite resistance prediction for young grazing cattle and suggest new candidate gene targets for development of disease-modifying therapies or future studies of host response to GI parasite infection.

Measuring individual inbreeding in the age of genomics: marker-based measures are better than pedigrees

Inbreeding (mating between relatives) can dramatically reduce the fitness of offspring by causing parts of the genome to be identical by descent. Thus, measuring individual inbreeding is crucial for ecology, evolution and conservation biology. We used computer simulations to test whether the realized proportion of the genome that is identical by descent (IBDG) is predicted better by the pedigree inbreeding coefficient (FP) or by genomic (marker-based) measures of inbreeding. Genomic estimators of IBDG included the increase in individual homozygosity relative to mean Hardy-Weinberg expected homozygosity (FH), and two measures (FROH and FE) that use mapped genetic markers to estimate IBDG. IBDG was more strongly correlated with FH, FE and FROH than with FP across a broad range of simulated scenarios when thousands of SNPs were used. For example, IBDG was more strongly correlated with FROH, FH and FE (estimated with ⩾10 000 SNPs) than with FP (estimated with 20 generations of complete pedigree) in populations with a recent reduction in the effective populations size (from Ne=500 to Ne=75). FROH, FH and FE generally explained >90% of the variance in IBDG (among individuals) when 35 K or more SNPs were used. FP explained <80% of the variation in IBDG on average in all simulated scenarios, even when pedigrees included 20 generations. Our results demonstrate that IBDG can be more precisely estimated with large numbers of genetic markers than with pedigrees. We encourage researchers to adopt genomic marker-based measures of IBDG as thousands of loci can now be genotyped in any species.

Prader-Willi syndrome with a long-contiguous stretch of homozygosity not covering the critical region

Prader-Willi syndrome is a common and complex disorder affecting multiple systems. Its main manifestations are infantile hypotonia with a poor sucking reflex, a characteristic facial appearance, mild mental retardation, hypogonadism and early-onset obesity. Prader-Willi syndrome is due to the absence of paternally expressed imprinted genes at 15q11.2-13, and 3 main mechanisms are known to be involved in its pathogenesis: paternal microdeletions, maternal uniparental disomy events, and imprinting defects. DNA methylation analysis can detect almost all individuals with Prader-Willi syndrome but is unable to distinguish between the molecular classes of the disease. Thus, additional methods are necessary to identify the molecular classes. Here, we employed chromosomal microarray analysis-single nucleotide polymorphism for diagnosis and detected a long-contiguous stretch of homozygosity on chromosome 15, which is highly predictive of maternal uniparental disomy on chromosome 15. Other methods, including fluorescence in situ hybridization, chromosomal microarray analysis-comparative genomic hybridization, genotyping and family linkage analysis, were performed for further validation. In conclusion, our study highlights the use of long-contiguous stretch of homozygosity detection for the diagnosis of Prader-Willi syndrome.

Assessment of autozygosity in Nellore cows (Bos indicus) through high-density SNP genotypes

The use of relatively low numbers of sires in cattle breeding programs, particularly on those for carcass and weight traits in Nellore beef cattle (Bos indicus) in Brazil, has always raised concerns about inbreeding, which affects conservation of genetic resources and sustainability of this breed. Here, we investigated the distribution of autozygosity levels based on runs of homozygosity (ROH) in a sample of 1,278 Nellore cows, genotyped for over 777,000 SNPs. We found ROH segments larger than 10 Mb in over 70% of the samples, representing signatures most likely related to the recent massive use of few sires. However, the average genome coverage by ROH (>1 Mb) was lower than previously reported for other cattle breeds (4.58%). In spite of 99.98% of the SNPs being included within a ROH in at least one individual, only 19.37% of the markers were encompassed by common ROH, suggesting that the ongoing selection for weight, carcass and reproductive traits in this population is too recent to have produced selection signatures in the form of ROH. Three short-range highly prevalent ROH autosomal hotspots (occurring in over 50% of the samples) were observed, indicating candidate regions most likely under selection since before the foundation of Brazilian Nellore cattle. The putative signatures of selection on chromosomes 4, 7, and 12 may be involved in resistance to infectious diseases and fertility, and should be subject of future investigation.

Extensive genome-wide autozygosity in the population isolates of Daghestan

Isolated populations are valuable resources for mapping disease genes, as inbreeding increases genome-wide homozygosity and enhances the ability to map disease alleles on a genetically uniform background within a relatively homogenous environment. The populations of Daghestan are thought to have resided in the Caucasus Mountains for hundreds of generations and are characterized by a high prevalence of certain complex diseases. To explore the extent to which their unique population history led to increased levels of inbreeding, we genotyped >550 000 autosomal single-nucleotide polymorphisms (SNPs) in a set of 14 population isolates speaking Nakh-Daghestanian (ND) languages. The ND-speaking populations showed greatly elevated coefficients of inbreeding, very high numbers and long lengths of Runs of Homozygosity, and elevated linkage disequilibrium compared with surrounding groups from the Caucasus, the Near East, Europe, Central and South Asia. These results are consistent with the hypothesis that most ND-speaking groups descend from a common ancestral population that fragmented into a series of genetic isolates in the Daghestanian highlands. They have subsequently maintained a long-term small effective population size as a result of constant inbreeding and very low levels of gene flow. Given these findings, Daghestanian population isolates are likely to be useful for mapping genes associated with complex diseases.

Detecting inbreeding depression for reproductive traits in Iberian pigs using genome-wide data

BACKGROUND

The current availability of genotypes for very large numbers of single nucleotide polymorphisms (SNPs) is leading to more accurate estimates of inbreeding coefficients and more detailed approaches for detecting inbreeding depression. In the present study, genome-wide information was used to detect inbreeding depression for two reproductive traits (total number of piglets born and number of piglets born alive) in an ancient strain of Iberian pigs (the Guadyerbas strain) that is currently under serious danger of extinction.

METHODS

A total of 109 sows with phenotypic records were genotyped with the PorcineSNP60 BeadChip v1. Inbreeding depression was estimated using a bivariate animal model in which the inbreeding coefficient was included as a covariate. We used two different measures of genomic inbreeding to perform the analyses: inbreeding estimated on a SNP-by-SNP basis and inbreeding estimated from runs of homozygosity. We also performed the analyses using pedigree-based inbreeding.

RESULTS

Significant inbreeding depression was detected for both traits using all three measures of inbreeding. Genome-wide information allowed us to identify one region on chromosome 13 associated with inbreeding depression. This region spans from 27 to 54 Mb and overlaps with a previously detected quantitative trait locus and includes the inter-alpha-trypsin inhibitor gene cluster that is involved with embryo implantation.

CONCLUSIONS

Our results highlight the value of high-density SNP genotyping for providing new insights on where genes causing inbreeding depression are located in the genome. Genomic measures of inbreeding obtained on a SNP-by-SNP basis or those based on the presence/absence of runs of homozygosity represent a suitable alternative to pedigree-based measures to detect inbreeding depression, and a useful tool for mapping studies. To our knowledge, this is the first study in domesticated animals using the SNP-by-SNP inbreeding coefficient to map specific regions within chromosomes associated with inbreeding depression.

Quantifying realized inbreeding in wild and captive animal populations

Most molecular measures of inbreeding do not measure inbreeding at the scale that is most relevant for understanding inbreeding depression-namely the proportion of the genome that is identical-by-descent (IBD). The inbreeding coefficient FPed obtained from pedigrees is a valuable estimator of IBD, but pedigrees are not always available, and cannot capture inbreeding loops that reach back in time further than the pedigree. We here propose a molecular approach to quantify the realized proportion of the genome that is IBD (propIBD), and we apply this method to a wild and a captive population of zebra finches (Taeniopygia guttata). In each of 948 wild and 1057 captive individuals we analyzed available single-nucleotide polymorphism (SNP) data (260 SNPs) spread over four different genomic regions in each population. This allowed us to determine whether any of these four regions was completely homozygous within an individual, which indicates IBD with high confidence. In the highly nomadic wild population, we did not find a single case of IBD, implying that inbreeding must be extremely rare (propIBD=0-0.00094, 95% CI). In the captive population, a five-generation pedigree strongly underestimated the average amount of realized inbreeding (FPed=0.013<propIBD=0.064), as expected given that pedigree founders were already related. We suggest that this SNP-based technique is generally useful for quantifying inbreeding at the individual or population level, and we show analytically that it can capture inbreeding loops that reach back up to a few hundred generations.

Genome-wide inbreeding estimation within Lebanese communities using SNP arrays

Consanguineous marriages have been widely practiced in several global communities with varying rates depending on religion, culture, and geography. In consanguineous marriages, parents pass to their children autozygous segments known as homozygous by descent segments. In this study, single-nucleotide polymorphisms were analyzed in 165 unrelated Lebanese people from Greek Orthodox, Maronite, Shiite and Sunni communities. Runs of homozygosity, total inbreeding levels, remote consanguinity, and population admixture and structure were estimated. The inbreeding coefficient value was estimated to be 1.61% in offspring of unrelated parents over three generations and 8.33% in offspring of first cousins. From these values, remote consanguinity values, resulting from genetic drift or recurrent consanguineous unions, were estimated in offspring of unrelated and first-cousin parents to be 0.61 and 1.2%, respectively. This remote consanguinity value suggests that for any unrelated marriages in Lebanon, the mates could be related as third cousins or as second cousins once removed. Under the assumption that 25% of marriages occur between first cousins, the mean inbreeding value of 2.3% may explain the increased incidence of recessive disease in offspring. Our analysis reveals a common ancestral population in the four Lebanese communities we studied.

Pedigree structure and kinship measurements of a mid-Michigan community: a new North American population isolate identified

Previous studies identified a cluster of individuals with an autosomal recessive form of deafness that reside in a small region of mid-Michigan. We hypothesized that affected members from this community descend from a defined founder population. Using public records and personal interviews, we constructed a genealogical database that includes the affected individuals and their extended families as descendants of 461 settlers who emigrated from the Eifel region of Germany between 1836 and 1875. The genealogical database represents a 13-generation pedigree that includes 27,747 descendants of these settlers. Among these descendants, 13,784 are presumed living. Many of the extant descendants reside in a 90-square-mile area, and 52% were born to parents who share at least one common ancestor. Among those born to related parents, the median kinship coefficient is 3.7 × 10(-3). While the pedigree contains 2,510 founders, 344 of the 461 settlers accounted for 67% of the genome in the extant population. These data suggest that we identified a new population isolate in North America and that, as demonstrated for congenital hearing loss, this rural mid-Michigan community is a new resource to discover heritable factors that contribute to common health-related conditions.

High inbreeding in the Faroe Islands does not appear to constitute a risk factor for multiple sclerosis

Background:

Large population-based genome-wide association studies have identified several multiple sclerosis (MS) genetic risk variants, but the existing missing heritability warrants different strategies. Isolated populations offer an alternative way of searching for rare genetic variants and evaluating the possible role of consanguinity in the development of MS. Studies of consanguinity and MS risk have yielded conflicting results.

Objectives:

In this study we investigated the role of consanguinity on MS risk in the relatively isolated Faroe Islands, which have a presumed high level of inbreeding.

Methods:

A total of 29 cases and 28 matched controls were genotyped and assessed for inbreeding coefficients, number of runs of homozygosity (ROH) at different lengths and observed number of homozygotes as measures of relatedness. Parametric and non-parametric statistical models were applied.

Results:

Both cases and controls exhibited considerable relatedness demonstrated by very high inbreeding coefficients, large number of observed homozygotes and many long ROH. However, apart from the number of ROH ≥ 2.5 mega base pairs, no significant differences between the two groups were observed.

Conclusions:

Overall, no significant difference between cases and controls were found, indicating that consanguinity in itself does not appear to be an important risk factor for MS in the population of the Faroe Islands.

Genetic characterization of Greek population isolates reveals strong genetic drift at missense and trait-associated variants

Isolated populations are emerging as a powerful study design in the search for low-frequency and rare variant associations with complex phenotypes. Here we genotype 2,296 samples from two isolated Greek populations, the Pomak villages (HELIC-Pomak) in the North of Greece and the Mylopotamos villages (HELIC-MANOLIS) in Crete. We compare their genomic characteristics to the general Greek population and establish them as genetic isolates. In the MANOLIS cohort, we observe an enrichment of missense variants among the variants that have drifted up in frequency by more than fivefold. In the Pomak cohort, we find novel associations at variants on chr11p15.4 showing large allele frequency increases (from 0.2% in the general Greek population to 4.6% in the isolate) with haematological traits, for example, with mean corpuscular volume (rs7116019, P=2.3 × 10(-26)). We replicate this association in a second set of Pomak samples (combined P=2.0 × 10(-36)). We demonstrate significant power gains in detecting medical trait associations.

Methods for collapsing multiple rare variants in whole-genome sequence data

Genetic Analysis Workshop 18 provided whole-genome sequence data in a pedigree-based sample and longitudinal phenotype data for hypertension and related traits, presenting an excellent opportunity for evaluating analysis choices. We summarize the nine contributions to the working group on collapsing methods, which evaluated various approaches for the analysis of multiple rare variants. One contributor defined a variant prioritization scheme, whereas the remaining eight contributors evaluated statistical methods for association analysis. Six contributors chose the gene as the genomic region for collapsing variants, whereas three contributors chose nonoverlapping sliding windows across the entire genome. Statistical methods spanned most of the published methods, including well-established burden tests, variance-components-type tests, and recently developed hybrid approaches. Lesser known methods, such as functional principal components analysis, higher criticism, and homozygosity association, and some newly introduced methods were also used. We found that performance of these methods depended on the characteristics of the genomic region, such as effect size and direction of variants under consideration. Except for MAP4 and FLT3, the performance of all statistical methods to identify rare casual variants was disappointingly poor, providing overall power almost identical to the type I error. This poor performance may have arisen from a combination of (1) small sample size, (2) small effects of most of the causal variants, explaining a small fraction of variance, (3) use of incomplete annotation information, and (4) linkage disequilibrium between causal variants in a gene and noncausal variants in nearby genes. Our findings demonstrate challenges in analyzing rare variants identified from sequence data.

Evaluation of bioelectrical impedance analysis for identifying overweight individuals at increased cardiometabolic risk: a cross-sectional study

OBJECTIVE

To investigate whether bioelectrical impedance analysis could be used to identify overweight individuals at increased cardiometabolic risk, defined as the presence of metabolic syndrome and/or diabetes.

DESIGN AND METHODS

Cross-sectional study of a Scottish population including 1210 women and 788 men. The diagnostic performance of thresholds of percentage body fat measured by bioelectrical impedance analysis to identify people at increased cardiometabolic risk was assessed using receiver-operating characteristic curves. Odds ratios for increased cardiometabolic risk in body mass index categories associated with values above compared to below sex-specific percentage body fat thresholds with optimal diagnostic performance were calculated using multivariable logistic regression analyses. The validity of bioelectrical impedance analysis to measure percentage body fat in this population was tested by examining agreement between bioelectrical impedance analysis and dual-energy X-ray absorptiometry in a subgroup of individuals.

RESULTS

Participants were aged 16-91 years and the optimal bioelectrical impedance analysis cut-points for percentage body fat for identifying people at increased cardiometabolic risk were 25.9% for men and 37.1% for women. Stratifying by these percentage body fat cut-points, the prevalence of increased cardiometabolic risk was 48% and 38% above the threshold and 24% and 19% below these thresholds for men and women, respectively. By comparison, stratifying by percentage body fat category had little impact on identifying increased cardiometabolic risk in normal weight and obese individuals. Fully adjusted odds ratios of being at increased cardiometabolic risk among overweight people with percentage body fat ≥ 25.9/37.1% compared with percentage body fat <25.9/37.1% as a reference were 1.93 (95% confidence interval: 1.20-3.10) for men and 1.79 (1.10-2.92) for women.

CONCLUSION

Percentage body fat measured using bioelectrical impedance analysis above a sex-specific threshold could be used in overweight people to identify individuals at increased cardiometabolic risk, who could benefit from risk factor management.

Genomic Analyses of Patients With Unexplained Early-Onset Scoliosis

STUDY DESIGN

To test for rare genetic mutations, a cohort of patients with unexplained early-onset scoliosis (EOS) was screened using high-density microarray genotyping. A cohort of patients with adolescent idiopathic scoliosis (AIS) was similarly screened and the results were compared.

SUMMARY OF BACKGROUND DATA

Patients with scoliosis in infancy or early childhood (EOS) are at high risk for progressive deformity and associated problems including respiratory compromise. Early-onset scoliosis is frequently associated with genetic disorders but many patients present with nonspecific clinical features and without an associated diagnosis. The authors hypothesized that EOS in these patients may be caused by rare genetic mutations detectable by next-generation genomic methods.

METHODS

The researchers identified 24 patients with unexplained EOS from pediatric orthopedic clinics. They genotyped them, along with 39 connecting family members, using the Illumina OmniExpress-12, version 1.0 beadchip. Resulting genotypes were analyzed for chromosomal changes, specifically copy number variation and absence of heterozygosity. They screened 482 adolescent idiopathic scoliosis (AIS) patients and 744 healthy controls, who were similarly genotyped with the same beadchip, for chromosomal changes identified in the EOS cohort.

RESULTS

Copy number variation and absence of heterozygosity analyses revealed a genetic diagnosis of chromosome 15q24 microdeletion syndrome in 1 patient and maternal uniparental disomy of chromosome 14 in a second one. Prior genetic testing and clinical evaluations had been negative in both cases. A large novel chromosome 10 deletion was likely causal in a third EOS patient. These mutations identified in the EOS patients were absent in AIS patients and controls, and thus were not associated with AIS or found in asymptomatic individuals.

CONCLUSIONS

These data underscore the usefulness of updated genetic evaluations including high-density microarray-based genotyping and other next-generation methods in patients with unexplained EOS, even when prior genetic studies were negative. These data also suggest the intriguing possibility that other mutations detectable by whole genome sequencing, as well as epigenetic effects, await discovery in the EOS population.

Characterizing runs of homozygosity and their impact on risk for psychosis in a population isolate

An increased abundance of runs of homozygosity (ROH) has been associated with risk for various diseases, including schizophrenia. Here we investigate the characteristics of ROH in Palau, an Oceanic population, evaluating whether these characteristics are related to risk for psychotic disorders and the nature of this association. To accomplish these aims we evaluate a sample of 203 cases with schizophrenia and related psychotic disorders-representing almost complete ascertainment of affected individuals in the population-and contrast their ROH to that of 125 subjects chosen to function as controls. While Palauan diagnosed with psychotic disorders tend to have slightly more ROH regions than controls, the distinguishing features are that they have longer ROH regions, greater total length of ROH, and their ROH tends to co-occur more often at the same locus. The nature of the sample allows us to investigate whether rare, highly penetrant recessive variants generate such case-control differences in ROH. Neither rare, highly penetrant recessive variants nor individual common variants of large effect account for a substantial proportion of risk for psychosis in Palau. These results suggest a more nuanced model for risk is required to explain patterns of ROH for this population.

Regions of homozygosity identified by oligonucleotide SNP arrays: evaluating the incidence and clinical utility

Copy neutral segments with allelic homozygosity, also known as regions of homozygosity (ROHs), are frequently identified in cases interrogated by oligonucleotide single-nucleotide polymorphism (oligo-SNP) microarrays. Presence of ROHs may be because of parental relatedness, chromosomal recombination or rearrangements and provides important clues regarding ancestral homozygosity, consanguinity or uniparental disomy. In this study of 14 574 consecutive cases, 832 (6%) were found to harbor one or more ROHs over 10 Mb, of which 651 cases (78%) had multiple ROHs, likely because of identity by descent (IBD), and 181 cases (22%) with ROHs involving a single chromosome. Parental relatedness was predicted to be first degree or closer in 5%, second in 9% and third in 19%. Of the 181 cases, 19 had ROHs for a whole chromosome revealing uniparental isodisomy (isoUPD). In all, 25 cases had significant ROHs involving a single chromosome; 5 cases were molecularly confirmed to have a mixed iso- and heteroUPD15 and 1 case each with segmental UPD9pat and segmental UPD22mat; 17 cases were suspected to have a mixed iso- and heteroUPD including 2 cases with small supernumerary marker and 2 cases with mosaic trisomy. For chromosome 15, 12 (92%) of 13 molecularly studied cases had either Prader-Willi or Angelman syndrome. Autosomal recessive disorders were confirmed in seven of nine cases from eight families because of the finding of suspected gene within a ROH. This study demonstrates that ROHs are much more frequent than previously recognized and often reflect parental relatedness, ascertain autosomal recessive diseases or unravel UPD in many cases.

Evaluation of European coeliac disease risk variants in a north Indian population

Studies in European populations have contributed to a better understanding of the genetics of complex diseases, for example, in coeliac disease (CeD), studies of over 23 000 European samples have reported association to the HLA locus and another 39 loci. However, these associations have not been evaluated in detail in other ethnicities. We sought to better understand how disease-associated loci that have been mapped in Europeans translate to a disease risk for a population with a different ethnic background. We therefore performed a validation of European risk loci for CeD in 497 cases and 736 controls of north Indian origin. Using a dense-genotyping platform (Immunochip), we confirmed the strong association to the HLA region (rs2854275, P=8.2 × 10⁻⁴⁹). Three loci showed suggestive association (rs4948256, P=9.3 × 10⁻⁷, rs4758538, P=8.6 × 10⁻⁵ and rs17080877, P=2.7 × 10⁻⁵). We directly replicated five previously reported European variants (P<0.05; mapping to loci harbouring FASLG/TNFSF18, SCHIP1/IL12A, PFKFB3/PRKCQ, ZMIZ1 and ICOSLG). Using a transferability test, we further confirmed association at PFKFB3/PRKCQ (rs2387397, P=2.8 × 10⁻⁴) and PTPRK/THEMIS (rs55743914, P=3.4 × 10⁻⁴). The north Indian population has a higher degree of consanguinity than Europeans and we therefore explored the role of recessively acting variants, which replicated the HLA locus (rs9271850, P=3.7 × 10⁻²³) and suggested a role of additional four loci. To our knowledge, this is the first replication study of CeD variants in a non-European population.

Great ape genomics

The great ape families are the species most closely related to our own, comprising chimpanzees, bonobos, gorillas, and orangutans. They live exclusively in tropical rainforests in Central Africa and the islands of Southeast Asia. Due to their close evolutionary relationship with humans, great apes share many cognitive, physiological, and morphological similarities with humans. The members of the great ape family make obvious models to facilitate the further understanding about humans' biology and history. This review will discuss how the recent addition of genome-wide data from great apes has furthered humans' understanding of these species and humanity, especially in the realm of evolutionary genetics.

Random inbreeding, isonymy, and population isolates in Argentina

Population isolates are an important tool in identifying and mapping genes of Mendelian diseases and complex traits. The geographical identification of isolates represents a priority from a genetic and health care standpoint. The purpose of this study is to analyze the spatial distribution of consanguinity by random isonymy (F ST) in Argentina and its relationship with the isolates previously identified in the country. F ST was estimated from the surname distribution of 22.6 million electors registered for the year 2001 in the 24 provinces, 5 geographical regions, and 510 departments of the country. Statistically significant spatial clustering of F ST was determined using the SaTScan V5.1 software. F ST exhibited a marked regional and departamental variation, showing the highest values towards the North and West of Argentina. The clusters of high consanguinity by random isonymy followed the same distribution. Recognized Argentinean genetic isolates are mainly localized at the north of the country, in clusters of high inbreeding. Given the availability of listings of surnames in high-capacity storage devices for different countries, estimating F ST from them can provide information on inbreeding for all levels of administrative subdivisions, to be used as a demographic variable for the identification of isolates within the country for public health purposes.