Genotyping for genetic association studies: Methods and applications

In this thesis, two separate single nucleotide polymorphism (SNP) genotyping techniques were set up at the Finnish Genome Center, pooled genotyping was evaluated as a screening method for large-scale association studies, and finally, the former approaches were used to identify genetic factors predisposing to two distinct complex diseases by utilizing large epidemiological cohorts and also taking environmental factors into account. The first genotyping platform was based on traditional but improved restriction-fragment-length-polymorphism (RFLP) utilizing 384-microtiter well plates, multiplexing, small reaction volumes (5 µl), and automated genotype calling. We participated in the development of the second genotyping method, based on single nucleotide primer extension (SNuPeTM by Amersham Biosciences), by carrying out the alpha- and beta tests for the chemistry and the allele-calling software. Both techniques proved to be accurate, reliable, and suitable for projects with thousands of samples and tens of markers. Pooled genotyping (genotyping of pooled instead of individual DNA samples) was evaluated with Sequenom s MassArray MALDI-TOF, in addition to SNuPeTM and PCR-RFLP techniques. We used MassArray mainly as a point of comparison, because it is known to be well suited for pooled genotyping. All three methods were shown to be accurate, the standard deviations between measurements being 0.017 for the MassArray, 0.022 for the PCR-RFLP, and 0.026 for the SNuPeTM. The largest source of error in the process of pooled genotyping was shown to be the volumetric error, i.e., the preparation of pools. We also demonstrated that it would have been possible to narrow down the genetic locus underlying congenital chloride diarrhea (CLD), an autosomal recessive disorder, by using the pooling technique instead of genotyping individual samples. Although the approach seems to be well suited for traditional case-control studies, it is difficult to apply if any kind of stratification based on environmental factors is needed. Therefore we chose to continue with individual genotyping in the following association studies. Samples in the two separate large epidemiological cohorts were genotyped with the PCR-RFLP and SNuPeTM techniques. The first of these association studies concerned various pregnancy complications among 100,000 consecutive pregnancies in Finland, of which we genotyped 2292 patients and controls, in addition to a population sample of 644 blood donors, with 7 polymorphisms in the potentially thrombotic genes. In this thesis, the analysis of a sub-study of pregnancy-related venous thromboses was included. We showed that the impact of factor V Leiden polymorphism on pregnancy-related venous thrombosis, but not the other tested polymorphisms, was fairly large (odds ratio 11.6; 95% CI 3.6-33.6), and increased multiplicatively when combined with other risk factors such as obesity or advanced age. Owing to our study design, we were also able to estimate the risks at the population level. The second epidemiological cohort was the Helsinki Birth Cohort of men and women who were born during 1924-1933 in Helsinki. The aim was to identify genetic factors that might modify the well known link between small birth size and adult metabolic diseases, such as type 2 diabetes and impaired glucose tolerance. Among ~500 individuals with detailed birth measurements and current metabolic profile, we found that an insertion/deletion polymorphism of the angiotensin converting enzyme (ACE) gene was associated with the duration of gestation, and weight and length at birth. Interestingly, the ACE insertion allele was also associated with higher indices of insulin secretion (p=0.0004) in adult life, but only among individuals who were born small (those among the lowest third of birth weight). Likewise, low birth weight was associated with higher indices of insulin secretion (p=0.003), but only among carriers of the ACE insertion allele. The association with birth measurements was also found with a common haplotype of the glucocorticoid receptor (GR) gene. Furthermore, the association between short length at birth and adult impaired glucose tolerance was confined to carriers of this haplotype (p=0.007). These associations exemplify the interaction between environmental factors and genotype, which, possibly due to altered gene expression, predisposes to complex metabolic diseases. Indeed, we showed that the common GR gene haplotype associated with reduced mRNA expression in thymus of three individuals (p=0.0002).

Oncolytic Adenoviruses for Gynecologic Cancer

Gene therapy is a promising novel approach for treating cancers resistant to or escaping currently available modalities. Treatment approaches are based on taking advantage of molecular differences between normal and tumor cells. Various strategies are currently in clinical development with adenoviruses as the most popular vehicle. Recent developments include improving targeting strategies for gene delivery to tumor cells with tumor specific promoters or infectivity enhancement. A rapidly developing field is as well replication competent agents, which allow improved tumor penetration and local amplification of the anti-tumor effect. Adenoviral cancer gene therapy approaches lack cross-resistance with other treatment options and therefore synergistic effects are possible. This study focused on development of adenoviral vectors suitable for treatment of various gynecologic cancer types, describing the development of the field from non-replicating adenoviral vectors to multiple-modified conditional replicating viruses. Transcriptional targeting of gynecologic cancer cells by the use of the promoter of vascular endothelial growth factor receptor type 1 (flt-1) was evaluated. Flt-1 is not expressed in the liver and thus an ideal promoter for transcriptional targeting of adenoviruses. Our studies implied that the flt-1 promoter is active in teratocarcinomas.and therefore a good candidate for development of oncolytic adenoviruses for treatment of this often problematic disease with then poor outcome. A tropism modified conditionally replicating adenovirus (CRAd), Ad5-Δ24RGD, was studied in gynecologic cancers. Ad5-Δ24RGD is an adenovirus selectively replication competent in cells defective in the p16/Rb pathway, including many or most tumor cells. The fiber of Ad5-Δ24RGD contains an integrin binding arginine-glycine-aspartic acid motif (RGD-4C), allowing coxackie-adenovirus receptor independent infection of cancer cells. This approach is attractive because expression levels of CAR are highly variable and often low on primary gynecological cancer cells. Oncolysis could be shown for a wide variety of ovarian and cervical cancer cell lines as well as primary ovarian cancer cell spheroids, a novel system developed for in vitro analysis of CRAds on primary tumor substrates. Biodistribution was evaluated and preclinical safety data was obtained by demonstrating lack of replication in human peripheral blood mononuclear cells. The efficicacy of Ad5-Δ24RGD was shown in different orthotopic murine models including a highly aggressive intraperitoneal model of disseminated ovarian cancer cells, where Ad5-Δ24RGD resulted in complete eradication of intraperitoneal disease in half of the mice. To further improve the selectivity and specificity of CRAds, triple-targeted oncolytic adenoviruses were cloned, featuring the cyclo-oxygenase-2 (cox-2) promoter, E1A transcomplementation and serotype chimerism. Those viruses were evaluated on ovarian cancer cells for specificity and oncolytic potency with regard to two different cox2 versions and three different variants of E1A (wild type, delta24 and delta2delta24). Ad5/3cox2Ld24 emerged as the best combination due to enhanced selectivity without potency lost in vitro or in an aggressive intraperitoneal orthotopic ovarian tumor model. In summary, the preclinical therapeutic efficacy of the CRAds tested in this study, taken together with promising biodistribution and safety data, suggest that these CRAds are interesting candidates for translation into clinical trials for gynecologic cancer.

Hybridisation of Australian chickpea cultivars with wild Cicer spp. increases resistance to root-lesion nematodes (Pratylenchus thornei and P. neglectus).

Australian and international chickpea (Cicer arietinum) cultivars and germplasm accessions, and wild annual Cicer spp. in the primary and secondary gene pools, were assessed in glasshouse experiments for levels of resistance to the root-lesion nematodes Pratylenchus thornei and P. neglectus. Lines were grown in replicated experiments in pasteurised soil inoculated with a pure culture of either P. thornei or P. neglectus and the population density of the nematodes in the soil plus roots after 16 weeks growth was used as a measure of resistance. Combined statistical analyses of experiments (nine for P. thornei and four for P. neglectus) were conducted and genotypes were assessed using best linear unbiased predictions. Australian and international chickpea cultivars possessed a similar range of susceptibilities through to partial resistance. Wild relatives from both the primary (C. reticulatum and C. echinospermum) and secondary (C. bijugum) gene pools of chickpea were generally more resistant than commercial chickpea cultivars to either P. thornei or P. neglectus or both. Wild relatives of chickpea have probably evolved to have resistance to endemic root-lesion nematodes whereas modern chickpea cultivars constitute a narrower gene pool with respect to nematode resistance. Resistant accessions of C. reticulatum and C. echinospermum were crossed and topcrossed with desi chickpea cultivars and resistant F(4) lines were obtained. Development of commercial cultivars with the high levels of resistance to P. thornei and P. neglectus in these hybrids will be most valuable for areas of the Australian grain region and other parts of the world where alternating chickpea and wheat crops are the preferred rotation.

Tree mortality and deadwood dynamics in late-successional boreal forests

Here I aimed at quantifying the main components of deadwood dynamics, i.e. tree mortality, deadwood pools, and their decomposition, in late-successional boreal forests. I focused on standing dead trees in three stand types dominated by Picea mariana and Abies balsamea in eastern Canada, and on standing and down dead trees in Picea abies-dominated stands in three areas in Northern Europe. Dead and living trees were measured on five sample plots of 1.6-ha size in each study area and stand type. Stem disks from dead trees were sampled to determine wood density and year of death, using dendrochronological methods. The results were applied to reconstruct past tree mortality and to model deadwood decay class dynamics. Site productivity, stand developmental stage, and the occurrence of episodic tree mortality influenced deadwood volume and quality. In all study areas tree mortality was continuous, leading to continuity in deadwood decay stage distribution. Episodic tree mortality due to either autogenic or allogenic causes influenced deadwood volume and quality in all but one study area. However, regardless of productivity and disturbance history deadwood was abundant, accounting for 20 53% of total wood volume in European study areas, and 15 27% of total standing volume in eastern Canada. Deadwood was a persistent structural component, since its expected residence time in early- and midstages of decay was 18 yr even in the area with the most rapid decomposition. The results indicated that in the absence of episodic tree mortality, stands may eventually develop to a steady state, in which deadwood volume fluctuates around an equilibrium state. However, in many forests deadwood is naturally variable, due to recurrent moderate-severity disturbances. This variability, the continuous tree mortality, and variation in rates of wood decomposition determine the dynamics and availability of deadwood as a habitat and carbon storage medium in boreal coniferous forest ecosystems.

Effect of fire and tree-grass patches on soil nitrogen in Australian tropical savannas

Fire is an important driver of nutrient cycling in savannas. Here, we determined the impact of fire frequency on total and soluble soil nitrogen (N) pools in tropical savanna. The study sites consisted of 1-ha experimental plots near Darwin, Australia, which remained unburnt for at least 14 years or were burnt at 1-, 2- or 5-year intervals over the past 6 years. Soil was analysed from patches underneath tree canopies and in inter-canopy patches at 1, 12, 28, 55 and 152 days after fire. Patch type had a significant effect on all soil N pools, with greater concentrations of total and soluble (nitrate, ammonium, amino acids) N under tree canopies than inter-canopy patches. The time since the last fire had no significant effect on N pools. Fire frequency similarly did not affect total soil N but it did influence soluble soil N. Soil amino acids were most prominent in burnt savanna, ammonium was highest in infrequently burnt (5-year interval) savanna and nitrate was highest in unburnt savanna. We suggest that the main effect of fire on soil N relations occurs indirectly through altered tree-grass dynamics. Previous studies have shown that high fire frequencies reduce tree cover by lowering recruitment and increasing mortality. Our findings suggest that these changes in tree cover could result in a 30% reduction in total soil N and 1060% reductions in soluble N pools. This finding is consistent with studies from savannas globally, providing further evidence for a general theory of patchiness as a key driver of nutrient cycling in the savanna biome.

Phosphorus in agricultural soils of Finland - characterization of reserves and retention in mineral soil profiles

An overwhelming majority of all the research on soil phosphorus (P) has been carried out with soil samples taken from the surface soils only, and our understanding of the forms and the reactions of P at a soil profile scale is based on few observations. In Finland, the interest in studying the P in complete soil profiles has been particularly small because of the lack of tradition in studying soil genesis, morphology, or classification. In this thesis, the P reserves and the retention of orthophosphate phosphorus (PO4-P) were examined in four cultivated mineral soil profiles in Finland (three Inceptisols and one Spodosol). The soils were classified according to the U.S. Soil Taxonomy and soil samples were taken from the genetic horizons in the profiles. The samples were analyzed for total P concentration, Chang and Jackson P fractions, P sorption properties, concentrations of water-extractable P, and for concentrations of oxalate-extractable Al and Fe. Theoretical P sorption capacities and degrees of P saturation were calculated with the data from the oxalate-extractions and the P fractionations. The studied profiles can be divided into sections with clearly differing P characteristics by their master horizons Ap, B and C. The C (or transitional BC) horizons below an approximate depth of 70 cm were dominated by, assumingly apatitic, H2SO4-soluble P. The concentration of total P in the C horizons ranged from 729 to 810 mg kg-1. In the B horizons between the depths of 30 and 70 cm, a significant part of the primary acid-soluble P has been weathered and transformed to secondary P forms. A mean weathering rate of the primary P in the soils was estimated to vary between 230 and 290 g ha-1 year-1. The degrees of P saturation in the B and C horizons were smaller than 7%, and the solubility of PO4-P was negligible. The P conditions in the Ap horizons differed drastically from those in the subsurface horizons. The high concentrations of total P (689-1870 mg kg-1) in the Ap horizons are most likely attributable to long-term cultivation with positive P balances. A significant proportion of the P in the Ap horizons occurred in the NH4F- and NaOH-extractable forms and as organic P. These three P pools, together with the concentrations of oxalate-extractable Al and Fe, seem to control the dynamics of PO4-P in the soils. The degrees of P saturation in the Ap horizons were greater (8-36%) than in the subsurface horizons. This was also reflected in the sorption experiments: Only the Ap horizons were able to maintain elevated PO4-P concentrations in the solution phase − all the subsoil horizons acted as sinks for PO4-P. Most of the available sorption capacity in the soils is located in the B horizons. The results suggest that this capacity could be utilized in reducing the losses of soluble P from excessively fertilized soils by mixing highly sorptive material from the B horizons with the P-enriched surface soil. The drastic differences in the P characteristics observed between adjoining horizons have to be taken into consideration when conducting soil sampling. Sampling of subsoils has to be made according to the genetic horizons or at small depth increments. Otherwise, contrasting materials are likely to be mixed in the same sample; and the results of such samples are not representative of any material present in the studied profile. Air-drying of soil samples was found to alter the results of the sorption experiments and the water extractions. This indicates that the studies on the most labile P forms in soil should be carried out with moist samples.

Estimating mineralisation of organic nitrogen from biosolids and other organic wastes applied to soils in subtropical Australia

One major benefit of land application of biosolids is to supply nitrogen (N) for agricultural crops, and understanding mineralisation processes is the key for better N-management strategies. Field studies were conducted to investigate the process of mineralisation of three biosolids products (aerobic, anaerobic, and thermally dried biosolids) incorporated into four different soils at rates of 7-90 wet t/ha in subtropical Queensland. Two of these studies also examined mineralisation rates of commonly used organic amendments (composts, manures, and sugarcane mill muds). Organic N in all biosolids products mineralised very rapidly under ambient conditions in subtropical Queensland, with rates much faster than from other common amendments. Biosolids mineralisation rates ranged from 30 to 80% of applied N during periods ranging from 3.5 to 18 months after biosolids application; these rates were much higher than those suggested in the biosolids land application guidelines established by the NSW EPA (15% for anaerobic and 25% for aerobic biosolids). There was no consistently significant difference in mineralisation rate between aerobic and anaerobic biosolids in our studies. When applied at similar rates of N addition, other organic amendments supplied much less N to the soil mineral N and plant N pools during the crop season. A significant proportion of the applied biosolids total N (up to 60%) was unaccounted for at the end of the observation period. High rates of N addition in calculated Nitrogen Limited Biosolids Application Rates (850-1250 kg N/ha) resulted in excessive accumulation of mineral N in the soil profile, which increases the environmental risks due to leaching, runoff, or gaseous N losses. Moreover, the rapid mineralisation of the biosolids organic N in these subtropical environments suggests that biosolids should be applied at lower rates than in temperate areas, and that care must be taken with the timing to maximise plant uptake and minimise possible leaching, runoff, or denitrification losses of mineralised N.

Concentrations of phthalates and DINCH metabolites in pooled urine from Queensland, Australia

Dialkyl phthalate esters (phthalates) are ubiquitous chemicals used extensively as plasticizers, solvents and adhesives in a range of industrial and consumer products. 1,2-Cyclohexane dicarboxylic acid, diisononyl ester (DINCH) is a phthalate alternative introduced due to a more favourable toxicological profile, but exposure is largely uncharacterised. The aim of this study was to provide the first assessment of exposure to phthalates and DINCH in the general Australian population. De-identified urine specimens stratified by age and sex were obtained from a community-based pathology laboratory and pooled (n = 24 pools of 100). Concentrations of free and total species were measured using online solid phase extraction isotope dilution high performance liquid chromatography tandem mass spectrometry. Concentrations ranged from 2.4 to 71.9 ng/mL for metabolites of di(2-ethylhexyl)phthalate, and from < 0.5 to 775 ng/mL for all other metabolites. Our data suggest that phthalate metabolites concentrations in Australia were at least two times higher than in the United States and Germany; and may be related to legislative differences among countries. DINCH metabolite concentrations were comparatively low and consistent with the limited data available. Ongoing biomonitoring among the general Australian population may help assess temporal trends in exposure and assess the effectiveness of actions aimed at reducing exposures.

Monitoring exposure to polycyclic aromatic hydrocarbons in an Australian population using pooled urine samples

Integrated exposure to polycyclic aromatic hydrocarbons (PAHs) can be assessed through monitoring of urinary mono-hydroxylated PAHs (OH-PAHs). The aim of this study was to provide the first assessment of exposure to PAHs in a large sample of the population in Queensland, Australia including exposure to infant (0-4. years). De-identified urine specimens, obtained from a pathology laboratory, were stratified by age and sex, and pooled (n. =. 24 pools of 100) and OH-PAHs were measured by gas chromatography-isotope dilution-tandem mass spectrometry. Geometric mean (GM) concentrations ranged from 30. ng/L (4-hydroxyphenanthrene) to 9221. ng/L (1-naphthol). GM of 1-hydroxypyrene, the most commonly used PAH exposure biomarker, was 142. ng/L. The concentrations of OH-PAHs found in this study are consistent with those in developed countries and lower than those in developing countries. We observed no association between sex and OH-PAH concentrations. However, we observed lower urinary concentrations of all OH-PAHs in samples from infants (0-4. years), children (5-14. years) and the elderly (>. 60. year old) compared with samples from other age groups (15-29, 30-44 and 45-59. years) which may be attributed to age-dependent behaviour-specific exposure sources.

A comparison of self-nominated and actual speeds in work zones

Despite significant research on drivers’ speeding behavior in work zones, little is known about how well drivers’ judgments of appropriate speeds match their actual speeds and what factors influence their judgments. This study aims to fill these two important gaps in the literature by comparing observed speeds in two work zones with drivers’ self-nominated speeds for the same work zones. In an online survey, drivers nominated speeds for the two work zones based on photographs in which the actual posted speed limits were not revealed. A simultaneous equation modelling approach was employed to examine the effects of driver characteristics on their self-nominated speeds. The results showed that survey participants nominated lower speeds (corresponding to higher compliance rates) than those which were observed. Higher speeds were nominated by males than females, young and middle aged drivers than older drivers, and drivers with truck driving experience than those who drive only cars. Larger differences between nominated and observed speeds were found among car drivers than truck drivers. These differences suggest that self-nominated speeds might not be valid indicators of the observed work zone speeds and therefore should not be used as an alternative to observed speed data.