Gene expression and genome-wide association studies analysis to select heavy pigs for protected designation of origin cured meats production

The high quality of protected designation of origin (PDO) dry-cured pork products depends largely on the chemical and physical parameters of the fresh meat and their variation during the production process of the final product. The discovery of the mechanisms that regulate the variability of these parameters was aided by the reference genome of swine adjuvant to genetic analysis methods. This thesis can contribute to the discovery of genetic mechanisms that regulate the variability of some quality parameters of fresh meat for PDO dry-cured pork production. The first study is of gene expression and showed that between low and high glycolytic potential (GP) samples of Semimembranosus muscle of Italian Large White (ILW) pigs in early postmortem, the differentially expressed genes were all but one over expressed in low GP. These were involved in ATP biosynthesis processes, calcium homeostasis, and lipid metabolism including the potential master regulator gene Peroxisome Proliferator-Activated Receptor Alpha (PPARA). The second is a study in commercial hybrid pigs to evaluate correlations between carcass and fresh ham traits, including carcass and fresh ham lean meat percentages, the former, a potential predictor of the latter. In addition, a genome-wide association study allowed the identification of chromosome-wide associations with phenotypic traits for 19 SNPs, and genome-wide associations for 14 SNPs for ferrochelatase activity. The latter could be a determinant for color variation in nitrite-free dry-cured ham. The third study showed gene expression differences in the Longissimus thoracis muscle of ILW pigs by feeding diets with extruded linseed (source of polyunsaturated fatty acids) and vitamin E and selenium (diet three) or natural (diet four) antioxidants. The diet three promoted a more rapid and massive immune system response possibly determined by improvement in muscle tissue function, while the diet four promoted oxidative stability and increased the anti-inflammatory potential of muscle tissue.

Applied genomics in livestock: genome-wide association studies and population genomics analysis in pig and rabbit breeds

The domestication and selection processes in pigs and rabbits have resulted in the constitution of multiple breeds with broad phenotypic diversity. Population genomics analysis and Genome-wide association study analysis can be utilized to gain insights into the ancestral origins, genetic diversity, and the presence of lethal mutations across these diverse breeds. In this thesis, we analysed the dataset obtained from three Italian Pig breeds to detect deleterious alleles. We screened the dataset for genetic markers showing homozygous deficiency using two approaches single marker and haplotype-based approach. Moreover, Genome-wide association study analyses were performed to detect genetic markers associated with pigs' reproductive traits. In rabbits, we investigated the application of SNP bead chip for detection signatures of selection in rabbits using different methods. This analysis was implemented for the first time in different fancy and meet rabbit breeds. Multiple approaches were utilized for the detection of the selection of signatures including Fst analysis, ROH analysis, PCAdapt analysis, and haplotype-based analysis. The analysis in pigs was able to identify five putative deleterious SNPs and nine putative deleterious haplotypes in the analysed Italian Pig breeds. The genomic regions of the detected putative deleterious genomic markers harboring loss of function variants such as the Frameshift variant, start lost, and splice donor variant. Those variants are close to important candidate genes such as IGF2BP1, ADGRL4, and HGF. In rabbits, multiple genomic regions were detected to be under selection of signature. These genomic regions harbor candidate genes associated with coat color phenotype (MC1R, TYR, and ASIP), hair structure (LIPH), and body size (HMGA2 and COL2A1). The described results in rabbits and pigs could be used to improve breeding programs by excluding the deleterious genetic markers carriers and incorporating candidate genes for coat color, body size, and meat production in rabbit breeding programs to enhance desired traits

High-throughput SNP genotyping in the highly heterozygous genome of Eucalyptus: assay success, polymorphism and transferability across species

Background: High-throughput SNP genotyping has become an essential requirement for molecular breeding and population genomics studies in plant species. Large scale SNP developments have been reported for several mainstream crops. A growing interest now exists to expand the speed and resolution of genetic analysis to outbred species with highly heterozygous genomes. When nucleotide diversity is high, a refined diagnosis of the target SNP sequence context is needed to convert queried SNPs into high-quality genotypes using the Golden Gate Genotyping Technology (GGGT). This issue becomes exacerbated when attempting to transfer SNPs across species, a scarcely explored topic in plants, and likely to become significant for population genomics and inter specific breeding applications in less domesticated and less funded plant genera. Results: We have successfully developed the first set of 768 SNPs assayed by the GGGT for the highly heterozygous genome of Eucalyptus from a mixed Sanger/454 database with 1,164,695 ESTs and the preliminary 4.5X draft genome sequence for E. grandis. A systematic assessment of in silico SNP filtering requirements showed that stringent constraints on the SNP surrounding sequences have a significant impact on SNP genotyping performance and polymorphism. SNP assay success was high for the 288 SNPs selected with more rigorous in silico constraints; 93% of them provided high quality genotype calls and 71% of them were polymorphic in a diverse panel of 96 individuals of five different species. SNP reliability was high across nine Eucalyptus species belonging to three sections within subgenus Symphomyrtus and still satisfactory across species of two additional subgenera, although polymorphism declined as phylogenetic distance increased. Conclusions: This study indicates that the GGGT performs well both within and across species of Eucalyptus notwithstanding its nucleotide diversity >= 2%. The development of a much larger array of informative SNPs across multiple Eucalyptus species is feasible, although strongly dependent on having a representative and sufficiently deep collection of sequences from many individuals of each target species. A higher density SNP platform will be instrumental to undertake genome-wide phylogenetic and population genomics studies and to implement molecular breeding by Genomic Selection in Eucalyptus.

Near full-length genome analysis of low prevalent human immunodeficiency virus type 1 subclade F1 in Sao Paulo, Brazil

Background: The genetic diversity of the human immunodeficiency virus type 1 (HIV-1) is critical to lay the groundwork for the design of successful drugs or vaccine. In this study we aimed to characterize and define the molecular prevalence of HIV-1 subclade F1 currently circulating in Sao Paulo, Brazil. Methods: A total of 36 samples were selected from 888 adult patients residing in Sao Paulo who had previously been diagnosed in two independent studies in our laboratory as being infected with subclade F1 based on pol subgenomic fragment sequencing. Proviral DNA was amplified from the purified genomic DNA of all 36 blood samples by 5 fragments overlapping PCR followed by direct sequencing. Sequence data were obtained from the 5 fragments of pure subclade F1 and phylogenetic trees were constructed and compared with previously published sequences. Subclades F1 that exhibited mosaic structure with other subtypes were omitted from any further analysis Results: Our methods of fragment amplification and sequencing confirmed that only 5 sequences inferred from pol region as subclade F1 also holds true for the genome as a whole and, thus, estimated the true prevalence at 0.56%. The results also showed a single phylogenetic cluster of the Brazilian subclade F1 along with non-Brazilian South American isolates in both subgenomic and the full-length genomes analysis with an overall intrasubtype nucleotide divergence of 6.9%. The nucleotide differences within the South American and Central African F1 strains, in the C2-C3 env, were 8.5% and 12.3%, respectively. Conclusion: All together, our findings showed a surprisingly low prevalence rate of subclade F1 in Brazil and suggest that these isolates originated in Central Africa and subsequently introduced to South America.

An assessment of population structure in eight breeds of cattle using a whole genome SNP panel

Background: Analyses of population structure and breed diversity have provided insight into the origin and evolution of cattle. Previously, these studies have used a low density of microsatellite markers, however, with the large number of single nucleotide polymorphism markers that are now available, it is possible to perform genome wide population genetic analyses in cattle. In this study, we used a high-density panel of SNP markers to examine population structure and diversity among eight cattle breeds sampled from Bos indicus and Bos taurus. Results: Two thousand six hundred and forty one single nucleotide polymorphisms ( SNPs) spanning all of the bovine autosomal genome were genotyped in Angus, Brahman, Charolais, Dutch Black and White Dairy, Holstein, Japanese Black, Limousin and Nelore cattle. Population structure was examined using the linkage model in the program STRUCTURE and Fst estimates were used to construct a neighbor-joining tree to represent the phylogenetic relationship among these breeds. Conclusion: The whole-genome SNP panel identified several levels of population substructure in the set of examined cattle breeds. The greatest level of genetic differentiation was detected between the Bos taurus and Bos indicus breeds. When the Bos indicus breeds were excluded from the analysis, genetic differences among beef versus dairy and European versus Asian breeds were detected among the Bos taurus breeds. Exploration of the number of SNP loci required to differentiate between breeds showed that for 100 SNP loci, individuals could only be correctly clustered into breeds 50% of the time, thus a large number of SNP markers are required to replace the 30 microsatellite markers that are currently commonly used in genetic diversity studies.

Integration and visualization of systems biology data in context of the genome

Background: High-density tiling arrays and new sequencing technologies are generating rapidly increasing volumes of transcriptome and protein-DNA interaction data. Visualization and exploration of this data is critical to understanding the regulatory logic encoded in the genome by which the cell dynamically affects its physiology and interacts with its environment. Results: The Gaggle Genome Browser is a cross-platform desktop program for interactively visualizing high-throughput data in the context of the genome. Important features include dynamic panning and zooming, keyword search and open interoperability through the Gaggle framework. Users may bookmark locations on the genome with descriptive annotations and share these bookmarks with other users. The program handles large sets of user-generated data using an in-process database and leverages the facilities of SQL and the R environment for importing and manipulating data. A key aspect of the Gaggle Genome Browser is interoperability. By connecting to the Gaggle framework, the genome browser joins a suite of interconnected bioinformatics tools for analysis and visualization with connectivity to major public repositories of sequences, interactions and pathways. To this flexible environment for exploring and combining data, the Gaggle Genome Browser adds the ability to visualize diverse types of data in relation to its coordinates on the genome. Conclusions: Genomic coordinates function as a common key by which disparate biological data types can be related to one another. In the Gaggle Genome Browser, heterogeneous data are joined by their location on the genome to create information-rich visualizations yielding insight into genome organization, transcription and its regulation and, ultimately, a better understanding of the mechanisms that enable the cell to dynamically respond to its environment.

Deep sequencing of New World screw-worm transcripts to discover genes involved in insecticide resistance

Background: The New World screw-worm (NWS), Cochliomyia hominivorax, is one of the most important myiasis-causing flies, causing severe losses to the livestock industry. In its current geographical distribution, this species has been controlled by the application of insecticides, mainly organophosphate (OP) compounds, but a number of lineages have been identified that are resistant to such chemicals. Despite its economic importance, only limited genetic information is available for the NWS. Here, as a part of an effort to characterize the C. hominivorax genome and identify putative genes involved in insecticide resistance, we sampled its transcriptome by deep sequencing of polyadenylated transcripts using the 454 sequencing technology. Results: Deep sequencing on the 454 platform of three normalized libraries (larval, adult male and adult female) generated a total of 548,940 reads. Eighteen candidate genes coding for three metabolic detoxification enzyme families, cytochrome P450 monooxygenases, glutathione S transferases and carboxyl/cholinesterases were selected and gene expression levels were measured using quantitative real-time polymerase chain reaction (qRT-PCR). Of the investigated candidates, only one gene was expressed differently between control and resistant larvae with, at least, a 10-fold down-regulation in the resistant larvae. The presence of mutations in the acetylcholinesterase (target site) and carboxylesterase E3 genes was investigated and all of the resistant flies presented E3 mutations previously associated with insecticide resistance. Conclusions: Here, we provided the largest database of NWS expressed sequence tags that is an important resource, not only for further studies on the molecular basis of the OP resistance in NWS fly, but also for functional and comparative studies among Calliphoridae flies. Among our candidates, only one gene was found differentially expressed in resistant individuals, and its role on insecticide resistance should be further investigated. Furthermore, the absence of mutations in the OP target site and the high frequency of mutant carboxylesterase E3 indicate that metabolic resistance mechanisms have evolved predominantly in this species.

Evaluation of the genetic variability of orchid fleck virus by single-strand conformational polymorphism analysis and nucleotide sequencing of a fragment from the nucleocapsid gene

The variability of a fragment of the nucleocapsid gene of orchid fleck virus (OFV) was investigated by single-strand conformational polymorphism (SSCP) analysis and nucleotide sequencing. Forty-eight samples of 18 genera of orchids were collected from Brazil, Costa Rica and Australia. The SSCP analysis yielded six different band patterns, and phylogenetic analysis based on the nucleotide fragment sequence obtained in this work and six available in GenBank showed two different groups, one with isolates 023Germany and So-Japan, and other with the rest of the isolates. None of the analyses showed geographic correlation among the Brazilian strains. The data obtained in this study showed a low genetic variation in this region of the genome; the d(N)/d(S) ratio of 0.251-0.405 demonstrated a negative selective pressure that maintains the stability of the analyzed fragments.

Distribution, expression, and motif variability of ankyrin domain genes in Wolbachia pipientis

The endosymbiotic bacterium Wolbachia pipientis infects a wide range of arthropods, in which it induces a variety of reproductive phenotypes, including cytoplasmic incompatibility (CI), parthenogenesis, male killing, and reversal of genetic sex determination. The recent sequencing and annotation of the first Wolbachia genome revealed an unusually high number of genes encoding ankyrin domain (ANK) repeats. These ANK genes are likely to be important in mediating the Wolbachia-host interaction. In this work we determined the distribution and expression of the different ANK genes found in the sequenced Wolbachia wMel genome in nine Wolbachia strains that induce different phenotypic effects in their hosts. A comparison of the ANK genes of wMel and the non-CI-inducing wAu Wolbachia strain revealed significant differences between the strains. This was reflected in sequence variability in shared genes that could result in alterations in the encoded proteins, such as motif deletions, amino acid insertions, and in some cases disruptions due to insertion of transposable elements and premature stops. In addition, one wMel ANK gene, which is part of an operon, was absent in the wAu genome. These variations are likely to affect the affinity, function, and cellular location of the predicted proteins encoded by these genes.

Wolbachia: why these bacteria are important to genome research

Intracellular bacteria of the genus Wolbachia were first discovered in mosquitoes in the 1920s. Their superficial similarity to pathogenic rickettsia initially raised interest in them as potential human pathogens. However, injection experiments with mice showed that they were non-pathogenic, and they were subsequently classified as symbionts of insects. Until the 1970s, Wolbachia was considered to infect a limited number of species of mosquitoes. It is now clear that Wolbachia is an extremely common intracellular agent of invertebrates, infecting nearly all the major groups of arthropods and other terrestrial invertebrates. Its wide host range and abundance can be attributed partly to the unusual phenotypes it exerts on the host it infects. These include the induction of parthenogenesis (the production of female offspring from unmated mothers) in certain insects, the feminization of genetic male crustaceans to functional phenotypic females, and the failure of fertilization in hosts when males and females have a different infection status (cytoplasmic incompatibility). All of these phenotypes favor maternal transmission of the intracellular Wolbachia. In the last year, Wolbachia has also been shown to be a widespread symbiont of filarial nematodes. It appears that Wolbachia is needed by the adult worm for normal fertility, indicating that Wolbachia is behaving like a classic mutualist in this case. This discovery exemplifies that the extent of the host range of Wolbachia and its associated phenotypes is still far from fully understood.

Algorithms for sequence analysis via mutagenesis

Despite many successes of conventional DNA sequencing methods, some DNAs remain difficult or impossible to sequence. Unsequenceable regions occur in the genomes of many biologically important organisms, including the human genome. Such regions range in length from tens to millions of bases, and may contain valuable information such as the sequences of important genes. The authors have recently developed a technique that renders a wide range of problematic DNAs amenable to sequencing. The technique is known as sequence analysis via mutagenesis (SAM). This paper presents a number of algorithms for analysing and interpreting data generated by this technique.

The p.A2215D Thyroglobulin Gene Mutation Leads to Deficient Synthesis and Secretion of the Mutated Protein and Congenital Hypothyroidism with Wide Phenotype Variation

Context: Thyroglobulin (TG) is a large glycoprotein and functions as a matrix for thyroid hormone synthesis. TG gene mutations give rise to goitrous congenital hypothyroidism (CH) with considerable phenotype variation. Objectives: The aim of the study was to report the genetic screening of 15 patients with CH due to TG gene mutations and to perform functional analysis of the p. A2215D mutation. Design: Clinical evaluation and DNA sequencing of the TG gene were performed in all patients. TG expression was analyzed in the goitrous tissue of one patient. Human cells were transfected with expression vectors containing mutated and wild-type human TG cDNA. Results: All patients had an absent rise of serum TG after stimulation with recombinant human TSH. Sequence analysis revealed three previously described mutations (p. A2215D, p. R277X, and g. IVS30 + 1G > T), and two novel mutations (p. Q2142X and g. IVS46-1G > A). Two known (g. IVS30 + 1G/p. A2215D and p. A2215D/p. R277X) and one novel (p. R277X/g. IVS46-1G > A) compound heterozygous constellations were also identified. Functional analysis indicated deficiency in TG synthesis, reduction of TG secretion, and retention of the mutant TG within the cell, leading to an endoplasmic reticulum storage disease, whereas small amounts of mutant TG were still secreted within the cell system. Conclusion: All studied patients were either homozygous or heterozygous for TG gene mutations. Two novel mutations have been detected, and we show that TG mutation p. A2215D promotes the retention of TG within the endoplasmic reticulum and reduces TG synthesis and secretion, causing mild hypothyroidism. In the presence of sufficient iodine supply, some patients with TG mutations are able to compensate the impaired hormonogenesis and generate thyroid hormone. (J Clin Endocrinol Metab 94: 2938-2944, 2009)

High frequency of submicroscopic chromosomal imbalances in patients with syndromic craniosynostosis detected by a combined approach of microsatellite segregation analysis, multiplex ligation-dependent probe amplification and array-based comparative genome hybridisation

We present the first comprehensive study, to our knowledge, on genomic chromosomal analysis in syndromic craniosynostosis. In total, 45 patients with craniosynostotic disorders were screened with a variety of methods including conventional karyotype, microsatellite segregation analysis, subtelomeric multiplex ligation-dependent probe amplification) and whole-genome array-based comparative genome hybridisation. Causative abnormalities were present in 42.2% (19/45) of the samples, and 27.8% (10/36) of the patients with normal conventional karyotype carried submicroscopic imbalances. Our results include a wide variety of imbalances and point to novel chromosomal regions associated with craniosynostosis. The high incidence of pure duplications or trisomies suggests that these are important mechanisms in craniosynostosis, particularly in cases involving the metopic suture.

Constitutional Telomerase Mutations Are Genetic Risk Factors for Cirrhosis

Some patients with liver disease progress to cirrhosis, but the risk factors for cirrhosis development are unknown. Dyskeratosis congenita, an inherited bone marrow failure syndrome associated with mucocutaneous anomalies, pulmonary fibrosis, and cirrhosis, is caused by germline mutations of genes in the telomerase complex. We examined whether telomerase mutations also occurred in sporadic cirrhosis. In all, 134 patients with cirrhosis of common etiologies treated at the Liver Research Institute, University of Arizona, between May 2008 and July 2009, and 528 healthy subjects were screened for variation in the TERT and TERC genes by direct sequencing; an additional 1,472 controls were examined for the most common genetic variation observed in patients. Telomere length of leukocytes was measured by quantitative polymerase chain reaction. Functional effects of genetic changes were assessed by transfection of mutation-containing vectors into telomerase-deficient cell lines, and telomerase activity was measured in cell lysates. Nine of the 134 patients with cirrhosis (7%) carried a missense variant in TERT, resulting in a cumulative carrier frequency significantly higher than in controls (P = 0.0009). One patient was homozygous and eight were heterozygous. The allele frequency for the most common missense TERT variant was significantly higher in patients with cirrhosis (2.6%) than in 2,000 controls (0.7%; P = 0.0011). One additional patient carried a TERC mutation. The mean telomere length of leukocytes in patients with cirrhosis, including six mutant cases, was shorter than in age-matched controls (P = 0.0004). Conclusion: Most TERT gene variants reduced telomerase enzymatic activity in vitro. Loss-of-function telomerase gene variants associated with short telomeres are risk factors for sporadic cirrhosis. (HEPATOLOGY 2011;53:1600-1607)

Linkage to chromosome 2q36.1 in autosomal dominant Dandy-Walker malformation with occipital cephalocele and evidence for genetic heterogeneity

We previously reported a Vietnamese-American family with isolated autosomal dominant occipital cephalocele. Upon further neuroimaging studies, we have recharacterized this condition as autosomal dominant Dandy-Walker with occipital cephalocele (ADDWOC). A similar ADDWOC family from Brazil was also recently described. To determine the genetic etiology of ADDWOC, we performed genome-wide linkage analysis on members of the Vietnamese-American and Brazilian pedigrees. Linkage analysis of the Vietnamese-American family identified the ADDWOC causative locus on chromosome 2q36.1 with a multipoint parametric LOD score of 3.3, while haplotype analysis refined the locus to 1.1 Mb. Sequencing of the five known genes in this locus did not identify any protein-altering mutations. However, a terminal deletion of chromosome 2 in a patient with an isolated case of Dandy-Walker malformation also encompassed the 2q36.1 chromosomal region. The Brazilian pedigree did not show linkage to this 2q36.1 region. Taken together, these results demonstrate a locus for ADDWOC on 2q36.1 and also suggest locus heterogeneity for ADDWOC.