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.

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.

Estimation of (co)variance components and genetic parameters for weights of red-winged tinamou using random regression models

Com este trabalho objetivou-se determinar parâmetros genéticos para peso corporal de perdizes em cativeiro. Foram utilizados modelos de regressão aleatória na análise dos dados considerando os efeitos genéticos aditivos diretos (AD) e de ambiente permanente de animal (AP) como aleatórios. As variâncias residuais foram modeladas utilizando-se funções de variância de ordem 5. A curva média da população foi ajustada por polinômios ortogonais de Legendre de ordem 6. Os efeitos genéticos aditivos diretos e de ambiente permanente de animal foram modelados utilizando-se polinômios de Legendre de segunda a nona ordem. Os melhores resultados foram obtidos pelos modelos de ordem 6 de ajuste para os efeitos genéticos aditivos diretos e de ordem 3 para os de ambiente permanente pelo Critério de Informação de Akaike e ordem 3 para ambos os efeitos pelos Critério de Informação Bayesiano de Schwartz e Teste de Razão de Verossimilhança. As herdabilidades estimadas variaram de 0,02 a 0,57. O primeiro autovalor respondeu por 94 e 90% da variação decorrente de efeitos aditivos diretos e de ambiente permanente, respectivamente. A seleção de perdizes para peso corporal é mais efetiva a partir de 112 dias de idade.

Female fertility in a Guzerat dairy subpopulation: Heterogeneity of variance components for calving intervals

The objectives of the present study were to determine if variance components of calving intervals varied with age at calving and if considering calving intervals as a longitudinal trait would be a useful approach for fertility analysis of Zebu dairy herds. With these purposes, calving records from females born from 1940 to 2006 in a Guzerat dairy subpopulation in Brazil were analyzed. The fixed effects of contemporary groups, formed by year and farm at birth or at calving, and the regressions of age at calving, equivalent inbreeding coefficient and day of the year on the studied traits were considered in the statistical models. In one approach, calving intervals (Cl) were analyzed as a single trait, by fitting a statistical model on which both animal and permanent environment effects were adjusted for the effect of age at calving by random regression. In a second approach, a four-trait analysis was conducted, including age at first calving (AFC) and three different female categories for the calving intervals: first calving females; young females (less than 80 months old, but not first calving); or mature females (80 months old or more). Finally, a two-trait analysis was performed, also including AFC and Cl, but calving intervals were regarded as a single trait in a repeatability model. Additionally, the ranking of sires was compared among approaches. Calving intervals decreased with age until females were about 80 months old, remaining nearly constant after that age. A quasi-linear increase of 11.5 days on the calving intervals was observed for each 10% increase in the female's equivalent inbreeding coefficient. The heritability of AFC was 0.37. For Cl. the genetic-phenotypic variance ratios ranged from 0.064 to 0.141, depending on the approach and on ages at calving. Differences among genetic variance components for calving intervals were observed along the animal's lifetime. Those differences confirmed the longitudinal aspect of that trait, indicating the importance of such consideration when accessing fertility of Zebu dairy females, especially in situations where the available information relies on their calving intervals. Spearman rank correlations among approaches ranged from 0.90 to 0.95, and changes observed in the ranking of sires suggested that the genetic progress of the population could be affected by the approach chosen for the analysis of calving intervals. (C) 2012 Elsevier ay. All rights reserved.

Leber Congenital Amaurosis: Comprehensive Survey of the Genetic Heterogeneity, Refinement of the Clinical Definition, and Genotype-Phenotype Correlations as a Strategy for Molecular Diagnosis

Leber congenital amaurosis (LCA) is the earliest and most severe form of all inherited retinal dystrophies, responsible for congenital blindness. Disease-associated mutations have been hitherto reported in seven genes. These genes are all expressed preferentially in the photoreceptor cells or the retinal pigment epithelium but they are involved in strikingly different physiologic pathways resulting in an unforeseeable physiopathologic variety. This wide genetic and physiologic heterogeneity that could largely increase in the coming years, hinders the molecular diagnosis in LCA patients. The genotyping is, however, required to establish genetically defined subgroups of patients ready for therapy. Here, we report a comprehensive mutational analysis of the all known genes in 179 unrelated LCA patients, including 52 familial and 127 sporadic (27/127 consanguineous) cases. Mutations were identified in 47.5% patients. GUCY2D appeared to account for most LCA cases of our series (21.2%), followed by CRB1 (10%), RPE65 (6.1%), RPGRIP1 (4.5%), AIPL1 (3.4%), TULP1 (1.7%), and CRX (0.6%). The clinical history of all patients with mutations was carefully revisited to search for phenotype variations. Sound genotype-phenotype correlations were found that allowed us to divide patients into two main groups. The first one includes patients whose symptoms fit the traditional definition of LCA, i.e., congenital or very early cone-rod dystrophy, while the second group gathers patients affected with severe yet progressive rod-cone dystrophy. Besides, objective ophthalmologic data allowed us to subdivide each group into two subtypes. Based on these findings, we have drawn decisional flowcharts directing the molecular analysis of LCA genes in a given case. These flowcharts will hopefully lighten the heavy task of genotyping new patients but only if one has access to the most precise clinical history since birth.

Identification of C7orf11 (TTDN1) gene mutations and genetic heterogeneity in nonphotosensitive trichothiodystrophy.

We have identified C7orf11, which localizes to the nucleus and is expressed in fetal hair follicles, as the first disease gene for nonphotosensitive trichothiodystrophy (TTD). C7orf11 maps to chromosome 7p14, and the disease locus has been designated "TTDN1" (TTD nonphotosensitive 1). Mutations were found in patients with Amish brittle-hair syndrome and in other nonphotosensititive TTD cases with mental retardation and decreased fertility but not in patients with Sabinas syndrome or Pollitt syndrome. Therefore, genetic heterogeneity in nonphotosensitive TTD is a feature similar to that observed in photosensitive TTD, which is caused by mutations in transcription factor II H (TFIIH) subunit genes. Comparative immunofluorescence analysis, however, suggests that C7orf11 does not influence TFIIH directly. Given the absence of cutaneous photosensitivity in the patients with C7orf11 mutations, together with the protein's nuclear localization, C7orf11 may be involved in transcription but not DNA repair.

European collaborative study of early-onset bipolar disorder: Evidence for genetic heterogeneity on 2q14 according to age at onset.

Bipolar disorder has a genetic component, but the mode of inheritance remains unclear. A previous genome scan conducted in 70 European families led to detect eight regions linked to bipolar disease. Here, we present an investigation of whether the phenotypic heterogeneity of the disorder corresponds to genetic heterogeneity in these regions using additional markers and an extended sample of families. The MLS statistic was used for linkage analyses. The predivided sample test and the maximum likelihood binomial methods were used to test genetic homogeneity between early-onset bipolar type I (cut-off of 22 years) and other types of the disorder (later onset of bipolar type I and early-onset bipolar type II), using a total of 138 independent bipolar-affected sib-pairs. Analysis of the extended sample of families supports linkage in four regions (2q14, 3p14, 16p23, and 20p12) of the eight regions of linkage suggested by our previous genome scan. Heterogeneity testing revealed genetic heterogeneity between early and late-onset bipolar type I in the 2q14 region (P = 0.0001). Only the early form of the bipolar disorder but not the late form appeared to be linked to this region. This region may therefore include a genetic factor either specifically involved in the early-onset bipolar type I or only influencing the age at onset (AAO). Our findings illustrate that stratification according to AAO may be valuable for the identification of genetic vulnerability polymorphisms.

Genetic evidence of heterogeneity in intrahepatic cholestasis of pregnancy

Background and aims: The aim of this study was to investigate the genetic aetiology of intrahepatic cholestasis of pregnancy (ICP) and the impact of known cholestasis genes (BSEP, FIC1, and MDR3) on the development of this disease. Patients and methods: Sixty nine Finnish ICP patients were prospectively interviewed for a family history of ICP, and clinical features were compared in patients with familial ICP (patients with a positive family history, n=11) and sporadic patients (patients with no known family history of ICP, n=58). For molecular genetic analysis, 16 individuals from two independently ascertained Finnish ICP families were genotyped for the flanking markers for BSEP, FIC1, and MDR3. Results: The pedigree structures in 16% (11/69) of patients suggested dominant inheritance. Patients with familial ICP had higher serum aminotransferase levels and a higher recurrence risk (92% v 40%). Both segregation of haplotypes and multipoint linkage analysis excluded BSEP, FIC1, and MDR3 genes in the studied pedigrees. Additionally, the MDR3 gene, previously shown to harbour mutations in ICP patients, was negative for mutations when sequenced in four affected individuals from the two families. Conclusions: These results support the hypothesis that the aetiology of ICP is heterogeneous and that ICP is due to a genetic predisposition in a proportion of patients. The results of molecular genetic analysis further suggest that the previously identified three cholestasis genes are not likely to be implicated in these Finnish ICP families with dominant inheritance.

The interplay of spatial and climatic landscapes in the genetic distribution of a South American parrot

AimOur aim was to understand the interplay of heterogeneous climatic and spatial landscapes in shaping the distribution of nuclear microsatellite variation in burrowing parrots, Cyanoliseus patagonus. Given the marked phenotypic differences between populations of burrowing parrots we hypothesized an important role of geographical as well climatic heterogeneity in the population structure of this species. LocationSouthern South America. MethodsWe applied a landscape genetics approach to investigate the explicit patterns of genetic spatial autocorrelation based on both geography and climate using spatial principal component analysis (sPCA). This necessitated a novel statistical estimation of the species climatic landscape, considering temperature- and precipitation-based variables separately to evaluate their weight in shaping the distribution of genetic variation in our model system. ResultsGeographical and climatic heterogeneity successfully explained molecular variance in burrowing parrots. sPCA divided the species distribution into two main areas, Patagonia and the pre-Andes, which were connected by an area of geographical and climatic transition. Moreover, sPCA revealed cryptic and conservation-relevant genetic structure: the pre-Andean populations and the transition localities were each divided into two groups, each management units for conservation. Main conclusionssPCA, a method originally developed for spatial genetics, allowed us to unravel the genetic structure related to spatial and climatic landscapes and to visualize these patterns in landscape space. These novel climatic inferences underscore the importance of our modified sPCA approach in revealing how climatic variables can drive cryptic patterns of genetic structure, making the approach potentially useful in the study of any species distributed over a climatically heterogeneous landscape.

Genetic determinants of clinical heterogeneity in sickle cell disease

L’anémie falciforme est une maladie monogénique causée par une mutation dans le locus de la β-globine. Malgré le fait que l’anémie falciforme soit une maladie monogénique, cette maladie présente une grande hétérogénéité clinique. On présume que des facteurs environnementaux et génétiques contribuent à cette hétérogénéité. Il a été observé qu’un haut taux d’hémoglobine fœtale (HbF) diminuait la sévérité et la mortalité des patients atteints de l’anémie falciforme. Le but de mon projet était d’identifier des variations génétiques modifiant la sévérité clinique de l’anémie falciforme. Dans un premier temps, nous avons effectué la cartographie-fine de trois régions précédemment associées avec le taux d’hémoglobine fœtale. Nous avons ensuite effectué des études d’association pan-génomiques avec deux complications cliniques de l’anémie falciforme ainsi qu’avec le taux d’hémoglobine fœtale. Hormis les régions déjà identifiées comme étant associées au taux d’hémoglobine fœtale, aucun locus n’a atteint le niveau significatif de la puce de génotypage. Pour identifier des groupes de gènes modérément associés au taux d’hémoglobine fœtale qui seraient impliqués dans de mêmes voies biologiques, nous avons effectué une étude des processus biologiques. Finalement, nous avons effectué l’analyse de 19 exomes de patients Jamaïcains ayant des complications cliniques mineures de l’anémie falciforme. Compte tenu de la taille des cohortes de réplication disponibles, nous n’avons pas les moyens de valider statistiquement les variations identifiées par notre étude. Cependant, nos résultats fournissent de bons gènes candidats pour des études fonctionnelles et pour les réplications futures. Nos résultats suggèrent aussi que le β-hydroxybutyrate en concentration endogène pourraient influencer le taux d’hémoglobine fœtale. De plus, nous montrons que la cartographie-fine des régions associées par des études pan-génomiques peut identifier des signaux d’association additionnels et augmenter la variation héritable expliquée par cette région.

Genetic heterogeneity in erythropoietic protoporphyria: a study of the enzymatic defect in nine affected families

Erythropoietic protoporphyria (EPP) is associated with a deficiency of protohaem ferrolyase. We have used a novel assay for this enzyme based on its ability to utilize zinc as a substrate to investigate the inheritance of EPP in nine affected families. Zinc chelatase activity was markedly reduced in peripheral blood mononuclear cells from 14 EPP patients (mean, 3.3 nmol Zn protohaem/h/mg protein; range, 0.3-8.0) when compared with 41 controls (16.8 +/- 3.6) p less than 0.01. In three families with parent-to-child transmission of disease, the asymptomatic parent had an enzymatic activity within the normal range. In three pedigrees where the parents were asymptomatic, enzymatic activities were below the 95% confidence limits in both. Zinc chelatase activity was below the mean control value in 17 of the 18 parents in nine affected pedigrees, and six of seven asymptomatic offspring of patients with protoporphyria. The findings suggest that EPP is not transmitted as a simple dominant trait and that inheritance of more than one gene may be required for disease expression.

Auriculo-condylar syndrome: mapping of a first locus and evidence for genetic heterogeneity

Auriculo-condylar syndrome (ACS), an autosomal dominant disorder of first and second pharyngeal arches, is characterized by malformed ears (`question mark ears`), prominent cheeks, microstomia, abnormal temporomandibular joint, and mandibular condyle hypoplasia. Penetrance seems to be complete, but there is high inter-and intra-familial phenotypic variation, with no evidence of genetic heterogeneity. We herein describe a new multigeneration family with 11 affected individuals (F1), in whom we confirm intra-familial clinical variability. Facial asymmetry, a clinical feature not highlighted in other ACS reports, was highly prevalent among the patients reported here. The gene responsible for ACS is still unknown and its identification will certainly contribute to the understanding of human craniofacial development. No chromosomal rearrangements have been associated with ACS, thus mapping and positional cloning is the best approach to identify this disease gene. To map the ACS gene, we conducted linkage analysis in two large ACS families, F1 and F2 (F2; reported elsewhere). Through segregation analysis, we first excluded three known loci associated with disorders of first and second pharyngeal arches (Treacher Collins syndrome, oculo-auriculo-vertebral spectrum, and Townes-Brocks syndrome). Next, we performed a wide genome search and we observed evidence of linkage to 1p21.1-q23.3 in F2 (LOD max 3.01 at theta = 0). Interestingly, this locus was not linked to the phenotype segregating in F1. Therefore, our results led to the mapping of a first locus of ACS (ACS1) and also showed evidence for genetic heterogeneity, suggesting that there are at least two loci responsible for this phenotype.

Random regressions models to describe the genetic variation of milk yield over multiple parities in Buffaloes

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estimation of genetic parameters for milk yield in Murrah buffaloes by Bayesian inference

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Genetic parameters of total milk yield and factors describing the shape of lactation curve in dairy buffaloes

The objective of this study was to apply factor analysis to describe lactation curves in dairy buffaloes in order to estimate the phenotypic and genetic association between common latent factors and cumulative milk yield. A total of 31 257 monthly test-day milk yield records from buffaloes belonging to herds located in the state of São Paulo were used to estimate two common latent factors, which were then analysed in a multi-trait animal model for estimating genetic parameters. Estimates of (co)variance components for the two common latent factors and cumulated 270-d milk yield were obtained by Bayesian inference using a multiple trait animal model. Contemporary group, number of milkings per day (two levels) and age of buffalo cow at calving (linear and quadratic) as covariate were included in the model as fixed effects. The additive genetic, permanent environmental and residual effects were included as random effects. The first common latent factor (F1) was associated with persistency of lactation and the second common latent factor (F2) with the level of production in early lactation. Heritability estimates for Fl and F2 were 0.12 and 0.07, respectively. Genetic correlation estimates between El and F2 with cumulative milk yield were positive and moderate (0.63 and 0.52). Multivariate statistics employing factor analysis allowed the extraction of two variables (latent factors) that described the shape of the lactation curve. It is expected that the response to selection to increase lactation persistency is higher than the response obtained from selecting animals to increase lactation peak. Selection for higher total milk yield would result in a favourable correlated response to increase the level of production in early lactation and the lactation persistency.