Genetic control of resistance to experimental infection with virulent Mycobacterium tuberculosis

Over 2 billion people are estimated to be infected with virulent Mycobacterium tuberculosis, yet fewer than 10% progress to clinical tuberculosis within their lifetime. Twin studies and variations in the outcome of tuberculosis infection after exposure to similar environmental risks suggest genetic heterogeneity among individuals in their susceptibility to disease. In a mouse model of tuberculosis, we have established that resistance and susceptibility to virulent M. tuberculosis is a complex genetic trait. A new locus with a major effect on tuberculosis susceptibility, designated sst1 (susceptibility to tuberculosis 1), was mapped to a 9-centimorgan (cM) interval on mouse chromosome 1. It is located 10–19 cM distal to a previously identified gene, Nramp1, that controls the innate resistance of mice to the attenuated bacillus Calmette–Guérin vaccine strain. The phenotypic expression of the newly identified locus is distinct from that of Nramp1 in that sst1 controls progression of tuberculosis infection in a lung-specific manner. Mice segregating at the sst1 locus exhibit marked differences in the growth rates of virulent tubercle bacilli in the lungs. Lung lesions in congenic sst1-susceptible mice are characterized by extensive necrosis and unrestricted extracellular multiplication of virulent mycobacteria, whereas sst1-resistant mice develop interstitial granulomas and effectively control multiplication of the bacilli. The resistant allele of sst1, although powerful in controlling infection, is not sufficient to confer full protection against virulent M. tuberculosis, indicating that other genes located outside of the sst1 locus are likely also to be important for controlling tuberculosis infection.

Lack of association between genetic markers on chromosome 16q22-Q24 and type 1 diabetes in Russian affected families

Aim To evaluate whether the T1D susceptibility locus on chromosome 16q contributes to the genetic susceptibility to T1D in Russian patients. Method Thirteen microsatellite markers, spanning a 47-centimorgan genomic region on 16q22-q24 were evaluated for linkage to T1D in 98 Russian multiplex families. Multipoint logarithm of odds (LOD) ratio (MLS) and nonparametric LOD (NPL) values were computed for each marker, using GENEHUNTER 2.1 software. Four microsatellites (D16S422, D16S504, D16S3037, and D16S3098) and 6 biallelic markers in 2 positional candidate genes, ICSBP1 and NQO1, were additionally tested for association with T1D in 114 simplex families, using transmission disequilibrium test (TDT). Results A peak of linkage (MLS = 1.35, NPL = 0.91) was shown for marker D16S750, but this was not significant (P = 0.18). The subsequent linkage analysis in the subset of 46 multiplex families carrying a common risk HLA-DR4 haplotype increased peak MLS and NPL values to 1.77 and 1.22, respectively, but showed no significant linkage (P = 0.11) to T1D in the 16q22-q24 genomic region. TDT analysis failed to find significant association between these markers and disease, even after the conditioning for the predisposing HLA-DR4 haplotype. Conclusion Our results did not support the evidence for the susceptibility locus to T1D on chromosome 16q22-24 in the Russian family data set. The lack of association could reflect genetic heterogeneity of type 1 diabetes in diverse ethnic groups.

Characterizing Genetic Drivers of Lymphoma through High-Throughput Sequencing

The advent of next-generation sequencing, now nearing a decade in age, has enabled, among other capabilities, measurement of genome-wide sequence features at unprecedented scale and resolution.

In this dissertation, I describe work to understand the genetic underpinnings of non-Hodgkin’s lymphoma through exploration of the epigenetics of its cell of origin, initial characterization and interpretation of driver mutations, and finally, a larger-scale, population-level study that incorporates mutation interpretation with clinical outcome.

In the first research chapter, I describe genomic characteristics of lymphomas through the lens of their cells of origin. Just as many other cancers, such as breast cancer or lung cancer, are categorized based on their cell of origin, lymphoma subtypes can be examined through the context of their normal B Cells of origin, Naïve, Germinal Center, and post-Germinal Center. By applying integrative analysis of the epigenetics of normal B Cells of origin through chromatin-immunoprecipitation sequencing, we find that differences in normal B Cell subtypes are reflected in the mutational landscapes of the cancers that arise from them, namely Mantle Cell, Burkitt, and Diffuse Large B-Cell Lymphoma.

In the next research chapter, I describe our first endeavor into understanding the genetic heterogeneity of Diffuse Large B Cell Lymphoma, the most common form of non-Hodgkin’s lymphoma, which affects 100,000 patients in the world. Through whole-genome sequencing of 1 case as well as whole-exome sequencing of 94 cases, we characterize the most recurrent genetic features of DLBCL and lay the groundwork for a larger study.

In the last research chapter, I describe work to characterize and interpret the whole exomes of 1001 cases of DLBCL in the largest single-cancer study to date. This highly-powered study enabled sub-gene, gene-level, and gene-network level understanding of driver mutations within DLBCL. Moreover, matched genomic and clinical data enabled the connection of these driver mutations to clinical features such as treatment response or overall survival. As sequencing costs continue to drop, whole-exome sequencing will become a routine clinical assay, and another diagnostic dimension in addition to existing methods such as histology. However, to unlock the full utility of sequencing data, we must be able to interpret it. This study undertakes a first step in developing the understanding necessary to uncover the genomic signals of DLBCL hidden within its exomes. However, beyond the scope of this one disease, the experimental and analytical methods can be readily applied to other cancer sequencing studies.

Thus, this dissertation leverages next-generation sequencing analysis to understand the genetic underpinnings of lymphoma, both by examining its normal cells of origin as well as through a large-scale study to sensitively identify recurrently mutated genes and their relationship to clinical outcome.

Peroxisome Proliferator-Activated Receptor-γ Coactivator 1-α (PPARGC1A) Genetic Associations with Type 2 Diabetes in Three Ethnicities

Genetic heterogeneity, lifestyle factors, gene-gene or gene-environment interactions are the determinants of T2D which puts Hispanics and populations with African ancestry at higher risk of developing T2D. In this dissertation, the genetic associations of PPARGC1A polymorphisms with T2D and its related phenotypes (metabolic markers) in Haitian Americans (cases=110, controls=116), African Americans (cases=120, controls=124) and Cuban Americans (cases=160, controls=181) of South Florida were explored. Five single nucleotide polymorphisms of gene PPARGC1A were evaluated in each ethnicity for their disease association. In Haitian Americans, rs7656250 (OR= 0.22, pp=0.03) had significant protective association with T2D but had risk association in African Americans for rs7656250 (OR=1.02, p=0.96) and rs4235308 (OR=2.53, p=0.03). We found that in Haitian American females, both rs7656250 (OR=0.23, pp=0.03) had protective association with T2D. In African American females, rs7656250 (OR=1.14, p=0.78) had risk association whereas in males, it had significant protective effect (OR=0.37, p=0.04). However, the risk association exhibited by rs4235308 was stronger in African American females (OR=2.69, p=0.03) than males (OR=1.16, p=0.72). In Cuban Americans, only rs7656250 showed significant risk association with T2D (OR=6.87, p=0.02) which was stronger in females alone (OR=7.67, p=0.01). We also observed significant differences among correlations of PPARGC1A SNPs and T2D phenotypes. Positive correlation was observed for log Hs-CRP with rs3774907 (pp=0.03) in Cuban Americans respectively. Correlation of log A1C with rs7656250 (p=0.02) was positive in Cuban Americans while it was negative for rs3774907 in Haitian Americans (ppPPARGC1A correlations with T2D and its phenotypes among the three ethnicities studied (ii) the associations of PPARGC1A SNPs showed significant effect modification by sex. The findings suggest that variations in effects of PPARGC1A gene polymorphisms among three ethnicities and between sexes may have biomedical implications for the development of T2D as well as the phenotypes related to T2D.

Intraclonal heterogeneity and distinct molecular mechanisms characterize the development of t(4;14) and t(11;14) myeloma.

We have used whole exome sequencing to compare a group of presentation t(4;14) with t(11;14) cases of myeloma to define the mutational landscape. Each case was characterized by a median of 24.5 exonic nonsynonymous single-nucleotide variations, and there was a consistently higher number of mutations in the t(4;14) group, but this number did not reach statistical significance. We show that the transition and transversion rates in the 2 subgroups are similar, suggesting that there was no specific mechanism leading to mutation differentiating the 2 groups. Only 3% of mutations were seen in both groups, and recurrently mutated genes include NRAS, KRAS, BRAF, and DIS3 as well as DNAH5, a member of the axonemal dynein family. The pattern of mutation in each group was distinct, with the t(4;14) group being characterized by deregulation of chromatin organization, actin filament, and microfilament movement. Recurrent RAS pathway mutations identified subclonal heterogeneity at a mutational level in both groups, with mutations being present as either dominant or minor subclones. The presence of subclonal diversity was confirmed at a single-cell level using other tumor-acquired mutations. These results are consistent with a distinct molecular pathogenesis underlying each subgroup and have important impacts on targeted treatment strategies. The Medical Research Council Myeloma IX trial is registered under ISRCTN68454111.

Understanding the genetic basis of phenotype variability in individuals with neurocognitive disorders

Thesis (Ph.D.)--University of Washington, 2016-06

Genome-scale analysis of the non-cultivable Treponema pallidum reveals extensive within-patient genetic variation

Insights into the genomic adaptive traits of Treponema pallidum, the causative bacterium of syphilis, have long been hampered due to the absence of in vitro culture models and the constraints associated with its propagation in rabbits. Here, we have bypassed the culture bottleneck by means of a targeted strategy never applied to uncultivable bacterial human pathogens to directly capture whole-genome T. pallidum data in the context of human infection. This strategy has unveiled a scenario of discreet T. pallidum interstrain single-nucleotide-polymorphism-based microevolution, contrasting with a rampant within-patient genetic heterogeneity mainly targeting multiple phase-variable loci and a major antigen-coding gene (tprK). TprK demonstrated remarkable variability and redundancy, intra- and interpatient, suggesting ongoing parallel adaptive diversification during human infection. Some bacterial functions (for example, flagella- and chemotaxis-associated) were systematically targeted by both inter- and intrastrain single nucleotide polymorphisms, as well as by ongoing within-patient phase variation events. Finally, patient-derived genomes possess mutations targeting a penicillin-binding protein coding gene (mrcA) that had never been reported, unveiling it as a candidate target to investigate the impact on the susceptibility to penicillin. Our findings decode the major genetic mechanisms by which T. pallidum promotes immune evasion and survival, and demonstrate the exceptional power of characterizing evolving pathogen subpopulations during human infection.

Assessment of genetic diversity of zoonotic Brucella spp. recovered from livestock in Egypt using multiple locus VNTR analysis

Brucellosis is endemic in most parts of Egypt, where it is caused mainly by Brucella melitensis biovar 3, and affects cattle and small ruminants in spite of ongoing efforts devoted to its control. Knowledge of the predominant Brucella species/strains circulating in a region is a prerequisite of a brucellosis control strategy. For this reason a study aiming at the evaluation of the phenotypic and genetic heterogeneity of a panel of 17 Brucella spp. isolates recovered from domestic ruminants (cattle, buffalo, sheep, and goat) from four governorates during a period of five years (2002-2007) was carried out using microbiological tests and molecular biology techniques (PCR, MLVA-15, and sequencing). Thirteen strains were identified as B. melitensis biovar 3 while all phenotypic and genetic techniques classified the remaining isolates as B. abortus (n = 2) and B. suis biovar 1 (n = 2). MLVA-15 yielded a high discriminatory power (h = 0.801), indicating a high genetic diversity among the B. melitensis strains circulating among domestic ruminants in Egypt. This is the first report of the isolation of B. suis from cattle in Egypt which, coupled with the finding of B. abortus, suggests a potential role of livestock as reservoirs of several zoonotic Brucella species in the region.

Linkage to chromosome 2q32.2-q33.3 in familial serrated neoplasia (Jass syndrome)

Abstract Causative genetic variants have to date been identified for only a small proportion of familial colorectal cancer (CRC). While conditions such as Familial Adenomatous Polyposis and Lynch syndrome have well defined genetic causes, the search for variants underlying the remainder of familial CRC is plagued by genetic heterogeneity. The recent identification of families with a heritable predisposition to malignancies arising through the serrated pathway (familial serrated neoplasia or Jass syndrome) provides an opportunity to study a subset of familial CRC in which heterogeneity may be greatly reduced. A genome-wide linkage screen was performed on a large family displaying a dominantly-inherited predisposition to serrated neoplasia genotyped using the Affymetrix GeneChip Human Mapping 10 K SNP Array. Parametric and nonparametric analyses were performed and resulting regions of interest, as well as previously reported CRC susceptibility loci at 3q22, 7q31 and 9q22, were followed up by finemapping in 10 serrated neoplasia families. Genome-wide linkage analysis revealed regions of interest at 2p25.2-p25.1, 2q24.3-q37.1 and 8p21.2-q12.1. Finemapping linkage and haplotype analyses identified 2q32.2-q33.3 as the region most likely to harbour linkage, with heterogeneity logarithm of the odds (HLOD) 2.09 and nonparametric linkage (NPL) score 2.36 (P = 0.004). Five primary candidate genes (CFLAR, CASP10, CASP8, FZD7 and BMPR2) were sequenced and no segregating variants identified. There was no evidence of linkage to previously reported loci on chromosomes 3, 7 and 9.

Minor head trauma – induced sporadic hemiplegic migraine coma

Familial hemiplegic migraine is a severe, rare subtype of migraine. Gene mutations on chromosome 19 have been identified in the calcium channel, voltage-dependent, P/Q type, alpha-1A subunit gene (chromosome 19p13) for familial hemiplegic migraine. Recently a gene mutation (Serine-218-Leucine) for a dramatic syndrome associated with familial hemiplegic migraine, commonly named “migraine coma”, has implicated exon 5 of this gene. The occurrence of trivial head trauma, in such familial hemiplegic migraine patients, may also be complicated by severe, sometimes even fatal, cerebral edema and coma occurring after a lucid interval. Sporadic hemiplegic migraine shares a similar spectrum of clinical presentation and genetic heterogeneity. The case report presented in this article implicates the involvement of the Serine-218-Leucine mutation in the extremely rare disorder of minor head trauma–induced migraine coma. We conclude that the Serine-218-Leucine mutation in the calcium channel, voltage-dependent, P/Q type, alpha-1A subunit gene is involved in sporadic hemiplegic migraine, delayed cerebral edema and coma after minor head trauma.

DNA probes in Charcot-Marie-Tooth neuropathy

Results of Duffy (Fy) linkage confirm genetic heterogeneity in Charcot-Marie-Tooth disease type 1 (CMT1). Of 11 families informative for Fy, four showed probable linkage with CMT1, seven showed probable non-linkage and two showed definite non-linkage. These results suggest that Fy linked CMT1 may be less common than previously thought. These results combined with those of another DNA probe for the antithrombin III gene confirm that there are at least two gene loci for CMT1, termed 1A and 1B.

Genome-wide association study identifies a variant in HDAC9 associated with large vessel ischemic stroke

Genetic factors have been implicated in stroke risk, but few replicated associations have been reported. We conducted a genome-wide association study (GWAS) for ischemic stroke and its subtypes in 3,548 affected individuals and 5,972 controls, all of European ancestry. Replication of potential signals was performed in 5,859 affected individuals and 6,281 controls. We replicated previous associations for cardioembolic stroke near PITX2 and ZFHX3 and for large vessel stroke at a 9p21 locus. We identified a new association for large vessel stroke within HDAC9 (encoding histone deacetylase 9) on chromosome 7p21.1 (including further replication in an additional 735 affected individuals and 28,583 controls) (rs11984041; combined P = 1.87 × 10 -11; odds ratio (OR) = 1.42, 95% confidence interval (CI) = 1.28-1.57). All four loci exhibited evidence for heterogeneity of effect across the stroke subtypes, with some and possibly all affecting risk for only one subtype. This suggests distinct genetic architectures for different stroke subtypes.

The IFITM5 mutation c.-14C > T results in an elongated transcript expressed in human bone; And causes varying phenotypic severity of osteogenesis imperfecta type V

Background The genetic mutation resulting in osteogenesis imperfecta (OI) type V was recently characterised as a single point mutation (c.-14C > T) in the 5' untranslated region (UTR) of IFITM5, a gene encoding a transmembrane protein with expression restricted to skeletal tissue. This mutation creates an alternative start codon and has been shown in a eukaryotic cell line to result in a longer variant of IFITM5, but its expression has not previously been demonstrated in bone from a patient with OI type V. Methods Sanger sequencing of the IFITM5 5' UTR was performed in our cohort of subjects with a clinical diagnosis of OI type V. Clinical data was collated from referring clinicians. RNA was extracted from a bone sample from one patient and Sanger sequenced to determine expression of wild-type and mutant IFITM5. Results: All nine subjects with OI type V were heterozygous for the c.-14C > T IFITM5 mutation. Clinically, there was heterogeneity in phenotype, particularly in the manifestation of bone fragility amongst subjects. Both wild-type and mutant IFITM5 mRNA transcripts were present in bone. Conclusions The c.-14C > T IFITM5 mutation does not result in an RNA-null allele but is expressed in bone. Individuals with identical mutations in IFITM5 have highly variable phenotypic expression, even within the same family.

The role of germline polymorphisms in the T-cell receptor in susceptibility to ankylosing spondylitis

The role of germline polymorphisms of the T-cell receptor A/D and B loci in susceptibility to ankylosing spondylitis was investigated by linkage studies using microsatellite markers in 215 affected sibling pairs. The presence of a significant susceptibility gene (lambda ≤ 1.6) at the TCRA/D locus was excluded (LOD score < -2.0). At the TCRB locus, there was weak evidence of the presence of a susceptibility gene (P = 0.01, LOD score 1.1). Further family studies will be required to determine whether this is a true or false-positive finding. It is unlikely that either the TCRA/D or TCRB loci contain genes responsible for more than a moderate proportion of the non-MHC genetic susceptibility to ankylosing spondylitis.

Whole exome sequencing is an efficient and sensitive method for detection of germline mutations in patients with phaeochromcytomas and paragangliomas

Background Genetic testing is recommended when the probability of a disease-associated germline mutation exceeds 10%. Germline mutations are found in approximately 25% of individuals with phaeochromcytoma (PCC) or paraganglioma (PGL); however, genetic heterogeneity for PCC/PGL means many genes may require sequencing. A phenotype-directed iterative approach may limit costs but may also delay diagnosis, and will not detect mutations in genes not previously associated with PCC/PGL. Objective To assess whether whole exome sequencing (WES) was efficient and sensitive for mutation detection in PCC/PGL. Methods Whole exome sequencing was performed on blinded samples from eleven individuals with PCC/PGL and known mutations. Illumina TruSeq™ (Illumina Inc, San Diego, CA, USA) was used for exome capture of seven samples, and NimbleGen SeqCap EZ v3.0 (Roche NimbleGen Inc, Basel, Switzerland) for five samples (one sample was repeated). Massive parallel sequencing was performed on multiplexed samples. Sequencing data were called using Genome Analysis Toolkit and annotated using annovar. Data were assessed for coding variants in RET, NF1, VHL, SDHD, SDHB, SDHC, SDHA, SDHAF2, KIF1B, TMEM127, EGLN1 and MAX. Target capture of five exome capture platforms was compared. Results Six of seven mutations were detected using Illumina TruSeq™ exome capture. All five mutations were detected using NimbleGen SeqCap EZ v3.0 platform, including the mutation missed using Illumina TruSeq™ capture. Target capture for exons in known PCC/PGL genes differs substantially between platforms. Exome sequencing was inexpensive (<$A800 per sample for reagents) and rapid (results <5 weeks from sample reception). Conclusion Whole exome sequencing is sensitive, rapid and efficient for detection of PCC/PGL germline mutations. However, capture platform selection is critical to maximize sensitivity.