VNtyper enables accurate alignment-free genotyping of MUC1 coding VNTR using short-read sequencing data in autosomal dominant tubulointerstitial kidney disease

The human genome comprises approximately 3% of tandem repeats with variable length (VNTR), a few of which have been linked to human rare diseases. Autosomal dominant tubulointerstitial kidney disease-MUC1 (ADTKD-MUC1) is caused by specific frameshift variants in the coding VNTR of the MUC1 gene. Calling variants from VNTR using short-read sequencing (SRS) is challenging due to poor read mappability. We developed a computational pipeline, VNtyper, for reliable detection of MUC1 VNTR pathogenic variants and demonstrated its clinical utility in two distinct cohorts: (1) a historical cohort including 108 families with ADTKD and (2) a replication naive cohort comprising 2,910 patients previously tested on a panel of genes involved in monogenic renal diseases. In the historical cohort all cases known to carry pathogenic MUC1 variants were re-identified, and a new 25bp-frameshift insertion in an additional mislaid family was detected. In the replication cohort, we discovered and validated 30 new patients.

Identification of Germline Non-coding Deletions in XIAP Gene Causing XIAP Deficiency Reveals a Key Promoter Sequence

PURPOSE

X-linked inhibitor of apoptosis protein (XIAP) deficiency, also known as the X-linked lymphoproliferative syndrome of type 2 (XLP-2), is a rare immunodeficiency characterized by recurrent hemophagocytic lymphohistiocytosis, splenomegaly, and inflammatory bowel disease. Variants in XIAP including missense, non-sense, frameshift, and deletions of coding exons have been reported to cause XIAP deficiency. We studied three young boys with immunodeficiency displaying XLP-2-like clinical features. No genetic variation in the coding exons of XIAP was identified by whole-exome sequencing (WES), although the patients exhibited a complete loss of XIAP expression.

METHODS

Targeted next-generation sequencing (NGS) of the entire locus of XIAP was performed on DNA samples from the three patients. Molecular investigations were assessed by gene reporter expression assays in HEK cells and CRISPR-Cas9 genome editing in primary T cells.

RESULTS

NGS of XIAP identified three distinct non-coding deletions in the patients that were predicted to be driven by repetitive DNA sequences. These deletions share a common region of 839 bp that encompassed the first non-coding exon of XIAP and contained regulatory elements and marks specific of an active promoter. Moreover, we showed that among the 839 bp, the exon was transcriptionally active. Finally, deletion of the exon by CRISPR-Cas9 in primary cells reduced XIAP protein expression.

CONCLUSIONS

These results identify a key promoter sequence contained in the first non-coding exon of XIAP. Importantly, this study highlights that sequencing of the non-coding exons that are not currently captured by WES should be considered in the genetic diagnosis when no variation is found in coding exons.

Battle of the Sex Chromosomes: Competition between X and Y Chromosome-Encoded Proteins for Partner Interaction and Chromatin Occupancy Drives Multicopy Gene Expression and Evolution in Muroid Rodents

Transmission distorters (TDs) are genetic elements that favor their own transmission to the detriments of others. Slx/Slxl1 (Sycp3-like-X-linked and Slx-like1) and Sly (Sycp3-like-Y-linked) are TDs, which have been coamplified on the X and Y chromosomes of Mus species. They are involved in an intragenomic conflict in which each favors its own transmission, resulting in sex ratio distortion of the progeny when Slx/Slxl1 versus Sly copy number is unbalanced. They are specifically expressed in male postmeiotic gametes (spermatids) and have opposite effects on gene expression: Sly knockdown leads to the upregulation of hundreds of spermatid-expressed genes, whereas Slx/Slxl1-deficiency downregulates them. When both Slx/Slxl1 and Sly are knocked down, sex ratio distortion and gene deregulation are corrected. Slx/Slxl1 and Sly are, therefore, in competition but the molecular mechanism remains unknown. By comparing their chromatin-binding profiles and protein partners, we show that SLX/SLXL1 and SLY proteins compete for interaction with H3K4me3-reader SSTY1 (Spermiogenesis-specific-transcript-on-the-Y1) at the promoter of thousands of genes to drive their expression, and that the opposite effect they have on gene expression is mediated by different abilities to recruit SMRT/N-Cor transcriptional complex. Their target genes are predominantly spermatid-specific multicopy genes encoded by the sex chromosomes and the autosomal Speer/Takusan. Many of them have coamplified with not only Slx/Slxl1/Sly but also Ssty during muroid rodent evolution. Overall, we identify Ssty as a key element of the X versus Y intragenomic conflict, which may have influenced gene content and hybrid sterility beyond Mus lineage since Ssty amplification on the Y predated that of Slx/Slxl1/Sly.

A Variety of Alu-Mediated Copy Number Variations Can Underlie IL-12Rβ1 Deficiency

PURPOSE

Inborn errors of IFN-γ immunity underlie Mendelian susceptibility to mycobacterial disease (MSMD). Autosomal recessive complete IL-12Rβ1 deficiency is the most frequent genetic etiology of MSMD. Only two of the 84 known mutations are copy number variations (CNVs), identified in two of the 213 IL-12Rβ1-deficient patients and two of the 164 kindreds reported. These two CNVs are large deletions found in the heterozygous or homozygous state. We searched for novel families with IL-12Rβ1 deficiency due to CNVs.

METHODS

We studied six MSMD patients from five unrelated kindreds displaying adverse reactions to BCG vaccination. Three of the patients also presented systemic salmonellosis, two had mucocutaneous candidiasis, and one had disseminated histoplasmosis. We searched for CNVs and other variations by IL12RB1-targeted next-generation sequencing (NGS).

RESULTS

We identified six new IL-12Rβ1-deficient patients with a complete loss of IL-12Rβ1 expression on phytohemagglutinin-activated T cells and/or EBV-transformed B cells. The cells of these patients did not respond to IL-12 and IL-23. Five different CNVs encompassing IL12RB1 (four deletions and one duplication) were identified in these patients by NGS coverage analysis, either in the homozygous state (n = 1) or in trans (n = 4) with a single-nucleotide variation (n = 3) or a small indel (n = 1). Seven of the nine mutations are novel. Interestingly, four of the five CNVs were predicted to be driven by nearby Alu elements, as well as the two previously reported large deletions. The IL12RB1 locus is actually enriched in Alu elements (44.7%), when compared with the rest of the genome (10.5%).

CONCLUSION

The IL12RB1 locus is Alu-enriched and therefore prone to rearrangements at various positions. CNVs should be considered in the genetic diagnosis of IL-12Rβ1 deficiency.

Expression and epigenomic landscape of the sex chromosomes in mouse post-meiotic male germ cells

BACKGROUND

During meiosis, the X and Y chromosomes are transcriptionally silenced. The persistence of repressive chromatin marks on the sex chromatin after meiosis initially led to the assumption that XY gene silencing persists to some extent in spermatids. Considering the many reports of XY-linked genes expressed and needed in the post-meiotic phase of mouse spermatogenesis, it is still unclear whether or not the mouse sex chromatin is a repressive or permissive environment, after meiosis.

RESULTS

To determine the transcriptional and chromatin state of the sex chromosomes after meiosis, we re-analyzed ten ChIP-Seq datasets performed on mouse round spermatids and four RNA-seq datasets from male germ cells purified at different stages of spermatogenesis. For this, we used the last version of the genome (mm10/GRCm38) and included reads that map to several genomic locations in order to properly interpret the high proportion of sex chromosome-encoded multicopy genes. Our study shows that coverage of active epigenetic marks H3K4me3 and Kcr is similar on the sex chromosomes and on autosomes. The post-meiotic sex chromatin nevertheless differs from autosomal chromatin in its enrichment in H3K9me3 and its depletion in H3K27me3 and H4 acetylation. We also identified a posttranslational modification, H3K27ac, which specifically accumulates on the Y chromosome. In parallel, we found that the X and Y chromosomes are enriched in genes expressed post-meiotically and display a higher proportion of spermatid-specific genes compared to autosomes. Finally, we observed that portions of chromosome 14 and of the sex chromosomes share specific features, such as enrichment in H3K9me3 and the presence of multicopy genes that are specifically expressed in round spermatids, suggesting that parts of chromosome 14 are under the same evolutionary constraints than the sex chromosomes.

CONCLUSIONS

Based on our expression and epigenomic studies, we conclude that, after meiosis, the mouse sex chromosomes are no longer silenced but are nevertheless regulated differently than autosomes and accumulate different chromatin marks. We propose that post-meiotic selective constraints are at the basis of the enrichment of spermatid-specific genes and of the peculiar chromatin composition of the sex chromosomes and of parts of chromosome 14.

Association of REL polymorphisms and outcome of patients with septic shock

Background

cRel, a subunit of NF-κB, is implicated in the inflammatory response observed in autoimmune disease. Hence, knocked-out mice for cRel had a significantly higher mortality, providing new and important functions of cRel in the physiopathology of septic shock. Whether genetic variants in the human REL gene are associated with severity of septic shock is unknown.

Methods

We genotyped a population of 1040 ICU patients with septic shock and 855 ICU controls for two known polymorphisms of REL; REL rs842647 and REL rs13031237. Outcome of patients according to the presence of REL variant alleles was compared.

Results

The distribution of REL variant alleles was not significantly different between patients and controls. Among the septic shock group, REL rs13031237*T minor allele was not associated with worse outcome. In contrast, REL rs842647*G minor allele was significantly associated with more multi-organ failure and early death [OR 1.4; 95 % CI (1.02–1.8)].

Conclusion

In a large ICU population, we report a significant clinical association between a variation in the human REL gene and severity and mortality of septic shock, suggesting for the first time a new insight into the role of cRel in response to infection in humans.

Assessing the impact of a combined analysis of four common low-risk genetic variants on autism risk

Background

Autism is a complex disorder characterized by deficits involving communication, social interaction, and repetitive and restrictive patterns of behavior. Twin studies have shown that autism is strongly heritable, suggesting a strong genetic component. In other disease states with a complex etiology, such as type 2 diabetes, cancer and cardiovascular disease, combined analysis of multiple genetic variants in a genetic score has helped to identify individuals at high risk of disease. Genetic scores are designed to test for association of genetic markers with disease.

Method

The accumulation of multiple risk alleles markedly increases the risk of being affected, and compared with studying polymorphisms individually, it improves the identification of subgroups of individuals at greater risk. In the present study, we show that this approach can be applied to autism by specifically looking at a high-risk population of children who have siblings with autism. A two-sample study design and the generation of a genetic score using multiple independent genes were used to assess the risk of autism in a high-risk population.

Results

In both samples, odds ratios (ORs) increased significantly as a function of the number of risk alleles, with a genetic score of 8 being associated with an OR of 5.54 (95% confidence interval [CI] 2.45 to 12.49). The sensitivities and specificities for each genetic score were similar in both analyses, and the resultant area under the receiver operating characteristic curves were identical (0.59).

Conclusions

These results suggest that the accumulation of multiple risk alleles in a genetic score is a useful strategy for assessing the risk of autism in siblings of affected individuals, and may be better than studying single polymorphisms for identifying subgroups of individuals with significantly greater risk.

Estimating haplotype relative risks in complex disease from unphased SNPs data in families using a likelihood adjusted for ascertainment

The understanding of complex diseases and insights to improve their medical management may be achieved through the deduction of how specific haplotypes may play a joint effect to change relative risk information. In this paper we describe an ascertainment adjusted likelihood-based method to estimate haplotype relative risks using pooled family data coming from association and/or linkage studies that were used to identify specific haplotypes. Haplotype-based analysis tends to require a large amount of parameters to capture all the information that leads to efficiency problems. An adaptation of the Stochastic Expectation Maximization algorithm is used for haplotypes inference from genotypic data and to reduce the number of nuisance parameters for risk estimation. Using different simulations, we show that this method provides unbiased relative risk estimates even in case of departure from Hardy-Weinberg equilibrium.

Haplotypes in the gene encoding protein kinase c-beta (PRKCB1) on chromosome 16 are associated with autism

Autism is a developmental disorder characterized by impairments in social interaction and communication associated with repetitive patterns of interest or behavior. Autism is highly influenced by genetic factors. Genome-wide linkage and candidate gene association approaches have been used to try and identify autism genes. A few loci have repeatedly been reported linked to autism. Several groups reported evidence for linkage to a region on chromosome 16p. We have applied a direct physical identity-by-descent (IBD) mapping approach to perform a high-density (0.85 megabases) genome-wide linkage scan in 116 families from the AGRE collection. Our results confirm linkage to a region on chromosome 16p with autism. High-resolution single-nucleotide polymorphism (SNP) genotyping and analysis of this region show that haplotypes in the protein kinase c-beta gene are strongly associated with autism. An independent replication of the association in a second set of 167 trio families with autism confirmed our initial findings. Overall, our data provide evidence that the PRKCB1 gene on chromosome 16p may be involved in the etiology of autism.

The art of pedigree drawing: algorithmic aspects

MOTIVATION

Giving a meaningful representation of a pedigree is not obvious when it includes consanguinity loops, individuals with multiple mates or several related families.

RESULTS

We show that finding a perfectly meaningful representation of a pedigree is equivalent to the interval graph sandwich problem and we propose an algorithm for drawing pedigrees.

Sib-pair linkage analysis of alcohol dependence taking into account covariates and age-of-onset variability: evaluation of the residual approach

The construction of residuals allowed us to consider a phenotype related to alcohol dependence (ALDX1) adjusted for age, gender and current smoking status. Using the Haseman and Elston regression test, we compared genome scan results for detection of loci for the unadjusted and adjusted ALDX1 phenotypes in the Collaborative Study on the Genetics of Alcoholism (COGA) data. More markers were significant at the 1% level for the binary than for the residual analyses. Overall, our residual analyses do not show substantial improvement in their ability to detect linkage for alcohol dependence in the COGA data. This may be due to the relatively small sample size of discordant sib pairs.