Case report: Novel compound heterozygous IL1RN mutations as the likely cause of a lethal form of deficiency of interleukin-1 receptor antagonist

Undiagnosed monogenic diseases represent a challenging group of human conditions highly suspicious to have a genetic origin, but without conclusive evidences about it. We identified two brothers born prematurely from a non-consanguineous healthy couple, with a neonatal-onset, chronic disease characterized by severe skin and bone inflammatory manifestations and a fatal outcome in infancy. We conducted DNA and mRNA analyses in the patients' healthy relatives to identify the genetic cause of the patients' disease. DNA analyses were performed by both Sanger and next-generation sequencing, which identified two novel heterozygous IL1RN variants: the intronic c.318 + 2T>G variant in the father and a ≈2,600-bp intragenic deletion in the mother. IL1RN mRNA production was markedly decreased in both progenitors when compared with healthy subjects. The mRNA sequencing performed in each parent identified two novel, truncated IL1RN transcripts. Additional experiments revealed a perfect intrafamilial phenotype-genotype segregation following an autosomal recessive inheritance pattern. The evidences shown here supported for the presence of two novel loss-of-function (LoF) IL1RN pathogenic variants in the analyzed family. Biallelic LoF variants at the IL1RN gene cause the deficiency of interleukin-1 receptor antagonist (DIRA), a monogenic autoinflammatory disease with marked similarities with the patients described here. Despite the non-availability of the patients' samples representing the main limitation of this study, the collected evidences strongly suggest that the patients described here suffered from a lethal form of DIRA likely due to a compound heterozygous genotype at IL1RN, thus providing a reliable genetic diagnosis based on the integration of old medical information with currently obtained genetic data.

Topoisomerase 1 facilitates nucleosome reassembly at stress genes during recovery

Chromatin remodeling is essential to allow full development of alternative gene expression programs in response to environmental changes. In fission yeast, oxidative stress triggers massive transcriptional changes including the activation of hundreds of genes, with the participation of histone modifying complexes and chromatin remodelers. DNA transcription is associated to alterations in DNA topology, and DNA topoisomerases facilitate elongation along gene bodies. Here, we test whether the DNA topoisomerase Top1 participates in the RNA polymerase II-dependent activation of the cellular response to oxidative stress. Cells lacking Top1 are resistant to H2O2 stress. The transcriptome of Δtop1 strain was not greatly affected in the absence of stress, but activation of the anti-stress gene expression program was more sustained than in wild-type cells. Top1 associated to stress open reading frames. While the nucleosomes of stress genes are partially and transiently evicted during stress, the chromatin configuration remains open for longer times in cells lacking Top1, facilitating RNA polymerase II progression. We propose that, by removing DNA tension arising from transcription, Top1 facilitates nucleosome reassembly and works in synergy with the chromatin remodeler Hrp1 as opposing forces to transcription and to Snf22 / Hrp3 opening remodelers.

Spanish cohort of VEXAS syndrome: clinical manifestations, outcome of treatments and novel evidences about UBA1 mosaicism

BACKGROUND

The vacuoles, E1-enzyme, X linked, autoinflammatory and somatic (VEXAS) syndrome is an adult-onset autoinflammatory disease (AID) due to postzygotic UBA1 variants.

OBJECTIVES

To investigate the presence of VEXAS syndrome among patients with adult-onset undiagnosed AID. Additional studies evaluated the mosaicism distribution and the circulating cytokines.

METHODS

Gene analyses were performed by both Sanger and amplicon-based deep sequencing. Patients' data were collected from their medical charts. Cytokines were quantified by Luminex.

RESULTS

Genetic analyses of enrolled patients (n=42) identified 30 patients carrying UBA1 pathogenic variants, with frequencies compatible for postzygotic variants. All patients were male individuals who presented with a late-onset disease (mean 67.5 years; median 67.0 years) characterised by cutaneous lesions (90%), fever (66.7%), pulmonary manifestations (66.7%) and arthritis (53.3%). Macrocytic anaemia and increased erythrocyte sedimentation rate and ferritin were the most relevant analytical abnormalities. Glucocorticoids ameliorated the inflammatory manifestations, but most patients became glucocorticoid-dependent. Positive responses were obtained when targeting the haematopoietic component of the disease with either decitabine or allogeneic haematopoietic stem cell transplantation. Additional analyses detected the UBA1 variants in both haematopoietic and non-haematopoietic tissues. Finally, analysis of circulating cytokines did not identify inflammatory mediators of the disease.

CONCLUSION

Thirty patients with adult-onset AID were definitively diagnosed with VEXAS syndrome through genetic analyses. Despite minor interindividual differences, their main characteristics were in concordance with previous reports. We detected for the first time the UBA1 mosaicism in non-haematopoietic tissue, which questions the previous concept of myeloid-restricted mosaicism and may have conceptual consequences for the disease mechanisms.

Population analysis of complete mitogenomes for 334 samples from El Salvador

The use of mitochondrial DNA (mtDNA) in the field of forensic genetics is widely spread mainly due to its advantages when identifying highly degraded samples. In this sense, massive parallel sequencing has made the analysis of the whole mitogenome more accessible, noticeably increasing the informativeness of mtDNA haplotypes. The civil war (1980-1992) in El Salvador caused many deaths and disappearances (including children) all across the country and the economic and social instability after the war forced many people to emigration. For this reason, different organizations have collected DNA samples from relatives with the aim of identifying missing people. Thus, we present a dataset containing 334 complete mitogenomes from the Salvadoran general population. To the best of our knowledge, this is the first publication of a nationwide forensic-quality complete mitogenome database of any Latin American country. We found 293 different haplotypes, with a random match probability of 0.0041 and 26.6 mean pairwise differences, which is similar to other Latin American populations, and which represent a marked improvement from the values obtained with just control region sequences. These haplotypes belong to 54 different haplogroups, being 91% of them of Native American origin. Over a third (35.9%) of the individuals carried at least a heteroplasmic site (excluding length heteroplasmies). Ultimately, the present database aims to represent mtDNA haplotype diversity in the general Salvadoran populations as a basis for the identification of people that disappeared during or after the civil war.

Somatic genetic variation in healthy tissue and non-cancer diseases

Somatic genetic variants have been studied for several years mostly concerning cancer, where they contribute to its origin and development. It is also clear that the somatic variants load is greater in aged individuals in comparison to younger ones, pointing to a cause/consequence of the senescence process. More recently, researchers have focused on the role of this type of variation in healthy tissue and its dynamics in cell lineages and different organs. In addition, somatic variants have been described to contribute to monogenic diseases, and the number of evidences of their role in complex disorders is also increasing. Thanks to recent advances in next-generation sequencing technologies, this type of genetic variation can be now more easily studied than in the past, although we still face some important limitations. Novel strategies for sampling, sequencing and filtering are being investigated to detect these variants, although validating them with an orthogonal approach will most likely still be needed. In this review, we aim to update our knowledge of somatic variation detection and its relation to healthy tissue and non-cancer diseases.

Healthcare workers' perception of the usefulness of a healthcare ethics consultation service

BACKGROUND AND OBJECTIVE

Clinical ethics consultation services (CEC) have arisen from Healthcare Ethics Committees (HEC) to address ethical conflicts in real-time. Our aim was to determine the perception of usefulness of a CEC service among healthcare workers one year after its creation as well as to assess changes in trends in the use of the CEC and HEC between 2015 and 2021.

MATERIALS AND METHODS

This observational, cross-sectional study was based on a standardized survey of healthcare workers at an urban tertiary care hospital. The results were also compared to those from an identical survey conducted in the same population in 2015.

RESULTS

A total of 213 professionals participated (mean age 44 ± 11 years, 69% women). The professionals were more familiar with the HEC than the CEC service (94% vs 61%; p < 0.001). Forty-five individuals (21%) had consulted the CEC since its implementation; 95% of them found the consultation useful. Physicians knew about and used the CEC more than other groups of professionals. The degree of knowledge of the HEC increased significantly by 2021 compared to 2015 (94% v. 76%; p < 0.001). Some areas for improvement identified were the need for greater dissemination of the service, guaranteeing institutional resources to maintain the service, and encouraging greater participation from different professional groups.

CONCLUSIONS

Knowledge of the institutional HEC and CEC services has increased in recent years among healthcare workers, who considered the CEC service to be useful for addressing ethical conflicts in daily practice.

Assessing the digenic model in rare disorders using population sequencing data

An important fraction of patients with rare disorders remains with no clear genetic diagnostic, even after whole-exome or whole-genome sequencing, posing a difficulty in giving adequate treatment and genetic counseling. The analysis of genomic data in rare disorders mostly considers the presence of single gene variants in coding regions that follow a concrete monogenic mode of inheritance. A digenic inheritance, with variants in two functionally-related genes in the same individual, is a plausible alternative that might explain the genetic basis of the disease in some cases. In this case, digenic disease combinations should be absent or underrepresented in healthy individuals. We develop a framework to evaluate the significance of digenic combinations and test its statistical power in different scenarios. We suggest that this approach will be relevant with the advent of new sequencing efforts including hundreds of thousands of samples.

A transient inflammatory response contributes to oxaliplatin neurotoxicity in mice

OBJECTIVES

Peripheral neuropathy is a relevant dose-limiting adverse event that can affect up to 90% of oncologic patients with colorectal cancer receiving oxaliplatin treatment. The severity of neurotoxicity often leads to dose reduction or even premature cessation of chemotherapy. Unfortunately, the limited knowledge about the molecular mechanisms related to oxaliplatin neurotoxicity leads to a lack of effective treatments to prevent the development of this clinical condition. In this context, the present work aimed to determine the exact molecular mechanisms involved in the development of oxaliplatin neurotoxicity in a murine model to try to find new therapeutical targets.

METHODS

By single-cell RNA sequencing (scRNA-seq), we studied the transcriptomic profile of sensory neurons and satellite glial cells (SGC) of the Dorsal Root Ganglia (DRG) from a well-characterized mouse model of oxaliplatin neurotoxicity.

RESULTS

Analysis of scRNA-seq data pointed to modulation of inflammatory processes in response to oxaliplatin treatment. In this line, we observed increased levels of NF-kB p65 protein, pro-inflammatory cytokines, and immune cell infiltration in DRGs and peripheral nerves of oxaliplatin-treated mice, which was accompanied by mechanical allodynia and decrease in sensory nerve amplitudes.

INTERPRETATION

Our data show that, in addition to the well-described DNA damage, oxaliplatin neurotoxicity is related to an exacerbated pro-inflammatory response in DRG and peripheral nerves, and open new insights in the development of anti-inflammatory strategies as a treatment for preventing peripheral neuropathy induced by oxaliplatin.

Somatic Mutations Detected in Parkinson Disease Could Affect Genes With a Role in Synaptic and Neuronal Processes

The role of somatic mutations in complex diseases, including neurodevelopmental and neurodegenerative disorders, is becoming increasingly clear. However, to date, no study has shown their relation to Parkinson disease’s phenotype. To explore the relevance of embryonic somatic mutations in sporadic Parkinson disease, we performed whole-exome sequencing in blood and four brain regions of ten patients. We identified 59 candidate somatic single nucleotide variants (sSNVs) through sensitive calling and a careful filtering strategy (COSMOS). We validated 27 of them with amplicon-based ultra-deep sequencing, with a 70% validation rate for the highest-confidence variants. The identified sSNVs are in genes with synaptic functions that are co-expressed with genes previously associated with Parkinson disease. Most of the sSNVs were only called in blood but were also found in the brain tissues with ultra-deep amplicon sequencing, demonstrating the strength of multi-tissue sampling designs.

A forensic population database in El Salvador: 58 STRs and 94 SNPs

We have genotyped the 58 STRs (27 autosomal, 24 Y-STRs and 7 X-STRs) and 94 autosomal SNPs in Illumina ForenSeq™ Primer Mix A in a sample of 248 men and 143 women from El Salvador, Central America. Regional division (Centro, Oriente, Occidente) showed in almost all cases F_ST values not significantly different from 0, and further analyses were applied only to the undivided, country-wide population. The overall random match probability (RMP) decreased from 6.79 × 10^-31 in length-based genotypes in the 27 autosomal STRs to 1.47 × 10^-34 in repeat-sequence based genotypes. Combining the autosomal loci in this set, RMP reaches 2.97 × 10^-70. In a population genetic analysis, El Salvador showed the lowest F_ST values with US Hispanics both for autosomal and X-STRs; however, it was much closer to Native Americans for the latter than for the former, in accordance with the well-known gender-biased admixture that created most Latin American populations.

First Description of Late-Onset Autoinflammatory Disease Due to Somatic NLRC4 Mosaicism

INTRODUCTION

Autoinflammatory diseases (AIDs) are inherited disorders of innate immunity that usually start during childhood. However, several studies have recently reported an increasing number of patients with AID starting in adulthood.

OBJECTIVE

This study was undertaken to characterize the cause underlying a patient with late-onset uncharacterized AID.

METHODS

Genetic studies were performed using Sanger sequencing and next-generation sequencing (NGS) methods. In silico, in vitro and ex vivo analyses were performed to determine the functional consequences of the detected variant.

RESULTS

We studied a 57 years-old woman who started at the age of 47 years with recurrent episodes of fever, myalgias, arthralgias, diffuse abdominal pain, diarrhea, adenopathies and systemic inflammation, which were relatively well controlled with anti-IL-1 drugs. NGS analyses did not detect germline variants in any of the known AID-associated genes, but they identified the p.Ser171Phe NLRC4 variant in unfractionated blood, with an allele fraction (2-4%) compatible with gene mosaicism. Structural modeling analyses suggest that this missense variant might favor the open, active conformation of the NLRC4 protein, and in vitro and ex vivo analyses confirmed its trend to oligomerize and activate the NLRC4-inflammasome, with a subsequent IL-18 overproduction.

CONCLUSION

We have identified the post-zygotic p.Ser171Phe NLRC4 variant as a plausible cause of the disease in the enrolled patient. Functional and structural studies clearly support for the first time its gain-of-function behavior consistent with previously reported NLRC4 pathogenic variants. These novel evidences should be considered in the diagnostic workup of patients with uncharacterized AID starting during adulthood.

Second GHEP-ISFG exercise for DVI: "DNA-led" victims' identification in a simulated air crash

The Spanish and Portuguese-Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) has organized a second collaborative exercise on a simulated case of Disaster Victim Identification (DVI), with the participation of eighteen laboratories. The exercise focused on the analysis of a simulated plane crash case of medium-size resulting in 66 victims with varying degrees of fragmentation of the bodies (with commingled remains). As an additional difficulty, this second exercise included 21 related victims belonging to 6 families among the 66 missings to be identified. A total number of 228 post-mortem samples were represented with aSTR and mtDNA profiles, with a proportion of partial aSTR profiles simulating charred remains. To perform the exercise, participants were provided with aSTR and mtDNA data of 51 reference pedigrees -some of which deficient-including 128 donors for identification purposes. The exercise consisted firstly in the comparison of the post-mortem genetic profiles in order to re-associate fragmented remains to the same individual and secondly in the identification of the re-associated remains by comparing aSTR and mtDNA profiles with reference pedigrees using pre-established thresholds to report a positive identification. Regarding the results of the post-mortem samples re-associations, only a small number of discrepancies among participants were detected, all of which were from just a few labs. However, in the identification process by kinship analysis with family references, there were more discrepancies in comparison to the correct results. The identification results of single victims yielded fewer problems than the identification of multiple related victims within the same family groups. Several reasons for the discrepant results were detected: a) the identity/non-identity hypotheses were sometimes wrongly expressed in the likelihood ratio calculations, b) some laboratories failed to use all family references to report the DNA match, c) In families with several related victims, some laboratories firstly identified some victims and then unnecessarily used their genetic information to identify the remaining victims within the family, d) some laboratories did not correctly use "prior odds" values for the Bayesian treatment of the episode for both post-mortem/post-mortem re-associations as well as the ante-mortem/post-mortem comparisons to evaluate the probability of identity. For some of the above reasons, certain laboratories failed to identify some victims. This simulated "DNA-led" identification exercise may help forensic genetic laboratories to gain experience and expertize for DVI or MPI in using genetic data and comparing their own results with the ones in this collaborative exercise.

Genome-wide postnatal changes in immunity following fetal inflammatory response

The fetal inflammatory response (FIR) increases the risk of perinatal brain injury, particularly in extremely low gestational age newborns (ELGANs, < 28 weeks of gestation). One of the mechanisms contributing to such a risk is a postnatal intermittent or sustained systemic inflammation (ISSI) following FIR. The link between prenatal and postnatal systemic inflammation is supported by the presence of well‐established inflammatory biomarkers in the umbilical cord and peripheral blood. However, the extent of molecular changes contributing to this association is unknown. Using RNA sequencing and mass spectrometry proteomics, we profiled the transcriptome and proteome of archived neonatal dried blood spot (DBS) specimens from 21 ELGANs. Comparing FIR‐affected and unaffected ELGANs, we identified 782 gene and 27 protein expression changes of 50% magnitude or more, and an experiment‐wide significance level below 5% false discovery rate. These expression changes confirm the robust postnatal activation of the innate immune system in FIR‐affected ELGANs and reveal for the first time an impairment of their adaptive immunity. In turn, the altered pathways provide clues about the molecular mechanisms triggering ISSI after FIR, and the onset of perinatal brain injury.

Databases

EGAS00001003635 (EGA); PXD011626 (PRIDE).

Severe Autoinflammatory Manifestations and Antibody Deficiency Due to Novel Hypermorphic PLCG2 Mutations

Autoinflammatory diseases (AIDs) were first described as clinical disorders characterized by recurrent episodes of seemingly unprovoked sterile inflammation. In the past few years, the identification of novel AIDs expanded their phenotypes toward more complex clinical pictures associating vasculopathy, autoimmunity, or immunodeficiency. Herein, we describe two unrelated patients suffering since the neonatal period from a complex disease mainly characterized by severe sterile inflammation, recurrent bacterial infections, and marked humoral immunodeficiency. Whole-exome sequencing detected a novel, de novo heterozygous PLCG2 variant in each patient (p.Ala708Pro and p.Leu845_Leu848del). A clear enhanced PLCγ2 activity for both variants was demonstrated by both ex vivo calcium responses of the patient's B cells to IgM stimulation and in vitro assessment of PLC activity. These data supported the autoinflammation and PLCγ2-associated antibody deficiency and immune dysregulation (APLAID) diagnosis in both patients. Immunological evaluation revealed a severe decrease of immunoglobulins and B cells, especially class-switched memory B cells, with normal T and NK cell counts. Analysis of bone marrow of one patient revealed a reduced immature B cell fraction compared with controls. Additional investigations showed that both PLCG2 variants activate the NLRP3-inflammasome through the alternative pathway instead of the canonical pathway. Collectively, the evidences here shown expand APLAID diversity toward more severe phenotypes than previously reported including dominantly inherited agammaglobulinemia, add novel data about its genetic basis, and implicate the alternative NLRP3-inflammasome activation pathway in the basis of sterile inflammation.

The first GHEP-ISFG collaborative exercise on forensic applications of massively parallel sequencing

One of the main goals of the Spanish and Portuguese-Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) is to promote and contribute to the development and dissemination of scientific knowledge in the field of forensic genetics. The GHEP-ISFG supports several Working Commissions which develop different scientific activities. One of them, the Working Commission on "Massively Parallel Sequencing (MPS): Forensic Applications", organized its first collaborative exercise on forensic applications of MPS technology in 2019. The aim of this exercise was to assess the concordance between the MPS results and those obtained with conventional technologies (capillary electrophoresis and Sanger sequencing), as well as to compare the results obtained within the different MPS platforms and/or the different kits/panels and analysis software packages (commercial and open-access) available on the market. The seven participating laboratories analyzed some samples of the annual GHEP-ISFG proficiency test (EIADN No. 27 (2019)), using Ion Torrent™ or MiSeq FGx® platforms. Six of them sent autosomal STR sequence data, five laboratories performed MPS analysis of individual identification SNPs, four laboratories reported MPS data of Y-chromosomal STRs, and X-chromosomal STRs, three laboratories performed MPS analysis of ancestry informative SNPs and phenotype informative SNPs, two labs performed MPS analysis of the mitochondrial DNA control region, and only one lab produced MPS data of lineage informative SNPs. Autosomal STR sequencing results were highly concordant to the consensus obtained by capillary electrophoresis in the EIADN No. 27 (2019) exercise. Furthermore, in general, a high level of concordance was observed between the results of the participating laboratories, regardless of the platform used. The main discordances were due to errors during the analysis process or from sequence data obtained with low depth of coverage. In this paper we highlight some issues that still arise, such as standardization of the nomenclature for STRs analyzed by sequencing with MPS, the universal uptake of a nomenclature framework by the analysis software, and well established validation and accreditation of the new MPS platforms for use in routine forensic case-work.

Historical human remains identification through maternal and paternal genetic signatures in a founder population with extensive genealogical record

OBJECTIVES

We describe a method to identify human remains excavated from unmarked graves in historical Québec cemeteries by combining parental-lineage genetic markers with the whole-population genealogy of Québec contained in the BALSAC database.

MATERIALS AND METHODS

The remains of six men were exhumed from four historical cemeteries in the province of Québec, Canada. DNA was extracted from the remains and genotyped to reveal their mitochondrial and Y-chromosome haplotypes, which were compared to a collection of haplotypes of genealogically-anchored modern volunteers. Maternal and paternal genealogies were searched in the BALSAC genealogical record for parental couples matching the mitochondrial and the Y-chromosome haplotypic signatures, to identify candidate sons from whom the remains could have originated.

RESULTS

Analysis of the matching genealogies identified the parents of one man inhumed in the cemetery of the investigated parish during its operating time. The candidate individual died in 1833 at the age of 58, a plausible age at death in light of osteological analysis of the remains.

DISCUSSION

This study demonstrates the promising potential of coupling genetic information from living individuals to genealogical data in BALSAC to identify historical human remains. If genetic coverage is increased, the genealogical information in BALSAC could enable the identification of 87% of the men (n = 178,435) married in Québec before 1850, with high discriminatory power in most cases since >75% of the parental couples have unique biparental signatures in most regions. Genotyping and identifying Québec's historical human remains are a key to reconstructing the genomes of the founders of Québec and reinhuming archeological remains with a marked grave.

FHLdb: A Comprehensive Database on the Molecular Basis of Familial Hemophagocytic Lymphohistiocytosis

Background: Primary immunodeficiencies (PIDs) are a heterogeneous group of disorders. The lack of comprehensive disease-specific mutation databases may hinder or delay classification of the genetic variants found in samples from these patients. This is especially true for familial hemophagocytic lymphohistiocytosis (FHL), a life-threatening PID classically considered an autosomal recessive condition, but with increasingly demonstrated genetic heterogeneity.

Objective: The aim of this study was to build an open-access repository to collect detailed information on the known genetic variants reported in FHL.

Methods: We manually reviewed more than 120 articles to identify all reported variants related to FHL. We retrieved relevant information about the allelic status, the number of patients with the same variant, and whether functional assays were done. We stored all the data retrieved in a PostgreSQL database and then built a website on top of it, using the Django framework.

Results: The database designed (FHLdb) (https://www.biotoclin.org/FHLdb) contains comprehensive information on reported variants in the 4 genes related to FHL (PRF1, UNC13D, STXBP2, STX11). It comprises 240 missense, 69 frameshift, 51 nonsense, 51 splicing, 10 in-frame indel, 7 deep intronic, and 5 large rearrangement variants together with their allelic status, carrier(s) information, and functional evidence. All genetic variants have been classified as pathogenic, likely pathogenic, uncertain significance, likely benign or benign, according to the American College of Medical Genetics guidelines. Additionally, it integrates information from other relevant databases: clinical evidence from ClinVar and UniProt, population allele frequency from ExAC and gnomAD, and pathogenicity predictions from well-recognized tools (e.g., PolyPhen-2, SIFT). Finally, a diagram depicts the location of the variant relative to the gene exon and protein domain structures.

Conclusion: FHLdb includes a broad range of data on the reported genetic variants in familial HLH genes. It is a free-access and easy-to-use resource that will facilitate the interpretation of molecular results of FHL patients, and it illustrates the potential value of disease-specific databases for other PIDs.

Flow Sorting Enrichment and Nanopore Sequencing of Chromosome 1 From a Chinese Individual

Sorting of individual chromosomes by Flow Cytometry (flow-sorting) is an enrichment method to potentially simplify genome assembly by isolating chromosomes from the context of the genome. We have recently developed a workflow to sequence native, unamplified DNA and applied it to the smallest human chromosome, the Y chromosome. Here, we modify improve upon that workflow to increase DNA recovery from chromosome sorting as well as sequencing yield. We apply it to sequence and assemble the largest human chromosome - chromosome 1 - of a Chinese individual using a single Oxford Nanopore MinION flow cell. We generate a selective and highly continuous assembly whose continuity reaches into the order of magnitude of the human reference GRCh38. We then use this assembly to call candidate structural variants against the reference and find 685 putative novel SV candidates. We propose this workflow as a potential solution to assemble structurally complex chromosomes, or the study of very large plant or animal genomes that might challenge traditional assembly strategies.

Phenotype, penetrance, and treatment of 133 cytotoxic T-lymphocyte antigen 4-insufficient subjects

BACKGROUND

Cytotoxic T-lymphocyte antigen 4 (CTLA-4) is a negative immune regulator. Heterozygous CTLA4 germline mutations can cause a complex immune dysregulation syndrome in human subjects.

OBJECTIVE

We sought to characterize the penetrance, clinical features, and best treatment options in 133 CTLA4 mutation carriers.

METHODS

Genetics, clinical features, laboratory values, and outcomes of treatment options were assessed in a worldwide cohort of CTLA4 mutation carriers.

RESULTS

We identified 133 subjects from 54 unrelated families carrying 45 different heterozygous CTLA4 mutations, including 28 previously undescribed mutations. Ninety mutation carriers were considered affected, suggesting a clinical penetrance of at least 67%; median age of onset was 11 years, and the mortality rate within affected mutation carriers was 16% (n = 15). Main clinical manifestations included hypogammaglobulinemia (84%), lymphoproliferation (73%), autoimmune cytopenia (62%), and respiratory (68%), gastrointestinal (59%), or neurological features (29%). Eight affected mutation carriers had lymphoma, and 3 had gastric cancer. An EBV association was found in 6 patients with malignancies. CTLA4 mutations were associated with lymphopenia and decreased T-, B-, and natural killer (NK) cell counts. Successful targeted therapies included application of CTLA-4 fusion proteins, mechanistic target of rapamycin inhibitors, and hematopoietic stem cell transplantation. EBV reactivation occurred in 2 affected mutation carriers after immunosuppression.

CONCLUSIONS

Affected mutation carriers with CTLA-4 insufficiency can present in any medical specialty. Family members should be counseled because disease manifestation can occur as late as 50 years of age. EBV- and cytomegalovirus-associated complications must be closely monitored. Treatment interventions should be coordinated in clinical trials.

Genomes reveal marked differences in the adaptive evolution between orangutan species

BACKGROUND

Integrating demography and adaptive evolution is pivotal to understanding the evolutionary history and conservation of great apes. However, little is known about the adaptive evolution of our closest relatives, in particular if and to what extent adaptions to environmental differences have occurred. Here, we used whole-genome sequencing data from critically endangered orangutans from North Sumatra (Pongo abelii) and Borneo (P. pygmaeus) to investigate adaptive responses of each species to environmental differences during the Pleistocene.

RESULTS

Taking into account the markedly disparate demographic histories of each species after their split ~ 1 Ma ago, we show that persistent environmental differences on each island had a strong impact on the adaptive evolution of the genus Pongo. Across a range of tests for positive selection, we find a consistent pattern of between-island and species differences. In the more productive Sumatran environment, the most notable signals of positive selection involve genes linked to brain and neuronal development, learning, and glucose metabolism. On Borneo, however, positive selection comprised genes involved in lipid metabolism, as well as cardiac and muscle activities.

CONCLUSIONS

We find strikingly different sets of genes appearing to have evolved under strong positive selection in each species. In Sumatran orangutans, selection patterns were congruent with well-documented cognitive and behavioral differences between the species, such as a larger and more complex cultural repertoire and higher degrees of sociality. However, in Bornean orangutans, selective responses to fluctuating environmental conditions appear to have produced physiological adaptations to generally lower and temporally more unpredictable food supplies.

Unexpected relevant role of gene mosaicism in patients with primary immunodeficiency diseases

BACKGROUND

Postzygotic de novo mutations lead to the phenomenon of gene mosaicism. The 3 main types are called somatic, gonadal, and gonosomal mosaicism, which differ in terms of the body distribution of postzygotic mutations. Mosaicism has been reported occasionally in patients with primary immunodeficiency diseases (PIDs) since the early 1990s, but its real involvement has not been systematically addressed.

OBJECTIVE

We sought to investigate the incidence of gene mosaicism in patients with PIDs.

METHODS

The amplicon-based deep sequencing method was used in the 3 parts of the study that establish (1) the allele frequency of germline variants (n = 100), (2) the incidence of parental gonosomal mosaicism in families with PIDs with de novo mutations (n = 92), and (3) the incidence of mosaicism in families with PIDs with moderate-to-high suspicion of gene mosaicism (n = 36). Additional investigations evaluated body distribution of postzygotic mutations, their stability over time, and their characteristics.

RESULTS

The range of allele frequency (44.1% to 55.6%) was established for germline variants. Those with minor allele frequencies of less than 44.1% were assumed to be postzygotic. Mosaicism was detected in 30 (23.4%) of 128 families with PIDs, with a variable minor allele frequency (0.8% to 40.5%). Parental gonosomal mosaicism was detected in 6 (6.5%) of 92 families with de novo mutations, and a high incidence of mosaicism (63.9%) was detected among families with moderate-to-high suspicion of gene mosaicism. In most analyzed cases mosaicism was found to be both uniformly distributed and stable over time.

CONCLUSION

This study represents the largest performed to date to investigate mosaicism in patients with PIDs, revealing that it affects approximately 25% of enrolled families. Our results might have serious consequences regarding treatment and genetic counseling and reinforce the use of next-generation sequencing-based methods in the routine analyses of PIDs.

Evaluating the Genetics of Common Variable Immunodeficiency: Monogenetic Model and Beyond

Common variable immunodeficiency (CVID) is the most frequent symptomatic primary immunodeficiency characterized by recurrent infections, hypogammaglobulinemia and poor response to vaccines. Its diagnosis is made based on clinical and immunological criteria, after exclusion of other diseases that can cause similar phenotypes. Currently, less than 20% of cases of CVID have a known underlying genetic cause. We have analyzed whole-exome sequencing and copy number variants data of 36 children and adolescents diagnosed with CVID and healthy relatives to estimate the proportion of monogenic cases. We have replicated an association of CVID to p.C104R in TNFRSF13B and reported the second case of homozygous patient to date. Our results also identify five causative genetic variants in LRBA, CTLA4, NFKB1, and PIK3R1, as well as other very likely causative variants in PRKCD, MAPK8, or DOCK8 among others. We experimentally validate the effect of the LRBA stop-gain mutation which abolishes protein production and downregulates the expression of CTLA4, and of the frameshift indel in CTLA4 producing expression downregulation of the protein. Our results indicate a monogenic origin of at least 15–24% of the CVID cases included in the study. The proportion of monogenic patients seems to be lower in CVID than in other PID that have also been analyzed by whole exome or targeted gene panels sequencing. Regardless of the exact proportion of CVID monogenic cases, other genetic models have to be considered for CVID. We propose that because of its prevalence and other features as intermediate penetrancies and phenotypic variation within families, CVID could fit with other more complex genetic scenarios. In particular, in this work, we explore the possibility of CVID being originated by an oligogenic model with the presence of heterozygous mutations in interacting proteins or by the accumulation of detrimental variants in particular immunological pathways, as well as perform association tests to detect association with rare genetic functional variation in the CVID cohort compared to healthy controls.

The impact of endogenous content, replicates and pooling on genome capture from faecal samples

Target-capture approach has improved over the past years, proving to be very efficient tool for selectively sequencing genetic regions of interest. These methods have also allowed the use of noninvasive samples such as faeces (characterized by their low quantity and quality of endogenous DNA) to be used in conservation genomic, evolution and population genetic studies. Here we aim to test different protocols and strategies for exome capture using the Roche SeqCap EZ Developer kit (57.5 Mb). First, we captured a complex pool of DNA libraries. Second, we assessed the influence of using more than one faecal sample, extract and/or library from the same individual, to evaluate its effect on the molecular complexity of the experiment. We validated our experiments with 18 chimpanzee faecal samples collected from two field sites as a part of the Pan African Programme: The Cultured Chimpanzee. Those two field sites are in Kibale National Park, Uganda (N = 9) and Loango National Park, Gabon (N = 9). We demonstrate that at least 16 libraries can be pooled, target enriched through hybridization, and sequenced allowing for the genotyping of 951,949 exome markers for population genetic analyses. Further, we observe that molecule richness, and thus, data acquisition, increase when using multiple libraries from the same extract or multiple extracts from the same sample. Finally, repeated captures significantly decrease the proportion of off-target reads from 34.15% after one capture round to 7.83% after two capture rounds, supporting our conclusion that two rounds of target enrichment are advisable when using complex faecal samples.

Morphometric, Behavioral, and Genomic Evidence for a New Orangutan Species

Six extant species of non-human great apes are currently recognized: Sumatran and Bornean orangutans, eastern and western gorillas, and chimpanzees and bonobos [1]. However, large gaps remain in our knowledge of fine-scale variation in hominoid morphology, behavior, and genetics, and aspects of great ape taxonomy remain in flux. This is particularly true for orangutans (genus: Pongo), the only Asian great apes and phylogenetically our most distant relatives among extant hominids [1]. Designation of Bornean and Sumatran orangutans, P. pygmaeus (Linnaeus 1760) and P. abelii (Lesson 1827), as distinct species occurred in 2001 [1, 2]. Here, we show that an isolated population from Batang Toru, at the southernmost range limit of extant Sumatran orangutans south of Lake Toba, is distinct from other northern Sumatran and Bornean populations. By comparing cranio-mandibular and dental characters of an orangutan killed in a human-animal conflict to those of 33 adult male orangutans of a similar developmental stage, we found consistent differences between the Batang Toru individual and other extant Ponginae. Our analyses of 37 orangutan genomes provided a second line of evidence. Model-based approaches revealed that the deepest split in the evolutionary history of extant orangutans occurred ∼3.38 mya between the Batang Toru population and those to the north of Lake Toba, whereas both currently recognized species separated much later, about 674 kya. Our combined analyses support a new classification of orangutans into three extant species. The new species, Pongo tapanuliensis, encompasses the Batang Toru population, of which fewer than 800 individuals survive. VIDEO ABSTRACT.

Y-chromosomal sequences of diverse Indian populations and the ancestry of the Andamanese

We present 42 new Y-chromosomal sequences from diverse Indian tribal and non-tribal populations, including the Jarawa and Onge from the Andaman Islands, which are analysed within a calibrated Y-chromosomal phylogeny incorporating South Asian (in total 305 individuals) and worldwide (in total 1286 individuals) data from the 1000 Genomes Project. In contrast to the more ancient ancestry in the South than in the North that has been claimed, we detected very similar coalescence times within Northern and Southern non-tribal Indian populations. A closest neighbour analysis in the phylogeny showed that Indian populations have an affinity towards Southern European populations and that the time of divergence from these populations substantially predated the Indo-European migration into India, probably reflecting ancient shared ancestry rather than the Indo-European migration, which had little effect on Indian male lineages. Among the tribal populations, the Birhor (Austro-Asiatic-speaking) and Irula (Dravidian-speaking) are the nearest neighbours of South Asian non-tribal populations, with a common origin in the last few millennia. In contrast, the Riang (Tibeto-Burman-speaking) and Andamanese have their nearest neighbour lineages in East Asia. The Jarawa and Onge shared haplogroup D lineages with each other within the last ~7000 years, but had diverged from Japanese haplogroup D Y-chromosomes ~53000 years ago, most likely by a split from a shared ancestral population. This analysis suggests that Indian populations have complex ancestry which cannot be explained by a single expansion model.

Natural Selection in the Great Apes

Natural selection is crucial for the adaptation of populations to their environments. Here, we present the first global study of natural selection in the Hominidae (humans and great apes) based on genome-wide information from population samples representing all extant species (including most subspecies). Combining several neutrality tests we create a multi-species map of signatures of natural selection covering all major types of natural selection. We find that the estimated efficiency of both purifying and positive selection varies between species and is significantly correlated with their long-term effective population size. Thus, even the modest differences in population size among the closely related Hominidae lineages have resulted in differences in their ability to remove deleterious alleles and to adapt to changing environments. Most signatures of balancing and positive selection are species-specific, with signatures of balancing selection more often being shared among species. We also identify loci with evidence of positive selection across several lineages. Notably, we detect signatures of positive selection in several genes related to brain function, anatomy, diet and immune processes. Our results contribute to a better understanding of human evolution by putting the evidence of natural selection in humans within its larger evolutionary context. The global map of natural selection in our closest living relatives is available as an interactive browser at http://tinyurl.com/nf8qmzh.

Detection of genomic rearrangements from targeted resequencing data in Parkinson's disease patients

Background

The analysis of coverage depth in next‐generation sequencing data allows the detection of gene dose alterations. We explore the frequency of such structural events in a Spanish cohort of sporadic PD cases.

Methods

Gene dose alterations were detected with the eXome‐Hidden Markov Model (XHMM) software from depth of coverage in resequencing data available for 38 Mendelian and other risk PD loci in 394 individuals (249 cases and 145 controls) and subsequently validated by quantitative PCR.

Results

We identified 10 PD patients with exon dosage alterations in PARK2, GBA‐GBAP1, and DJ1. Additional functional variants, including 2 novel nonsense mutations (p.Arg1552Ter in LRRK2 and p.Trp90Ter in PINK1), were confirmed by Sanger sequencing. This combined approach disclosed the genetic cause of 12 PD cases.

Conclusions

Gene dose alterations related to PD can be correctly identified from targeting resequencing data. This approach substantially improves the detection rate of cases with causal genetic alterations. © 2016 The Authors. Movement Disorders published by Wiley Periodicals, Inc. on behalf of International Parkinson and Movement Disorder Society.

Chimpanzee genomic diversity reveals ancient admixture with bonobos

Our closest living relatives, chimpanzees and bonobos, have a complex demographic history. We analyzed the high-coverage whole genomes of 75 wild-born chimpanzees and bonobos from 10 countries in Africa. We found that chimpanzee population substructure makes genetic information a good predictor of geographic origin at country and regional scales. Multiple lines of evidence suggest that gene flow occurred from bonobos into the ancestors of central and eastern chimpanzees between 200,000 and 550,000 years ago, probably with subsequent spread into Nigeria-Cameroon chimpanzees. Together with another, possibly more recent contact (after 200,000 years ago), bonobos contributed less than 1% to the central chimpanzee genomes. Admixture thus appears to have been widespread during hominid evolution.

Genetic Load of Loss-of-Function Polymorphic Variants in Great Apes

Loss of function (LoF) genetic variants are predicted to disrupt gene function, and are therefore expected to substantially reduce individual's viability. Knowing the genetic burden of LoF variants in endangered species is of interest for a better understanding of the effects of declining population sizes on species viability. In this study, we have estimated the number of LoF polymorphic variants in six great ape populations, based on whole-genome sequencing data in 79 individuals. Our results show that although the number of functional variants per individual is conditioned by the effective population size, the number of variants with a drastic phenotypic effect is very similar across species. We hypothesize that for those variants with high selection coefficients, differences in effective population size are not important enough to affect the efficiency of natural selection to remove them. We also describe that mostly CpG LoF mutations are shared across species, and an accumulation of LoF variants at olfactory receptor genes in agreement with its pseudogenization in humans and other primate species.

Population and genomic lessons from genetic analysis of two Indian populations

Indian demographic history includes special features such as founder effects, interpopulation segregation, complex social structure with a caste system and elevated frequency of consanguineous marriages. It also presents a higher frequency for some rare mendelian disorders and in the last two decades increased prevalence of some complex disorders. Despite the fact that India represents about one-sixth of the human population, deep genetic studies from this terrain have been scarce. In this study, we analyzed high-density genotyping and whole-exome sequencing data of a North and a South Indian population. Indian populations show higher differentiation levels than those reported between populations of other continents. In this work, we have analyzed its consequences, by specifically assessing the transferability of genetic markers from or to Indian populations. We show that there is limited genetic marker portability from available genetic resources such as HapMap or the 1,000 Genomes Project to Indian populations, which also present an excess of private rare variants. Conversely, tagSNPs show a high level of portability between the two Indian populations, in contrast to the common belief that North and South Indian populations are genetically very different. By estimating kinship from mates and consanguinity in our data from trios, we also describe different patterns of assortative mating and inbreeding in the two populations, in agreement with distinct mating preferences and social structures. In addition, this analysis has allowed us to describe genomic regions under recent adaptive selection, indicating differential adaptive histories for North and South Indian populations. Our findings highlight the importance of considering demography for design and analysis of genetic studies, as well as the need for extending human genetic variation catalogs to new populations and particularly to those with particular demographic histories.

Whole-exome sequencing reveals a rapid change in the frequency of rare functional variants in a founding population of humans

Whole-exome or gene targeted resequencing in hundreds to thousands of individuals has shown that the majority of genetic variants are at low frequency in human populations. Rare variants are enriched for functional mutations and are expected to explain an important fraction of the genetic etiology of human disease, therefore having a potential medical interest. In this work, we analyze the whole-exome sequences of French-Canadian individuals, a founder population with a unique demographic history that includes an original population bottleneck less than 20 generations ago, followed by a demographic explosion, and the whole exomes of French individuals sampled from France. We show that in less than 20 generations of genetic isolation from the French population, the genetic pool of French-Canadians shows reduced levels of diversity, higher homozygosity, and an excess of rare variants with low variant sharing with Europeans. Furthermore, the French-Canadian population contains a larger proportion of putatively damaging functional variants, which could partially explain the increased incidence of genetic disease in the province. Our results highlight the impact of population demography on genetic fitness and the contribution of rare variants to the human genetic variation landscape, emphasizing the need for deep cataloguing of genetic variants by resequencing worldwide human populations in order to truly assess disease risk.

Exome sequencing identifies mutations in the gene TTC7A in French-Canadian cases with hereditary multiple intestinal atresia

BACKGROUND

Congenital multiple intestinal atresia (MIA) is a severe, fatal neonatal disorder, involving the occurrence of obstructions in the small and large intestines ultimately leading to organ failure. Surgical interventions are palliative but do not provide long-term survival. Severe immunodeficiency may be associated with the phenotype. A genetic basis for MIA is likely. We had previously ascertained a cohort of patients of French-Canadian origin, most of whom were deceased as infants or in utero. The goal of the study was to identify the molecular basis for the disease in the patients of this cohort.

METHODS

We performed whole exome sequencing on samples from five patients of four families. Validation of mutations and familial segregation was performed using standard Sanger sequencing in these and three additional families with deceased cases. Exon skipping was assessed by reverse transcription-PCR and Sanger sequencing.

RESULTS

Five patients from four different families were each homozygous for a four base intronic deletion in the gene TTC7A, immediately adjacent to a consensus GT splice donor site. The deletion was demonstrated to have deleterious effects on splicing causing the skipping of the attendant upstream coding exon, thereby leading to a predicted severe protein truncation. Parents were heterozygous carriers of the deletion in these families and in two additional families segregating affected cases. In a seventh family, an affected case was compound heterozygous for the same 4bp deletion and a second missense mutation p.L823P, also predicted as pathogenic. No other sequenced genes possessed deleterious variants explanatory for all patients in the cohort. Neither mutation was seen in a large set of control chromosomes.

CONCLUSIONS

Based on our genetic results, TTC7A is the likely causal gene for MIA.

Long-standing balancing selection in the THBS4 gene: influence on sex-specific brain expression and gray matter volumes in Alzheimer disease

The THBS4 gene encodes a glycoprotein involved in inflammatory responses and synaptogenesis. THBS4 is expressed at higher levels in the brain of humans compared with nonhuman primates, and the protein accumulates in β-amyloid plaques. We analyzed THBS4 genetic variability in humans and show that two haplotypes (hap1 and hap2) are maintained by balancing selection and modulate THBS4 expression in lymphocytes. Indeed, the balancing selection region covers a predicted transcriptional enhancer. In humans, but not in macaques and chimpanzees, THBS4 brain expression increases with age, and variants in the balancing selection region interact with sex in influencing THBS4 expression (pinteraction = 0.038), with hap1 homozygous females showing lowest expression. In Alzheimer disease (AD) patients, significant interactions between sex and THBS4 genotype were detected for peripheral gray matter (pinteraction = 0.014) and total gray matter (pinteraction = 0.012) volumes. Similarly to the gene expression results, the interaction is mainly mediated by hap1 homozygous AD females, who show reduced volumes. Thus, the balancing selection target in THBS4 is likely represented by one or more variants that regulate tissue-specific and sex-specific gene expression. The selection signature associated with THBS4 might not be related to AD pathogenesis, but rather to inflammatory responses.

Rare allelic forms of PRDM9 associated with childhood leukemogenesis

One of the most rapidly evolving genes in humans, PRDM9, is a key determinant of the distribution of meiotic recombination events. Mutations in this meiotic-specific gene have previously been associated with male infertility in humans and recent studies suggest that PRDM9 may be involved in pathological genomic rearrangements. In studying genomes from families with children affected by B-cell precursor acute lymphoblastic leukemia (B-ALL), we characterized meiotic recombination patterns within a family with two siblings having hyperdiploid childhood B-ALL and observed unusual localization of maternal recombination events. The mother of the family carries a rare PRDM9 allele, potentially explaining the unusual patterns found. From exomes sequenced in 44 additional parents of children affected with B-ALL, we discovered a substantial and significant excess of rare allelic forms of PRDM9. The rare PRDM9 alleles are transmitted to the affected children in half the cases; nonetheless there remains a significant excess of rare alleles among patients relative to controls. We successfully replicated this latter observation in an independent cohort of 50 children with B-ALL, where we found an excess of rare PRDM9 alleles in aneuploid and infant B-ALL patients. PRDM9 variability in humans is thought to influence genomic instability, and these data support a potential role for PRDM9 variation in risk of acquiring aneuploidies or genomic rearrangements associated with childhood leukemogenesis.

Next-generation sequencing approaches for genetic mapping of complex diseases

The advent of next generation sequencing technologies has opened new possibilities in the analysis of human disease. In this review we present the main next-generation sequencing technologies, with their major contributions and possible applications to the study of the genetic etiology of complex diseases.

Genetic adaptation of the antibacterial human innate immunity network

BACKGROUND

Pathogens have represented an important selective force during the adaptation of modern human populations to changing social and other environmental conditions. The evolution of the immune system has therefore been influenced by these pressures. Genomic scans have revealed that immune system is one of the functions enriched with genes under adaptive selection.

RESULTS

Here, we describe how the innate immune system has responded to these challenges, through the analysis of resequencing data for 132 innate immunity genes in two human populations. Results are interpreted in the context of the functional and interaction networks defined by these genes. Nucleotide diversity is lower in the adaptors and modulators functional classes, and is negatively correlated with the centrality of the proteins within the interaction network. We also produced a list of candidate genes under positive or balancing selection in each population detected by neutrality tests and showed that some functional classes are preferential targets for selection.

CONCLUSIONS

We found evidence that the role of each gene in the network conditions the capacity to evolve or their evolvability: genes at the core of the network are more constrained, while adaptation mostly occurred at particular positions at the network edges. Interestingly, the functional classes containing most of the genes with signatures of balancing selection are involved in autoinflammatory and autoimmune diseases, suggesting a counterbalance between the beneficial and deleterious effects of the immune response.

Similarity in recombination rate estimates highly correlates with genetic differentiation in humans

Recombination varies greatly among species, as illustrated by the poor conservation of the recombination landscape between humans and chimpanzees. Thus, shorter evolutionary time frames are needed to understand the evolution of recombination. Here, we analyze its recent evolution in humans. We calculated the recombination rates between adjacent pairs of 636,933 common single-nucleotide polymorphism loci in 28 worldwide human populations and analyzed them in relation to genetic distances between populations. We found a strong and highly significant correlation between similarity in the recombination rates corrected for effective population size and genetic differentiation between populations. This correlation is observed at the genome-wide level, but also for each chromosome and when genetic distances and recombination similarities are calculated independently from different parts of the genome. Moreover, and more relevant, this relationship is robustly maintained when considering presence/absence of recombination hotspots. Simulations show that this correlation cannot be explained by biases in the inference of recombination rates caused by haplotype sharing among similar populations. This result indicates a rapid pace of evolution of recombination, within the time span of differentiation of modern humans.

A population genetic approach to mapping neurological disorder genes using deep resequencing

Deep resequencing of functional regions in human genomes is key to identifying potentially causal rare variants for complex disorders. Here, we present the results from a large-sample resequencing (n = 285 patients) study of candidate genes coupled with population genetics and statistical methods to identify rare variants associated with Autism Spectrum Disorder and Schizophrenia. Three genes, MAP1A, GRIN2B, and CACNA1F, were consistently identified by different methods as having significant excess of rare missense mutations in either one or both disease cohorts. In a broader context, we also found that the overall site frequency spectrum of variation in these cases is best explained by population models of both selection and complex demography rather than neutral models or models accounting for complex demography alone. Mutations in the three disease-associated genes explained much of the difference in the overall site frequency spectrum among the cases versus controls. This study demonstrates that genes associated with complex disorders can be mapped using resequencing and analytical methods with sample sizes far smaller than those required by genome-wide association studies. Additionally, our findings support the hypothesis that rare mutations account for a proportion of the phenotypic variance of these complex disorders.

It is widely accepted that genetic factors play important roles in the etiology of neurological diseases. However, the nature of the underlying genetic variation remains unclear. Critical questions in the field of human genetics relate to the frequency and size effects of genetic variants associated with disease. For instance, the common disease–common variant model is based on the idea that sets of common variants explain a significant fraction of the variance found in common disease phenotypes. On the other hand, rare variants may have strong effects and therefore largely contribute to disease phenotypes. Due to their high penetrance and reduced fitness, such variants are maintained in the population at low frequencies, thus limiting their detection in genome-wide association studies. Here, we use a resequencing approach on a cohort of 285 Autism Spectrum Disorder and Schizophrenia patients and preformed several analyses, enhanced with population genetic approaches, to identify variants associated with both diseases. Our results demonstrate an excess of rare variants in these disease cohorts and identify genes with negative (deleterious) selection coefficients, suggesting an accumulation of variants of detrimental effects. Our results present further evidence for rare variants explaining a component of the genetic etiology of autism and schizophrenia.

Recent human evolution has shaped geographical differences in susceptibility to disease

Background

Searching for associations between genetic variants and complex diseases has been a very active area of research for over two decades. More than 51,000 potential associations have been studied and published, a figure that keeps increasing, especially with the recent explosion of array-based Genome-Wide Association Studies. Even if the number of true associations described so far is high, many of the putative risk variants detected so far have failed to be consistently replicated and are widely considered false positives. Here, we focus on the world-wide patterns of replicability of published association studies.

Results

We report three main findings. First, contrary to previous results, genes associated to complex diseases present lower degrees of genetic differentiation among human populations than average genome-wide levels. Second, also contrary to previous results, the differences in replicability of disease associated-loci between Europeans and East Asians are highly correlated with genetic differentiation between these populations. Finally, highly replicated genes present increased levels of high-frequency derived alleles in European and Asian populations when compared to African populations.

Conclusions

Our findings highlight the heterogeneous nature of the genetic etiology of complex disease, confirm the importance of the recent evolutionary history of our species in current patterns of disease susceptibility and could cast doubts on the status as false positives of some associations that have failed to replicate across populations.

Life-history traits of the stone loach Barbatula barbatula

The life-history tactics of the stone loach Barbatula barbatula were studied in a Mediterranean-type climate stream (Matarranya River) located in the Ebro River basin (north-east Spain). Maximum observed ages were 2+ years in both sexes (1% of individuals), although only 0+ and 1+ year age groups were well represented. It is the lowest longevity reported for this species in its entire distribution. The seasonal growth period started in June and continued until November, but the pattern observed was different to northern populations. Barbatula barbatula in the Matarranya River was a multiple spawner, releasing small batches of oocytes between April and June. The fecundity of females was higher and the size of oocytes smaller in 1984 than in 1985. The relative fecundity (number of ripening and ripe oocytes g(-1) of fish) was lower than in northern European populations. The role of the particular environmental conditions of a Mediterranean stream was discussed in relation to the life-history tactics of B. barbatula.