Epidemiology and genetic architecture of blood pressure: a family based study of Generation Scotland

Hypertension is a major risk factor for cardiovascular disease and mortality, and a growing global public health concern, with up to one-third of the world’s population affected. Despite the vast amount of evidence for the benefits of blood pressure (BP) lowering accumulated to date, elevated BP is still the leading risk factor for disease and disability worldwide. It is well established that hypertension and BP are common complex traits, where multiple genetic and environmental factors contribute to BP variation. Furthermore, family and twin studies confirmed the genetic component of BP, with a heritability estimate in the range of 30-50%. Contemporary genomic tools enabling the genotyping of millions of genetic variants across the human genome in an efficient, reliable, and cost-effective manner, has transformed hypertension genetics research. This is accompanied by the presence of international consortia that have offered unprecedentedly large sample sizes for genome-wide association studies (GWASs). While GWAS for hypertension and BP have identified more than 60 loci, variants in these loci are associated with modest effects on BP and in aggregate can explain less than 3% of the variance in BP. The aims of this thesis are to study the genetic and environmental factors that influence BP and hypertension traits in the Scottish population, by performing several genetic epidemiological analyses. In the first part of this thesis, it aims to study the burden of hypertension in the Scottish population, along with assessing the familial aggregation and heritialbity of BP and hypertension traits. In the second part, it aims to validate the association of common SNPs reported in the large GWAS and to estimate the variance explained by these variants. In this thesis, comprehensive genetic epidemiology analyses were performed on Generation Scotland: Scottish Family Health Study (GS:SFHS), one of the largest population-based family design studies. The availability of clinical, biological samples, self-reported information, and medical records for study participants has allowed several assessments to be performed to evaluate factors that influence BP variation in the Scottish population. Of the 20,753 subjects genotyped in the study, a total of 18,470 individuals (grouped into 7,025 extended families) passed the stringent quality control (QC) criteria and were available for all subsequent analysis. Based on the BP-lowering treatment exposure sources, subjects were further classified into two groups. First, subjects with both a self-reported medications (SRMs) history and electronic-prescription records (EPRs; n =12,347); second, all the subjects with at least one medication history source (n =18,470). In the first group, the analysis showed a good concordance between SRMs and EPRs (kappa =71%), indicating that SRMs can be used as a surrogate to assess the exposure to BP-lowering medication in GS:SFHS participants. Although both sources suffer from some limitations, SRMs can be considered the best available source to estimate the drug exposure history in those without EPRs. The prevalence of hypertension was 40.8% with higher prevalence in men (46.3%) compared to women (35.8%). The prevalence of awareness, treatment and controlled hypertension as defined by the study definition were 25.3%, 31.2%, and 54.3%, respectively. These findings are lower than similar reported studies in other populations, with the exception of controlled hypertension prevalence, which can be considered better than other populations. Odds of hypertension were higher in men, obese or overweight individuals, people with a parental history of hypertension, and those living in the most deprived area of Scotland. On the other hand, deprivation was associated with higher odds of treatment, awareness and controlled hypertension, suggesting that people living in the most deprived area may have been receiving better quality of care, or have higher comorbidity levels requiring greater engagement with doctors. These findings highlight the need for further work to improve hypertension management in Scotland. The family design of GS:SFHS has allowed family-based analysis to be performed to assess the familial aggregation and heritability of BP and hypertension traits. The familial correlation of BP traits ranged from 0.07 to 0.20, and from 0.18 to 0.34 for parent-offspring pairs and sibling pairs, respectively. A higher correlation of BP traits was observed among first-degree relatives than other types of relative pairs. A variance-component model that was adjusted for sex, body mass index (BMI), age, and age-squared was used to estimate heritability of BP traits, which ranged from 24% to 32% with pulse pressure (PP) having the lowest estimates. The genetic correlation between BP traits showed a high correlation between systolic (SBP), diastolic (DBP) and mean arterial pressure (MAP) (G: 81% to 94%), but lower correlations with PP (G: 22% to 78%). The sibling recurrence risk ratio (λS) for hypertension and treatment were calculated as 1.60 and 2.04 respectively. These findings confirm the genetic components of BP traits in GS:SFHS, and justify further work to investigate genetic determinants of BP. Genetic variants reported in the recent large GWAS of BP traits were selected for genotyping in GS:SFHS using a custom designed TaqMan® OpenArray®. The genotyping plate included 44 single nucleotide polymorphisms (SNPs) that have been previously reported to be associated with BP or hypertension at genome-wide significance level. A linear mixed model that is adjusted for age, age-squared, sex, and BMI was used to test for the association between the genetic variants and BP traits. Of the 43 variants that passed the QC, 11 variants showed statistically significant association with at least one BP trait. The phenotypic variance explained by these variant for the four BP traits were 1.4%, 1.5%, 1.6%, and 0.8% for SBP, DBP, MAP, and PP, respectively. The association of genetic risk score (GRS) that were constructed from selected variants has showed a positive association with BP level and hypertension prevalence, with an average effect of one mmHg increase with each 0.80 unit increases in the GRS across the different BP traits. The impact of BP-lowering medication on the genetic association study for BP traits has been established, with typical practice of adding a fixed value (i.e. 15/10 mmHg) to the measured BP values to adjust for BP treatment. Using the subset of participants with the two treatment exposure sources (i.e. SRMs and EPRs), the influence of using either source to justify the addition of fixed values in SNP association signal was analysed. BP phenotypes derived from EPRs were considered the true phenotypes, and those derived from SRMs were considered less accurate, with some phenotypic noise. Comparing SNPs association signals between the four BP traits in the two model derived from the different adjustments showed that MAP was the least impacted by the phenotypic noise. This was suggested by identifying the same overlapped significant SNPs for the two models in the case of MAP, while other BP traits had some discrepancy between the two sources

Small-scale genetic structure of Cerastoderma glaucum in a lagoonal environment: potential significance of habitat discontinuity and unstable population dynamics

Environmental heterogeneity in coastal lagoons is expected to facilitate local adaptation in response to different ecological conditions, causing significant genetic structuring within lagoon populations at a small scale and also differentiation between lagoons. However, these patterns and processes of genetic structuring are still poorly understood. The aims of our study were (1) to seek genetic structure at a small scale in Cerastoderma glaucum inside the Mar Menor coastal lagoon using a mitochondrial DNA marker (COI) that has previously detected genetic differentiation inside the lagoon in other species and (2) to evaluate the influence of extreme environmental conditions and habitat discontinuity on its genetic composition. The results indicate high levels of haplotype diversity and low values of nucleotide diversity. COI data provide evidence of significant population differentiation among some localities within the lagoon. Limited gene flow and unstable population dynamics (i.e. fluctuations in population size caused by local extinction and recolonization), probably due to the high environmental heterogeneity, could generate the small-scale genetic divergence detected between populations within the lagoon.

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 Structure of the Florida Key Tree Cactus, Pilosocereus robinii, using Restriction Site associated DNA (RAD) markers

Rare plant conservation efforts must utilize current genetic methods to ensure the evolutionary potential of populations is preserved. One such effort involves the Key Tree Cactus, Pilosocereus robinii, which is an endangered columnar cactus native to the Florida Keys. The populations have precipitously declined over the past decade because of habitat loss and increasing soil salinity from rising sea levels and storm surge. Next-generation DNA sequencing was used to assess the genetic structure of the populations. Twenty individuals representative of both wild and extirpated cacti were chosen for Restriction Site Associated DNA (RAD) analysis. Samples processed using the HindIII and NotIII restriction enzymes produced 82,382,440 high quality reads used for genetic mapping, from which 5,265 Single Nucleotide Polymorphisms (SNPs) were discovered. The analysis revealed that the Keys’ populations are closely related with little population differentiation. In addition, the populations display evidence of inbreeding and low genetic diversity.

Comparative Genomics of Field Isolates of Mycobacterium bovis and M. caprae Provides Evidence for Possible Correlates with Bacterial Viability and Virulence

Mycobacteria of the Mycobacterium tuberculosis complex (MTBC) greatly affect humans and animals worldwide. The life cycle of mycobacteria is complex and the mechanisms resulting in pathogen infection and survival in host cells are not fully understood. Recently, comparative genomics analyses have provided new insights into the evolution and adaptation of the MTBC to survive inside the host. However, most of this information has been obtained using M. tuberculosis but not other members of the MTBC such as M. bovis and M. caprae. In this study, the genome of three M. bovis (MB1, MB3, MB4) and one M. caprae (MB2) field isolates with different lesion score, prevalence and host distribution phenotypes were sequenced. Genome sequence information was used for whole-genome and protein-targeted comparative genomics analysis with the aim of finding correlates with phenotypic variation with potential implications for tuberculosis (TB) disease risk assessment and control. At the whole-genome level the results of the first comparative genomics study of field isolates of M. bovis including M. caprae showed that as previously reported for M. tuberculosis, sequential chromosomal nucleotide substitutions were the main driver of the M. bovis genome evolution. The phylogenetic analysis provided a strong support for the M. bovis/M. caprae clade, but supported M. caprae as a separate species. The comparison of the MB1 and MB4 isolates revealed differences in genome sequence, including gene families that are important for bacterial infection and transmission, thus highlighting differences with functional implications between isolates otherwise classified with the same spoligotype. Strategic protein-targeted analysis using the ESX or type VII secretion system, proteins linking stress response with lipid metabolism, host T cell epitopes of mycobacteria, antigens and peptidoglycan assembly protein identified new genetic markers and candidate vaccine antigens that warrant further study to develop tools to evaluate risks for TB disease caused by M. bovis/M.caprae and for TB control in humans and animals.

Occurrence of Mycobacterium avium subspecies paratuberculosis across host species and European countries with evidence for transmission between wildlife and domestic ruminants

BACKGROUND Mycobacterium avium subspecies paratuberculosis (Map) causes an infectious chronic enteritis (paratuberculosis or Johne's disease) principally of ruminants. The epidemiology of Map is poorly understood, particularly with respect to the role of wildlife reservoirs and the controversial issue of zoonotic potential (Crohn's disease). Genotypic discrimination of Map isolates is pivotal to descriptive epidemiology and resolving these issues. This study was undertaken to determine the genetic diversity of Map, enhance our understanding of the host range and distribution and assess the potential for interspecies transmission. RESULTS 164 Map isolates from seven European countries representing 19 different host species were genotyped by standardized IS900--restriction fragment length polymorphism (IS900-RFLP), pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphisms (AFLP) and mycobacterial interspersed repeat unit-variable number tandem repeat (MIRU-VNTR) analyses. Six PstI and 17 BstEII IS900-RFLP, 31 multiplex [SnaBI-SpeI] PFGE profiles and 23 MIRU-VNTR profiles were detected. AFLP gave insufficient discrimination of isolates for meaningful genetic analysis. Point estimates for Simpson's index of diversity calculated for the individual typing techniques were in the range of 0.636 to 0.664 but a combination of all three methods increased the discriminating power to 0.879, sufficient for investigating transmission dynamics. Two predominant strain types were detected across Europe with all three typing techniques. Evidence for interspecies transmission between wildlife and domestic ruminants on the same property was demonstrated in four cases, between wildlife species on the same property in two cases and between different species of domestic livestock on one property. CONCLUSION The results of this study showed that it is necessary to use multiple genotyping techniques targeting different sources of genetic variation to obtain the level of discrimination necessary to investigate transmission dynamics and trace the source of Map infections. Furthermore, the combination of genotyping techniques may depend on the geographical location of the population to be tested. Identical genotypes were obtained from Map isolated from different host species co-habiting on the same property strongly suggesting that interspecies transmission occurs. Interspecies transmission of Map between wildlife species and domestic livestock on the same property provides further evidence to support a role for wildlife reservoirs of infection.

Genetic architecture and sexual conflict in the life history of Drosophila

Because males and females of a species express many homologous traits, sex-specific selection on these traits can shift the opposite sex away from its phenotypic optimum. This mode of sexually antagonistic selection, known as intralocus sexual conflict (IaSC), arises when the evolution of sexual dimorphism is constrained by the two sexes sharing a common gene pool. As IaSC has been historically overlooked, many outstanding questions remain. For example, what is its contribution in maintaining genetic variation for fitness in populations? What characters underlie this variation in fitness? How does the selection history of the population influence the standing genetic variation? I used the model organism Drosophila melanogaster to attempt to resolve some of these questions. The first part of my Master’s project involved assessing the detectability of sexually antagonistic alleles in populations at different stages of adaptation to the laboratory. For the second part of my Master’s project, I looked for evidence of conflict during the development of body size, a well-known sexually dimorphic trait. While the first part of my thesis proved inconclusive, the second part revealed a surprising source of sexual conflict in pre-adult stages of D. melanogaster.

Novel ciliate genetic code variants including the reassignment of all three stop codons to sense codons in Condylostoma magnum

mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum. We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3' end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum. We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3' end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.

Prevalence and genetic diversity of Blastocystis in family units living in the United States

The human gut is host to a diversity of microorganisms including the single-celled microbial eukaryote Blastocystis. Although Blastocystis has a global distribution, there is dearth of information relating to its prevalence and diversity in many human populations. The mode of Blastocystis transmission to humans is also insufficiently characterised, however, it is speculated to vary between different populations. Here we investigated the incidence and genetic diversity of Blastocystis in a US population and also the possibility of Blastocystis human-human transmission between healthy individuals using family units (N = 50) living in Boulder, Colorado as our sample-set. Ten of the 139 (~ 7%) individuals in our dataset were positive for Blastocystis, nine of whom were adults and one individual belonging to the children/adolescents group. All positive cases were present in different family units. A number of different Blastocystis subtypes (species) were detected with no evidence of mixed infections. The prevalence of Blastocystis in this subset of the US population is comparatively low relative to other industrialised populations investigated to date; however, subtype diversity was largely consistent with that previously reported in studies of European populations. The distribution of Blastocystis within family units indicates that human-human transmission is unlikely to have occurred within families that participated in this study. It is not unexpected that given the world-wide variation in human living conditions and lifestyles between different populations, both the prevalence of Blastocystis and its mode of transmission to humans may vary considerably.

Genetic signature of the northward expansion of the Egyptian mongoose Herpestes ichneumon (Herpestidae) in the Iberian Peninsula

In the last three decades, the range of the Egyptian mongoose (Herpestes ichneumon) has increased in the Iberian Peninsula. A panel of microsatellites was used to confront the patterns of genetic diversity of the species with the scenario of its recent northward expansion in its Iberian range. Evidence of substructure and significant genetic differentiation within the studied population were recorded, with a central-northern subpopulation (CNorth) and a southern subpopulation (S). Northward range expansion was supported by the observed allelic frequencies, diversity parameters, and observed heterozygosity of the studied loci, with S showing a higher allelic diversity and a higher number of private alleles than CNorth. Patterns of isolation-by-distance and isolation-by-barrier as a result of the Tagus River were demonstrated, suggesting that the river acted as a semi-permeable barrier, possibly leading to genetic differentiation of the studied population. The observed individuals from CNorth in southern locations and individuals from S in central/northern areas might comprise evidence for long range dispersals across the studied range. A bottleneck event after population expansion was supported by a significant heterozygosity deficiency in CNorth, which is in agreement with a scenario of founder events occurring in recently colonized areas after the crossing of the Tagus River.