The intestinal mucosa is made up of complex secreted mucus layer consist of mainly mucin 2 (MUC2) and antimicrobial components that defend the underlining cellular barrier from intrusion by luminal microbiota and toxins. In inflammatory bowel diseases (IBD), the mucosal integrity is compromised. This can result from a combination of altered host genetics, gut immune responses and environment factors. However, it is the presence of intestinal bacteria that is central to the pathogenesis of IBD. As part of the dynamic gut microbial flora, mucolytic bacteria produce a wide range of glycosidases that are able to remove the outer oligosaccharide chains of MUC2, which allow other luminal bacteria to further degrade the mucin. We hypothesised that increased mucolytic bacteria will cause excessive degradation of the mucus layer, which in turn, allow more luminal bacteria to be in close proximity to the underlining epithelial cells resulting in inflammation. Consistent with our group’s previous semi-quantitative bacterial 16S rRNA gene clone library analysis, we found increased Ruminococcus gnavus in non-inflamed ulcerative colitis (UC) mucosa. R. gnavus was previously isolated by others based on its mucolytic property. In this study, we quantify total mucosa-associated bacteria and mucolytic bacteria, namely, R. gnavus, R. torques, Akkermansia muciniphila and bifidobacteria. We were able to show quantitatively that total mucosa-associated bacteria were increased in IBD. There was also a population shift in the mucosa-associated mucolytic bacteria, which were increased overall. There was significantly more R. gnavus in non-inflamed IBD biopsies. For the first time, we were also able to demonstrate that R. gnavus can degrade human MUC2 in vitro. To examine whether the numerical association of R. gnavus in IBD does have functional influence on intestinal inflammation and Paneth cell antimicrobial peptide gene expression, we fed mice with R. gnavus. Interestingly, R. gnavus feeding did not result in histological or molecular evidence of gut inflammation; however, it was able to specifically induce Paneth cell cryptdins and lysozyme P genes expression in 3 week old, antibiotic pre-treated C57BL/6 mice. This demonstrated that R. gnavus is not a pathogenic bacterium, which will directly cause colitis. However, the increased Paneth cell response suggested the need for host innate defence when R. gnavus is increased. Other than bacterial degradation, altered host genetics will also influence the mucus barrier. There is evidence to suggest that the MUC2 gene is highly unstable and is susceptible to gene copy number variation (CNV). Therefore, we hypothesised that MUC2 CNV is present, which may result in altered oligomerisation of the MUC2 glycoprotein causing endoplasmic reticulum stress of the goblet cells that appears to be characteristic of UC. Currently, our data partly support the presence of MUC2 CNV. However, further investigation is required to verify the MUC2 CNV identity. Only then can a high throughput methodology be designed to screen a large population for any association with IBD.

The human genome is a dynamic structure, recently recognized to present with significant large-scale structural variation. DNA-copy number changes represent one common type of such variation and is found both between individuals and within the somatic cells of the same individual, especially in disease states like cancer.  Apart from DNA-rearrangements, epigenomic changes are increasingly acknowledged as important events in the maintenance of genomic integrity. In this thesis, different array-based methods have been applied for global genomic and epigenomic profiling of both normal and cancer cells. In paper I, a genomic microarray was established and used to determine DNA-copy number variants (CNVs) in a cohort of 76 healthy individuals from three ethnic populations. We identified 315 CNV regions that in total encompassed ~3,5% of the genome. In paper II, the array was utilized to discover CNVs within several differentiated tissues from the same subject. Six variants were identified providing evidence for somatic mosaicism. In paper III and IV we studied pheochromocytomas and paragangliomas, rare endocrine tumours that most often present as benign and sporadic with unclear genetic/epigenetic cause. Genome-wide DNA-copy number analysis of 53 benign and malignant samples in paper III revealed numerous common and novel chromosomal regions of losses and gains. High frequencies of relatively small overlapping regions of deletions were detected on chromosome 1p arm, encompassing several candidate tumour suppressor genes. In paper IV, an epigenomic map for two histone modifications associated with silent (H3K27me3) or active (H3K4me3) gene transcription, was generated for one malignant pheochromocytoma. Integrated analysis of global histone methylation, copy number alterations and gene expression data aided in the identification of candidate tumour genes. In conclusion, the performed studies have contributed to gain knowledge of CNVs in healthy individuals, and identified regions and genes which are likely associated with the development and progression of pheochromocytoma/paraganglioma.

P1 prophage is stably maintained in E. coli as a low-copy-number plasmid. Stable maintenance of P1 is dependent on the function of the plasmid encoded partition system, parABS. ParA is the partition ATPase, ParB is the partition-site binding protein, and parS is the partition site. The concerted action of these proteins results in dynamic movement of the plasmid over the bacterial nucleoid, which results in its stable maintenance. Plasmid movement has been proposed to be caused by interactions between parS bound ParB and nucleoid bound ParA. In this thesis, I have identified a complex of ParA, ParB, and DNA that is capable of promoting plasmid stability. ParA, ParB, DNA interactions required the ATP bound conformation of ParA. The ParA-ParB-DNA complex was dynamically regulated by nucleotide hydrolysis, which promoted complex disassembly. Complex formation resulted from the cooperative binding of ParA and ParB to DNA. ParA-ParB and ParB-DNA interactions were both necessary for complex formation. ParA-ParB-DNA complex size was regulated by ParB stimulation of ParA-ATP hydrolysis. Microscopy demonstrated that complexes resulted in the association of multiple DNA molecules due to protein binding. The properties of complex assembly, dynamics, and DNA grouping lead me to propose a model where associations between ParA bound to the bacterial nucleoid and the partition complex mediated plasmid movement and localization.

P1 prophage is stably maintained in E. coli as a low-copy-number plasmid. Stable maintenance of P1 is dependent on the function of the plasmid encoded partition system, parABS. ParA is the partition ATPase, ParB is the partition-site binding protein, and parS is the partition site. The concerted action of these proteins results in dynamic movement of the plasmid over the bacterial nucleoid, which results in its stable maintenance. Plasmid movement has been proposed to be caused by interactions between parS bound ParB and nucleoid bound ParA. In this thesis, I have identified a complex of ParA, ParB, and DNA that is capable of promoting plasmid stability. ParA, ParB, DNA interactions required the ATP bound conformation of ParA. The ParA-ParB-DNA complex was dynamically regulated by nucleotide hydrolysis, which promoted complex disassembly. Complex formation resulted from the cooperative binding of ParA and ParB to DNA. ParA-ParB and ParB-DNA interactions were both necessary for complex formation. ParA-ParB-DNA complex size was regulated by ParB stimulation of ParA-ATP hydrolysis. Microscopy demonstrated that complexes resulted in the association of multiple DNA molecules due to protein binding. The properties of complex assembly, dynamics, and DNA grouping lead me to propose a model where associations between ParA bound to the bacterial nucleoid and the partition complex mediated plasmid movement and localization.

Eukaryotic Translation Elongation Factor 1 Alpha (EEF1A) exists as two forms with different tissue specificities and encoded by separate loci: eEF1A1 on 6q13 and eEF1A2 on 20q13.3. eEF1A1 is ubiquitously expressed whereas eEF1A2 expression is normally limited to the heart, brain and skeletal muscles. eEF1A proteins are GTP-binding proteins that recruit an amino-acylated tRNA to the ribosome during the elongation phase of protein translation. eEF1A2 mRNA and protein are highly expressed in 50–60% of primary human breast tumors and metastases but not in normal breast epithelium. Since it is also overexpressed in 30% of primary human ovarian tumors, transforms rodent fibroblasts and increases their tumorigenicity in nude mice, eEF1A2 is considered to be a potential human oncogene. The mechanism of eEF1A2 expression is yet to be determined. Studies showed no gene mutation and no correlation between locus amplification or methylation and gene expression. Phosphate Tyrosine Kinase-6 (PTK6) is also located on 20q13.3. It is 48kb upstream of EEF1A2. PTK6 is a non-receptor tyrosine-kinase that is normally expressed in epithelial linings, prostate, skin and oral epithelium but it is not detected in the normal human mammary epithelium. PTK6 has been found to be expressed in many breast cancer cell lines and in approximately 60% of primary human breast tumors but it has not been detected in normal human breast tissue nor in fibroadenomas. Like other tyrosine kinases, PTK6 phosphorylates and activates downstream substrates that would be predicted to lead to increased transcriptional activity and therefore mediates proliferation of breast cancer cells. PTK6 is considered a prognostic marker of metastasis-free survival in breast cancer independent of the classical markers of tumor size, lymph node involvement and HER2 status. The aim of this project was to characterize for the first time the genomic region containing the two mentioned breast cancer oncogenes and understand their various roleswhether they act in tandem or independently in breast tumorigenesis. Immunohistochemistry was performed on tissue microarrays from 300 breast cancer patients to detect the expression levels of eEF1A2 and PTK6. Tumors that showed a high co-expression were analyzed for the genes’ copy number. An increased copy number of PTK6 was detected but not of eEF1A2 nor of adjacent genes on the 20q13.3 amplicon. To understand the effect of EEF1A2 expression on other genes, microarray analysis was performed on NIH-3T3 cells stably transfected with EEF1A2. Many upregulated genes were associated with different types of cancers. This was further confirmed by real-time PCR. However, when the NIH-3T3 cells were transiently transfected with EEF1A2, the genes that were upregulated in the microarray study showed no change in expression. In conclusion, EEF1A2 and PTK6 act independently and each acts through a different mechanism in breast tumorigenesis.

Aberrationen in der DNA-Kopienzahl sind häufige genetische Veränderungen bei malignen Lymphomerkrankungen. Zugewinne sowie Deletionen stellen dabei Mechanismen zur Onkogen-Aktivierung sowie Tumorsuppressorgen-Inaktivierung dar und tragen somit zur Pathogenese der Erkrankung bei. Array-CGH und SNP-Array sind Messplattformen, die die genomweite Bestimmung von Kopienzahlaberrationen in einem Experiment ermöglichen. Die bei der Analyse entstehenden Datensätze sind komplex und erfordern automatische Methoden zur Unterstützung der Analyse und Interpretation der Messergebnisse. In dieser Promotionsarbeit wurden Methoden entwickelt, welche die Analyse von Array-CGH- und SNP-Array-Messungen ermöglichen. Diese Methoden wurden für die Auswertung umfangreicher Datensätze von malignen Non-Hodgkin-Lymphomen verwendet. Dabei wurden Lymphome der Entitäten Burkitt-Lymphom, diffus großzelliges B-Zell-Lymphom, Mantelzelllymphom, primäres ZNS-Lymphom und peripheres T-Zell-Lymphom – nicht anderweitig spezifiziert – analysiert. Für die untersuchten Lymphom-Entitäten konnten hierbei zahlreiche neue rekurrente Kopienzahlaberrationen sowie uniparentale Disomien gezeigt werden, die neue Einblicke in die Pathogenese der jeweiligen Erkrankungen erlauben. Darüber hinaus erfolgte ein Vergleich beider Messplattformen anhand eines Datensatzes mit gepaarten Array-CGH- und SNP-Array-Daten. Für die eingesetzten Plattformen (2800k-BAC-Array vs. Affymetrix 250k-Sty-SNP-Array) konnte eine circa zwölffach höhere effektive Auflösung der SNP-Array-Plattform gezeigt werden. Die wesentlichen Ergebnisse dieser Arbeit sind in sieben Publikationen eingeflossen.

This dissertation brings to the fore the “Characterization of the SULT1A1 polymorphism in a South Africa Tswana population group.” The primary experimental group studied came from South African homogeneous Tswana individuals who participated voluntarily in an ongoing large-scale epidemiological Prospective Urban and Rural Epidemiological (PURE) study the North-West University (Potchefstroom Campus) participates in, as one of the 16 low- middleand high-income countries across the world. The primary aspect investigated was the comprehensive profile of the single nucleotide polymorphism (SNP) and copy number variation (CNP) of the SULT1A1 gene. Using the PCRbased RFLP method, SULT1A1 genotypes, and allele frequency distributions in an experimental group of 1 867 individuals were determined. According to the literature this is by far the largest and most homogeneous group from which such information has been acquired to date. The SULT1A1*1, SULT1A1*1/*2 and SULT1A1*2 genotypes were found to be present at a percentage of 43.76, 47.12 and 9.11 respectively. In comparison to similar studies in other population groups, results from this study indicate that there are ethnic differences in the SULT1A1 genotypes incidence. Asian group differs from Caucasian and Tswana groups because of its exceptionally high prevalence of individuals with the SULT1A1*1 genotype and a very low incidence of the SULT1A1*2 genotype. The SULT1A1*1 genotype profiles of Caucasian and Tswana groups were comparable, but notable differences were observed for the SULT1A1*2 genotype. Using a quantitative multiplex PCR method for the CNV study, the numbers of copies of the SULT1A1 gene in the Tswana population were determined, and the results showed 1 to ~5 copies: only 0.65% of the subjects had a single copy, whereas 59.69% of the subjects had 3 or more copies. This result shows a significant discrepancy between the Caucasian-American samples, which showed that only 26% from that group had more than three copies. However, there is a significant relationship with the African-American population, which presented 63% with 3 or more copies. This finding confirms results from a much smaller African-American study, and suggests a possible genetic link between the African Tswana and the heritage of the African-Americans. These findings were submitted for publication to the South African Journal of Science, as that journal specializes in publication of new knowledge that has a regional focus on Africa. Simultaneous phenotypic consequences of the SNP and CNP of the SULT1A1 gene, as well as the thermo-stable and thermo-labile forms of the sulfotransferases were determined. For this, the formation of [35S]-4-nitrophenyl sulphate from 4-nitrophenol and [35S]-3’-phosphoadenosine- 5’-phosphosulfate ([35S]-PAPS) in platelet homogenates were measured, with the data normalized to a common platelet count. This investigation required fresh blood for enzyme activity. These samples came from 98 Caucasian subjects who voluntarily participated in this part of the study. The experimental data presented a unique challenge to develop a statistical model to accommodate the complexity of the distribution of the data in the phenotype and genotype components, which could be achieved by the development of a mixed model. The model indicated that product formation increased through increasing copy number, but did not differ for SULT1A1*1 and SULT1A1*1/*2. However, the rate of increase in product for the thermo-stable forms of the SULTs was greater than that of thermo-labile forms. In contrast, copy number effect for SULT1A1*2 differed considerably from that of the other two genotypes. Since genotype is also a significant factor, it was concluded from Tukey post-hoc tests that the population group means for product formation differ significantly (for all levels). These results are presently being prepared for publication in an accredited international journal. Finally, perturbations in 23 biochemical parameters measured in the PURE study were analyzed as a function of the SULT1A1 SNP and CNP were evaluated. No group separation in this regard could be found. It could be shown however, that sulfonation of the iodothyronines, which are endogenous substrates for the SULTs, was influenced by the SULT1A1 genotype. The relative concentrations in plasma of the sulphonated iodothyronines may be expressed as T2S > T3S >> T4S, which coincides with the substrate preference of the SULT1A1 enzymes. This observation may, however, only be qualitatively interpreted as (1) the targeted metabolomics mass spectrometric method used for the quantitative analysis of these substances needs further development, and (2) the influence of deiodonation was not taken into account in these studies. In conclusion, three perspectives are given at the end of the thesis which might be considered for further investigations. / Thesis (Ph.D. (Biochemistry))--North-West University, Potchefstroom Campus, 2010.

This dissertation brings to the fore the “Characterization of the SULT1A1 polymorphism in a South Africa Tswana population group.” The primary experimental group studied came from South African homogeneous Tswana individuals who participated voluntarily in an ongoing large-scale epidemiological Prospective Urban and Rural Epidemiological (PURE) study the North-West University (Potchefstroom Campus) participates in, as one of the 16 low- middleand high-income countries across the world. The primary aspect investigated was the comprehensive profile of the single nucleotide polymorphism (SNP) and copy number variation (CNP) of the SULT1A1 gene. Using the PCRbased RFLP method, SULT1A1 genotypes, and allele frequency distributions in an experimental group of 1 867 individuals were determined. According to the literature this is by far the largest and most homogeneous group from which such information has been acquired to date. The SULT1A1*1, SULT1A1*1/*2 and SULT1A1*2 genotypes were found to be present at a percentage of 43.76, 47.12 and 9.11 respectively. In comparison to similar studies in other population groups, results from this study indicate that there are ethnic differences in the SULT1A1 genotypes incidence. Asian group differs from Caucasian and Tswana groups because of its exceptionally high prevalence of individuals with the SULT1A1*1 genotype and a very low incidence of the SULT1A1*2 genotype. The SULT1A1*1 genotype profiles of Caucasian and Tswana groups were comparable, but notable differences were observed for the SULT1A1*2 genotype. Using a quantitative multiplex PCR method for the CNV study, the numbers of copies of the SULT1A1 gene in the Tswana population were determined, and the results showed 1 to ~5 copies: only 0.65% of the subjects had a single copy, whereas 59.69% of the subjects had 3 or more copies. This result shows a significant discrepancy between the Caucasian-American samples, which showed that only 26% from that group had more than three copies. However, there is a significant relationship with the African-American population, which presented 63% with 3 or more copies. This finding confirms results from a much smaller African-American study, and suggests a possible genetic link between the African Tswana and the heritage of the African-Americans. These findings were submitted for publication to the South African Journal of Science, as that journal specializes in publication of new knowledge that has a regional focus on Africa. Simultaneous phenotypic consequences of the SNP and CNP of the SULT1A1 gene, as well as the thermo-stable and thermo-labile forms of the sulfotransferases were determined. For this, the formation of [35S]-4-nitrophenyl sulphate from 4-nitrophenol and [35S]-3’-phosphoadenosine- 5’-phosphosulfate ([35S]-PAPS) in platelet homogenates were measured, with the data normalized to a common platelet count. This investigation required fresh blood for enzyme activity. These samples came from 98 Caucasian subjects who voluntarily participated in this part of the study. The experimental data presented a unique challenge to develop a statistical model to accommodate the complexity of the distribution of the data in the phenotype and genotype components, which could be achieved by the development of a mixed model. The model indicated that product formation increased through increasing copy number, but did not differ for SULT1A1*1 and SULT1A1*1/*2. However, the rate of increase in product for the thermo-stable forms of the SULTs was greater than that of thermo-labile forms. In contrast, copy number effect for SULT1A1*2 differed considerably from that of the other two genotypes. Since genotype is also a significant factor, it was concluded from Tukey post-hoc tests that the population group means for product formation differ significantly (for all levels). These results are presently being prepared for publication in an accredited international journal. Finally, perturbations in 23 biochemical parameters measured in the PURE study were analyzed as a function of the SULT1A1 SNP and CNP were evaluated. No group separation in this regard could be found. It could be shown however, that sulfonation of the iodothyronines, which are endogenous substrates for the SULTs, was influenced by the SULT1A1 genotype. The relative concentrations in plasma of the sulphonated iodothyronines may be expressed as T2S > T3S >> T4S, which coincides with the substrate preference of the SULT1A1 enzymes. This observation may, however, only be qualitatively interpreted as (1) the targeted metabolomics mass spectrometric method used for the quantitative analysis of these substances needs further development, and (2) the influence of deiodonation was not taken into account in these studies. In conclusion, three perspectives are given at the end of the thesis which might be considered for further investigations. / Thesis (Ph.D. (Biochemistry))--North-West University, Potchefstroom Campus, 2010.

Ces dix dernières années, l’investigation des maladies génétiques a été bouleversée par l’émergence des techniques de séquençage haut-débit. Celles-ci permettent désormais de ne plus séquencer les gènes un par un, mais d’avoir accès à l’intégralité de la séquence génomique ou transcriptomique d’un individu. La difficulté devient alors d’identifier les variants causaux parmi une multitude d’artefacts techniques et de variants bénins, pour ensuite comprendre la physiopathologie des gènes identifiés.L’application du séquençage haut débit est particulièrement prometteuse dans le champ de la génétique de l’infertilité masculine car il s’agit d’une pathologie dont l’étiologie est souvent génétique, qui est génétiquement très hétérogène et pour laquelle peu de gènes ont été identifiés. Mon travail de thèse est donc centré sur la l’infertilité et comporte deux parties majeures : l’analyse des données issues du séquençage haut débit d’homme infertiles et de modèles animaux et la caractérisation moléculaire d’un phénotype spécifique d’infertilité, laglobozoospermie.Le nombre de variants identifiés dans le cadre d’un séquençage exomique pouvant s’élever à plusieurs dizaines de milliers, l’utilisation d’un outil informatique performant est indispensable. Pour arriver à une liste de variants suffisamment restreinte pour pouvoir être interprétée, plusieurs traitements sont nécessaires. Ainsi, j’ai développé un pipeline d’analyse de données issues de séquençage haut-débit effectuant de manière successive l’intégralité des étapes de l’analyse bio-informatique, c’est-à-dire l’alignement des reads sur un génome de référence, l’appel des génotypes, l’annotation des variants obtenus ainsi que le filtrage de ceux considérés comme non pertinents dans le contexte de l’analyse. L’ensemble de ces étapes étant interdépendantes,les réaliser au sein du même pipeline permet de mieux les calibrer pour ainsi réduire le nombre d’erreurs générées. Ce pipeline a été utilisé dans cinq études au sein du laboratoire, et a permis l’identification de variants impactant des gènes candidats prometteurs pouvant expliquer le phénotype d’infertilité des patients.L’ensemble des variants retenus ont ensuite pu être validés expérimentalement.J’ai également pris part aux investigations génétiques et moléculaires permettant la caractérisation du gène DPY19L2, identifié au laboratoire et dont la délétion homozygote entraine une globozoospermie, caractériséepar la présence dans l’éjaculât de spermatozoïdes à tête ronde dépourvus d’acrosome. Pour cela, j’ai contribué à caractériser les mécanismes responsables de cette délétion récurrente, puis, en utilisant le modèle murin Dpy19l2 knock out (KO) mimant le phénotype humain, j’ai réalisé une étude comparative des transcriptomes testiculaires de souris sauvages et de souris KO Dpy19l2-/-. Cette étude a ainsi permis de mettre en évidence la dérégulation de 76 gènes chez la souris KO. Parmi ceux-ci, 23 sont impliqués dans la liaison d’acides nucléiques et de protéines, pouvant ainsi expliquer les défauts d’ancrage de l’acrosome au noyau chez les spermatozoïdes globozoocéphales.Mon travail a donc permis de mieux comprendre la globozoospermie et de développer un pipeline d’analyse bioinformatique qui a déjà permis l’identification de plus de 15 gènes de la gamétogenèse humaine impliqués dans différents phénotypes d’infertilité. / In the last decade, the investigations of genetic diseases have been revolutionized by the rise of high throughput sequencing (HTS). Thanks to these new techniques it is now possible to analyze the totality of the coding sequences of an individual (exome sequencing) or even the sequences of his entire genome or transcriptome.The understanding of a pathology and of the genes associated with it now depends on our ability to identify causal variants within a plethora of technical artifact and benign variants.HTS is expected to be particularly useful in the field infertility as this pathology is expected to be highly genetically heterogeneous and only a few genes have so far been associated with it. My thesis focuses on male infertility and is divided into two main parts: HTS data analysis of infertile men and the molecular characterization of a specific phenotype, globozoospermia.Several thousands of distinct variants can be identified in a single exome, thereby using effective informatics is essential in order to obtain a short and actionable list of variants. It is for this purpose that I developed a HTS data analysis pipeline performing successively all bioinformatics analysis steps: 1) reads mapping along a reference genome, 2) genotype calling, 3) variant annotation and 4) the filtering of the variants considered as non-relevant for the analysis. Performing all these independent steps within a single pipeline is a good way to calibrate them and therefore to reduce the number of erroneous calls. This pipeline has been used in five studies and allowed the identification of variants impacting candidate genes that may explain the patients’ infertility phenotype. All these variants have been experimentally validated using Sanger sequencing.I also took part in the genetic and molecular investigations which permitted to demonstrate that the absence of the DPY192 gene induces male infertility due to globozoospermia, the presence in the ejaculate of only round-headed and acrosomeless spermatozoa. Most patients with globozoospermia have a homozygous deletion of the whole gene. I contributed to the characterization of the mechanisms responsible for this recurrent deletion, then, using Dpy19l2 knockout (KO) mice, I realized the comparative study of testicular transcriptome of wild type and Dpy19l2 -/- KO mice. This study highlighted a dysregulation of 76 genes in KO mice. Among them, 23 are involved in nucleic acid and protein binding, which may explain acrosome anchoring defaults observed in the sperm of globozoospermic patients.My work allowed a better understanding of globozoospermia and the development of a HTS data analysis pipeline. The latter allowed the identification of more than 15 human gametogenesis genes involved in different infertility phenotypes.

Submitted by MARCOS VINICIUS ANTUNES DE LEMOS (marcoslemoszootec@gmail.com) on 2017-06-02T15:32:45Z No. of bitstreams: 1 Tese_Marcos_VA_Lemos1.pdf: 3226396 bytes, checksum: 1f97ad6fd98b95977fbdc4ba573a83bb (MD5) / Approved for entry into archive by Luiz Galeffi (luizgaleffi@gmail.com) on 2017-06-02T17:51:15Z (GMT) No. of bitstreams: 1 lemos_mva_dr_jabo.pdf: 3226396 bytes, checksum: 1f97ad6fd98b95977fbdc4ba573a83bb (MD5) / Made available in DSpace on 2017-06-02T17:51:15Z (GMT). No. of bitstreams: 1 lemos_mva_dr_jabo.pdf: 3226396 bytes, checksum: 1f97ad6fd98b95977fbdc4ba573a83bb (MD5) Previous issue date: 2017-05-05 / Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) / Objetivou-se identificar regiões no genôma de bovinos da raça Nelore que apresentam variações no número de cópias (CNV) e, associar estes CNV com o perfil de ácidos graxos da carne. Além disso, objetivou-se realizar associação genica ampla utilizando os método de single step (GWASss) a fim de detectar regiões genômicas associadas aos ácidos graxos dos grupos saturados, mono e poliinsaturados, assim como os omegas 3, 6 e sua relação. O estudo de caracterização e distribuição dos CNVs ao longo do genoma de bovinos Nelore, foi realizado através do software PennCNV utizando dados genotípicos de 3.794 aniamais, resultando em 399.361 CNVs identificados. Após controle de qualidade, 2.902 foram mantidos nas analises, resultando em 195.873 CNVs, com tamanho medio de 54,744 pb, maximo de 8.7 Mb e minimo 3 kb. As regiões de CNV foram geradas pela sobreposição dos CNVs através do software CNVRuler. Os cromossomos que mostraram maior incidencia de CNVR foram BTA19 (24,26%), BTA23 (18,68%) e BTA25 (18,05%). Ja os que mostraram menor incidencia foram BTA29 (1,63%), BTA13 (9,72%) and BTA8 (9,72%). As 9.805 regiões da CNV estimadas no presente estudo cobrem aproximadamente 13.05% do genoma bovino e sobrepõem-se a 5.495 genes conhecidos que envolvem processos biológicos que poderiam estar envolvidos na adaptação ambiental da subespécie a áreas tropicais. O estudo de GWASss identificou 115 janelas que explicaram mais de 1% da variação genética aditiva para os 22 ácidos graxos estudados. A identificação destas regiões e seus genes genes, tais qual ELOVL5, ESRRG, PCYT1A e os genes do grupo ABC (ABCA5, ABCA6 e ABCA10) são genes que estão relacionados direto e inderamente ao metabolismo lipídico. O GWAS entres os fenótipos de AG e os CNVs resultaram em um total de 186 CNVR siginificantivos para os grupos dos ácidos graxos saturados (43), monosaturados (42), poliinsaturados (66) e omegas (35), nas quais foram identificados 278 genes com funçao descrita. Estes resultados apontaram genes associados a AG de varias saturações, podendo ser destacados os genes SAMD8 e BSCL2, os quais estão relacionados ao metabolism lipidico; e o gene RAPGEF6, relacionado ao metabolism energetico. Assim os inúmeras regiões genômicas encontradas neste estudo, bem como os genes identificados nas mesmas, devem contribuir para a formação de uma base genética do perfil de ácidos graxos da carne de bovino Nelore (Bos indicus), podendo contribuir para uma melhor seleção das características associadas à melhora da saúde humana. saúde. O conhecimento desses CNVs deverá melhorar a compreensão dos mecanismos genéticos e fisiológicos que contribuem para as características produtivas, bem como na seleção de animais mais produtivos e eficientes, contribuindo para o melhoramento genético das características produtivas. / The objective of this study was to identify genomic regions of Nellore cattle that present variations in the number of copies (CNV) and to associate these CNV with the fatty acid profile of the meat. In addition, the objective of this study was to carry out a genome-wid association using the single step method (GWASss) to detect genomic regions associated with saturated, mono and polyunsaturated fatty acids, as well as omega 3, 6 and their relationship. The study of the characterization and distribution of CNVs along the Nellore genome was performed by PennCNV software using genotypic data of 3,794 animals, resulting in 399,361 CNVs identified. After quality control, 2,902 were maintained in the analyzes, resulting in 195,873 CNVs, with an average size of 54,744 pb, maximum of 8.7 Mb and minimum 3 kb. The CNV regions were generated by the overlap of the CNVs by CNVRuler software. The chromosomes that showed the highest incidence of CNVR were BTA19 (24.26%), BTA23 (18.68%) and BTA25 (18.05%). Those with the lowest incidence were BTA29 (1.63%), BTA13 (9.72%) and BTA8 (9.72%). The 9,805 CNV regions estimated in the present study covered approximately 13.05% of the bovine genome and overlaped 5,495 genes known to envolve in biological processes that could be involved in the environmental adaptation of the subspecies to tropical areas. The GWASss study showed 115 windows that explained more than 1% of the additive genetic variation for the 22 fatty acids studied. The identification of these regions and their genes, such as ELOVL5, ESRRG, PCYT1A and the ABC group genes (ABCA5, ABCA6 and ABCA10) are genes directly and indirectly related to lipid metabolism. The GWAS between AG phenotypes and CNVs resulted in a total of 186 CNVRs that were significant for the saturated (43), monosaturated (42), polyunsaturated (66) and omegas (35) fatty acids groups, in which 278 genes with described function. These results pointed genes associated to AG of various saturations, and the genes SAMD8 and BSCL2, which are related to lipid metabolism; and the RAPGEF6 gene, related to energetic metabolism. Thus, the numerous genomic regions found in this study, as well as the genes identified in them, should contribute to the formation of a genetic basis of the Nellore beef (Bos indicus) fatty acid profile, contributing to a better selection of the traits associated with improvement of human health. The knowledge of these CNVs could improve the understanding of the genetic and physiological mechanisms that contribute to the productive traits, as well as the selection of more productive and efficient animals. / FAPESP: 2013/11853-4