Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Reliable Means of Diagnosis and Serovar Determination of Blood-Borne Salmonella Strains: Quick PCR Amplification of Unique Genomic Loci by Novel Primer Sets

Typhoid fever is becoming an ever increasing threat in the developing countries. We have improved considerably upon the existing PCR-based diagnosis method by designing primers against a region that is unique to Salmonella enterica subsp. enterica serovar Typhi and Salmonella enterica subsp. enterica serovar Paratyphi A, corresponding to the STY0312 gene in S. Typhi and its homolog SPA2476 in S. Paratyphi A. An additional set of primers amplify another region in S. Typhi CT18 and S. Typhi Ty2 corresponding to the region between genes STY0313 to STY0316 but which is absent in S. Paratyphi A. The possibility of a false-negative result arising due to mutation in hypervariable genes has been reduced by targeting a gene unique to typhoidal Salmonella serovars as a diagnostic marker. The amplified region has been tested for genomic stability by amplifying the region from clinical isolates of patients from various geographical locations in India, thereby showing that this region is potentially stable. These set of primers can also differentiate between S. Typhi CT18, S. Typhi Ty2, and S. Paratyphi A, which have stable deletions in this specific locus. The PCR assay designed in this study has a sensitivity of 95% compared to the Widal test which has a sensitivity of only 63%. As observed, in certain cases, the PCR assay was more sensitive than the blood culture test was, as the PCR-based detection could also detect dead bacteria.

Genomic and functional profiling of gastric cancer

Helicobacter pylori infection is a risk factor for gastric cancer, which is a major health issue worldwide. Gastric cancer has a poor prognosis due to the unnoticeable progression of the disease and surgery is the only available treatment in gastric cancer. Therefore, gastric cancer patients would greatly benefit from identifying biomarker genes that would improve diagnostic and prognostic prediction and provide targets for molecular therapies. DNA copy number amplifications are the hallmarks of cancers in various anatomical locations. Mechanisms of amplification predict that DNA double-strand breaks occur at the margins of the amplified region. The first objective of this thesis was to identify the genes that were differentially expressed in H. pylori infection as well as the transcription factors and signal transduction pathways that were associated with the gene expression changes. The second objective was to identify putative biomarker genes in gastric cancer with correlated expression and copy number, and the last objective was to characterize cancers based on DNA copy number amplifications. DNA microarrays, an in vitro model and real-time polymerase chain reaction were used to measure gene expression changes in H. pylori infected AGS cells. In order to identify the transcription factors and signal transduction pathways that were activated after H. pylori infection, gene expression profiling data from the H. pylori experiments and a bioinformatics approach accompanied by experimental validation were used. Genome-wide expression and copy number microarray analysis of clinical gastric cancer samples and immunohistochemistry on tissue microarray were used to identify putative gastric cancer genes. Data mining and machine learning techniques were applied to study amplifications in a cross-section of cancers. FOS and various stress response genes were regulated by H. pylori infection. H. pylori regulated genes were enriched in the chromosomal regions that are frequently changed in gastric cancer, suggesting that molecular pathways of gastric cancer and premalignant H. pylori infection that induces gastritis are interconnected. 16 transcription factors were identified as being associated with H. pylori infection induced changes in gene expression. NF-κB transcription factor and p50 and p65 subunits were verified using elecrophoretic mobility shift assays. ERBB2 and other genes located in 17q12- q21 were found to be up-regulated in association with copy number amplification in gastric cancer. Cancers with similar cell type and origin clustered together based on the genomic localization of the amplifications. Cancer genes and large genes were co-localized with amplified regions and fragile sites, telomeres, centromeres and light chromosome bands were enriched at the amplification boundaries. H. pylori activated transcription factors and signal transduction pathways function in cellular mechanisms that might be capable of promoting carcinogenesis of the stomach. Intestinal and diffuse type gastric cancers showed distinct molecular genetic profiles. Integration of gene expression and copy number microarray data allowed the identification of genes that might be involved in gastric carcinogenesis and have clinical relevance. Gene amplifications were demonstrated to be non-random genomic instabilities. Cell lineage, properties of precursor stem cells, tissue microenvironment and genomic map localization of specific oncogenes define the site specificity of DNA amplifications, whereas labile genomic features define the structures of amplicons. These conclusions suggest that the definition of genomic changes in cancer is based on the interplay between the cancer cell and the tumor microenvironment.

Multi-trait assessment of early-in-life female, male and genomic measures for use in genetic selection to improve female reproductive performance of Brahman cattle

Early-in-life female and male measures with potential to be practical genetic indicators were chosen from earlier analyses and examined together with genomic measures for multi-trait use to improve female reproduction of Brahman cattle. Combinations of measures were evaluated on the genetic gains expected from selection of sires and dams for each of age at puberty (AGECL, i.e. first observation of a corpus luteum), lactation anoestrous interval in 3-year-old cows (LAI), and lifetime annual weaning rate (LAWR, i.e. the weaning rate of cows based on the number of annual matings they experienced over six possible matings). Selection was on an index of comparable records for each combination. Selection intensities were less than theoretically possible but assumed a concerted selection effort was able to be made across the Brahman breed. The results suggested that substantial genetic gains could be possible but need to be confirmed in other data. The estimated increase in LAWR in 10 years, for combinations without or with genomic measures, ranged from 8 to 12 calves weaned per 100 cows from selection of sires, and from 12 to 15 calves weaned per 100 cows from selection of sires and dams. Corresponding reductions in LAI were 60-103 days or 94-136 days, and those for AGECL were 95-125 or 141-176 days, respectively. Coat score (a measure of the sleekness or wooliness of the coat) and hip height in females, and preputial eversion and liveweight in males, were measures that may warrant wider recording for Brahman female reproduction genetic evaluation. Pregnancy-test outcomes from Matings 1 and 2 also should be recorded. Percentage normal sperm may be important to record for reducing LAI and scrotal size and serum insulin-like growth factor-I concentration in heifers at 18 months for reducing AGECL. Use of a genomic estimated breeding value (EBV) in combination with other measures added to genetic gains, especially at genomic EBV accuracies of 40%. Accuracies of genomic EBVs needed to approach 60% for the genomic EBV to be the most important contributor to gains in the combinations of measures studied.

Genomic deletions and mutations resulting in the loss of eight genes reduce the in vivo replication capacity of Meleagrid herpesvirus 1

Meleagrid herpesvirus 1 (MeHV-1 or turkey herpesvirus) has been widely used as a vaccine in commercial poultry. Initially, these vaccine applications were for the prevention of Marek’s disease resulting from Gallid herpesvirus 2 infections, while more recently MeHV-1 has been used as recombinant vector for other poultry infections. The construction of herpesvirus infectious clones that permit propagation and manipulation of the viral genome in bacterial hosts has advanced the studies of herpesviral genetics. The current study reports the construction of five MeHV-1 infectious clones. The in vitro properties of viruses recovered from these clones were indistinguishable from the parental MeHV-1. In contrast, the rescued MeHV-1 viruses were significantly attenuated when used in vivo. Complete sequencing of the infectious clones identified the absence of two regions of the MeHV-1 genome compared to the MeHV-1 reference sequence. These analyses determined the rescued viruses have seven genes, UL43, UL44, UL45, UL56, HVT071, sorf3 and US2 either partially or completely deleted. In addition, single nucleotide polymorphisms were identified in all clones compared with the MeHV-1 reference sequence. As a consequence of one of the polymorphisms identified in the UL13 gene, four of the rescued viruses were predicted to encode a serine/threonine protein kinase lacking two of three domains required for activity. Thus four of the recovered viruses have a total of eight missing or defective genes. The implications of these findings in the context of herpesvirus biology and infectious clone construction are discussed.

Genomic deletions and mutations resulting in the loss of eight genes reduce the in vivo replication capacity of Meleagrid herpesvirus 1

Meleagrid herpesvirus 1 (MeHV-1 or turkey herpesvirus) has been widely used as a vaccine in commercial poultry. Initially, these vaccine applications were for the prevention of Marek’s disease resulting from Gallid herpesvirus 2 infections, while more recently MeHV-1 has been used as recombinant vector for other poultry infections. The construction of herpesvirus infectious clones that permit propagation and manipulation of the viral genome in bacterial hosts has advanced the studies of herpesviral genetics. The current study reports the construction of five MeHV-1 infectious clones. The in vitro properties of viruses recovered from these clones were indistinguishable from the parental MeHV-1. In contrast, the rescued MeHV-1 viruses were significantly attenuated when used in vivo. Complete sequencing of the infectious clones identified the absence of two regions of the MeHV-1 genome compared to the MeHV-1 reference sequence. These analyses determined the rescued viruses have seven genes, UL43, UL44, UL45, UL56, HVT071, sorf3 and US2 either partially or completely deleted. In addition, single nucleotide polymorphisms were identified in all clones compared with the MeHV-1 reference sequence. As a consequence of one of the polymorphisms identified in the UL13 gene, four of the rescued viruses were predicted to encode a serine/threonine protein kinase lacking two of three domains required for activity. Thus four of the recovered viruses have a total of eight missing or defective genes. The implications of these findings in the context of herpesvirus biology and infectious clone construction are discussed.

Computational Methods for Detecting Large-Scale Chromosome Rearrangements in SNP Data

Large-scale chromosome rearrangements such as copy number variants (CNVs) and inversions encompass a considerable proportion of the genetic variation between human individuals. In a number of cases, they have been closely linked with various inheritable diseases. Single-nucleotide polymorphisms (SNPs) are another large part of the genetic variance between individuals. They are also typically abundant and their measuring is straightforward and cheap. This thesis presents computational means of using SNPs to detect the presence of inversions and deletions, a particular variety of CNVs. Technically, the inversion-detection algorithm detects the suppressed recombination rate between inverted and non-inverted haplotype populations whereas the deletion-detection algorithm uses the EM-algorithm to estimate the haplotype frequencies of a window with and without a deletion haplotype. As a contribution to population biology, a coalescent simulator for simulating inversion polymorphisms has been developed. Coalescent simulation is a backward-in-time method of modelling population ancestry. Technically, the simulator also models multiple crossovers by using the Counting model as the chiasma interference model. Finally, this thesis includes an experimental section. The aforementioned methods were tested on synthetic data to evaluate their power and specificity. They were also applied to the HapMap Phase II and Phase III data sets, yielding a number of candidates for previously unknown inversions, deletions and also correctly detecting known such rearrangements.

Genomic regions influencing seminal root traits in barley

Water availability is a major limiting factor for crop production, making drought adaptation and its many component traits a desirable attribute of plant cultivars. Previous studies in cereal crops indicate that root traits expressed at early plant developmental stages, such as seminal root angle and root number, are associated with water extraction at different depths. Here, we conducted the first study to map seminal root traits in barley (Hordeum vulgare L.). Using a recently developed high-throughput phenotyping method, a panel of 30 barley genotypes and a doubled-haploid (DH) population (ND24260 × 'Flagship') comprising 330 lines genotyped with diversity array technology (DArT) markers were evaluated for seminal root angle (deviation from vertical) and root number under controlled environmental conditions. A high degree of phenotypic variation was observed in the panel of 30 genotypes: 13.5 to 82.2 and 3.6 to 6.9° for root angle and root number, respectively. A similar range was observed in the DH population: 16.4 to 70.5 and 3.6 to 6.5° for root angle and number, respectively. Seven quantitative trait loci (QTL) for seminal root traits (root angle, two QTL; root number, five QTL) were detected in the DH population. A major QTL influencing both root angle and root number (RAQ2/RNQ4) was positioned on chromosome 5HL. Across-species analysis identified 10 common genes underlying root trait QTL in barley, wheat (Triticum aestivum L.), and sorghum [Sorghum bicolor (L.) Moench]. Here, we provide insight into seminal root phenotypes and provide a first look at the genetics controlling these traits in barley.

Genomic approach to organ-specific growth control

Growth is a fundamental aspect of life cycle of all organisms. Body size varies highly in most animal groups, such as mammals. Moreover, growth of a multicellular organism is not uniform enlargement of size, but different body parts and organs grow to their characteristic sizes at different times. Currently very little is known about the molecular mechanisms governing this organ-specific growth. The genome sequencing projects have provided complete genomic DNA sequences of several species over the past decade. The amount of genomic sequence information, including sequence variants within species, is constantly increasing. Based on the universal genetic code, we can make sense of this sequence information as far as it codes proteins. However, less is known about the molecular mechanisms that control expression of genes, and about the variations in gene expression that underlie many pathological states in humans. This is caused in part by lack of information about the second genetic code that consists of the binding specificities of transcription factors and the combinatorial code by which transcription factor binding sites are assembled to form tissue-specific and/or ligand-regulated enhancer elements. This thesis presents a high-throughput assay for identification of transcription factor binding specificities, which were then used to measure the DNA binding profiles of transcription factors involved in growth control. We developed ‘enhancer element locator’, a computational tool, which can be used to predict functional enhancer elements. A genome-wide prediction of human and mouse enhancer elements generated a large database of enhancer elements. This database can be used to identify target genes of signaling pathways, and to predict activated transcription factors based on changes in gene expression. Predictions validated in transgenic mouse embryos revealed the presence of multiple tissue-specific enhancers in mouse c- and N-Myc genes, which has implications to organ specific growth control and tumor type specificity of oncogenes. Furthermore, we were able to locate a variation in a single nucleotide, which carries a susceptibility to colorectal cancer, to an enhancer element and propose a mechanism by which this SNP might be involved in generation of colorectal cancer.

Genomic Profiling of Gastric Cancer

Gastric cancer is the fourth most common cancer and the second most common cause of cancer-related death worldwide. Due to lack of early symptoms, gastric cancer is characterized by late stage diagnosis and unsatisfactory options for curative treatment. Several genomic alterations have been identified in gastric cancer, but the major factors contributing to initiation and progression of gastric cancer remain poorly known. Gene copy number alterations play a key role in the development of gastric cancer, and a change in gene copy number is one of the fundamental mechanisms for a cancer cell to control the expression of potential oncogenes and tumor suppressor genes. This thesis aims at clarifying the complex genomic alterations of gastric cancer to identify novel molecular biomarkers for diagnostic purposes as well as for targeted treatment. To highlight genes of potential biological and clinical relevance, we carried out a systematic microarray-based survey of gene expression and copy number levels in primary gastric tumors and gastric cancer cell lines. Results were validated using immunohistochemistry, real-time qRT-PCR, and affinity capture-based transcript (TRAC) assay. Altogether 192 clinical gastric tissue samples and 7 gastric cancer cell lines were included in this study. Multiple chromosomal regions with recurrent copy number alterations were detected. The most frequent chromosomal alterations included gains at 7q, 8q, 17q, 19q, and 20q and losses at 9p, 18q, and 21q. Distinctive patterns of copy number alterations were detected for different histological subtypes (intestinal and diffuse) and for cancers located in different parts of the stomach. The impact of copy number alterations on gene expression was significant, as 6-10% of genes located in the regions of gains and losses also showed concomitant alterations in their expression. By combining the information from the DNA- and RNA-level analyses many novel gastric cancer-related genes, such as ALPK2, ENAH, HHIPL2, and OSMR, were identified. Independent genome-wide gene expression analysis of Finnish and Japanese gastric tumors revealed an additional set of genes that was differentially expressed in cancerous gastric tissues compared with normal tissue. Overexpression of one of these genes, CXCL1, was associated with an improved survival of gastric cancer. Thus, using an integrative microarray analysis, several novel genes were identified that may be critically important for gastric carcinogenesis. Further studies of these genes may lead to novel biomarkers for gastric cancer diagnosis and targeted therapy.

Chromosomal rearrangements, phenotypic variation and modularity: A case study from a contact zone between house mouse Robertsonian races in Central Italy

The Western European house mouse, Mus musculus domesticus, is well-known for the high frequency of Robertsonian fusions that have rapidly produced more than 50 karyotipic races, making it an ideal model for studying the mechanisms of chromosomal speciation. The mouse mandible is one of the traits studied most intensively to investigate the effect of Robertsonian fusions on phenotypic variation within and between populations. This complex bone structure has also been widely used to study the level of integration between different morphogenetic units. Here, with the aim of testing the effect of different karyotypic assets on the morphology of the mouse mandible and on its level of modularity, we performed morphometric analyses of mice from a contact area between two highly metacentric races in Central Italy. We found no difference in size, while the mandible shape was found to be different between the two Robertsonian races, even after accounting for the genetic relationships among individuals and geographic proximity. Our results support the existence of two modules that indicate a certain degree of evolutionary independence, but no difference in the strength of modularity between chromosomal races. Moreover, the ascending ramus showed more pronounced interpopulation/race phenotypic differences than the alveolar region, an effect that could be associated to their different polygenic architecture. This study suggests that chromosomal rearrangements play a role in the house mouse phenotypic divergence, and that the two modules of the mouse mandible are differentially affected by environmental factors and genetic makeup.

Exploring the post-genomic world: Differing explanatory and manipulatory functions of post-genomic sciences

Richard Lewontin proposed that the ability of a scientific field to create a narrative for public understanding garners it social relevance. This article applies Lewontin's conceptual framework of the functions of science (manipulatory and explanatory) to compare and explain the current differences in perceived societal relevance of genetics/genomics and proteomics. We provide three examples to illustrate the social relevance and strong cultural narrative of genetics/genomics for which no counterpart exists for proteomics. We argue that the major difference between genetics/genomics and proteomics is that genomics has a strong explanatory function, due to the strong cultural narrative of heredity. Based on qualitative interviews and observations of proteomics conferences, we suggest that the nature of proteins, lack of public understanding, and theoretical complexity exacerbates this difference for proteomics. Lewontin's framework suggests that social scientists may find that omics sciences affect social relations in different ways than past analyses of genetics.