Interpreting Human Genetic Variation through Genomics and In Vivo Models

Improvements in genomic technology, both in the increased speed and reduced cost of sequencing, have expanded the appreciation of the abundance of human genetic variation. However the sheer amount of variation, as well as the varying type and genomic content of variation, poses a challenge in understanding the clinical consequence of a single mutation. This work uses several methodologies to interpret the observed variation in the human genome, and presents novel strategies for the prediction of allele pathogenicity.

Using the zebrafish model system as an in vivo assay of allele function, we identified a novel driver of Bardet-Biedl Syndrome (BBS) in CEP76. A combination of targeted sequencing of 785 cilia-associated genes in a cohort of BBS patients and subsequent in vivo functional assays recapitulating the human phenotype gave strong evidence for the role of CEP76 mutations in the pathology of an affected family. This portion of the work demonstrated the necessity of functional testing in validating disease-associated mutations, and added to the catalogue of known BBS disease genes.

Further study into the role of copy-number variations (CNVs) in a cohort of BBS patients showed the significant contribution of CNVs to disease pathology. Using high-density array comparative genomic hybridization (aCGH) we were able to identify pathogenic CNVs as small as several hundred bp. Dissection of constituent gene and in vivo experiments investigating epistatic interactions between affected genes allowed for an appreciation of several paradigms by which CNVs can contribute to disease. This study revealed that the contribution of CNVs to disease in BBS patients is much higher than previously expected, and demonstrated the necessity of consideration of CNV contribution in future (and retrospective) investigations of human genetic disease.

Finally, we used a combination of comparative genomics and in vivo complementation assays to identify second-site compensatory modification of pathogenic alleles. These pathogenic alleles, which are found compensated in other species (termed compensated pathogenic deviations [CPDs]), represent a significant fraction (from 3 – 10%) of human disease-associated alleles. In silico pathogenicity prediction algorithms, a valuable method of allele prioritization, often misrepresent these alleles as benign, leading to omission of possibly informative variants in studies of human genetic disease. We created a mathematical model that was able to predict CPDs and putative compensatory sites, and functionally showed in vivo that second-site mutation can mitigate the pathogenicity of disease alleles. Additionally, we made publically available an in silico module for the prediction of CPDs and modifier sites.

These studies have advanced the ability to interpret the pathogenicity of multiple types of human variation, as well as made available tools for others to do so as well.

Dissecting the Functional Impacts of Non-Coding Genetic Variation

A large proportion of the variation in traits between individuals can be attributed to variation in the nucleotide sequence of the genome. The most commonly studied traits in human genetics are related to disease and disease susceptibility. Although scientists have identified genetic causes for over 4,000 monogenic diseases, the underlying mechanisms of many highly prevalent multifactorial inheritance disorders such as diabetes, obesity, and cardiovascular disease remain largely unknown. Identifying genetic mechanisms for complex traits has been challenging because most of the variants are located outside of protein-coding regions, and determining the effects of such non-coding variants remains difficult. In this dissertation, I evaluate the hypothesis that such non-coding variants contribute to human traits and diseases by altering the regulation of genes rather than the sequence of those genes. I will specifically focus on studies to determine the functional impacts of genetic variation associated with two related complex traits: gestational hyperglycemia and fetal adiposity. At the genomic locus associated with maternal hyperglycemia, we found that genetic variation in regulatory elements altered the expression of the HKDC1 gene. Furthermore, we demonstrated that HKDC1 phosphorylates glucose in vitro and in vivo, thus demonstrating that HKDC1 is a fifth human hexokinase gene. At the fetal-adiposity associated locus, we identified variants that likely alter VEPH1 expression in preadipocytes during differentiation. To make such studies of regulatory variation high-throughput and routine, we developed POP-STARR, a novel high throughput reporter assay that can empirically measure the effects of regulatory variants directly from patient DNA. By combining targeted genome capture technologies with STARR-seq, we assayed thousands of haplotypes from 760 individuals in a single experiment. We subsequently used POP-STARR to identify three key features of regulatory variants: that regulatory variants typically have weak effects on gene expression; that the effects of regulatory variants are often coordinated with respect to disease-risk, suggesting a general mechanism by which the weak effects can together have phenotypic impact; and that nucleotide transversions have larger impacts on enhancer activity than transitions. Together, the findings presented here demonstrate successful strategies for determining the regulatory mechanisms underlying genetic associations with human traits and diseases, and value of doing so for driving novel biological discovery.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

Genetic variation in growth and wood-quality traits of Corymbia citriodora subsp. variegata across three sites in south-east Queensland, Australia

Ten growth or wood-quality traits were assessed in three nearby Corymbia citriodora subsp. variegata (CCV) open-pollinated family-within-provenance trials (18 provenances represented by a total of 374 families) to provide information for the development of a breeding program targeting both pulp and solid-wood products. Growth traits (diameter at breast high over bark [DBH], height and conical volume) were assessed at 3 and 7 years of age. Wood-quality traits (density [DEN], Kraft pulp yield [KPY], modulus of elasticity [MoE] and microfibril angle [MfA]) were predicted using near-infrared spectroscopy on wood samples collected from these trials when aged between 10 and 12 years. The high average KPY, DEN and MoE, and low average MfA observed indicates CCV is very suitable for both pulp and timber products. All traits were under moderate to strong genetic control. In across- trials analyses, high (>0.4) heritability estimates were observed for height, DEN, MoE and MfA, while moderate heritability estimates (0.24 to 0.34) were observed for DBH, volume and KPY. Most traits showed very low levels of genotype × site interaction. Estimated age–age genetic correlations for growth traits were strong at both the family (0.97) and provenance (0.99) levels. Relationships among traits (additive genetic correlation estimates) were favourable, with strong and positive estimates between growth traits (0.84 to 0.98), moderate and positive values between growth and wood-quality traits (0.32 to 0.68), moderate and positive between KPY and MoE (0.64), and high and positive between DEN and MoE (0.82). However, negative (but favourable) correlations were detected between MfA and all other evaluated traits (−0.31 to −0.96). The genetic correlation between the same trait expressed on two different sites, at family level, ranged from 0.24 to 0.42 for growth traits, and from 0.29 to 0.53 for wood traits. Therefore simultaneous genetic improvement of growth and wood property traits in CCV for the target environment in south-east Queensland should be possible, given the moderate to high estimates of heritability and favourable correlations amongst all traits studied, unless genotype × site interactions are greater than was evident. © 2016 NISC (Pty) Ltd

Low-frequency and common genetic variation in ischemic stroke: The METASTROKE collaboration

Erratum in: Low-frequency and common genetic variation in ischemic stroke: The METASTROKE collaboration. [Neurology. 2016]

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

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

GENETIC VARIATION, LOCAL ADAPTATION AND POPULATION STRUCTURE IN NORTH AMERICAN RED OAK SPECIES, QUERCUS RUBRA L. AND Q. ELLIPSOIDALIS E. J. HILL

Forest trees, like oaks, rely on high levels of genetic variation to adapt to varying environmental conditions. Thus, genetic variation and its distribution are important for the long-term survival and adaptability of oak populations. Climate change is projected to lead to increased drought and fire events as well as a northward migration of tree species, including oaks. Additionally, decline in oak regeneration has become increasingly concerning since it may lead to decreased gene flow and increased inbreeding levels. This will in turn lead to lowered levels of genetic diversity, negatively affecting the growth and survival of populations. At the same time, populations at the species’ distribution edge, like those in this study, could possess important stores of genetic diversity and adaptive potential, while also being vulnerable to climatic or anthropogenic changes. A survey of the level and distribution of genetic variation and identification of potentially adaptive genes is needed since adaptive genetic variation is essential for their long-term survival. Oaks possess a remarkable characteristic in that they maintain their species identity and specific environmental adaptations despite their propensity to hybridize. Thus, in the face of interspecific gene flow, some areas of the genome remain differentiated due to selection. This characteristic allows the study of local environmental adaptation through genetic variation analyses. Furthermore, using genic markers with known putative functions makes it possible to link those differentiated markers to potential adaptive traits (e.g., flowering time, drought stress tolerance). Demographic processes like gene flow and genetic drift also play an important role in how genes (including adaptive genes) are maintained or spread. These processes are influenced by disturbances, both natural and anthropogenic. An examination of how genetic variation is geographically distributed can display how these genetic processes and geographical disturbances influence genetic variation patterns. For example, the spatial clustering of closely related trees could promote inbreeding with associated negative effects (inbreeding depression), if gene flow is limited. In turn this can have negative consequences for a species’ ability to adapt to changing environmental conditions. In contrast, interspecific hybridization may also allow the transfer of genes between species that increase their adaptive potential in a changing environment. I have studied the ecologically divergent, interfertile red oaks, Quercus rubra and Q. ellipsoidalis, to identify genes with potential roles in adaptation to abiotic stress through traits such as drought tolerance and flowering time, and to assess the level and distribution of genetic variation. I found evidence for moderate gene flow between the two species and low interspecific genetic differences at most genetic markers (Lind and Gailing 2013). However, the screening of genic markers with potential roles in phenology and drought tolerance led to the identification of a CONSTANS-like (COL) gene, a candidate gene for flowering time and growth. This marker, located in the coding region of the gene, was highly differentiated between the two species in multiple geographical areas, despite interspecific gene flow, and may play a role in reproductive isolation and adaptive divergence between the two species (Lind-Riehl et al. 2014). Since climate change could result in a northward migration of trees species like oaks, this gene could be important in maintaining species identity despite increased contact zones between species (e.g., increased gene flow). Finally I examined differences in spatial genetic structure (SGS) and genetic variation between species and populations subjected to different management strategies and natural disturbances. Diverse management activities combined with various natural disturbances as well as species specific life history traits influenced SGS patterns and inbreeding levels (Lind-Riehl and Gailing submitted).

Genetic and environmental components of phenotypic variation in immune response and body size of a colonial bird, Delichon urbica (the house martin).

Directional selection for parasite resistance is often intense in highly social host species. Using a partial cross-fostering experiment we studied environmental and genetic variation in immune response and morphology in a highly colonial bird species, the house martin (Delichon urbica). We manipulated intensity of infestation of house martin nests by the haematophagous parasitic house martin bug Oeciacus hirundinis either by spraying nests with a weak pesticide or by inoculating them with 50 bugs. Parasitism significantly affected tarsus length, T cell response, immunoglobulin and leucocyte concentrations. We found evidence of strong environmental effects on nestling body mass, body condition, wing length and tarsus length, and evidence of significant additive genetic variance for wing length and haematocrit. We found significant environmental variance, but no significant additive genetic variance in immune response parameters such as T cell response to the antigenic phytohemagglutinin, immunoglobulins, and relative and absolute numbers of leucocytes. Environmental variances were generally greater than additive genetic variances, and the low heritabilities of phenotypic traits were mainly a consequence of large environmental variances and small additive genetic variances. Hence, highly social bird species such as the house martin, which are subject to intense selection by parasites, have a limited scope for immediate microevolutionary response to selection because of low heritabilities, but also a limited scope for long-term response to selection because evolvability as indicated by small additive genetic coefficients of variation is weak.

Genetic variability in a population of arbuscular mycorrhizal fungi causes variation in plant growth.

Different species of arbuscular mycorrhizal fungi (AMF) alter plant growth and affect plant coexistence and diversity. Effects of within-AMF species or within-population variation on plant growth have received less attention. High genetic variation exists within AMF populations. However, it is unknown whether genetic variation contributes to differences in plant growth. In our study, a population of AMF was cultivated under identical conditions for several generations prior to the experiments thus avoiding environmental maternal effects. We show that genetically different Glomus intraradices isolates from one AMF population significantly alter plant growth in an axenic system and in greenhouse experiments. Isolates increased or reduced plant growth meaning that plants potentially receive benefits or are subject to costs by forming associations with different individuals in the AMF population. This shows that genetic variability in AMF populations could affect host-plant fitness and should be considered in future research to understand these important soil organisms.

Progressive erosion of genetic and epigenetic variation in callus-derived cocoa (Theobroma cacao) plants

Relatively little is known about the timing of genetic and epigenetic forms of somaclonal variation arising from callus growth. We surveyed for both types of change in cocoa (Theobroma cacao) plants regenerated from calli of various ages, and also between tissues from the source trees. For genetic change, we used 15 single sequence repeat (SSR) markers from four source trees and from 233 regenerated plants. For epigenetic change, we used 386 methylation-sensitive amplified polymorphism (MSAP) markers on leaf and explant (staminode) DNA from two source trees and on leaf DNA from 114 regenerants. Genetic variation within source trees was limited to one slippage mutation in one leaf. Regenerants were far more variable, with 35% exhibiting at least one mutation. Genetic variation initially accumulated with culture age but subsequently declined. MSAP (epigenetic) profiles diverged between leaf and staminode samples from source trees. Multivariate analysis revealed that leaves from regenerants occupied intermediate eigenspace between leaves and staminodes of source plants but became progressively more similar to source tree leaves with culture age. Statistical analysis confirmed this rather counterintuitive finding that leaves of ‘late regenerants’ exhibited significantly less genetic and epigenetic divergence from source leaves than those exposed to short periods of callus growth.

Genome-wide association mapping in a wild avian population identifies a link between genetic and phenotypic variation in a life-history trait

Understanding the genetic basis of traits involved in adaptation is a major challenge in evolutionary biology but remains poorly understood. Here, we use genome-wide association mapping using a custom 50 k single nucleotide polymorphism (SNP) array in a natural population of collared flycatchers to examine the genetic basis of clutch size, an important life-history trait in many animal species. We found evidence for an association on chromosome 18 where one SNP significant at the genome-wide level explained 3.9% of the phenotypic variance. We also detected two suggestive quantitative trait loci (QTLs) on chromosomes 9 and 26. Fitness differences among genotypes were generally weak and not significant, although there was some indication of a sex-by-genotype interaction for lifetime reproductive success at the suggestive QTL on chromosome 26. This implies that sexual antagonism may play a role in maintaining genetic variation at this QTL. Our findings provide candidate regions for a classic avian life-history trait that will be useful for future studies examining the molecular and cellular function of, as well as evolutionary mechanisms operating at, these loci.

Genetic and morphological variation and differentiation of Raja clavata populations in the European seas: marker development and preliminary data

In this study the population structure and connectivity of the Mediterranean and Atlantic Raja clavata (L., 1758) were investigated by analyzing the genetic variation of six population samples (N = 144) at seven nuclear microsatellite loci. The genetic dataset was generated by selecting population samples available in the tissue databases of the GenoDREAM laboratory (University of Bologna) and of the Department of Life Sciences and Environment (University of Cagliari), all collected during past scientific surveys (MEDITS, GRUND) from different geographical locations in the Mediterranean basin and North-east Atlantic sea, as North Sea, Sardinian coasts, Tuscany coasts and Cyprus Island. This thesis deals with to estimate the genetic diversity and differentiation among 6 geographical samples, in particular, to assess the presence of any barrier (geographic, hydrogeological or biological) to gene flow evaluating both the genetic diversity (nucleotide diversity, observed and expected heterozygosity, Hardy- Weinberg equilibrium analysis) and population differentiation (Fst estimates, population structure analysis). In addition to molecular analysis, quantitative representation and statistical analysis of morphological individuals shape are performed using geometric morphometrics methods and statistical tests. Geometric coordinates call landmarks are fixed in 158 individuals belonging to two population samples of Raja clavata and in population samples of closely related species, Raja straeleni (cryptic sibling) and Raja asterias, to assess significant morphological differences at multiple taxonomic levels. The results obtained from the analysis of the microsatellite dataset suggested a geographic and genetic separation between populations from Central-Western and Eastern Mediterranean basins. Furthermore, the analysis also showed that there was no separation between geographic samples from North Atlantic Ocean and central-Western Mediterranean, grouping them to a panmictic population. The Landmark-based geometric morphometry method results showed significant differences of body shape able to discriminate taxa at tested levels (from species to populations).

Genetic basis of variation for root traits and response to heat stress in durum wheat

Durum wheat is the second most important wheat species worldwide and the most important crop in several Mediterranean countries including Italy. Durum wheat is primarily grown under rainfed conditions where episodes of drought and heat stress are major factors limiting grain yield. The research presented in this thesis aimed at the identification of traits and genes that underlie root system architecture (RSA) and tolerance to heat stress in durum wheat, in order to eventually contribute to the genetic improvement of this species. In the first two experiments we aimed at the identification of QTLs for root trait architecture at the seedling level by studying a bi-parental population of 176 recombinant inbred lines (from the cross Meridiano x Claudio) and a collection of 183 durum elite accessions. Forty-eight novel QTLs for RSA traits were identified in each of the two experiments, by means of linkage- and association mapping-based QTL analysis, respectively. Important QTLs controlling the angle of root growth in the seedling were identified. In a third experiment, we investigated the phenotypic variation of root anatomical traits by means of microscope-based analysis of root cross sections in 10 elite durum cultivars. The results showed the presence of sizeable genetic variation in aerenchyma-related traits, prompting for additional studies aimed at mapping the QTLs governing such variation and to test the role of aerenchyma in the adaptive response to abiotic stresses. In the fourth experiment, an association mapping experiment for cell membrane stability at the seedling stage (as a proxy trait for heat tolerance) was carried out by means of association mapping. A total of 34 QTLs (including five major ones), were detected. Our study provides information on QTLs for root architecture and heat tolerance which could potentially be considered in durum wheat breeding programs.

Are there genetic breaks between Atlantic and Pacific Yellowfin tuna (Thunnus albacares) populations? A preliminary study based on microsatellites gene variation

Yellowfin tuna (Thunnus albacares, YFT, Bonnaterre 1788) is one of the most important market tuna species in the world. The high mortality of juveniles is in part caused by their bycatch. Indeed, if unregulated, it could permanently destabilize stocks health. For this reason investigating and better knowing the stock boundaries represent a crucial concern. Aim of this thesis was to preliminary investigate the YFT population structure within and between Atlantic and Pacific Oceans through the analysis of genetic variation at eight microsatellite loci and assess the occurrence of barriers to the gene flow between Oceans. For this propouse we collected 4 geographical samples coming from Atlantic and Pacific Ocean and selected a panel of 8 microsatellites loci developped by Antoni et al., (2014). Samples 71-2-Y and 77-2-Y, came from rispectively west central pacific ocean (WCPO) and east central pacific ocean (ECPO), instead samples 41-1-Y and 34-2-Y derive from west central atlantic ocean (WCAO) and east central atlantic ocean (ECAO). Total 160 specimens were analyzed (40 per sample) and were carried out several genetic information as allele frequencies, allele number, allelic richness, HWE (using He and Ho) and pairwise Fst genetic distance. Results obtained, may support the panmictic theory of this species, only one of pairwise Fst obtained is statistically significant (Fst= 0.00927; pV= 0.00218) between 41-1-Y and 71-2-Y samples. Results suggest low genetic differentiation and consequent high level of gene flow between Atlantic and Pacific populations. Furthermore, we performed an analysis of molecular taxonomy through the use of ATCO (the flaking region between ATPse6 and cytochrome oxidase subunit III genes mt DNA, to discriminate within the gener Thunnus two of the related species (Yellofin and bigeye tuna) according with their difficult recognition at certain size (<40 cm). ATCO analysis in this thesis, has provided strong discriminate evidence between the target species proving to be one of the most reliable genetic tools capable to indagate within the genus Thunnus. Thus, our study has provided useful information for possible use of this protocol for conservation plans and management of this fish stocks.

Diversification in the South American Pampas: the genetic and morphological variation of the widespread Petunia axillaris complex (Solanaceae)

Understanding the spatiotemporal distribution of genetic variation and the ways in which this distribution is connected to the ecological context of natural populations is fundamental for understanding the nature and mode of intraspecific and, ultimately, interspecific differentiation. The Petunia axillaris complex is endemic to the grasslands of southern South America and includes three subspecies: P.a.axillaris, P.a.parodii and P.a.subandina. These subspecies are traditionally delimited based on both geography and floral morphology, although the latter is highly variable. Here, we determined the patterns of genetic (nuclear and cpDNA), morphological and ecological (bioclimatic) variation of a large number of P.axillaris populations and found that they are mostly coincident with subspecies delimitation. The nuclear data suggest that the subspecies are likely independent evolutionary units, and their morphological differences may be associated with local adaptations to diverse climatic and/or edaphic conditions and population isolation. The demographic dynamics over time estimated by skyline plot analyses showed different patterns for each subspecies in the last 100000years, which is compatible with a divergence time between 35000 and 107000years ago between P.a.axillaris and P.a.parodii, as estimated with the IMa program. Coalescent simulation tests using Approximate Bayesian Computation do not support previous suggestions of extensive gene flow between P.a.axillaris and P.a.parodii in their contact zone.