Association between 15 known or potential breast cancer susceptibility genes and breast cancer risks in Chinese women.

OBJECTIVE: There are many hereditary breast cancer patients in China, and multigene panel testing has been a new paradigm of genetic testing for these patients and their relatives. However, the magnitude of breast cancer risks related to multiple breast cancer susceptibility genes are largely unknown in Chinese women. METHODS: We screened pathogenic variants in 15 established or potential breast cancer susceptibility genes from 8,067 consecutive Chinese female breast cancer patients and 13,129 Chinese cancer-free female controls. These breast cancer patients were unselected for age at diagnosis or family history. RESULTS: We found that pathogenic variants in TP53 [odds ratio (OR): 16.9, 95% confidence interval (CI): 5.2-55.2]; BRCA2 (OR: 10.4, 95% CI: 7.6-14.2); BRCA1 (OR: 9.7, 95% CI: 6.3-14.8); and PALB2 (OR: 5.2, 95% CI: 3.0-8.8) were associated with a high risk of breast cancer. ATM, BARD1, CHEK2, and RAD51D were associated with a moderate risk of breast cancer with ORs ranging from 2-fold to 4-fold. In contrast, pathogenic variants of NBN, RAD50, BRIP1, and RAD51C were not associated with increased risk of breast cancer in Chinese women. The pathogenic variants of PTEN, CDH1, and STK11 were very rare, so they had a limited contribution to Chinese breast cancer. Patients with pathogenic variants of TP53, BRCA1, BRCA2, and PALB2 more often had early-onset breast cancer, bilateral breast cancer, and a family history of breast cancer and/or any cancer. CONCLUSIONS: This study provided breast cancer risk assessment data for multiple genes in Chinese women, which is useful for genetic testing and clinical management of Chinese hereditary breast cancer.

Monoallelic Mutations in CC2D1A Suggest a Novel Role in Human Heterotaxy and Ciliary Dysfunction.

BACKGROUND: Human heterotaxy is a group of congenital disorders characterized by misplacement of one or more organs according to the left-right axis. The genetic causes of human heterotaxy are highly heterogeneous. METHODS: We performed exome sequencing in a cohort of 26 probands with heterotaxy followed by gene burden analysis for the enrichment of novel rare damaging mutations. Transcription activator-like effector nuclease was used to generate somatic loss-of-function mutants in a zebrafish model. Ciliary defects were examined by whole-mount immunostaining of acetylated α-tubulin. RESULTS: We identified a significant enrichment of novel rare damaging mutations in the CC2D1A gene. Seven occurrences of CC2D1A mutations were found to affect 4 highly conserved amino acid residues of the protein. Functional analyses in the transcription activator-like effector nuclease-mediated zebrafish knockout models were performed, and heterotaxy phenotypes of the cardiovascular and gastrointestinal systems in both somatic and germline mutants were observed. Defective cilia were demonstrated by whole-mount immunostaining of acetylated α-tubulin. These abnormalities were rescued by wild-type cc2d1a mRNA but not cc2d1a mutant mRNA, strongly suggesting a loss-of-function mechanism. On the other hand, overexpression of cc2d1a orthologous mutations cc2d1a P559L and cc2d1a G808V (orthologous to human CC2D1A P532L and CC2D1A G781V) did not affect embryonic development. CONCLUSIONS: Using a zebrafish model, we were able to establish a novel association of CC2D1A with heterotaxy and ciliary dysfunction in the F2 generation via a loss-of-function mechanism. Future mechanistic studies are needed for a better understanding of the role of CC2D1A in left-right patterning and ciliary dysfunction.

Development and validation of next generation sequencing based 35-gene hereditary cancer panel.

BACKGROUND: Understanding the genetic basis of cancer risk is a major international endeavor. The emergence of next-generation sequencing (NGS) in late 2000's has further accelerated the discovery of many cancer susceptibility genes. The use of targeted NGS-based multigene testing panels to provide comprehensive analysis of cancer susceptible genes has proven to be a viable option, with the accurate and robust detection of a wide range of clinically relevant variants in the targeted genes being crucial. METHODS: We have developed and validated a targeted NGS-based test for hereditary cancer risk assessment using Illumina's NGS platform by analyzing the protein-coding regions of 35 hereditary cancer genes with a bioinformatics pipeline that utilizes standard practices in the field. This 35-gene hereditary cancer panel is designed to identify germline cancer-causing mutations for 8 different cancers: breast, ovarian, prostate, uterine, colorectal, pancreatic, stomach cancers and melanoma. The panel was validated using well-characterized DNA specimens [NIGMS Human Genetic Cell Repository], where DNA had been extracted using blood of individuals whose genetic variants had been previously characterized by the 1000 Genome Project and the Coriell Catalog. RESULTS: The 35-gene hereditary cancer panel shows high sensitivity (99.9%) and specificity (100%) across 4820 variants including single nucleotide variants (SNVs) and small insertions and deletions (indel; up to 25 bp). The reproducibility and repeatability are 99.8 and 100%, respectively. CONCLUSIONS: The use of targeted NGS-based multigene testing panels to provide comprehensive analysis of cancer susceptible genes has been considered a viable option. In the present study, we developed and validated a 35-gene panel for testing 8 common cancers using next-generation sequencing (NGS). The performance of our hereditary cancer panel is assessed across a board range of variants in the 35 genes to support clinical use.

A powerful approach reveals numerous expression quantitative trait haplotypes in multiple tissues.

Motivation: Recently many studies showed single nucleotide polymorphisms (SNPs) affect gene expression and contribute to development of complex traits/diseases in a tissue context-dependent manner. However, little is known about haplotype's influence on gene expression and complex traits, which reflects the interaction effect between SNPs. Results: In the present study, we firstly proposed a regulatory region guided eQTL haplotype association analysis approach, and then systematically investigate the expression quantitative trait loci (eQTL) haplotypes in 20 different tissues by the approach. The approach has a powerful design of reducing computational burden by the utilization of regulatory predictions for candidate SNP selection and multiple testing corrections on non-independent haplotypes. The application results in multiple tissues showed that haplotype-based eQTLs not only increased the number of eQTL genes in a tissue specific manner, but were also enriched in loci that associated with complex traits in a tissue-matched manner. In addition, we found that tag SNPs of eQTL haplotypes from whole blood were selectively enriched in certain combination of regulatory elements (e.g. promoters and enhancers) according to predicted chromatin states. In summary, this eQTL haplotype detection approach, together with the application results, shed insights into synergistic effect of sequence variants on gene expression and their susceptibility to complex diseases. Availability and implementation: The executable application 'eHaplo' is implemented in Java and is publicly available at http://grass.cgs.hku.hk/limx/ehaplo/. Supplementary information: Supplementary data are available at Bioinformatics online.

cepip: context-dependent epigenomic weighting for prioritization of regulatory variants and disease-associated genes.

It remains challenging to predict regulatory variants in particular tissues or cell types due to highly context-specific gene regulation. By connecting large-scale epigenomic profiles to expression quantitative trait loci (eQTLs) in a wide range of human tissues/cell types, we identify critical chromatin features that predict variant regulatory potential. We present cepip, a joint likelihood framework, for estimating a variant's regulatory probability in a context-dependent manner. Our method exhibits significant GWAS signal enrichment and is superior to existing cell type-specific methods. Furthermore, using phenotypically relevant epigenomes to weight the GWAS single-nucleotide polymorphisms, we improve the statistical power of the gene-based association test.

CFTR founder mutation causes protein trafficking defects in Chinese patients with cystic fibrosis.

BACKGROUND: Cystic fibrosis (CF) is a rare condition in Asians. Since 1985, only about 30 Chinese patients have been reported with molecular confirmation. METHOD: Using our in-house next-generation sequencing (NGS) pipeline for childhood bronchiectasis, we identified disease-causing CFTR mutations in CF patients in Hong Kong. After identifying p.I1023R in multiple patients, haplotype analysis was performed with genome-wide microarray to ascertain the likelihood of this being a founder mutation. We also assessed the processing and gating activity of the mutant protein by Western hybridization and patch-clamp test. RESULTS: Molecular diagnoses were confirmed in four patients, three of whom shared a missense mutation: CFTR:c.3068T>G:p.I1023R. The results suggested that p.I1023R is a founder mutation in southern Han Chinese. In addition, the processing and gating activity of the mutant protein was assessed by gel electrophoresis and a patch-clamp test. The mutant protein exhibited trafficking defects, suggesting that the dysfunction is caused by reduced cell surface expression of the fully glycosylated proteins. CONCLUSION: Together with other previously reported mutations, the specific founder mutation presented herein suggests a unique CFTR mutation spectrum in the southern Chinese populations, and this finding has vital implications for improving molecular testing and mutation-specific treatments for Chinese patients with CF.

Identification of mutations in the PI3K-AKT-mTOR signalling pathway in patients with macrocephaly and developmental delay and/or autism.

Background: Macrocephaly, which is defined as a head circumference greater than or equal to + 2 standard deviations, is a feature commonly observed in children with developmental delay and/or autism spectrum disorder. Although PTEN is a well-known gene identified in patients with this syndromic presentation, other genes in the PI3K-AKT-mTOR signalling pathway have also recently been suggested to have important roles. The aim of this study is to characterise the mutation spectrum of this group of patients. Methods: We performed whole-exome sequencing of 21 patients with macrocephaly and developmental delay/autism spectrum disorder. Sources of genomic DNA included blood, buccal mucosa and saliva. Germline mutations were validated by Sanger sequencing, whereas somatic mutations were validated by droplet digital PCR. Results: We identified ten pathogenic/likely pathogenic mutations in PTEN (n = 4), PIK3CA (n = 3), MTOR (n = 1) and PPP2R5D (n = 2) in ten patients. An additional PTEN mutation, which was classified as variant of unknown significance, was identified in a patient with a pathogenic PTEN mutation, making him harbour bi-allelic germline PTEN mutations. Two patients harboured somatic PIK3CA mutations, and the level of somatic mosaicism in blood DNA was low. Patients who tested positive for mutations in the PI3K-AKT-mTOR pathway had a lower developmental quotient than the rest of the cohort (DQ = 62.8 vs. 76.1, p = 0.021). Their dysmorphic features were non-specific, except for macrocephaly. Among the ten patients with identified mutations, brain magnetic resonance imaging was performed in nine, all of whom showed megalencephaly. Conclusion: We identified mutations in the PI3K-AKT-mTOR signalling pathway in nearly half of our patients with macrocephaly and developmental delay/autism spectrum disorder. These patients have subtle dysmorphic features and mild developmental issues. Clinically, patients with germline mutations are difficult to distinguish from patients with somatic mutations, and therefore, sequencing of buccal or saliva DNA is important to identify somatic mosaicism. Given the high diagnostic yield and the management implications, we suggest implementing comprehensive genetic testing in the PI3K-AKT-mTOR pathway in the clinical evaluation of patients with macrocephaly and developmental delay and/or autism spectrum disorder.

Prenatal Tobacco Exposure Shortens Telomere Length in Children.

INTRODUCTION: Preliminary evidence suggests a possible association between prenatal tobacco exposure and telomere length in children. This study was conducted to investigate whether maternal smoking during pregnancy was associated with telomere shortening in their children and whether prenatal and childhood exposure to environmental tobacco had any impact on this association. METHODS: This is a population-representative study on the association between prenatal tobacco exposure and telomere length in children. Ninety-eight Hong Kong Chinese children aged under 15 years with prenatal tobacco exposure and 98 age- and gender-matched controls were recruited from a population health study with stratified random sampling. RESULTS: Telomere length in children with prenatal tobacco exposure was significantly shorter than in those with no exposure (mean T/S ratio = 24.9 [SD = 8.58] in exposed vs. 28.97 [14.15] in control groups; P = 0.02). A negative dose-response relationship was observed between the T/S ratio and tobacco exposure duration: the longer the duration of maternal smoking in pregnancy, the shorter the child's telomere length. The association between the child's telomere length and prenatal tobacco exposure remained significant after considering the influence of family socioeconomic status and exposure to environmental tobacco smoke during pregnancy and childhood. CONCLUSIONS: Prenatal tobacco exposure was associated with telomere shortening in children. As this may impose significant health impacts through fetal genetic programming, more efforts should be made to reduce fetal tobacco exposure by educating pregnant women to not smoke and motivating smokers to quit in early pregnancy. IMPLICATIONS: As reflected by telomere shortening, prenatal tobacco exposure in children can cause premature aging and increased health risks, which we suggest is entirely preventable. Not smoking during pregnancy or quitting smoking is critical to improving the health outcome of our future generations as prenatal tobacco exposure may affect children's biological programming.

Gene-based meta-analysis of genome-wide association study data identifies independent single-nucleotide polymorphisms in ANXA6 as being associated with systemic lupus erythematosus in Asian populations.

OBJECTIVE: Previous genome-wide association studies (GWAS), which were mainly based on single-variant analysis, have identified many systemic lupus erythematosus (SLE) susceptibility loci. However, the genetic architecture of this complex disease is far from being understood. The aim of this study was to investigate whether using a gene-based analysis may help to identify novel loci, by considering global evidence of association from a gene or a genomic region rather than focusing on evidence for individual variants. METHODS: Based on the results of a meta-analysis of 2 GWAS of SLE conducted in 2 Asian cohorts, we performed an in-depth gene-based analysis followed by replication in a total of 4,626 patients and 7,466 control subjects of Asian ancestry. Differential allelic expression was measured by pyrosequencing. RESULTS: More than one-half of the reported SLE susceptibility loci showed evidence of independent effects, and this finding is important for understanding the mechanisms of association and explaining disease heritability. ANXA6 was detected as a novel SLE susceptibility gene, with several single-nucleotide polymorphisms (SNPs) contributing independently to the association with disease. The risk allele of rs11960458 correlated significantly with increased expression of ANXA6 in peripheral blood mononuclear cells from heterozygous healthy control subjects. Several other associated SNPs may also regulate ANXA6 expression, according to data obtained from public databases. Higher expression of ANXA6 in patients with SLE was also reported previously. CONCLUSION: Our study demonstrated the merit of using gene-based analysis to identify novel susceptibility loci, especially those with independent effects, and also demonstrated the widespread presence of loci with independent effects in SLE susceptibility genes.

HaploShare: identification of extended haplotypes shared by cases and evaluation against controls.

Recent founder mutations may play important roles in complex diseases and Mendelian disorders. Detecting shared haplotypes that are identical by descent (IBD) could facilitate discovery of these mutations. Several programs address this, but are usually limited to detecting pair-wise shared haplotypes and not providing a comparison of cases and controls. We present a novel algorithm and software package, HaploShare, which detects extended haplotypes that are shared by multiple individuals, and allows comparisons between cases and controls. Testing on simulated and real cases demonstrated significant improvements in detection power and reduction of false positive rate by HaploShare relative to other programs.

HLAreporter: a tool for HLA typing from next generation sequencing data.

Human leukocyte antigen (HLA) typing from next generation sequencing (NGS) data has the potential for widespread applications. Here we introduce a novel tool (HLAreporter) for HLA typing from NGS data based on read-mapping using a comprehensive reference panel containing all known HLA alleles, followed by de novo assembly of the gene-specific short reads. Accurate HLA typing at high-digit resolution was achieved when it was tested on publicly available NGS data, outperforming other newly developed tools such as HLAminer and PHLAT. HLAreporter can be downloaded from http://paed.hku.hk/genome/.

Meta-analysis of GWAS on two Chinese populations followed by replication identifies novel genetic variants on the X chromosome associated with systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a prototypic autoimmune disease that affects mainly females. What role the X chromosome plays in the disease has always been an intriguing question. In this study, we examined the genetic variants on the X chromosome through meta-analysis of two genome-wide association studies (GWAS) on SLE on Chinese Han populations. Prominent association signals from the meta-analysis were replicated in 4 additional Asian cohorts, with a total of 5373 cases and 9166 matched controls. We identified a novel variant in PRPS2 on Xp22.3 as associated with SLE with genome-wide significance (rs7062536, OR = 0.84, P = 1.00E-08). Association of the L1CAM-MECP2 region with SLE was reported previously. In this study, we identified independent contributors in this region in NAA10 (rs2071128, OR = 0.81, P = 2.19E-13) and TMEM187 (rs17422, OR = 0.75, P = 1.47E-15), in addition to replicating the association from IRAK1-MECP2 region (rs1059702, OR = 0.71, P = 2.40E-18) in Asian cohorts. The X-linked susceptibility variants showed higher effect size in males than that in females, similar to results from a genome-wide survey of associated SNPs on the autosomes. These results suggest that susceptibility genes identified on the X chromosome, while contributing to disease predisposition, might not contribute significantly to the female predominance of this prototype autoimmune disease.

Three SNPs in chromosome 11q23.3 are independently associated with systemic lupus erythematosus in Asians.

Systemic lupus erythematosus (SLE) has a complex etiology and is affected by both genetic and environmental factors. Although more than 40 loci have shown robust association with SLE, the details of these loci, such as the independent contributors and the genes involved, are still unclear. In this study, we performed meta-analysis of two existing genome-wide association studies (GWASs) on Chinese Han populations from Hong Kong and Anhui, China, and followed the findings by further replication on three additional Chinese and Thailand cohorts with a total of 4254 cases and 6262 controls matched geographically and ethnically. We discovered multiple susceptibility variants for SLE in the 11q23.3 region, including variants in/near PHLDB1 (rs11603023, P(_combined) = 1.25E-08, OR = 1.20), DDX6 (rs638893, P(_combined) = 5.19E-07, OR = 1.22) and CXCR5 (rs10892301, P(_combined) = 2.51E-08, OR = 0.85). Genetic contributions from the newly identified variants were all independent of SNP rs4639966, whose association was reported from the previous GWAS. In addition, the three newly identified variants all showed independent association with the disease through modeling by both stepwise and conditional logistic regression. The presence of multiple independent variants in this region emphasizes its role in SLE susceptibility, and also hints the possibility that distinct biological mechanisms might be involved in the disease involving this genomic region.

Epistatic interaction between genetic variants in susceptibility gene ETS1 correlates with IL-17 levels in SLE patients.

T-helper cells that produce IL-17 (Th17 cells) are a subset of CD4(+) T-cells with pathological roles in autoimmune diseases including systemic lupus erythematosus (SLE), and ETS1 is a negative regulator of Th17 cell differentiation. Our previous work on genome-wide association study (GWAS) identified two variants in the ETS1 gene (rs10893872 and rs1128334) as being associated with SLE. However, like many other risk alleles for complex diseases, little is known on how these genetic variants might affect disease pathogenesis. In this study, we examined serum IL-17 levels from 283 SLE cases and observed a significant correlation between risk variants in ETS1 and serum IL-17 concentration in patients, which suggests a potential mechanistic link between these variants and the disease. Furthermore, we found that the two variants act synergistically in influencing IL-17 production, with evidence of significant genetic interaction between them as well as higher correlation between the haplotype formed by the risk alleles and IL-17 level in patient serum. In addition, the correlation between ETS1 variants and IL-17 level seems to be more significant in SLE patients manifesting renal involvement, dsDNA autoantibody production or early-onset.

Meta-analysis followed by replication identifies loci in or near CDKN1B, TET3, CD80, DRAM1, and ARID5B as associated with systemic lupus erythematosus in Asians.

Systemic lupus erythematosus (SLE) is a prototype autoimmune disease with a strong genetic involvement and ethnic differences. Susceptibility genes identified so far only explain a small portion of the genetic heritability of SLE, suggesting that many more loci are yet to be uncovered for this disease. In this study, we performed a meta-analysis of genome-wide association studies on SLE in Chinese Han populations and followed up the findings by replication in four additional Asian cohorts with a total of 5,365 cases and 10,054 corresponding controls. We identified genetic variants in or near CDKN1B, TET3, CD80, DRAM1, and ARID5B as associated with the disease. These findings point to potential roles of cell-cycle regulation, autophagy, and DNA demethylation in SLE pathogenesis. For the region involving TET3 and that involving CDKN1B, multiple independent SNPs were identified, highlighting a phenomenon that might partially explain the missing heritability of complex diseases.

PriVar: a toolkit for prioritizing SNVs and indels from next-generation sequencing data.

UNLABELLED: Next-generation sequencing has become a valuable tool for detecting mutations involved in Mendelian diseases. However, it is a challenge to identify the small subset of functionally important mutations from tens of thousands of rare variants in a whole exome/genome. Therefore, we developed a toolkit called PriVar, a systematic prioritization pipeline that takes into consideration calling quality of the variants, their predicted functional impact, known connection of the gene to the disease and the number of mutations in a gene, and inference from linkage analysis. AVAILABILITY: Executable jar package is available at http://paed.hku.hk/uploadarea/yangwl/html/software.html.

Homozygosity mapping on a single patient: identification of homozygous regions of recent common ancestry by using population data.

Homozygosity mapping has played an important role in detecting recessive mutations using families of consanguineous marriages. However, detection of regions identical and homozygosity by descent (HBD) when family data are not available, or when relationships are unknown, is still a challenge. Making use of population data from high-density SNP genotyping may allow detection of regions HBD from recent common founders in singleton patients without genealogy information. We report a novel algorithm that detects such regions by estimating the population haplotype frequencies (HF) for an entire homozygous region. We also developed a simulation method to evaluate the probability of HBD and linkage to disease for a homozygous region by examining the best regions in unaffected controls from the host population. The method can be applied to diseases of Mendelian inheritance but can also be extended to complex diseases to detect rare founder mutations that affect a very small number of patients using either multiplex families or sporadic cases. Testing of the method on both real cases (singleton affected) and simulated data demonstrated its superb sensitivity and robustness under genetic heterogeneity.

Genome-wide association study in Asian populations identifies variants in ETS1 and WDFY4 associated with systemic lupus erythematosus.

Systemic lupus erythematosus is a complex and potentially fatal autoimmune disease, characterized by autoantibody production and multi-organ damage. By a genome-wide association study (320 patients and 1,500 controls) and subsequent replication altogether involving a total of 3,300 Asian SLE patients from Hong Kong, Mainland China, and Thailand, as well as 4,200 ethnically and geographically matched controls, genetic variants in ETS1 and WDFY4 were found to be associated with SLE (ETS1: rs1128334, P = 2.33x10(-11), OR = 1.29; WDFY4: rs7097397, P = 8.15x10(-12), OR = 1.30). ETS1 encodes for a transcription factor known to be involved in a wide range of immune functions, including Th17 cell development and terminal differentiation of B lymphocytes. SNP rs1128334 is located in the 3'-UTR of ETS1, and allelic expression analysis from peripheral blood mononuclear cells showed significantly lower expression level from the risk allele. WDFY4 is a conserved protein with unknown function, but is predominantly expressed in primary and secondary immune tissues, and rs7097397 in WDFY4 changes an arginine residue to glutamine (R1816Q) in this protein. Our study also confirmed association of the HLA locus, STAT4, TNFSF4, BLK, BANK1, IRF5, and TNFAIP3 with SLE in Asians. These new genetic findings may help us to gain a better understanding of the disease and the functions of the genes involved.

In-depth cDNA library sequencing provides quantitative gene expression profiling in cancer biomarker discovery.

Quantitative gene expression analysis plays an important role in identifying differentially expressed genes in various pathological states, gene expression regulation and co-regulation, shedding light on gene functions. Although microarray is widely used as a powerful tool in this regard, it is suboptimal quantitatively and unable to detect unknown gene variants. Here we demonstrated effective detection of differential expression and co-regulation of certain genes by expressed sequence tag analysis using a selected subset of cDNA libraries. We discussed the issues of sequencing depth and library preparation, and propose that increased sequencing depth and improved preparation procedures may allow detection of many expression features for less abundant gene variants. With the reduction of sequencing cost and the emerging of new generation sequencing technology, in-depth sequencing of cDNA pools or libraries may represent a better and powerful tool in gene expression profiling and cancer biomarker detection. We also propose using sequence-specific subtraction to remove hundreds of the most abundant housekeeping genes to increase sequencing depth without affecting relative expression ratio of other genes, as transcripts from as few as 300 most abundantly expressed genes constitute about 20% of the total transcriptome. In-depth sequencing also represents a unique advantage of detecting unknown forms of transcripts, such as alternative splicing variants, fusion genes, and regulatory RNAs, as well as detecting mutations and polymorphisms that may play important roles in disease pathogenesis.