Association mapping

Association mapping seeks to identify marker alleles present at significantly different frequencies in cases carrying a particular disease or trait compared with controls. Genome-wide association studies are increasingly replacing candidate gene-based association studies for complex diseases, where a number of loci are likely to contribute to disease risk and the effect size of each particular risk allele is typically modest or low. Good study design is essential to the success of an association study, and factors such as the heritability of the disease under investigation, the choice of controls, statistical power, multiple testing and whether the association can be replicated need to be considered before beginning. Likewise, thorough quality control of the genotype data needs to be undertaken prior to running any association analyses. Finally, it should be kept in mind that a significant genetic association is not proof positive that a particular genetic locus causes a disease, but rather an important first step in discovering the genetic variants underlying a complex disease.

A genome wide linkage scan for dizygotic twinning in 525 families of mothers of dizygotic twins

BACKGROUND: The tendency to conceive dizygotic (DZ) twins is a complex trait influenced by genetic and environmental factors. To search for new candidate loci for twinning, we conducted a genome-wide linkage scan in 525 families using microsatellite and single nucleotide polymorphism marker panels. METHODS AND RESULTS: Non-parametric linkage analyses, including 523 families containing a total of 1115 mothers of DZ twins (MODZT) from Australia and New Zealand (ANZ) and The Netherlands (NL), produced four linkage peaks above the threshold for suggestive linkage, including a highly suggestive peak at the extreme telomeric end of chromosome 6 with an exponential logarithm of odds \[(exp)LOD] score of 2.813 (P = 0.0002). Since the DZ twinning rate increases steeply with maternal age independent of genetic effects, we also investigated linkage including only families where at least one MODZT gave birth to her first set of twins before the age of 30. These analyses produced a maximum expLOD score of 2.718 (P = 0.0002), largely due to linkage signal from the ANZ cohort, however, ordered subset analyses indicated this result is most likely a chance finding in the combined dataset. Linkage analyses were also performed for two large DZ twinning families from the USA, one of which produced a peak on chromosome 2 in the region of two potential candidate genes. Sequencing of FSHR and FIGLA, along with INHBB in MODZTs from two large NL families with family specific linkage peaks directly over this gene, revealed a potentially functional variant in the 5' untranslated region of FSHR that segregated with the DZ twinning phenotype in the Utah family. CONCLUSION: Our data provide further evidence for complex inheritance of familial DZ twinning.

Identification of the semaphorin receptor PLXNA2 as a candidate for susceptibility to schizophrenia

The discovery of genetic factors that contribute to schizophrenia susceptibility is a key challenge in understanding the etiology of this disease. Here, we report the identification of a novel schizophrenia candidate gene on chromosome 1q32, plexin A2 (PLXNA2), in a genome-wide association study using 320 patients with schizophrenia of European descent and 325 matched controls. Over 25,000 single-nucleotide polymorphisms (SNPs) located within approximately 14,000 genes were tested. Out of 62 markers found to be associated with disease status, the most consistent finding was observed for a candidate locus on chromosome 1q32. The marker SNP rs752016 showed suggestive association with schizophrenia (odds ratio (OR) = 1.49, P = 0.006). This result was confirmed in an independent case-control sample of European Americans (combined OR = 1.38, P = 0.035) and similar genetic effects were observed in smaller subsets of Latin Americans (OR = 1.26) and Asian Americans (OR = 1.37). Supporting evidence was also obtained from two family-based collections, one of which reached statistical significance (OR = 2.2, P = 0.02). High-density SNP mapping showed that the region of association spans approximately 60 kb of the PLXNA2 gene. Eight out of 14 SNPs genotyped showed statistically significant differences between cases and controls. These results are in accordance with previous genetic findings that identified chromosome 1q32 as a candidate region for schizophrenia. PLXNA2 is a member of the transmembrane semaphorin receptor family that is involved in axonal guidance during development and may modulate neuronal plasticity and regeneration. The PLXNA2 ligand semaphorin 3A has been shown to be upregulated in the cerebellum of individuals with schizophrenia. These observations, together with the genetic results, make PLXNA2 a likely candidate for the 1q32 schizophrenia susceptibility locus.

Association study of the dystrobrevin-binding gene with schizophrenia in Australian and Indian samples

Numerous studies have reported association between variants in the dystrobrevin binding protein 1 (dysbindin) gene (DTNBP1) and schizophrenia. However, the pattern of results is complex and to date, no specific risk marker or haplotype has been consistently identified. The number of single nucleotide polymorphisms (SNPs) tested in these studies has ranged from 5 to 20. We attempted to replicate previous findings by testing 16 SNPs in samples of 41 Australian pedigrees, 194 Australian cases and 180 controls, and 197 Indian pedigrees. No globally significant evidence for association was observed in any sample, despite power calculations indicating sufficient power to replicate several previous findings. Possible explanations for our results include sample differences in background linkage disequilibrium and/or risk allele effect size, the presence of multiple risk alleles upon different haplotypes, or the presence of a single risk allele upon multiple haplotypes. Some previous associations may also represent false positives. Examination of Caucasian HapMap phase II genotype data spanning the DTNBP1 region indicates upwards of 40 SNPs are required to satisfactorily assess all nonredundant variation within DTNBP1 and its potential regulatory regions for association with schizophrenia. More comprehensive studies in multiple samples will be required to determine whether specific DTNBP1 variants function as risk factors for schizophrenia.

Separate and interacting effects within the catechol-O-methyltransferase (COMT) are associated with schizophrenia

Several lines of evidence have implicated the catechol-O-methyltransferase (COMT) gene as a candidate for schizophrenia (SZ) susceptibility, not only because it encodes a key dopamine catabolic enzyme but also because it maps to the velocardiofacial syndrome region of chromosome 22q11 which has long been associated with SZ predisposition. The interest in COMT as a candidate SZ risk factor has led to numerous case-control and family-based studies, with the majority placing emphasis on examining a functional Val/Met polymorphism within this enzyme. Unfortunately, these studies have continually produced conflicting results. To assess the genetic contribution of other COMT variants to SZ susceptibility, we investigated three single-nucleotide polymorphisms (SNPs) (rs737865, rs4633, rs165599) in addition to the Val/Met variant (rs4680) in a highly selected sample of Australian Caucasian families containing 107 patients with SZ. The Val/Met and rs4633 variants showed nominally significant associations with SZ (P<0.05), although neither of the individual SNPs remained significant after adjusting for multiple testing (most significant P=0.1174). However, haplotype analyses showed strong evidence of an association; the most significant being the three-marker haplotype rs737865-rs4680-rs165599 (global P=0.0022), which spans more than 26 kb. Importantly, conditional analyses indicated the presence of two separate and interacting effects within this haplotype, irrespective of gender. In addition, our results indicate the Val/Met polymorphism is not disease-causing and is simply in strong linkage disequilibrium with a causative effect, which interacts with another as yet unidentified variant approximately 20 kb away. These results may help explain the inconsistent results reported on the Val/Met polymorphism and have important implications for future investigations into the role of COMT in SZ susceptibility.

Further evidence for linkage of Gilles de la Tourette syndrome (GTS) susceptibility loci on chromosomes 2p11, 8q22 and 11q23-24 in South African Afrikaners

Utilizing DNA samples from 91 Afrikaner nuclear families with one or more affected children, five genomic regions on chromosomes 2p, 8q, 11q, 20q, and 21q that gave evidence for association with GTS in previous case-control association studies were investigated for linkage and association with GTS. Highly polymorphic markers with mean heterozygosity of 0.77 were typed and resulting genotypes evaluated using single marker transmission disequilibrium (TDT), single marker haplotype relative risk (HRR), and multi-marker "extended" TDT and HRR methods. Single marker TDT analysis showed evidence for linkage or association, with p-values near 0.05, for markers D2S139, GATA28F12, and D11S1377 on chromosomes 2p11, 8q22 and 11q23-24, respectively. Extended, two-locus TDT and HRR analysis provided further evidence for linkage or association on chromosome 2 with p-values of 0.007 and 0.025, and chromosome 8 with p-values of 0.059 and 0.013, respectively. These results provide important additional evidence for the location of GTS susceptibility loci.

A genomewide screen for autism susceptibility loci

We report the analysis of 335 microsatellite markers genotyped in 110 multiplex families with autism. All families include at least two "affected" siblings, at least one of whom has autism; the remaining affected sibs carry diagnoses of either Asperger syndrome or pervasive developmental disorder. Affected sib-pair analysis yielded multipoint maximum LOD scores (MLS) that reach the accepted threshold for suggestive linkage on chromosomes 5, X, and 19. Nominal evidence for linkage (point-wise P<.05) was obtained on chromosomes 2, 3, 4, 8, 10, 11, 12, 15, 16, 18, and 20, and secondary loci were found on chromosomes 5 and 19. Analysis of families sharing alleles at the putative X chromosomal linked locus and one or more other putative linked loci produced an MLS of 3.56 for the DXS470-D19S174 marker combination. In an effort to increase power to detect linkage, scan statistics were used to evaluate the significance of peak LOD scores based on statistical evidence at adjacent marker loci. This analysis yielded impressive evidence for linkage to autism and autism-spectrum disorders with significant genomewide P values <.05 for markers on chromosomes 5 and 8 and with suggestive linkage evidence for a marker on chromosome 19.

Exclusion of angiotensinogen gene in molecular basis of human hypertension: Sibpair linkage and association analyses in Australian anglo-caucasians

Linkage with essential hypertension has been claimed for a microsatellite marker near the angiotensinogen gene (AGT; chromosome 1q42), as has association for the AGT variants M235T, G(-6)A and A(-20)C. To more rigorously evaluate AGT as a candidate gene for hypertension we performed sibpair analysis with multiple microsatellite markers surrounding this locus and using more sophisticated analysis programs. We also performed an association study of the AGT variants in unrelated subjects with a strong family history (two affected parents). For the linkage study, single and multiplex polymerase chain reaction (PCRs) and automated genescan analysis were conducted on DNA from 175 Australian Anglo-Celtic Caucasian hypertensives for the following markers: D1S2880-(2.1 cM)-D1S213-(2.8 cM)-D1S251-(6.5 cM)-AGT-(2.0 cM) -D1S235. Statistical evaluation of genotype data by nonparametric methods resulted in the following scores: Single-point analysis - SPLINK, P > 0.18; APM method, P > 0.25; ASPEX, MLOD < 0.28; SIB-PAIR, P > 0. 24; Multipoint analysis - MAPMAKER/SIBS, MLOD < 0.24; GENEHUNTER, P > 0.35. Exclusion scores of Lod -4.1 to -5.1 were obtained for these markers using MAPMAKER/SIBS for a lambda(s) of 1.6. The association study of G(-6)A, A(-20)C and M235T variants in 111 hypertensives with strong family history and 190 normotensives with no family history showed significant linkage disequilibrium between particular haplotypes, but we could find no association with hypertension. The present study therefore excludes AGT in the etiology of hypertension, at least in the population of Australian Anglo-Celtic Caucasians studied.

Corneal confocal microscopy detects neuropathy in patients with type 1 diabetes without retinopathy or microalbuminuria

Objective Corneal innervation is increasingly used as a surrogate marker of human diabetic peripheral neuropathy (DPN) however its temporal relationship with the other microvascular complications of diabetes is not fully established. In this cross-sectional, observational study we aimed to assess whether neuropathy occurred in patients with type 1 diabetes, without retinopathy or microalbuminuria. Materials and Methods All participants underwent detailed assessment of peripheral neuropathy [neuropathy disability score (NDS), vibration perception threshold (VPT), peroneal motor nerve conduction velocity (PMNCV), sural sensory nerve conduction velocity (SSNCV) and in vivo corneal confocal microscopy (IVCCM)], retinopathy (digital fundus photography) and albuminuria status [albumin: creatinine ratio (ACR)]. Results 53 patients with Type 1 diabetes with (n=37) and without retinopathy (n=16) were compared to control subjects (n=27). SSNCV, corneal nerve fibre (CNFD) and branch (CNBD) density and length (CNFL) were reduced significantly (p<0.001) in diabetic patients without retinopathy compared to control subjects. Furthermore, CNFD, CNBD and CNFL were also significantly (p<0.001) reduced in diabetic patients without microalbuminuria (n=39), compared to control subjects. Greater neuropathic severity was associated with established retinopathy and microalbuminuria. Conclusions IVCCM detects early small fibre damage in the absence of retinopathy or microalbuminuria in patients with Type 1 diabetes.

MPC controlled multirotor with suspended slung load: System architecture and visual load detection

There is an increased interest in the use of Unmanned Aerial Vehicles for load transportation from environmental remote sensing to construction and parcel delivery. One of the main challenges is accurate control of the load position and trajectory. This paper presents an assessment of real flight trials for the control of an autonomous multi-rotor with a suspended slung load using only visual feedback to determine the load position. This method uses an onboard camera to take advantage of a common visual marker detection algorithm to robustly detect the load location. The load position is calculated using an onboard processor, and transmitted over a wireless network to a ground station integrating MATLAB/SIMULINK and Robotic Operating System (ROS) and a Model Predictive Controller (MPC) to control both the load and the UAV. To evaluate the system performance, the position of the load determined by the visual detection system in real flight is compared with data received by a motion tracking system. The multi-rotor position tracking performance is also analyzed by conducting flight trials using perfect load position data and data obtained only from the visual system. Results show very accurate estimation of the load position (~5% Offset) using only the visual system and demonstrate that the need for an external motion tracking system is not needed for this task.

The relationship between forward head posture and cervical muscle performance in healthy individuals

Background Forward head postures (FHP) are proposed to adversely load cervical spine structures. Neck muscles provide support for the neck, and thus an imbalance in neck muscle performance could potentially contribute to the development of FHP. Previous studies have not considered the interaction of multiple muscle groups with regard to postural orientation. Given the interdependence of muscles along the cervical spine for optimal orientation and physical support of the vertebral column, the performance of a single muscle group may not accurately reflect the coordinated ability of the muscles to maintain a neutral neck posture. Purpose The purpose of this study was to investigate the relationship between FHP and the balance between the cervical extensor and flexor muscle groups in healthy individuals. We hypothesised that the magnitude of FHP would be associated with the strength and endurance performance ratios between the cervical extensor and flexor muscle groups. Methods Twenty male and 24 female volunteers were photographed in the sagittal plane wearing surface markers. The FHP of each participant was measured via the tragus-sternum marker distance over two conditions: (1)in relaxed standing and (2)during a sustained sitting task. Maximal strength (Nm) and endurance (s) performance of the extensor and flexor muscle groups were recorded at the upper (craniocervical flexion/extension (CCF/CCE)) and lower (cervicothoracic flexion/extension (CTF/CTE)) cervical regions. Muscle performance measures were expressed as extension:flexion ratios and their relation to FHP evaluated. A stepwise multiple regression analysis using backward elimination was utilised to examine the relationship between the postural measures and the muscle performance ratio measures. Separate models were used for the two different postural conditions (standing, sustained sitting). Gender was included as a constant correction factor in all regression models. Where gender was a significant variable in the model, analyses were repeated separately for males and females. Results Greater FHP in standing was significantly associated with reduced proportional CTE to CCF strength in females (R2 = 0.21, P = 0.03) and greater proportional CTE to CTF strength in males (R2 = 0.23, P = 0.03). A greater drift into FHP during sustained sitting was associated with a relative reduction in CCE endurance proportional to CTF endurance in females only (R2 = 0.27, P = 0.017). Conclusion(s) This initial study indicates that the balance in performance between the cervical flexor and extensor muscle groups may impact FHP in healthy individuals. However, the findings were inconsistent across different muscle performance ratios and gender. Larger scale studies are therefore now needed to further clarify the relationship between FHP and muscle performance. Implications The findings suggest that relative performance of the various cervical muscle groups needs to be accounted for when considering postural correction strategies in the clinical setting, as is often recommended.

Development of type I genetic markers from expressed sequence tags in highly polymorphic species

Expressed sequence tag (EST) databases provide a primary source of nuclear DNA sequences for genetic marker development in non-model organisms. To date, the process has been relatively inefficient for several reasons: - 1) priming site polymorphism in the template leads to inferior or erratic amplification; - 2) introns in the target amplicon are too large and/or numerous to allow effective amplification under standard screening conditions, and; - 3) at least occasionally, a PCR primer straddles an exon–intron junction and is unable to bind to genomic DNA template. The first is only a minor issue for species or strains with low heterozygosity but becomes a significant problem for species with high genomic variation, such as marine organisms with extremely large effective population sizes. Problems arising from unanticipated introns are unavoidable but are most pronounced in intron-rich species, such as vertebrates and lophotrochozoans. We present an approach to marker development in the Pacific oyster Crassostrea gigas, a highly polymorphic and intron-rich species, which minimizes these problems, and should be applicable to other non-model species for which EST databases are available. Placement of PCR primers in the 3′ end of coding sequence and 3′ UTR improved PCR success rate from 51% to 97%. Almost all (37 of 39) markers developed for the Pacific oyster were polymorphic in a small test panel of wild and domesticated oysters.

Diagnostic values for skin temperature assessment to detect diabetes-related foot complications

Background Skin temperature assessment is a promising modality for early detection of diabetic foot problems, but its diagnostic value has not been studied. Our aims were to investigate the diagnostic value of different cutoff skin temperature values for detecting diabetes-related foot complications such as ulceration, infection, and Charcot foot and to determine urgency of treatment in case of diagnosed infection or a red-hot swollen foot. Materials and Methods The plantar foot surfaces of 54 patients with diabetes visiting the outpatient foot clinic were imaged with an infrared camera. Nine patients had complications requiring immediate treatment, 25 patients had complications requiring non-immediate treatment, and 20 patients had no complications requiring treatment. Average pixel temperature was calculated for six predefined spots and for the whole foot. We calculated the area under the receiver operating characteristic curve for different cutoff skin temperature values using clinical assessment as reference and defined the sensitivity and specificity for the most optimal cutoff temperature value. Mean temperature difference between feet was analyzed using the Kruskal–Wallis tests. Results The most optimal cutoff skin temperature value for detection of diabetes-related foot complications was a 2.2°C difference between contralateral spots (sensitivity, 76%; specificity, 40%). The most optimal cutoff skin temperature value for determining urgency of treatment was a 1.35°C difference between the mean temperature of the left and right foot (sensitivity, 89%; specificity, 78%). Conclusions Detection of diabetes-related foot complications based on local skin temperature assessment is hindered by low diagnostic values. Mean temperature difference between two feet may be an adequate marker for determining urgency of treatment.

PolyPatEx : an R package for paternity exclusion in autopolyploids

Microsatellite markers have demonstrated their value for performing paternity exclusion and hence exploring mating patterns in plants and animals. Methodology is well established for diploid species, and several software packages exist for elucidating paternity in diploids; however, these issues are not so readily addressed in polyploids due to the increased complexity of the exclusion problem and a lack of available software. We introduce polypatex, an r package for paternity exclusion analysis using microsatellite data in autopolyploid, monoecious or dioecious/bisexual species with a ploidy of 4n, 6n or 8n. Given marker data for a set of offspring, their mothers and a set of candidate fathers, polypatex uses allele matching to exclude candidates whose marker alleles are incompatible with the alleles in each offspring–mother pair. polypatex can analyse marker data sets in which allele copy numbers are known (genotype data) or unknown (allelic phenotype data) – for data sets in which allele copy numbers are unknown, comparisons are made taking into account all possible genotypes that could arise from the compared allele sets. polypatex is a software tool that provides population geneticists with the ability to investigate the mating patterns of autopolyploids using paternity exclusion analysis on data from codominant markers having multiple alleles per locus.

Acute fluid shifts influence the assessment of serum vitamin D status in critically ill patients

Introduction Recent reports have highlighted the prevalence of vitamin D deficiency and suggested an association with excess mortality in critically ill patients. Serum vitamin D concentrations in these studies were measured following resuscitation. It is unclear whether aggressive fluid resuscitation independently influences serum vitamin D. Methods Nineteen patients undergoing cardiopulmonary bypass were studied. Serum 25(OH)D3, 1α,25(OH)2D3, parathyroid hormone, C-reactive protein (CRP), and ionised calcium were measured at five defined timepoints: T1 - baseline, T2 - 5 minutes after onset of cardiopulmonary bypass (CPB) (time of maximal fluid effect), T3 - on return to the intensive care unit, T4 - 24 hrs after surgery and T5 - 5 days after surgery. Linear mixed models were used to compare measures at T2-T5 with baseline measures. Results Acute fluid loading resulted in a 35% reduction in 25(OH)D3 (59 ± 16 to 38 ± 14 nmol/L, P < 0.0001) and a 45% reduction in 1α,25(OH)2D3 (99 ± 40 to 54 ± 22 pmol/L P < 0.0001) and i(Ca) (P < 0.01), with elevation in parathyroid hormone (P < 0.0001). Serum 25(OH)D3 returned to baseline only at T5 while 1α,25(OH)2D3 demonstrated an overshoot above baseline at T5 (P < 0.0001). There was a delayed rise in CRP at T4 and T5; this was not associated with a reduction in vitamin D levels at these time points. Conclusions Hemodilution significantly lowers serum 25(OH)D3 and 1α,25(OH)2D3, which may take up to 24 hours to resolve. Moreover, delayed overshoot of 1α,25(OH)2D3 needs consideration. We urge caution in interpreting serum vitamin D in critically ill patients in the context of major resuscitation, and would advocate repeating the measurement once the effects of the resuscitation have abated.