The heritability of a trait is the proportion of its variance explained by genetic factors; it has historically been estimated using familial data. However, new methods have appeared for estimating heritabilities using genomewide data from unrelated individuals. A drawback of this strategy is that population stratification can bias the estimates. Indeed, an environmental factor associated with the phenotype may differ among population subgroups. This factor being associated both with the phenotype and the genetic variation in the population would be a confounder. A common solution consists in adjusting on the first Principal Components (PCs) of the genomic data. We study this procedure on simulated data and on 6000 individuals from the Three-City Study. We analyse the geographical coordinates of the birth cities, which are not genetically determined, but the heritability of which should be overestimated due to population stratification. We also analyse various anthropometric traits. The procedure fails to correct the bias in geographical coordinates heritability estimates. The heritability estimates of the anthropometric traits are affected by the inclusion of the first PC, but not by the following PCs, contrarily to geographical coordinates. We recommend to be cautious with heritability estimates obtained from a large population.

The European Prospective Investigation into Cancer and nutrition (EPIC) is a long-term, multi-centric prospective study in Europe investigating the relationships between cancer and nutrition. This study has served as a basis for a number of Genome-Wide Association Studies (GWAS) and other types of genetic analyses. Over a period of 5 years, 52,256 EPIC DNA samples have been extracted using an automated DNA extraction platform. Here we have evaluated the pre-analytical factors affecting DNA yield, including anthropometric, epidemiological and technical factors such as center of subject recruitment, age, gender, body-mass index, disease case or control status, tobacco consumption, number of aliquots of buffy coat used for DNA extraction, extraction machine or procedure, DNA quantification method, degree of haemolysis and variations in the timing of sample processing. We show that the largest significant variations in DNA yield were observed with degree of haemolysis and with center of subject recruitment. Age, gender, body-mass index, cancer case or control status and tobacco consumption also significantly impacted DNA yield. Feedback from laboratories which have analyzed DNA with different SNP genotyping technologies demonstrate that the vast majority of samples (approximately 88%) performed adequately in different types of assays. To our knowledge this study is the largest to date to evaluate the sources of pre-analytical variations in DNA extracted from peripheral leucocytes. The results provide a strong evidence-based rationale for standardized recommendations on blood collection and processing protocols for large-scale genetic studies.

Fibroblast growth factor 23 is an emerging vascular biomarker, recently associated with cerebral small vessel disease and poor cognition in patients on dialysis. It also interacts with klotho, an anti-aging and cognition enhancing protein.

Most modern citrus cultivars have an interspecific origin. As a foundational step towards deciphering the interspecific genome structures, a reference whole genome sequence was produced by the International Citrus Genome Consortium from a haploid derived from Clementine mandarin. The availability of a saturated genetic map of Clementine was identified as an essential prerequisite to assist the whole genome sequence assembly. Clementine is believed to be a 'Mediterranean' mandarin × sweet orange hybrid, and sweet orange likely arose from interspecific hybridizations between mandarin and pummelo gene pools. The primary goals of the present study were to establish a Clementine reference map using codominant markers, and to perform comparative mapping of pummelo, sweet orange, and Clementine.

Mutations in RYR1 give rise to diverse skeletal muscle phenotypes, ranging from classical central core disease to susceptibility to malignant hyperthermia. Next-generation sequencing has recently shown that RYR1 is implicated in a wide variety of additional myopathies, including centronuclear myopathy. In this work, we established an international cohort of 21 patients from 18 families with autosomal recessive RYR1-related centronuclear myopathy, to better define the clinical, imaging, and histological spectrum of this disorder. Early onset of symptoms with hypotonia, motor developmental delay, proximal muscle weakness, and a stable course were common clinical features in the cohort. Ptosis and/or ophthalmoparesis, facial weakness, thoracic deformities, and spinal involvement were also frequent but variable. A common imaging pattern consisted of selective involvement of the vastus lateralis, adductor magnus, and biceps brachii in comparison to adjacent muscles. In addition to a variable prominence of central nuclei, muscle biopsy from 20 patients showed type 1 fiber predominance and a wide range of intermyofibrillary architecture abnormalities. All families harbored compound heterozygous mutations, most commonly a truncating mutation combined with a missense mutation. This work expands the phenotypic characterization of patients with recessive RYR1-related centronuclear myopathy by highlighting common and variable clinical, histological, and imaging findings in these patients.

At the dawn of the second millennium, the expansion of the Indian Ocean trading network aligned with the emergence of an outward-oriented community along the East African coast to create a cosmopolitan cultural and trading zone known as the Swahili Corridor. On the basis of analyses of new genome-wide genotyping data and uniparental data in 276 individuals from coastal Kenya and the Comoros islands, along with large-scale genetic datasets from the Indian Ocean rim, we reconstruct historical population dynamics to show that the Swahili Corridor is largely an eastern Bantu genetic continuum. Limited gene flows from the Middle East can be seen in Swahili and Comorian populations at dates corresponding to historically documented contacts. However, the main admixture event in southern insular populations, particularly Comorian and Malagasy groups, occurred with individuals from Island Southeast Asia as early as the 8th century, reflecting an earlier dispersal from this region. Remarkably, our results support recent archaeological and linguistic evidence-based suggestions that the Comoros archipelago was the earliest location of contact between Austronesian and African populations in the Swahili Corridor.

RORα, the RAR-related orphan nuclear receptor alpha, is essential for cerebellar development. The spontaneous mutant mouse staggerer, with an ataxic gait caused by neurodegeneration of cerebellar Purkinje cells, was discovered two decades ago to result from homozygous intragenic Rora deletions. However, RORA mutations were hitherto undocumented in humans. Through a multi-centric collaboration, we identified three copy-number variant deletions (two de novo and one dominantly inherited in three generations), one de novo disrupting duplication, and nine de novo point mutations (three truncating, one canonical splice site, and five missense mutations) involving RORA in 16 individuals from 13 families with variable neurodevelopmental delay and intellectual disability (ID)-associated autistic features, cerebellar ataxia, and epilepsy. Consistent with the human and mouse data, disruption of the D. rerio ortholog, roraa, causes significant reduction in the size of the developing cerebellum. Systematic in vivo complementation studies showed that, whereas wild-type human RORA mRNA could complement the cerebellar pathology, missense variants had two distinct pathogenic mechanisms of either haploinsufficiency or a dominant toxic effect according to their localization in the ligand-binding or DNA-binding domains, respectively. This dichotomous direction of effect is likely relevant to the phenotype in humans: individuals with loss-of-function variants leading to haploinsufficiency show ID with autistic features, while individuals with de novo dominant toxic variants present with ID, ataxia, and cerebellar atrophy. Our combined genetic and functional data highlight the complex mutational landscape at the human RORA locus and suggest that dual mutational effects likely determine phenotypic outcome.

Spermatogenesis defects concern millions of men worldwide, yet the vast majority remains undiagnosed. Here we report men with primary infertility due to multiple morphological abnormalities of the sperm flagella with severe disorganization of the sperm axoneme, a microtubule-based structure highly conserved throughout evolution. Whole-exome sequencing was performed on 78 patients allowing the identification of 22 men with bi-allelic mutations in DNAH1 (n = 6), CFAP43 (n = 10), and CFAP44 (n = 6). CRISPR/Cas9 created homozygous CFAP43/44 male mice that were infertile and presented severe flagellar defects confirming the human genetic results. Immunoelectron and stimulated-emission-depletion microscopy performed on CFAP43 and CFAP44 orthologs in Trypanosoma brucei evidenced that both proteins are located between the doublet microtubules 5 and 6 and the paraflagellar rod. Overall, we demonstrate that CFAP43 and CFAP44 have a similar structure with a unique axonemal localization and are necessary to produce functional flagella in species ranging from Trypanosoma to human.

Neurodevelopmental disorders with periventricular nodular heterotopia (PNH) are etiologically heterogeneous, and their genetic causes remain in many cases unknown. Here we show that missense mutations in NEDD4L mapping to the HECT domain of the encoded E3 ubiquitin ligase lead to PNH associated with toe syndactyly, cleft palate and neurodevelopmental delay. Cellular and expression data showed sensitivity of PNH-associated mutants to proteasome degradation. Moreover, an in utero electroporation approach showed that PNH-related mutants and excess wild-type NEDD4L affect neurogenesis, neuronal positioning and terminal translocation. Further investigations, including rapamycin-based experiments, found differential deregulation of pathways involved. Excess wild-type NEDD4L leads to disruption of Dab1 and mTORC1 pathways, while PNH-related mutations are associated with deregulation of mTORC1 and AKT activities. Altogether, these data provide insights into the critical role of NEDD4L in the regulation of mTOR pathways and their contributions in cortical development.

The senescence of mammalian cells is characterized by a proliferative arrest in response to stress and the expression of an inflammatory phenotype. Here we show that histone H2A.J, a poorly studied H2A variant found only in mammals, accumulates in human fibroblasts in senescence with persistent DNA damage. H2A.J also accumulates in mice with aging in a tissue-specific manner and in human skin. Knock-down of H2A.J inhibits the expression of inflammatory genes that contribute to the senescent-associated secretory phenotype (SASP), and over expression of H2A.J increases the expression of some of these genes in proliferating cells. H2A.J accumulation may thus promote the signalling of senescent cells to the immune system, and it may contribute to chronic inflammation and the development of aging-associated diseases.

Cerebral β-amyloid (Aβ) deposition and atrophy are central features of Alzheimer disease. Studies of Alzheimer disease biomarkers have largely focused on Aβ in cerebrospinal fluid (CSF), and there is uncertainty as to what plasma Aβ may be a marker. We examined the association of Aβ levels in the plasma with magnetic resonance imaging (MRI)-markers of brain aging, including longitudinal changes in global and regional brain volumes, in dementia-free persons. We studied 1530 participants of the Three-City-Dijon cohort, aged 65-80 years. Plasma Aβ measurement and magnetic resonance imaging were performed at baseline and after a 4-year follow up. Total brain, gray matter, and hippocampal volume were estimated using voxel-based morphometry, and annualized change in brain volumes was calculated. Increased plasma Aβ1-40 was associated with lower baseline hippocampal volume. Although baseline plasma Aβ levels were not associated with longitudinal change in brain volumes, consistently high plasma Aβ1-40 levels were associated with faster total brain atrophy and consistently low plasma Aβ1-42/Aβ1-40 ratio, with increased total brain atrophy and gray matter atrophy. In dementia-free older adults, high plasma Aβ1-40 and low plasma Aβ1-42/Aβ1-40 ratio were associated with smaller hippocampal volume and accelerated global and regional brain atrophy respectively.

Astrocyte reactivity and neuroinflammation are hallmarks of CNS pathological conditions such as Alzheimer's disease. However, the specific role of reactive astrocytes is still debated. This controversy may stem from the fact that most strategies used to modulate astrocyte reactivity and explore its contribution to disease outcomes have only limited specificity. Moreover, reactive astrocytes are now emerging as heterogeneous cells and all types of astrocyte reactivity may not be controlled efficiently by such strategies.Here, we used cell type-specific approaches in vivo and identified the JAK2-STAT3 pathway, as necessary and sufficient for the induction and maintenance of astrocyte reactivity. Modulation of this cascade by viral gene transfer in mouse astrocytes efficiently controlled several morphological and molecular features of reactivity. Inhibition of this pathway in mouse models of Alzheimer's disease improved three key pathological hallmarks by reducing amyloid deposition, improving spatial learning and restoring synaptic deficits.In conclusion, the JAK2-STAT3 cascade operates as a master regulator of astrocyte reactivity in vivo. Its inhibition offers new therapeutic opportunities for Alzheimer's disease.

To identify common genetic variants associated with the presence of brain microbleeds (BMBs).

Human longevity is heritable, but genome-wide association (GWA) studies have had limited success. Here, we perform two meta-analyses of GWA studies of a rigorous longevity phenotype definition including 11,262/3484 cases surviving at or beyond the age corresponding to the 90th/99th survival percentile, respectively, and 25,483 controls whose age at death or at last contact was at or below the age corresponding to the 60th survival percentile. Consistent with previous reports, rs429358 (apolipoprotein E (ApoE) ε4) is associated with lower odds of surviving to the 90th and 99th percentile age, while rs7412 (ApoE ε2) shows the opposite. Moreover, rs7676745, located near GPR78, associates with lower odds of surviving to the 90th percentile age. Gene-level association analysis reveals a role for tissue-specific expression of multiple genes in longevity. Finally, genetic correlation of the longevity GWA results with that of several disease-related phenotypes points to a shared genetic architecture between health and longevity.

The human RNA helicase DDX6 is an essential component of membrane-less organelles called processing bodies (PBs). PBs are involved in mRNA metabolic processes including translational repression via coordinated storage of mRNAs. Previous studies in human cell lines have implicated altered DDX6 in molecular and cellular dysfunction, but clinical consequences and pathogenesis in humans have yet to be described. Here, we report the identification of five rare de novo missense variants in DDX6 in probands presenting with intellectual disability, developmental delay, and similar dysmorphic features including telecanthus, epicanthus, arched eyebrows, and low-set ears. All five missense variants (p.His372Arg, p.Arg373Gln, p.Cys390Arg, p.Thr391Ile, and p.Thr391Pro) are located in two conserved motifs of the RecA-2 domain of DDX6 involved in RNA binding, helicase activity, and protein-partner binding. We use functional studies to demonstrate that the first variants identified (p.Arg373Gln and p.Cys390Arg) cause significant defects in PB assembly in primary fibroblast and model human cell lines. These variants' interactions with several protein partners were also disrupted in immunoprecipitation assays. Further investigation via complementation assays included the additional variants p.Thr391Ile and p.Thr391Pro, both of which, similarly to p.Arg373Gln and p.Cys390Arg, demonstrated significant defects in P-body assembly. Complementing these molecular findings, modeling of the variants on solved protein structures showed distinct spatial clustering near known protein binding regions. Collectively, our clinical and molecular data describe a neurodevelopmental syndrome associated with pathogenic missense variants in DDX6. Additionally, we suggest DDX6 join the DExD/H-box genes DDX3X and DHX30 in an emerging class of neurodevelopmental disorders involving RNA helicases.

Leprosy is a chronic infectious disease of the skin and peripheral nerves with a strong genetic predisposition. Recent genome-wide approaches have identified numerous common variants associated with leprosy, almost all in the Chinese population. We conducted the first family-based genome-wide association study of leprosy in 622 affected offspring from Vietnam, followed by replication in an independent sample of 1181 leprosy cases and 668 controls of the same ethnic origin. The most significant results were observed within the HLA region, in which six SNPs displayed genome-wide significant associations, all of which were replicated in the independent case/control sample. We investigated the signal in the HLA region in more detail, by conducting a multivariate analysis on the case/control sample of 319 GWAS-suggestive HLA hits for which evidence for replication was obtained. We identified three independently associated SNPs, two located in the HLA class I region (rs1265048: OR = 0.69 [0.58-0.80], combined p-value = 5.53x10-11; and rs114598080: OR = 1.47 [1.46-1.48], combined p-value = 8.77x10-13), and one located in the HLA class II region (rs3187964 (OR = 1.67 [1.55-1.80], combined p-value = 8.35x10-16). We also validated two previously identified risk factors for leprosy: the missense variant rs3764147 in the LACC1 gene (OR = 1.52 [1.41-1.63], combined p-value = 5.06x10-14), and the intergenic variant rs6871626 located close to the IL12B gene (OR = 0.73 [0.61-0.84], combined p-value = 6.44x10-8). These results shed new light on the genetic control of leprosy, by dissecting the influence of HLA SNPs, and validating the independent role of two additional variants in a large Vietnamese sample.

New Guineans represent one of the oldest locally continuous populations outside Africa, harboring among the greatest linguistic and genetic diversity on the planet. Archeological and genetic evidence suggest that their ancestors reached Sahul (present day New Guinea and Australia) by at least 55,000 years ago (kya). However, little is known about this early settlement phase or subsequent dispersal and population structuring over the subsequent period of time. Here we report 379 complete Papuan mitochondrial genomes from across Papua New Guinea, which allow us to reconstruct the phylogenetic and phylogeographic history of northern Sahul. Our results support the arrival of two groups of settlers in Sahul within the same broad time window (50-65 kya), each carrying a different set of maternal lineages and settling Northern and Southern Sahul separately. Strong geographic structure in northern Sahul remains visible today, indicating limited dispersal over time despite major climatic, cultural, and historical changes. However, following a period of isolation lasting nearly 20 ky after initial settlement, environmental changes postdating the Last Glacial Maximum stimulated diversification of mtDNA lineages and greater interactions within and beyond Northern Sahul, to Southern Sahul, Wallacea and beyond. Later, in the Holocene, populations from New Guinea, in contrast to those of Australia, participated in early interactions with incoming Asian populations from Island Southeast Asia and continuing into Oceania.

Mitochondrial DNA (mtDNA) is present in multiple copies in human cells. We evaluated cross-sectional associations of whole blood mtDNA copy number (CN) with several cardiometabolic disease traits in 408,361 participants of multiple ancestries in TOPMed and UK Biobank. Age showed a threshold association with mtDNA CN: among younger participants (<65 years of age), each additional 10 years of age was associated with 0.03 standard deviation (s.d.) higher level of mtDNA CN (P = 0.0014) versus a 0.14 s.d. lower level of mtDNA CN (P = 1.82 × 10-13) among older participants (≥65 years). At lower mtDNA CN levels, we found age-independent associations with increased odds of obesity (P = 5.6 × 10-238), hypertension (P = 2.8 × 10-50), diabetes (P = 3.6 × 10-7), and hyperlipidemia (P = 6.3 × 10-5). The observed decline in mtDNA CN after 65 years of age may be a key to understanding age-related diseases.

Cerebrovascular reactivity (CVR), which measures the ability of cerebral blood vessels to dilate or constrict in response to vasoactive stimuli such as CO2 inhalation, is an important index of the brain's vascular health. Quantification of CVR using BOLD MRI with hypercapnia challenge has shown great promises in research and clinical studies. However, in order for it to be used as a potential imaging biomarker in large-scale and multi-site studies, the reliability of CO2-CVR quantification across different MRI acquisition platforms and researchers/raters must be examined. The goal of this report from the MarkVCID small vessel disease biomarkers consortium is to evaluate the reliability of CO2-CVR quantification in three studies. First, the inter-rater reliability of CO2-CVR data processing was evaluated by having raters from 5 MarkVCID sites process the same 30 CVR datasets using a cloud-based CVR data processing pipeline. Second, the inter-scanner reproducibility of CO2-CVR quantification was assessed in 10 young subjects across two scanners of different vendors. Third, test-retest repeatability was evaluated in 20 elderly subjects from 4 sites with a scan interval of less than 2 weeks. In all studies, the CO2 CVR measurements were performed using the fixed inspiration method, where the subjects wore a nose clip and a mouthpiece and breathed room air and 5% CO2 air contained in a Douglas bag alternatively through their mouth. The results showed that the inter-rater CoV of CVR processing was 0.08 ± 0.08% for whole-brain CVR values and ranged from 0.16% to 0.88% in major brain regions, with ICC of absolute agreement above 0.9959 for all brain regions. Inter-scanner CoV was found to be 6.90 ± 5.08% for whole-brain CVR values, and ranged from 4.69% to 12.71% in major brain regions, which are comparable to intra-session CoVs obtained from the same scanners on the same day. ICC of consistency between the two scanners was 0.8498 for whole-brain CVR and ranged from 0.8052 to 0.9185 across major brain regions. In the test-retest evaluation, test-retest CoV across different days was found to be 18.29 ± 17.12% for whole-brain CVR values, and ranged from 16.58% to 19.52% in major brain regions, with ICC of absolute agreement ranged from 0.6480 to 0.7785. These results demonstrated good inter-rater, inter-scanner, and test-retest reliability in healthy volunteers, and suggested that CO2-CVR has suitable instrumental properties for use as an imaging biomarker of cerebrovascular function in multi-site and longitudinal observational studies and clinical trials.

Autosomal Dominant Hypercholesterolemia (ADH) is a genetic disorder caused by pathogenic variants in LDLR, APOB, PCSK9 and APOE genes. We sought to identify new candidate genes responsible for the ADH phenotype in patients without pathogenic variants in the known ADH-causing genes by focusing on a French family with affected and non-affected members who presented a high ADH polygenic risk score (wPRS). Linkage analysis, whole exome and whole genome sequencing resulted in the identification of variants p.(Pro398Ala) in CYP7A1, p.(Val1382Phe) in LRP6 and p.(Ser202His) in LDLRAP1. A total of 6 other variants were identified in 6 of 160 unrelated ADH probands: p.(Ala13Val) and p.(Aps347Asn) in CYP7A1; p.(Tyr972Cys), p.(Thr1479Ile) and p.(Ser1612Phe) in LRP6; and p.(Ser202LeufsTer19) in LDLRAP1. All six probands presented a moderate wPRS. Serum analyses of carriers of the p.(Pro398Ala) variant in CYP7A1 showed no differences in the synthesis of bile acids compared to the serums of non-carriers. Functional studies of the four LRP6 mutants in HEK293T cells resulted in contradictory results excluding a major effect of each variant alone. Within the family, none of the heterozygous for only the LDLRAP1 p.(Ser202His) variant presented ADH. Altogether, each variant individually does not result in elevated LDL-C; however, the oligogenic combination of two or three variants reveals the ADH phenotype.