Next-generation sequencing for identifying new genes in rare genetic diseases: many challenges and a pinch of luck

A report on the European Society of Human Genetics conference, held in Paris, France, June 8-11, 2013.     

A Weighted Polygenic Risk Score Using 14 Known Susceptibility Variants to Estimate Risk and Age Onset of Psoriasis in Han Chinese

With numbers of common variants identified mainly through genome-wide association studies (GWASs), there is great interest in incorporating the findings into screening individuals at high risk of psoriasis. The purpose of this study is to establish genetic prediction models and evaluate its discriminatory ability in psoriasis in Han Chinese population. We built the genetic prediction models through weighted polygenic risk score (PRS) using 14 susceptibility variants in 8,819 samples. We found the risk of psoriasis among individuals in the top quartile of PRS was significantly larger than those in the lowest quartile of PRS (OR = 28.20, P < 2.0×10^-16). We also observed statistically significant associations between the PRS, family history and early age onset of psoriasis. We also built a predictive model with all 14 known susceptibility variants and alcohol consumption, which achieved an area under the curve statistic of ~ 0.88. Our study suggests that 14 psoriasis known susceptibility loci have the discriminating potential, as is also associated with family history and age of onset. This is the genetic predictive model in psoriasis with the largest accuracy to date.     

Sequencing of 19,219 exomes identifies a low-frequency variant in FKBP5 promoter predisposing to high myopia in a Han Chinese population

1. Download : Download high-res image (134KB)
2. Download : Download full-size image

     

Genetic Variants at 20p11 Confer Risk to Androgenetic Alopecia in the Chinese Han Population

Background

Androgenetic alopecia (AGA) is a well-characterized type of progressive hair loss commonly seen in men, with different prevalences in different ethnic populations. It is generally considered to be a polygenic heritable trait. Several susceptibility genes/loci, such as AR/EDA2R, HDAC9 and 20p11, have been identified as being involved in its development in European populations. In this study, we aim to validate whether these loci are also associated with AGA in the Chinese Han population.

Methods

We genotyped 16 previously reported single nucleotide polymorphisms (SNPs) with 445 AGA cases and 546 healthy controls using the Sequenom iPlex platform. The trend test was used to evaluate the association between these loci and AGA in the Chinese Han population. Conservatively accounting for multiple testing by the Bonferroni correction, the threshold for statistical significance was P ≤3.13×10⁻³.

Results

We identified that 5 SNPs at 20p11 were significantly associated with AGA in the Chinese Han population (1.84×10⁻¹¹≤P≤2.10×10⁻⁶).

Conclusions

This study validated, for the first time, that 20p11 also confers risk for AGA in the Chinese Han population and implicated the potential common genetic factors for AGA shared by both Chinese and European populations.     

Polymorphisms of ECE1 May Contribute to Susceptibility to Ischemic Stroke in Han Chinese of Northern China

Endothelin (ET) converting enzyme 1 (ECE1) is well known for its critical role in the process of ET. Recent studies have demonstrated that two genetic variants of ECE1 gene, rs212528 and rs213045 (C338A), are associated with hypertension and atherosclerosis formation. To investigate the association between the ECE1 gene polymorphisms and ischemic stroke (IS) in Chinese population. With a candidate loci strategy, we conducted a case–control study involving 381 IS cases and 366 non-IS controls in Han population of the Northern China. Two single nucleotide polymorphisms of ECE1 were genotyped and assessed the association with the risk of IS. Furthermore, stratified analyses were also carried out to evaluate the association between the gender or two etiologic subtypes [small-artery occlusion (SAO) and large-artery atherosclerosis (LAA)] and IS. Compared with rs213045 G homozygote, rs213045 TG genotype and rs213045 TT/TG genotypes are in dominant model significantly increased the risk of IS [adjusted odds ratio (OR) = 1.47, 95 % confidence interval (CI) = 1.04–2.07, P = 0.03; adjusted OR = 1.43, 95 % CI = 1.04–1.99, P = 0.029, respectively]. However, we did not find the significant association between the ECE1 gene polymorphisms and SAO or LAA of IS in this study. Our results indicated that ECE1 gene polymorphisms may contribute to the susceptibility of IS in Han population of the Northern China.     

Association between TL1A gene polymorphisms and systemic lupus erythematosus in a Chinese Han population

Our previous studies showed elevated tumor necrosis factor-like ligand 1 aberrance (TL1A) expression in systemic lupus erythematosus (SLE). However, TL1A polymorphisms with SLE susceptibility remain to be elucidated. In addition, we made meta-analysis to evaluate the relationship of TL1A polymorphisms and autoimmune diseases owing to inconsistent results. The present research was carried out by 404 SLE, 150 primary Sjogren's syndrome (pSS) patients, and 574 healthy individuals. Three TL1A polymorphisms (rs3810936, rs6478109, rs7848647) were genotyped using TaqMan genotyping assay. Then, the meta-analysis was performed by collecting the present case-control study and previously published research. Results showed that genotypes of rs3810936, rs7848647 were different between SLE patients and healthy controls, whereas no significant association was observed in the three polymorphisms and pSS patients. Genotypes distribution of rs6478109, rs7848647 were strongly related to lupus nephritis within SLE (p = 0.004, p = 0.011), respectively. Moreover, combined meta-analysis consisted of ten comparative research involving 4,305 patients and 5,600 controls. An association between autoimmune diseases and rs6478109 polymorphism was found. Our findings indicate that gene polymorphisms (rs3810936, rs7848647) of TL1A might correlate with lupus.     

A comprehensive study of CD44 rs 187115 variant and cancer risk in a central Chinese population

The rs187115, an intronic variant of CD44 gene, has been previously reported to play a potential role in genetic susceptibility to cancer. Here, we comprehensively examined the association between CD44 rs187115 variant and cancer risk (breast cancer, cervical cancer, lung cancer, gastric cancer, liver cancer, colon cancer, and rectal cancer) in a central Chinese population. The rs187115 variant was genotyped with the polymerase chain reaction-restriction fragment length polymorphism assay. In this study, we observed that rs187115 was associated with the risk of cervical, lung, and liver cancer, but not with the risk of breast, gastric, colon, rectal or colorectal cancer. Of note, the G allele and G allele genotypes of rs187115 conferred an increased risk of cervical, lung, and liver cancer. To improve the statistical strength, a followed meta-analysis was conducted. The results demonstrated that rs187115 was significantly associated with cancer risk, and the significant association remained in the stratification analysis by ethnicity, genotyping method, and cancer type. Collectively, the CD44 rs187115 variant may be associated with the risk of cervical, lung, and liver cancer in the central Chinese population, and may be used as a potential biomarker for cancer predisposition in the Asian population, especially in the Chinese population.     

Bi-allelic variants in DNHD1 cause flagellar axoneme defects and asthenoteratozoospermia in humans and mice

1. Download : Download high-res image (285KB)
2. Download : Download full-size image

     

Genetic Variants in PNPLA3 and Risk of Non-Alcoholic Fatty Liver Disease in a Han Chinese Population

We investigated the possible association between genetic variants in the Patatin like phospholipase-3 (PNPLA3) gene and nonalcoholic fatty liver disease (NAFLD) in a Han Chinese population. We evaluated twelve tagging single-nucleotide polymorphisms (tSNPs) of the PNPLA3 gene in a frequency matched case–control study from Fuzhou city of China (553 cases, 553 controls). In the multivariate logistic regression analysis, the rs738409 GG or GC, and rs139051 TT genotypes were found to be associated with increased risk of NAFLD, and a significant trend of increased risk with increasing numbers of risk genotype was observed in the cumulative effect analysis of these single nucleotide polymorphisms. Furthermore, haplotype association analysis showed that, compared with the most common haplotype, the CAAGAATGCGTG and CGAAGGTGTCCG haplotypes conferred a statistically significant increased risk for NAFLD, while the CGGGAACCCGCG haplotype decreased the risk of NAFLD. Moreover, rs738409 C>G appeared to have a multiplicative joint effect with tea drinking (P<0.005) and an additive joint effect with obesity (Interaction contrast ratio (ICR) = 2.31, 95% CI: 0.7–8.86), hypertriglyceridemia (ICR = 3.07, 95% CI: 0.98–5.09) or hypertension (ICR = 1.74, 95% CI: 0.52–3.12). Our data suggests that PNPLA3 genetic polymorphisms might influence the susceptibility to NAFLD development independently or jointly in Han Chinese.     

A novel mutation in PCD-associated gene DNAAF3 causes male infertility due to asthenozoospermia

Primary ciliary dyskinesia (PCD) is a rare autosomal-recessive disease manifested with recurrent infections of respiratory tract and infertility. DNAAF3 is identified as a novel gene associated with PCD and different mutations in DNAAF3 results in different clinical features of PCD patients, such as situs inversus, sinusitis and bronchiectasis. However, the sperm phenotypic characteristics of PCD males are generally poorly investigated. Our reproductive medicine centre received a case of PCD patient with infertility, who presented with sinusitis, recurrent infections of the lower airway and severe asthenozoospermia; However, no situs inversus was found in the patient. A novel homozygous mutation in DNAAF3(c.551T>A; p.V184E) was identified in the PCD patient by whole-exome sequencing. Subsequent Sanger sequencing further confirmed that the DNAAF3 had a homozygous missense variant in the fifth exon. Transmission electron microscopy and immunostaining analysis of the sperms from the patient showed a complete absence of outer dynein arms and partial absence of inner dynein arms, which resulted in the reduction in sperm motility. However, this infertility was overcome by intracytoplasmic sperm injections, as his wife achieved successful pregnancy. These findings showed that the PCD-associated pathogenic mutation within DNAAF3 also causes severe asthenozoospermia and male infertility ultimately due to sperm flagella axoneme defect in humans. Our study not only contributes to understand the sperm phenotypic characteristics of patients with DNAAF3 mutations but also expands the spectrum of DNAAF3 mutations and may contribute to the genetic diagnosis and therapy for infertile patient with PCD.     

Influence of CTNNB1 rs2953 polymorphism on schizophrenia susceptibility in Chinese Han population through modifying miR‐485 binding to CTNNB1

Schizophrenia (SCZ) and bipolar disorder (BD) are two major neuropsychiatric diseases that are the most substantial causes of disability and mortality worldwide. CTNNB1 encodes beta‐catenin, an important protein in canonical Wnt signaling. We aimed to investigate the association between the rs2953 of CTNNB1 and the risk of SCZ and BD and to further explore the function of rs2953. A total of 1658 samples (548 SCZ cases, 512 BD cases, and 598 controls) were examined in terms of the genotype of CTNNB1 rs2953. The mRNA expression level of CTNNB1 significantly increased in the SCZ and BD groups compared with that in the control group. Significant association remained between CTNNB1 3′‐untranslated region (UTR) variant rs2953 and SCZ susceptibility (additive and dominant model) after gender and age were adjusted. rs2953 disrupted the binding of CTNNB1 and miR‐485. miR‐485 significantly suppressed the luciferase activity of CTNNB1‐T vector by binding to the CTNNB1 3′‐UTR containing the T allele of rs2953. The mRNA expression of CTNNB1 can be used as a biomarker for the diagnosis of SCZ and BD. The 3′‐UTR variant rs2953 in CTNNB1 influences the risk of SCZ in the Han Chinese population and modifies the binding of miR‐485 to CTNNB1.     

Genetic association between IL-27 rs153109 polymorphism and cancer risk in Chinese population: a meta-analysis

IL-27 plays an important role in anti-cancer activity. The -964A/G polymorphism in IL-27 gene has been implicated in susceptibility to cancer, but the results were conflicting. The aim of this study was to assess the association between this polymorphism and cancer risk. Pubmed and Wanfang database were searched for all publications concerning IL-27 -964A/G polymorphism and cancer risk. Odds ratio (OR) and 95% confidence interval (CI) were used to assess the strength of association. Statistical analysis was performed using Stata 11.0 software. A total of eight case–control studies including 2044 cancer cases and 2197 controls were identified. Overall, significant association between IL-27 -964A/G polymorphism and cancer risk was observed (GG versus AA: OR?=?1.26, 95% CI?=?1.03–1.52; GG versus AG?+?AA: OR?=?1.20, 95% CI?=?1.00–1.44). In subgroup analysis based on cancer type, significant association was found in colorectal cancer (GG versus AA: OR?=?1.55, 95% CI?=?1.07–2.27; AG versus AA: OR?=?1.31, 95% CI?=?1.02–1.67). The current meta-analysis suggests that IL-27 -964A/G polymorphism might enhance cancer risk. However, large-scale and well-designed studies are still needed to confirm the result of our meta-analysis. The association of IL-27 polymorphism with colorectal cancer may provide insight for future therapies.     

CTLA4 Variants and Haplotype Contribute Genetic Susceptibility to Myasthenia Gravis in Northern Chinese Population

Background

Cytotoxic T lymphocyte-associated antigen-4 (CTLA4), a critical negative regulator of the T-cell response, has been considered a candidate for many autoimmune diseases. Evidence from Caucasians supported a genetic predisposition of CTLA4 to myasthenia gravis (MG), but the contribution in East Asians has not been established.

Objectives

To investigate the role of CTLA4 variants in the susceptibility to MG and the contribution to subtypes of MG.

Methods

Six autoimmune disease-related risk alleles of CTLA4 (rs1863800, rs733618, rs4553808, rs5742909, rs231775, and rs3087243) were investigated for MG in northern Chinese. 168 patients with MG (mean age 37.1±20.5 years, 64 men and 104 women) and 233 healthy controls (mean age 53.3±8.7 years, 96 men and 137 women) were screened, and the contribution of CTLA4 to the general risk of MG and each subgroup was explored.

Results

rs1863800*C, rs733618*C, and rs231775*G were significantly associated with the whole cohort of patients with MG after permutation correction for multiple-testing adjustment (P = 0.027, 0.001, and 0.032, respectively). A risk haplotype (CCACG) [odds ratio (OR) = 1.535, range = 1.150–2.059, P = 0.004)] was also identified. The stratified subtype analysis indicated that the positive contribution was possibly derived from early onset MG (EOMG), seropositive MG (SPMG), female patients, and MG without thymoma. No association was observed in juvenile MG/LOMG, and MG coupled with thymoma.

Conclusion

A predisposing effect of rs1863800*C, rs733618*C, and rs231775*G of CTLA4 gene to general risk of MG in Chinese was demonstrated for the first time, which was likely derived from EOMG, SPMG, MG without thymoma and the female patients.     

Validation of Type 2 Diabetes Risk Variants Identified by Genome-Wide Association Studies in Han Chinese Population: A Replication Study and Meta-Analysis

Background

Several genome-wide association studies (GWAS) involving European populations have successfully identified risk genetic variants associated with type 2 diabetes mellitus (T2DM). However, the effects conferred by these variants in Han Chinese population have not yet been fully elucidated.

Methods

We analyzed the effects of 24 risk genetic variants with reported associations from European GWAS in 3,040 Han Chinese subjects in Taiwan (including 1,520 T2DM cases and 1,520 controls). The discriminative power of the prediction models with and without genotype scores was compared. We further meta-analyzed the association of these variants with T2DM by pooling all candidate-gene association studies conducted in Han Chinese.

Results

Five risk variants in IGF2BP2 (rs4402960, rs1470579), CDKAL1 (rs10946398), SLC30A8 (rs13266634), and HHEX (rs1111875) genes were nominally associated with T2DM in our samples. The odds ratio was 2.22 (95% confidence interval, 1.81-2.73, P<0.0001) for subjects with the highest genetic score quartile (score>34) as compared with subjects with the lowest quartile (score<29). The incoporation of genotype score into the predictive model increased the C-statistics from 0.627 to 0.657 (P<0.0001). These estimates are very close to those observed in European populations. Gene-environment interaction analysis showed a significant interaction between rs13266634 in SLC30A8 gene and age on T2DM risk (P<0.0001). Further meta-analysis pooling 20 studies in Han Chinese confirmed the association of 10 genetic variants in IGF2BP2, CDKAL1, JAZF1, SCL30A8, HHEX, TCF7L2, EXT2, and FTO genes with T2DM. The effect sizes conferred by these risk variants in Han Chinese were similar to those observed in Europeans but the allele frequencies differ substantially between two populations.

Conclusion

We confirmed the association of 10 variants identified by European GWAS with T2DM in Han Chinese population. The incorporation of genotype scores into the prediction model led to a small but significant improvement in T2DM prediction.     

Low Frequency Variants in the Exons Only Encoding Isoform A of HNF1A Do Not Contribute to Susceptibility to Type 2 Diabetes

Background

There is considerable interest in the hypothesis that low frequency, intermediate penetrance variants contribute to the proportion of Type 2 Diabetes (T2D) susceptibility not attributable to the common variants uncovered through genome-wide association approaches. Genes previously implicated in monogenic and multifactorial forms of diabetes are obvious candidates in this respect. In this study, we focussed on exons 8–10 of the HNF1A gene since rare, penetrant mutations in these exons (which are only transcribed in selected HNF1A isoforms) are associated with a later age of diagnosis of Maturity onset diabetes of the young (MODY) than mutations in exons 1–7. The age of diagnosis in the subgroup of HNF1A-MODY individuals with exon 8–10 mutations overlaps with that of early multifactorial T2D, and we set out to test the hypothesis that these exons might also harbour low-frequency coding variants of intermediate penetrance that contribute to risk of multifactorial T2D.

Methodology and Principal Findings

We performed targeted capillary resequencing of HNF1A exons 8–10 in 591 European T2D subjects enriched for genetic aetiology on the basis of an early age of diagnosis (≤45 years) and/or family history of T2D (≥1 affected sibling). PCR products were sequenced and compared to the published HNF1A sequence. We identified several variants (rs735396 [IVS9−24T>C], rs1169304 [IVS8+29T>C], c.1768+44C>T [IVS9+44C>T] and rs61953349 [c.1545G>A, p.T515T] but no novel non-synonymous coding variants were detected.

Conclusions and Significance

We conclude that low frequency, nonsynonymous coding variants in the terminal exons of HNF1A are unlikely to contribute to T2D-susceptibility in European samples. Nevertheless, the rationale for seeking low-frequency causal variants in genes known to contain rare, penetrant mutations remains strong and should motivate efforts to screen other genes in a similar fashion.     

A metagenomic approach to dissect the genetic composition of enterotypes in Han Chinese and two Muslim groups

Distinct enterotypes have been observed in the human gut but little is known about the genetic basis of the microbiome. Moreover, it is not clear how many genetic differences exist between enterotypes within or between populations. In this study, both the 16S rRNA gene and the metagenomes of the gut microbiota were sequenced from 48 Han Chinese, 48 Kazaks, and 96 Uyghurs, and taxonomies were assigned after de novo assembly. Single nucleotide polymorphisms were also identified by referring to data from the Human Microbiome Project. Systematic analysis of the gut communities in terms of their abundance and genetic composition was also performed, together with a genome-wide association study of the host genomes. The gut microbiota of 192 subjects was clearly classified into two enterotypes (Bacteroides and Prevotella). Interestingly, both enterotypes showed a clear genetic differentiation in terms of their functional catalogue of genes, especially for genes involved in amino acid and carbohydrate metabolism. In addition, several differentiated genera and genes were found among the three populations. Notably, one human variant (rs878394) was identified that showed significant association with the abundance of Prevotella, which is linked to LYPLAL1, a gene associated with body fat distribution, the waist-hip ratio and insulin sensitivity. Taken together, considerable differentiation was observed in gut microbes between enterotypes and among populations that was reflected in both the taxonomic composition and the genetic makeup of their functional genes, which could have been influenced by a variety of factors, such as diet and host genetic variation.