Alleviating misclassified germline variants in underrepresented populations: A strategy using popmax

PURPOSE

Germline variant interpretation often depends on population-matched control cohorts. This is not feasible for population groups that are underrepresented in current population reference databases.

METHODS

We classify germline variants with population-matched controls for 2 ancestrally diverse cohorts of patients: 132 early-onset or familial colorectal carcinoma patients from Singapore and 100 early-onset colorectal carcinoma patients from the United States. The effects of using a population-mismatched control cohort are simulated by swapping the control cohorts used for each patient cohort, with or without the popmax computational strategy.

RESULTS

Population-matched classifications revealed a combined 62 pathogenic or likely pathogenic (P/LP) variants in 34 genes across both cohorts. Using a population-mismatched control cohort resulted in misclassification of non-P/LP variants as P/LP, driven by the absence of ancestry-specific rare variants in the control cohort. Popmax was more effective in alleviating misclassifications for the Singapore cohort than the US cohort.

CONCLUSION

Underrepresented population groups can suffer from higher rates of false-positive P/LP results. Popmax can partially alleviate these misclassifications, but its efficacy still depends on the degree with which the population groups are represented in the control cohort.

Degradation of methylation signals in cryopreserved DNA

BACKGROUND

Blood-based DNA methylation has shown great promise as a biomarker in a wide variety of diseases. Studies of DNA methylation in blood often utilize samples which have been cryopreserved for years or even decades. Therefore, changes in DNA methylation associated with long-term cryopreservation can introduce biases or otherwise mislead methylation analyses of cryopreserved DNA. However, previous studies have presented conflicting results with studies reporting hypomethylation, no effect, or even hypermethylation of DNA following long-term cryopreservation. These studies may have been limited by insufficient sample sizes, or by their profiling of methylation only on an aggregate global scale, or profiling of only a few CpGs.

RESULTS

We analyzed two large prospective cohorts: a discovery (n = 126) and a validation (n = 136) cohort, where DNA was cryopreserved for up to four years. In both cohorts there was no detectable change in mean global methylation across increasing storage durations as DNA. However, when analysis was performed on the level of individual CpG methylation both cohorts exhibited a greater number of hypomethylated than hypermethylated CpGs at q-value < 0.05 (4049 hypomethylated but only 50 hypermethylated CpGs in discovery, and 63 hypomethylated but only 6 hypermethylated CpGs in validation). The results were the same even after controlling for age, storage duration as buffy coat prior to DNA extraction, and estimated cell type composition. Furthermore, we find that in both cohorts, CpGs have a greater likelihood to be hypomethylated the closer they are to a CpG island; except for CpGs at the CpG islands themselves which are less likely to be hypomethylated.

CONCLUSION

Cryopreservation of DNA after a few years results in a detectable bias toward hypomethylation at the level of individual CpG methylation, though when analyzed in aggregate there is no detectable change in mean global methylation. Studies profiling methylation in cryopreserved DNA should be mindful of this hypomethylation bias, and more attention should be directed at developing more stable methods of DNA cryopreservation for biomedical research or clinical use.

Landscape of germline pathogenic variants in patients with dual primary breast and lung cancer

BACKGROUND

Cancer predisposition is most often studied in the context of single cancers. However, inherited cancer predispositions can also give rise to multiple primary cancers. Yet, there is a paucity of studies on genetic predisposition in multiple primary cancers, especially those outside of well-defined cancer predisposition syndromes. This study aimed to identify germline variants associated with dual primary cancers of the breast and lung.

METHODS

Exome sequencing was performed on germline DNA from 55 Singapore patients (52 [95%] never-smokers) with dual primaries in the breast and lung, confirmed by histopathology. Using two large control cohorts: the local SG10K_Health (n = 9770) and gnomAD non-cancer East Asians (n = 9626); and two additional local case cohorts of early-onset or familial breast cancer (n = 290), and lung cancer (n = 209), variants were assessed for pathogenicity in accordance with ACMG/AMP guidelines. In particular, comparisons were made with known pathogenic or likely pathogenic variants in the ClinVar database, pathogenicity predictions were obtained from in silico prediction software, and case-control association analyses were performed.

RESULTS

Altogether, we identified 19 pathogenic or likely pathogenic variants from 16 genes, detected in 17 of 55 (31%) patients. Six of the 19 variants were identified using ClinVar, while 13 variants were classified pathogenic or likely pathogenic using ACMG/AMP guidelines. The 16 genes include well-known cancer predisposition genes such as BRCA2, TP53, and RAD51D; but also lesser known cancer genes EXT2, WWOX, GATA2, and GPC3. Most of these genes are involved in DNA damage repair, reaffirming the role of impaired DNA repair mechanisms in the development of multiple malignancies. These variants warrant further investigations in additional populations.

CONCLUSIONS

We have identified both known and novel variants significantly enriched in patients with primary breast and lung malignancies, expanding the body of known cancer predisposition variants for both breast and lung cancer. These variants are mostly from genes involved in DNA repair, affirming the role of impaired DNA repair in the predisposition and development of multiple cancers.

Whole-exome sequencing of BRCA-negative breast cancer patients and case-control analyses identify variants associated with breast cancer susceptibility

BACKGROUND

For the majority of individuals with early-onset or familial breast cancer referred for genetic testing, the genetic basis of their familial breast cancer remains unexplained. To identify novel germline variants associated with breast cancer predisposition, whole-exome sequencing (WES) was performed.

METHODS

WES on 290 BRCA1/BRCA2-negative Singaporeans with early-onset breast cancer and/or a family history of breast cancer was done. Case-control analysis against the East-Asian subpopulation (EAS) from the Genome Aggregation Database (gnomAD) identified variants enriched in cases, which were further selected by occurrence in cancer gene databases. Variants were further evaluated in repeated case-control analyses using a second case cohort from the database of Genotypes and Phenotypes (dbGaP) comprising 466 early-onset breast cancer patients from the United States, and a Singapore SG10K_Health control cohort.

RESULTS

Forty-nine breast cancer-associated germline pathogenic variants in 37 genes were identified in Singapore cases versus gnomAD (EAS). Compared against SG10K_Health controls, 13 of 49 variants remain significantly enriched (False Discovery Rate (FDR)-adjusted p < 0.05). Comparing these 49 variants in dbGaP cases against gnomAD (EAS) and SG10K_Health controls revealed 23 concordant variants that were significantly enriched (FDR-adjusted p < 0.05). Fourteen variants were consistently enriched in breast cancer cases across all comparisons (FDR-adjusted p < 0.05). Seven variants in GPRIN2, NRG1, MYO5A, CLIP1, CUX1, GNAS and MGA were confirmed by Sanger sequencing.

CONCLUSIONS

In conclusion, we have identified pathogenic variants in genes associated with breast cancer predisposition. Importantly, many of these variants were significant in a second case cohort from dbGaP, suggesting that the strategy of using case-control analysis to select variants could potentially be utilized for identifying variants associated with cancer susceptibility.

The dimer-dependent catalytic activity of RAF family kinases is revealed through characterizing their oncogenic mutants

Although extensively studied for three decades, the molecular mechanisms that regulate the RAF/MEK/ERK kinase cascade remain ambiguous. Recent studies identified the dimerization of RAF as a key event in the activation of this cascade. Here, we show that in-frame deletions in the β3-αC loop activate ARAF as well as BRAF and other oncogenic kinases by enforcing homodimerization. By characterizing these RAF mutants, we find that ARAF has less allosteric and catalytic activity than the other two RAF isoforms, which arises from its non-canonical APE motif. Further, these RAF mutants exhibit a strong oncogenic potential, and a differential inhibitor resistance that correlates with their dimer affinity. Using these unique mutants, we demonstrate that active RAFs, including the BRAF(V600E) mutant, phosphorylate MEK in a dimer-dependent manner. This study characterizes a special category of oncogenic kinase mutations, and elucidates the molecular basis that underlies the differential ability of RAF isoforms to stimulate MEK-ERK pathway. Further, this study reveals a unique catalytic feature of RAF family kinases that can be exploited to control their activities for cancer therapies.

Comprehensive spectrum of BRCA1 and BRCA2 deleterious mutations in breast cancer in Asian countries

Approximately 5%-10% of breast cancers are due to genetic predisposition caused by germline mutations; the most commonly tested genes are BRCA1 and BRCA2 mutations. Some mutations are unique to one family and others are recurrent; the spectrum of BRCA1/BRCA2 mutations varies depending on the geographical origins, populations or ethnic groups. In this review, we compiled data from 11 participating Asian countries (Bangladesh, Mainland China, Hong Kong SAR, Indonesia, Japan, Korea, Malaysia, Philippines, Singapore, Thailand and Vietnam), and from ethnic Asians residing in Canada and the USA. We have additionally conducted a literature review to include other Asian countries mainly in Central and Western Asia. We present the current pathogenic mutation spectrum of BRCA1/BRCA2 genes in patients with breast cancer in various Asian populations. Understanding BRCA1/BRCA2 mutations in Asians will help provide better risk assessment and clinical management of breast cancer.

Detecting novel genetic variants associated with isoniazid-resistant Mycobacterium tuberculosis

Background

Isoniazid (INH) is a highly effective antibiotic central for the treatment of Mycobacterium tuberculosis (MTB). INH-resistant MTB clinical isolates are frequently mutated in the katG gene and the inhA promoter region, but 10 to 37% of INH-resistant clinical isolates have no detectable alterations in currently known gene targets associated with INH-resistance. We aimed to identify novel genes associated with INH-resistance in these latter isolates.

Methodology/Principal Findings

INH-resistant clinical isolates of MTB were pre-screened for mutations in the katG, inhA, kasA and ndh genes and the regulatory regions of inhA and ahpC. Twelve INH-resistant isolates with no mutations, and 17 INH-susceptible MTB isolates were subjected to whole genome sequencing. Phylogenetically related variants and synonymous mutations were excluded and further analysis revealed mutations in 60 genes and 4 intergenic regions associated with INH-resistance. Sanger sequencing verification of 45 genes confirmed that mutations in 40 genes were observed only in INH-resistant isolates and not in INH-susceptible isolates. The ratios of non-synonymous to synonymous mutations (dN/dS ratio) for the INH-resistance associated mutations identified in this study were 1.234 for INH-resistant and 0.654 for INH-susceptible isolates, strongly suggesting that these mutations are indeed associated with INH-resistance.

Conclusion

The discovery of novel targets associated with INH-resistance described in this study may potentially be important for the development of improved molecular detection strategies.

LARG at chromosome 11q23 has functional characteristics of a tumor suppressor in human breast and colorectal cancer

Deletion of 11q23-q24 is frequent in a diverse variety of malignancies, including breast and colorectal carcinoma, implicating the presence of a tumor suppressor gene at that chromosomal region. We examined a 6-Mb region on 11q23 by high-resolution deletion mapping, using both loss of heterozygosity analysis and customized microarray comparative genomic hybridization. LARG (leukemia-associated Rho guanine-nucleotide exchange factor) (also called ARHGEF12), identified from the analysed region, is frequently underexpressed in breast and colorectal carcinomas with a reduced expression observed in all breast cancer cell lines (n=11), in 12 of 38 (32%) primary breast cancers, 5 of 10 (50%) colorectal cell lines and in 20 of 37 (54%) primary colorectal cancers. Underexpression of the LARG transcript was significantly associated with genomic loss (P=0.00334). Hypermethylation of the LARG promoter was not detected in either breast or colorectal cancer, and treatment of four breast and four colorectal cancer cell lines with 5-aza-2'-deoxycytidine and/or trichostatin A did not result in a reactivation of LARG. Enforced expression of LARG in breast and colorectal cancer cells by stable transfection resulted in reduced cell proliferation and colony formation, as well as in a markedly slower cell migration rate in colorectal cancer cells, providing functional evidence for LARG as a candidate tumor suppressor gene.

Analysis of the role of Mycobacterium tuberculosis kasA gene mutations in isoniazid resistance

Previous studies have suggested that Mycobacterium tuberculosis kasA G312S and G269S gene mutations may represent sequence polymorphisms of the M. tuberculosis East-African-Indian (EAI) and T families, respectively, rather than relating to isoniazid resistance. The present study examined polymorphisms of these two codons in 98 drug-susceptible M. tuberculosis isolates (68 EAI and 30 T isolates). Twenty-eight isolates belonging to a sub-lineage of the EAI family had the kasA G312S mutation, but none of the 30 T isolates had the G269S mutation. The data suggest that the kasA G312S mutation is not related to isoniazid resistance, but represents a sequence polymorphism in a sub-lineage of the EAI family.

Genotype and phenotype relationships and transmission analysis of drug-resistant tuberculosis in Singapore

SETTING

The small urban country of Singapore.

OBJECTIVES

To investigate the relationships between Mycobacterium tuberculosis genotypes and drug-resistant phenotypes and to analyse the transmission of drug-resistant tuberculosis (DR-TB).

DESIGN

A 29-month population-based study comparing drug-resistant and drug-susceptible M. tuberculosis isolates.

RESULTS

We found that multidrug-resistant (MDR) isolates (n = 41, OR 2.66, 95%CI 1.28-5.50), rifampicin-resistant isolates (n = 48, OR 2.88, 95%CI 1.44-5.76), and streptomycin (SM) resistant isolates (n = 103, OR 3.35, 95%CI 1.99-5.62) were more common among Beijing genotype strains than among non-Beijing strains, while SM-resistant isolates were less common in East-African-Indian (EAI) genotype strains than in non-EAI strains (OR 0.30, 95%CI 0.14-0.64). Based on clustering analysis and drug-resistant patterns, 22 of 230 drug-resistant isolates were found to have likely resulted from recent transmission. The estimated transmission rate of DR-TB was 9.6% and that of MDR-TB was 7.7%. The transmission rate of DR-TB was significantly higher among Beijing genotype strains than non-Beijing strains (12.9% vs. 4.4%; P = 0.034).

CONCLUSIONS

Compared to other genotypes, Beijing genotype strains are associated with a higher frequency of drug resistance, including multidrug resistance, and are more transmissible. However, the overall transmission rate of DR-TB in Singapore is low.

EGFR Intron 1 polymorphism in Asian Populations and its correlation with EGFR gene expression and amplification in breast tumor tissues

AIM

The purpose of this study was to test the hypothesis if longer CA dinucleotide repeats are more common in the Asian population and also to gain insights into the interplay between the CA dinucleotide repeats and the frequencies of EGFR gene expression and amplifications as this might have therapeutic implications with regards to treatment with tyrosine kinase inhibitors.

MATERIALS AND METHODS

The EGFR intron 1 polymorphism was analysed in three distinct healthy Asian subjects, namely, Chinese (N = 96), Malays (N = 98) and Indians (N = 100). Comparative genomic hybridisation was performed to investigate for changes in DNA copy number in relation to the polymorphic CA dinucleotide repeats in breast tumor tissues (N = 22).

RESULTS

The frequency of short alleles with 14 and 15 CA repeats were most common in the Asian populations and significantly higher than those reported for Caucasians. The frequency of 20 CA repeats was 5%, almost 13-fold lower than previous reports. EGFR amplifications were detected in 23% and 11% of breast tumor tissues harboring short and long CA repeats, respectively.

CONCLUSION

Our results show that the frequency of alleles encoding for short CA dinucleotide repeats is common in Asian populations. EGFR expression and amplification levels were also higher in Asian breast tumor tissues with short CA dinucleotide repeats. These findings suggest that the EGFR intron 1 polymorphism may influence response to treatment with tyrosine kinase inhibitors in breast cancer patients and further studies are warranted.

Association of Mycobacterium tuberculosis Beijing genotype with tuberculosis relapse in Singapore

The relationship of Mycobacterium tuberculosis Beijing genotype with tuberculosis relapse was examined. Beijing strains were detected from 32 out of 45 (71%) relapsed cases and 148 out of 290 (51%) non-relapsed cases. Multivariate logistic regression analysis revealed that Beijing genotype was significantly associated with tuberculosis relapse (OR 2.64, 95% CI 1.30-5.34, P=0.005).

Evidence for deletion of 9q as a two-step process in chronic myeloid leukemia

Evidence for deletion of 9q as a two-step process in chronic myeloid leukemia. Chronic myeloid leukemia (CML) is characterized by the Philadelphia translocation t(9;22)(q34;q11) resulting in the BCR/ABL fusion gene. Submicroscopic deletion of the derivative chromosome 9 occurs in a subset of these patients and is associated with poor prognosis. In the current study, we present two unusual cases of CML selected from a series of 54 consecutive cases. A detailed study using classical cytogenetics and fluorescence in situ hybridization (FISH) analysis was done using dual color extra signal FISH and whole chromosome paint in order to elucidate the mechanism of 9q deletion. One case had two clones on interphase FISH, one with and one without chromosome 9q deletion. The other case had two clones on both cytogenetic and FISH analyses, one with and one without a marker chromosome carrying chromosome 9q sequences. In this latter case, the clone with deletion of the derivative chromosome 9 comprised 21.1% at diagnosis, increasing to 36.8% after 11 months, suggesting a growth advantage. We report here evidence that deletions on 9q in CML may occur through breakage and rearrangement of chromosomes resulting in derivative chromosomes and either a marker chromosome or fragment/episome, followed by loss of chromosome material from the cell.

Characterization of ancestral Mycobacterium tuberculosis by multiple genetic markers and proposal of genotyping strategy

Sixty-eight ancestral Mycobacterium tuberculosis isolates were previously identified by using the tuberculosis-specific deletion 1 (TbD1) PCR and mycobacterial interspersed-repetitive-unit-variable-number tandem repeat (MIRU-VNTR) typing (Y. J. Sun, R. Bellamy, A. S. G. Lee, S. T. Ng, S. Ravindran, S.-Y. Wong, C. Locht, P. Supply, and N. I. Paton, J. Clin. Microbiol. 42:1986-1993, 2004). These TbD1(+) ancestral isolates were further characterized and typed in this study by IS6110 restriction fragment length polymorphism (RFLP) typing, VNTR typing using exact tandem repeats (VNTR-ETR), and spoligotyping of the direct-repeat region. To our knowledge, this is the first characterization of this genogroup by multiple genetic markers based on a fairly large sample size. In this genogroup, all spoligotypes were characterized by the absence of spacers 29 to 32 and 34. In addition, VNTR-ETR typing could add further resolution to the clustered isolates identified by MIRU-VNTR, and the combination of MIRU-VNTR and VNTR-ETR, called MIRU-ETR, showed the highest discriminatory power for these strains compared to IS6110 RFLP typing and spoligotyping alone. However, MIRU-ETR appeared to still cluster some probably epidemiologically unrelated strains, as judged by IS6110 RFLP divergence. Therefore, a typing strategy based on stepwise combination of MIRU-ETR and IS6110 RFLP is proposed to achieve maximal discrimination for unrelated TbD1(+) strains. This typing strategy may be useful in areas where TbD1(+) ancestral strains are prevalent.

Novel mutations within the embB gene in ethambutol-susceptible clinical isolates of Mycobacterium tuberculosis

Genetic analysis of the embB gene revealed mutations in 17 (68%) of 25 ethambutol (EMB) resistant isolates (M306I, M306V, M306L, Q497R) but also in 4 (20%) of 20 EMB-susceptible isolates of Mycobacterium tuberculosis, namely, an ATG-->ATM substitution resulting in M306I, G406N, and the novel alterations M423I and A659T.

Use of mycobacterial interspersed repetitive unit-variable-number tandem repeat typing to examine genetic diversity of Mycobacterium tuberculosis in Singapore

Strain typing using variable-number tandem repeats of mycobacterial interspersed repetitive units (MIRU-VNTR) is a powerful tool for studying the epidemiology and genetic relationships of Mycobacterium tuberculosis isolates. For this study, isolates from 291 patients in Singapore were genotyped by this method. One hundred sixty-six distinct MIRU-VNTR patterns were detected. One hundred sixty-two strains were grouped into 1 of 35 different MIRU-VNTR clusters and 131 isolates were unique. In this sample collection, 9 of the 12 MIRU-VNTR loci were moderately or highly discriminative according to their allelic diversities. The Hunter-Gaston discriminatory index was 0.975, indicating the high power of discrimination of MIRU-VNTR typing. By direct comparisons with previously typed MIRU-VNTR patterns and by genetic relationship analyses, we could identify and clearly define four epidemic groups of M. tuberculosis in our sample, corresponding to the W/Beijing, East-Africa-Indian, Haarlem, and Delhi genotype families. Furthermore, MIRU-VNTR typing was able to clearly distinguish ancestral and modern M. tuberculosis strains as defined by TbD1 genomic deletion analysis. These results indicate that MIRU-VNTR typing can be a useful first-line tool for studying the genetic diversity of M. tuberculosis isolates in a large urban setting such as Singapore.

Founder mutation in the BRCA1 gene in Malay breast cancer patients from Singapore

The mutation spectrum of the BRCA1 gene among ethnic groups from Asia has not been well studied. We investigated the frequency of mutations in the BRCA1 gene among Malay breast cancer patients from Singapore, independent of family history. By using the protein truncation test (PTT) and direct sequencing, BRCA1 mutations were detected in 6 of 49 (12.2%) unrelated patients. Four novel missense mutations in exon 11, T557A (1788A>G), T582A (1863A>G), N656S (2086A>G) and P684S (2169C>T) were identified in one patient. Two patients had missense mutations in exon 23, V1809A (5545T>C), which has been previously detected in individuals from Central and Eastern Europe. Three unrelated patients had the deleterious 2846insA frameshift mutation in exon 11. Methylation specific PCR (MSP) of the promoter region of the BRCA1 gene detected hypermethylation of tumor DNA in an additional 2 patients. Haplotype analysis using the microsatellite markers D17S855, D17S1323 and D17S1325 revealed a common haplotype for the three unrelated patients and their three relatives with the 2846insA mutation. These findings strongly suggest that the 2846insA mutation, the most common deleterious mutation in this study, may possibly be a founder mutation in breast cancer patients of Malay ethnic background.

Characterization of pyrazinamide and ofloxacin resistance among drug resistant Mycobacterium tuberculosis isolates from Singapore

OBJECTIVES

To evaluate rapid molecular approaches for the detection of pyrazinamide (PZA) and ofloxacin resistance, by screening 100 known drug-resistant Mycobacterium tuberculosis isolates.

METHODS

Mycobacterium tuberculosis isolates were tested for phenotypic resistance to pyrazinamide and ofloxacin using the BACTEC 460 radiometric method and the E-test, respectively. Mutation screening was done by amplifying the pncA, gyrA, and gyrB genes by the polymerase chain reaction (PCR) and direct automated sequencing.

RESULTS

Twelve isolates were PZA-resistant and 8 of 12 (66.7%) isolates had missense mutations or deletions at the pncA gene, suggesting that mutation or deletion at the pncA gene is the major molecular mechanism of PZA resistance among the Singaporean isolates. Using the E-test, 48 isolates were resistant to ofloxacin, with minimum inhibitory concentrations of 4 microg/mL or higher. No mutations were observed at the quinolone resistance-determining region (QRDR) of gyrA in all isolates. At the QRDR of gyrB, mutations were present in 1 of 48 ofloxacin-resistant isolates and 0 of 19 ofloxacin-susceptible isolates.

CONCLUSIONS

In Singapore, genotypic analysis of resistance to PZA and ofloxacin is inadequate and should be complemented by conventional methods.

Discrimination of single-copy IS6110 DNA fingerprints of Mycobacterium tuberculosis isolates by high-resolution minisatellite-based typing

Seven isoniazid-resistant isolates with mutations in the NADH dehydrogenase (ndh) gene were molecularly typed by IS6110-based restriction fragment length polymorphism analysis. All seven isolates with the R268H mutation had identical 1.4-kb IS6110 fingerprints. High-resolution minisatellite-based typing discriminated five of these isolates; two isolates were identical.

Prevalence of hepatitis G virus infection in patients with liver diseases in Singapore

The prevalence of hepatitis G virus (HGV) infection in patients with liver diseases in Singapore and its pathogenic role in these patients was studied. One hundred and forty-eight patients who had chronic hepatitis or acute non A-E hepatitis were studied. Presence of HGV RNA was determined by nested polymerase chain reaction of the 5'non-coding region of the virus in all the patients. Hepatitis G IgG antibody to the envelope (E2) antigen was tested with an enzyme immunoassay (Boehringer Mannheim, Singapore) in 76 of them. Most patients (93%) were ethnically Chinese, predominantly males (74%) and chronic hepatitis B (72%) patients. Others had chronic hepatitis C (19%) or cryptogenic cirrhosis (6%). Four patients had acute non A-E hepatitis. HGV RNA and anti-HGenv were present in 3.5% and 8.3% of those with chronic liver disease. HGV infection did not account for any of the acute non A-E hepatitis and most of the cryptogenic cirrhosis.