An AluYa5 Insertion in the 3'UTR of COL4A1 and Cerebral Small Vessel Disease.

Importance: Cerebral small vessel diseases (CSVDs) account for one-fifth of stroke cases. Numerous familial cases remain unresolved after routine screening of known CSVD genes. Objective: To identify novel genes and mechanisms associated with familial CSVD. Design, Setting, and Participants: This 2-stage study involved linkage analysis and a case-control study; linkage analysis and whole exome and genome sequencing were used to identify candidate gene variants in 2 large families with CSVD (9 patients with CSVD). Then, a case-control analysis was conducted on 246 unrelated probands, including probands from these 2 families and 244 additional probands. All probands (clinical onset

Tracking clusters of patients over time enables extracting information from medico-administrative databases.

CONTEXT: Identifying clusters (i.e., subgroups) of patients from the analysis of medico-administrative databases is particularly important to better understand disease heterogeneity. However, these databases contain different types of longitudinal variables which are measured over different follow-up periods, generating truncated data. It is therefore fundamental to develop clustering approaches that can handle this type of data. OBJECTIVE: We propose here cluster-tracking approaches to identify clusters of patients from truncated longitudinal data contained in medico-administrative databases. MATERIAL AND METHODS: We first cluster patients at each age. We then track the identified clusters over ages to construct cluster-trajectories. We compared our novel approaches with three classical longitudinal clustering approaches by calculating the silhouette score. As a use-case, we analyzed antithrombotic drugs used from 2008 to 2018 contained in the Échantillon Généraliste des Bénéficiaires (EGB), a French national cohort. RESULTS: Our cluster-tracking approaches allow us to identify several cluster-trajectories with clinical significance without any imputation of data. The comparison of the silhouette scores obtained with the different approaches highlights the better performances of the cluster-tracking approaches. CONCLUSION: The cluster-tracking approaches are a novel and efficient alternative to identify patient clusters from medico-administrative databases by taking into account their specificities.

Deficiency of the minor spliceosome component U4atac snRNA secondarily results in ciliary defects in human and zebrafish.

In the human genome, about 750 genes contain one intron excised by the minor spliceosome. This spliceosome comprises its own set of snRNAs, among which U4atac. Its noncoding gene, RNU4ATAC, has been found mutated in Taybi-Linder (TALS/microcephalic osteodysplastic primordial dwarfism type 1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes. These rare developmental disorders, whose physiopathological mechanisms remain unsolved, associate ante- and post-natal growth retardation, microcephaly, skeletal dysplasia, intellectual disability, retinal dystrophy, and immunodeficiency. Here, we report bi-allelic RNU4ATAC mutations in five patients presenting with traits suggestive of the Joubert syndrome (JBTS), a well-characterized ciliopathy. These patients also present with traits typical of TALS/RFMN/LWS, thus widening the clinical spectrum of RNU4ATAC-associated disorders and indicating ciliary dysfunction as a mechanism downstream of minor splicing defects. Intriguingly, all five patients carry the n.16G>A mutation, in the Stem II domain, either at the homozygous or compound heterozygous state. A gene ontology term enrichment analysis on minor intron-containing genes reveals that the cilium assembly process is over-represented, with no less than 86 cilium-related genes containing at least one minor intron, among which there are 23 ciliopathy-related genes. The link between RNU4ATAC mutations and ciliopathy traits is supported by alterations of primary cilium function in TALS and JBTS-like patient fibroblasts, as well as by u4atac zebrafish model, which exhibits ciliopathy-related phenotypes and ciliary defects. These phenotypes could be rescued by WT but not by pathogenic variants-carrying human U4atac. Altogether, our data indicate that alteration of cilium biogenesis is part of the physiopathological mechanisms of TALS/RFMN/LWS, secondarily to defects of minor intron splicing.

Mutations in the non-coding RNU4ATAC gene affect the homeostasis and function of the Integrator complex.

Various genetic diseases associated with microcephaly and developmental defects are due to pathogenic variants in the U4atac small nuclear RNA (snRNA), a component of the minor spliceosome essential for the removal of U12-type introns from eukaryotic mRNAs. While it has been shown that a few RNU4ATAC mutations result in impaired binding of essential protein components, the molecular defects of the vast majority of variants are still unknown. Here, we used lymphoblastoid cells derived from RNU4ATAC compound heterozygous (g.108_126del;g.111G>A) twin patients with MOPD1 phenotypes to analyze the molecular consequences of the mutations on small nuclear ribonucleoproteins (snRNPs) formation and on splicing. We found that the U4atac108_126del mutant is unstable and that the U4atac111G>A mutant as well as the minor di- and tri-snRNPs are present at reduced levels. Our results also reveal the existence of 3'-extended snRNA transcripts in patients' cells. Moreover, we show that the mutant cells have alterations in splicing of INTS7 and INTS10 minor introns, contain lower levels of the INTS7 and INTS10 proteins and display changes in the assembly of Integrator subunits. Altogether, our results show that compound heterozygous g.108_126del;g.111G>A mutations induce splicing defects and affect the homeostasis and function of the Integrator complex.

Tracking Temporal Clusters from Patient Networks.

Creating homogeneous groups (clusters) of patients from medico-administrative databases provides a better understanding of health determinants. But in these databases, patients have truncated care pathways. We developed an approach based on patient networks to construct care trajectories from such truncated data. We tested this approach on antithrombotic treatments prescribed from 2008 to 2018 contained in the échantillon généraliste des bénéficiaires (EGB). We constructed a patient network for each patients' age (years from birth). We then applied the Markov clustering algorithm in each network. The care trajectories were finally constructed by matching clusters identified in two consecutive networks. We calculated the silhouette score to assess the performance of this network approach compared to three existing approaches. We identified 12 care trajectories that we were able to associate with pathologies. The best silhouette score was obtained for the network approach. Our approach allowed to highlight care trajectories taking into account the longitudinal, multidimensional and truncated nature of data from medico-administrative databases.

COL4A1/COL4A2 and inherited platelet disorder gene variants in fetuses showing intracranial hemorrhage.

BACKGROUND: Variants of COL4A1/COL4A2 genes have been reported in fetal intracranial hemorrhage (ICH) cases but their prevalence and characteristics have not been established in a large series of fetuses. Fetal neonatal alloimmune thrombocytopenia is a major acquired ICH factor but the prevalence and characteristics of inherited platelet disorder (IPD) gene variants leading to thrombocytopenia are unknown. Herein, we screened COL4A1/COL4A2 and IPD genes in a large series of ICH fetuses. METHODS: A cohort of 194 consecutive ICH fetuses were first screened for COL4A1/COL4A2 variants. We manually curated a list of 64 genes involved in IPD and investigated them in COL4A1/COL4A2 negative fetuses, using exome sequencing data from 101 of these fetuses. RESULT: Pathogenic variants of COL4A1/COL4A2 genes were identified in 36 fetuses (19%). They occurred de novo in 70% of the 32 fetuses for whom parental DNA was available. Pathogenic variants in two megakaryopoiesis genes (MPL and MECOM genes) were identified in two families with recurrent and severe fetal ICH, with variable extraneurological pathological features. CONCLUSION: Our study emphasizes the genetic heterogeneity of fetal ICH and the need to screen both COL4A1/COL4A2 and IPD genes in the etiological investigation of fetal ICH to allow proper genetic counseling.

Genetics of PlGF plasma levels highlights a role of its receptors and supports the link between angiogenesis and immunity.

Placental growth factor (PlGF) is a member of the vascular endothelial growth factor family and is involved in bone marrow-derived cell activation, endothelial stimulation and pathological angiogenesis. High levels of PlGF have been observed in several pathological conditions especially in cancer, cardiovascular, autoimmune and inflammatory diseases. Little is known about the genetics of circulating PlGF levels. Indeed, although the heritability of circulating PlGF levels is around 40%, no studies have assessed the relation between PlGF plasma levels and genetic variants at a genome-wide level. In the current study, PlGF plasma levels were measured in a population-based sample of 2085 adult individuals from three isolated populations of South Italy. A GWAS was performed in a discovery cohort (N = 1600), followed by a de novo replication (N = 468) from the same populations. The meta-analysis of the discovery and replication samples revealed one signal significantly associated with PlGF circulating levels. This signal was mapped to the PlGF co-receptor coding gene NRP1, indicating its important role in modulating the PlGF plasma levels. Two additional signals, at the PlGF receptor coding gene FLT1 and RAPGEF5 gene, were identified at a suggestive level. Pathway and TWAS analyses highlighted genes known to be involved in angiogenesis and immune response, supporting the link between these processes and PlGF regulation. Overall, these data improve our understanding of the genetic variation underlying circulating PlGF levels. This in turn could lead to new preventive and therapeutic strategies for a wide variety of PlGF-related pathologies.

Multiethnic genome-wide association study of differentiated thyroid cancer in the EPITHYR consortium.

Incidence of differentiated thyroid carcinoma (DTC) varies considerably between ethnic groups, with particularly high incidence rates in Pacific Islanders. DTC is one of the cancers with the highest familial risk suggesting a major role of genetic risk factors, but only few susceptibility loci were identified so far. In order to assess the contribution of known DTC susceptibility loci and to identify new ones, we conducted a multiethnic genome-wide association study (GWAS) in individuals of European ancestry and of Oceanian ancestry from Pacific Islands. Our study included 1554 cases/1973 controls of European ancestry and 301 cases/348 controls of Oceanian ancestry from seven population-based case-control studies participating to the EPITHYR consortium. All participants were genotyped using the OncoArray-500K Beadchip (Illumina). We confirmed the association with the known DTC susceptibility loci at 2q35, 8p12, 9q22.33 and 14q13.3 in the European ancestry population and suggested two novel signals at 1p31.3 and 16q23.2, which were associated with thyroid-stimulating hormone levels in previous GWAS. We additionally replicated an association with 5p15.33 reported previously in Chinese and European populations. Except at 1p31.3, all associations were in the same direction in the population of Oceanian ancestry. We also observed that the frequencies of risk alleles at 2q35, 5p15.33 and 16q23.2 were significantly higher in Oceanians than in Europeans. However, additional GWAS and epidemiological studies in Oceanian populations are needed to fully understand the highest incidence observed in these populations.

Clinical interpretation of variants identified in RNU4ATAC, a non-coding spliceosomal gene.

Biallelic variants in RNU4ATAC, a non-coding gene transcribed into the minor spliceosome component U4atac snRNA, are responsible for three rare recessive developmental diseases, namely Taybi-Linder/MOPD1, Roifman and Lowry-Wood syndromes. Next-generation sequencing of clinically heterogeneous cohorts (children with either a suspected genetic disorder or a congenital microcephaly) recently identified mutations in this gene, illustrating how profoundly these technologies are modifying genetic testing and assessment. As RNU4ATAC has a single non-coding exon, the bioinformatic prediction algorithms assessing the effect of sequence variants on splicing or protein function are irrelevant, which makes variant interpretation challenging to molecular diagnostic laboratories. In order to facilitate and improve clinical diagnostic assessment and genetic counseling, we present i) an update of the previously reported RNU4ATAC mutations and an analysis of the genetic variations affecting this gene using the Genome Aggregation Database (gnomAD) resource; ii) the pathogenicity prediction performances of scores computed based on an RNA structure prediction tool and of those produced by the Combined Annotation Dependent Depletion tool for the 285 RNU4ATAC variants identified in patients or in large-scale sequencing projects; iii) a method, based on a cellular assay, that allows to measure the effect of RNU4ATAC variants on splicing efficiency of a minor (U12-type) reporter intron. Lastly, the concordance of bioinformatic predictions and cellular assay results was investigated.

Novel Chronic Mouse Model of Cerebral Cavernous Malformations.

Background and Purpose- Cerebral cavernous malformations (CCMs) are vascular malformations of the brain that lead to cerebral hemorrhages. A pharmacological treatment is needed especially for patients with nonoperable deep-seated lesions. We and others obtained CCM mouse models that were useful for mechanistic studies and rapid trials testing the preventive effects of candidate drugs. The shortened lifespan of acute mouse models hampered evaluation of compounds that would not only prevent lesion appearance but also cure preexisting lesions. Indirubin-3'-monoxime previously demonstrated its efficacy to reverse the cardiac phenotype of ccm2m201 zebrafish mutants and to prevent lesion development in an acute CCM2 mouse model. In the present article, we developed and characterized a novel chronic CCM2 mouse model and evaluated the curative therapeutic effect of indirubin-3'-monoxime after CCM lesion development. Methods- The chronic mouse model was obtained by a postnatal induction of brain-endothelial-cell-specific ablation of the Ccm2 gene using the inducible Slco1c1-CreERT2 mouse line. Results- We obtained a fully penetrant novel CCM chronic mouse model without any obvious off-target phenotypes and compatible with long-term survival. By 3 months of age, CCM lesions ranging in size from small isolated lesions to multiple caverns developed throughout the brain. Lesion burden was quantified in animals from 1 week to 5 months of age. Clear signs of intracerebral hemorrhages were noticed in brain-endothelial-cell-specific ablation of the Ccm2 gene. In contrast with its preventive effect in the acute CCM2 mouse model, a 20 mg/kg indirubin-3'-monoxime treatment for 3 weeks in 3-month old animals neither had any beneficial effect on the lesion burden nor alleviated cerebral hemorrhages. Conclusions- The brain-endothelial-cell-specific ablation of the Ccm2 gene chronic model is a strongly improved disease model for the CCM community whose challenge today is to decipher which candidate drugs might have a curative effect on patients' preexisting lesions. Visual Overview- An online visual overview is available for this article.

New insights into minor splicing-a transcriptomic analysis of cells derived from TALS patients.

Minor intron splicing plays a central role in human embryonic development and survival. Indeed, biallelic mutations in RNU4ATAC, transcribed into the minor spliceosomal U4atac snRNA, are responsible for three rare autosomal recessive multimalformation disorders named Taybi-Linder (TALS/MOPD1), Roifman (RFMN), and Lowry-Wood (LWS) syndromes, which associate numerous overlapping signs of varying severity. Although RNA-seq experiments have been conducted on a few RFMN patient cells, none have been performed in TALS, and more generally no in-depth transcriptomic analysis of the ∼700 human genes containing a minor (U12-type) intron had been published as yet. We thus sequenced RNA from cells derived from five skin, three amniotic fluid, and one blood biosamples obtained from seven unrelated TALS cases and from age- and sex-matched controls. This allowed us to describe for the first time the mRNA expression and splicing profile of genes containing U12-type introns, in the context of a functional minor spliceosome. Concerning RNU4ATAC-mutated patients, we show that as expected, they display distinct U12-type intron splicing profiles compared to controls, but that rather unexpectedly mRNA expression levels are mostly unchanged. Furthermore, although U12-type intron missplicing concerns most of the expressed U12 genes, the level of U12-type intron retention is surprisingly low in fibroblasts and amniocytes, and much more pronounced in blood cells. Interestingly, we found several occurrences of introns that can be spliced using either U2, U12, or a combination of both types of splice site consensus sequences, with a shift towards splicing using preferentially U2 sites in TALS patients' cells compared to controls.

Detecting the dominance component of heritability in isolated and outbred human populations.

Inconsistencies between published estimates of dominance heritability between studies of human genetic isolates and human outbred populations incite investigation into whether such differences result from particular trait architectures or specific population structures. We analyse simulated datasets, characteristic of genetic isolates and of unrelated individuals, before analysing the isolate of Cilento for various commonly studied traits. We show the strengths of using genetic relationship matrices for variance decomposition over identity-by-descent based methods in a population isolate and that heritability estimates in isolates will avoid the downward biases that may occur in studies of samples of unrelated individuals; irrespective of the simulated distribution of causal variants. Yet, we also show that precise estimates of dominance in isolates are demonstrably problematic in the presence of shared environmental effects and such effects should be accounted for. Nevertheless, we demonstrate how studying isolates can help determine the existence or non-existence of dominance for complex traits, and we find strong indications of non-zero dominance for low-density lipoprotein level in Cilento. Finally, we recommend future study designs to analyse trait variance decomposition from ensemble data across multiple population isolates.

Strategies for phasing and imputation in a population isolate.

In the search for genetic associations with complex traits, population isolates offer the advantage of reduced genetic and environmental heterogeneity. In addition, cost-efficient next-generation association approaches have been proposed in these populations where only a subsample of representative individuals is sequenced and then genotypes are imputed into the rest of the population. Gene mapping in such populations thus requires high-quality genetic imputation and preliminary phasing. To identify an effective study design, we compare by simulation a range of phasing and imputation software and strategies. We simulated 1,115,604 variants on chromosome 10 for 477 members of the large complex pedigree of Campora, a village within the established isolate of Cilento in southern Italy. We assessed the phasing performance of identical by descent based software ALPHAPHASE and SLRP, LD-based software SHAPEIT2, SHAPEIT3, and BEAGLE, and new software EAGLE that combines both methodologies. For imputation we compared IMPUTE2, IMPUTE4, MINIMAC3, BEAGLE, and new software PBWT. Genotyping errors and missing genotypes were simulated to observe their effects on the performance of each software. Highly accurate phased data were achieved by all software with SHAPEIT2, SHAPEIT3, and EAGLE2 providing the most accurate results. MINIMAC3, IMPUTE4, and IMPUTE2 all performed strongly as imputation software and our study highlights the considerable gain in imputation accuracy provided by a genome sequenced reference panel specific to the population isolate.

Variation in worldwide incidence of amyotrophic lateral sclerosis: a meta-analysis.

Background: To assess the worldwide variation of amyotrophic lateral sclerosis (ALS) incidence, we performed a systematic review and meta-analysis of population-based data published to date. Methods: We reviewed Medline and Embase up to June 2015 and included all population-based studies of newly diagnosed ALS cases, using multiple sources for case ascertainment. ALS crude and standardized incidence (on age and sex using the US 2010 population) were calculated. Random effect meta-analysis and meta-regression were performed using the subcontinent as the main study level covariate. Sources of heterogeneity related to the characteristics of the study population and the study methodology were investigated. Results: Among 3216 records, 44 studies were selected, covering 45 geographical areas in 11 sub-continents. A total of 13 146 ALS cases and 825 million person-years of follow-up (PYFU) were co-nsidered. The overall pooled worldwide crude ALS incidence was at 1.75 (1.55-1.96)/100 000 PYFU; 1.68 (1.50-1.85)/100 000 PYFU after standardization. Heterogeneity was identified in ALS standardized incidence between North Europe [1.89 (1.46-2.32)/100 000 PYFU] and East Asia [0.83 (0.42-1.24)/100 000 PYFU, China and Japan P = 0.001] or South Asia [0.73 (0.58-0.89)/100 000/PYFU Iran, P = 0.02]. Conversely, homogeneous rates have been reported in populations from Europe, North America and New Zealand [pooled ALS standardized incidence of 1.81 (1.66-1.97)/100 000 PYFU for those areas]. Conclusion: This review confirms a heterogeneous distribution worldwide of ALS, and sets the scene to sustain a collaborative study involving a wide international consortium to investigate the link between ancestry, environment and ALS incidence.

Loss of VPS13C Function in Autosomal-Recessive Parkinsonism Causes Mitochondrial Dysfunction and Increases PINK1/Parkin-Dependent Mitophagy.

Autosomal-recessive early-onset parkinsonism is clinically and genetically heterogeneous. The genetic causes of approximately 50% of autosomal-recessive early-onset forms of Parkinson disease (PD) remain to be elucidated. Homozygozity mapping and exome sequencing in 62 isolated individuals with early-onset parkinsonism and confirmed consanguinity followed by data mining in the exomes of 1,348 PD-affected individuals identified, in three isolated subjects, homozygous or compound heterozygous truncating mutations in vacuolar protein sorting 13C (VPS13C). VPS13C mutations are associated with a distinct form of early-onset parkinsonism characterized by rapid and severe disease progression and early cognitive decline; the pathological features were striking and reminiscent of diffuse Lewy body disease. In cell models, VPS13C partly localized to the outer membrane of mitochondria. Silencing of VPS13C was associated with lower mitochondrial membrane potential, mitochondrial fragmentation, increased respiration rates, exacerbated PINK1/Parkin-dependent mitophagy, and transcriptional upregulation of PARK2 in response to mitochondrial damage. This work suggests that loss of function of VPS13C is a cause of autosomal-recessive early-onset parkinsonism with a distinctive phenotype of rapid and severe progression.

Six Novel Loci Associated with Circulating VEGF Levels Identified by a Meta-analysis of Genome-Wide Association Studies.

Vascular endothelial growth factor (VEGF) is an angiogenic and neurotrophic factor, secreted by endothelial cells, known to impact various physiological and disease processes from cancer to cardiovascular disease and to be pharmacologically modifiable. We sought to identify novel loci associated with circulating VEGF levels through a genome-wide association meta-analysis combining data from European-ancestry individuals and using a dense variant map from 1000 genomes imputation panel. Six discovery cohorts including 13,312 samples were analyzed, followed by in-silico and de-novo replication studies including an additional 2,800 individuals. A total of 10 genome-wide significant variants were identified at 7 loci. Four were novel loci (5q14.3, 10q21.3, 16q24.2 and 18q22.3) and the leading variants at these loci were rs114694170 (MEF2C, P = 6.79 x 10(-13)), rs74506613 (JMJD1C, P = 1.17 x 10(-19)), rs4782371 (ZFPM1, P = 1.59 x 10(-9)) and rs2639990 (ZADH2, P = 1.72 x 10(-8)), respectively. We also identified two new independent variants (rs34528081, VEGFA, P = 1.52 x 10(-18); rs7043199, VLDLR-AS1, P = 5.12 x 10(-14)) at the 3 previously identified loci and strengthened the evidence for the four previously identified SNPs (rs6921438, LOC100132354, P = 7.39 x 10(-1467); rs1740073, C6orf223, P = 2.34 x 10(-17); rs6993770, ZFPM2, P = 2.44 x 10(-60); rs2375981, KCNV2, P = 1.48 x 10(-100)). These variants collectively explained up to 52% of the VEGF phenotypic variance. We explored biological links between genes in the associated loci using Ingenuity Pathway Analysis that emphasized their roles in embryonic development and function. Gene set enrichment analysis identified the ERK5 pathway as enriched in genes containing VEGF associated variants. eQTL analysis showed, in three of the identified regions, variants acting as both cis and trans eQTLs for multiple genes. Most of these genes, as well as some of those in the associated loci, were involved in platelet biogenesis and functionality, suggesting the importance of this process in regulation of VEGF levels. This work also provided new insights into the involvement of genes implicated in various angiogenesis related pathologies in determining circulating VEGF levels. The understanding of the molecular mechanisms by which the identified genes affect circulating VEGF levels could be important in the development of novel VEGF-related therapies for such diseases.

Clinical and demographic factors and outcome of amyotrophic lateral sclerosis in relation to population ancestral origin.

BACKGROUND: To review how the phenotype and outcome of amyotrophic lateral sclerosis (ALS) change with variations in population ancestral origin (PAO). Knowledge of how PAO modifies ALS phenotype may provide important insight into the risk factors and pathogenic mechanisms of the disease. METHODS: We performed a systematic review and meta-analysis of the literature concerning differences in phenotype and outcome of ALS that relate to PAO. RESULTS: A review of 3111 records identified 78 population-based studies. The 40 that were included covered 40 geographical areas in 10 subcontinents. Around 12,700 ALS cases were considered. The results highlight the phenotypic heterogeneity of ALS at time of onset [age, sex ratio (SR), bulbar onset], age at diagnosis, occurrence of comorbidities in the first year after diagnosis, and outcome (survival). Subcontinent is a major explanatory factor for the variability of the ALS phenotype in population-based studies. Some markers of ALS phenotype were homogeneously distributed in western countries (SR, mean age at onset/diagnosis) but their distributions in other subcontinents were remarkably different. Other markers presented variations in European subcontinents (familial ALS, bulbar onset) and in other continents. As a consequence, ALS outcome strongly varied, with a median survival time from onset ranging from 24 months (Northern Europe) to 48 months (Central Asia). DISCUSSION: This review sets the scene for a collaborative study involving a wide international consortium to investigate, using a standard methodology, the link between ancestry, environment, and ALS phenotype.

Relationship inference from the genetic data on parents or offspring: A comparative study.

Relationship inference in a population is of interest for many areas of research from anthropology to genetics. It is possible to directly infer the relationship between the two individuals in a couple from their genetic data or to indirectly infer it from the genetic data of one of their offspring. For this reason, one can wonder if it is more advantageous to sample couples or single individuals to study relationships of couples in a population. Indeed, sampling two individuals is more informative than sampling one as we are looking at four haplotypes instead of two, but it also doubles the cost of the study and is a more complex sampling scheme. To answer this question, we performed simulations of 1000 trios from 10 different relationships using real human haplotypes to have realistic genome-wide genetic data. Then, we compared the genome sharing coefficients and the relationship inference obtained from either a pair of individuals or one of their offspring using both single-point and multi-point approaches. We observed that for relationships closer than 1st cousin, pairs of individuals were more informative than one of their offspring for relationship inference, and kinship coefficients obtained from single-point methods gave more accurate or equivalent genome sharing estimations. For more remote relationships, offspring were more informative for relationship inference, and inbreeding coefficients obtained from multi-point methods gave more accurate genome sharing estimations. In conclusion, relationship inference on a parental pair or on one of their offspring provides complementary information. When possible, sampling trios should be encouraged as it could allow spanning a wider range of potential relationships.

High level of inbreeding in final phase of 1000 Genomes Project.

The 1000 Genomes Project provides a unique source of whole genome sequencing data for studies of human population genetics and human diseases. The last release of this project includes more than 2,500 sequenced individuals from 26 populations. Although relationships among individuals have been investigated in some of the populations, inbreeding has never been studied. In this article, we estimated the genomic inbreeding coefficient of each individual and found an unexpected high level of inbreeding in 1000 Genomes data: nearly a quarter of the individuals were inbred and around 4% of them had inbreeding coefficients similar or greater than the ones expected for first-cousin offspring. Inbred individuals were found in each of the 26 populations, with some populations showing proportions of inbred individuals above 50%. We also detected 227 previously unreported pairs of close relatives (up to and including first-cousins). Thus, we propose subsets of unrelated and outbred individuals, for use by the scientific community. In addition, because admixed populations are present in the 1000 Genomes Project, we performed simulations to study the robustness of inbreeding coefficient estimates in the presence of admixture. We found that our multi-point approach (FSuite) was quite robust to admixture, unlike single-point methods (PLINK).