GOLD--graphical overview of linkage disequilibrium

Joint Linkage and Association Analysis Using GENEHUNTER-MODSCORE with an Application to Familial Pancreatic Cancer

INTRODUCTION

Joint linkage and association (JLA) analysis combines two disease gene mapping strategies: linkage information contained in families and association information contained in populations. Such a JLA analysis can increase mapping power, especially when the evidence for both linkage and association is low to moderate. Similarly, an association analysis based on haplotypes instead of single markers can increase mapping power when the association pattern is complex.

METHODS

In this paper, we present an extension to the GENEHUNTER-MODSCORE software package that enables a JLA analysis based on haplotypes and uses information from arbitrary pedigree types and unrelated individuals. Our new JLA method is an extension of the MOD score approach for linkage analysis, which allows the estimation of trait-model and linkage disequilibrium (LD) parameters, i.e., penetrance, disease-allele frequency, and haplotype frequencies. LD is modeled between alleles at a single diallelic disease locus and up to three diallelic test markers. Linkage information is contributed by additional multi-allelic flanking markers. We investigated the statistical properties of our JLA implementation using extensive simulations, and we compared our approach to another commonly used single-marker JLA test. To demonstrate the applicability of our new method in practice, we analyzed pedigree data from the German National Case Collection for Familial Pancreatic Cancer (FaPaCa).

RESULTS

Based on the simulated data, we demonstrated the validity of our JLA-MOD score analysis implementation and identified scenarios in which haplotype-based tests outperformed the single-marker test. The estimated trait-model and LD parameters were in good accordance with the simulated values. Our method outperformed another commonly used JLA single-marker test when the LD pattern was complex. The exploratory analysis of the FaPaCa families led to the identification of a promising genetic region on chromosome 22q13.33, which can serve as a starting point for future mutation analysis and molecular research in pancreatic cancer.

CONCLUSION

Our newly proposed JLA-MOD score method proves to be a valuable gene mapping and characterization tool, especially when either linkage or association information alone provide insufficient power to identify the disease-causing genetic variants.

Cluster of Differentiation 24 Polymorphism Has No Significant Association with Chronic Hepatitis B Virus Infection in the Chinese Han Population: A Family-Based Association Study

Objective

Studies have shown that cluster of differentiation (CD) 24 gene polymorphism is associated with several diseases. Among these, chronic hepatitis B (CHB) infection has not been investigated. This study aimed to assess the function of CD24 in CHB.

Methods

The study included 478 cases of CHB and 318 cases without CHB from 230 families that underwent genotyping. Polymerase chain reaction-restriction fragment length polymorphism was performed to assess the single nucleotide polymorphism (SNP) P170 of the CD24 gene. The detected genotypes were TT, CT, and CC. Then, family based-association analysis was carried out to investigate the association between CD24 gene polymorphism and susceptibility to CHB.

Results

In the 478 patients with CHB, the frequencies of CD24 P170 T and C alleles were 35.5% and 64.5%, respectively, and the frequencies of CD24 P170 CC, CT, and TT genotypes were 39.3%, 50.4% and 10.3%, respectively. In a CD24 single-locus analysis by a family-based association test of P170 polymorphisms, T and C were not significantly associated with CHB in the additive (Z = 0.169, P = 0.866; Z = −0.169, P = 0.866, respectively), dominant (Z = 0.522, P = 0.602; Z = 0.428, P = 0.669, respectively), or recessive (Z = −0.428, P = 0.669; Z = −0.522, P = 0.602, respectively) models. Transmission-disequilibrium (TD) and sib-transmission disequilibrium (STD) tests revealed no excess of T or C alleles from heterozygous parents to their children with the disease or higher frequencies of these alleles in patients compared with their normal siblings (χ² = 0.06, P = 0.897).

Conclusion

The study findings suggest that the SNP P170 of CD24 has no significant association with susceptibility to the HB virus and related phenotypes in Chinese patients.

Associations between Omega-3 Index, Dopaminergic Genetic Variants and Aggressive and Metacognitive Traits: A Study in Adult Male Prisoners

Omega-3 long-chain polyunsaturated fatty acids (n-3 LCPUFA) are critical for cell membrane structure and function. Human beings have a limited ability to synthesise docosahexaenoic acid (DHA), the main n-3 LCPUFA required for neurological development. Inadequate levels of n-3 LCPUFA can affect the dopaminergic system in the brain and, when combined with genetic and other factors, increase the risk of developing aggression, inattention and impulse-control disorders. In this study, male prisoners were administered questionnaires assessing aggressive behaviour and executive functions. Participants also produced blood sampling for the measurement of the Omega-3 Index and the genotyping of dopaminergic genetic variants. Significant associations were found between functional genetic polymorphism in DBH rs1611115 and verbal aggression and between DRD2 rs4274224 and executive functions. However, the Omega-3 Index was not significantly associated with the tested dopaminergic polymorphisms. Although previous interactions between specific genotypes and n-3 LCPUFA were previously reported, they remain limited and poorly understood. We did not find any association between n-3 LCPUFA and dopaminergic polymorphisms in adult male prisoners; however, we confirmed the importance of genetic predisposition for dopaminergic genes (DBH and DRD2) in aggressive behaviour, memory dysfunction and attention-deficit disorder.

ARID5B rs10821936 and rs10994982 gene polymorphism and susceptibility to juvenile systemic lupus erythematosus and lupus nephritis

BACKGROUND

The prevalence of SLE and the spectrum of clinical manifestations vary widely in different races and geographical populations.

OBJECTIVE

To investigate the possible role of ARID5B rs10821936 and rs10994982 polymorphism as a risk factor for the development of SLE in children (jSLE) and to evaluate their role in relation to clinical manifestations especially lupus nephritis (LN).

METHODS

DNA extraction and Real-time PCR genotyping of ARID5B rs10821936 and rs10994982 were done for 104 jSLE and 282 healthy controls.

RESULTS

The C allele and C containing genotypes (CC, CT and CC+CT) of ARID5B rs10821936 were higher in children with SLE (p = 0.009, OR = 1.56, 0.037, OR = 2.35, 0.016, OR = 1.81 and 0.008 OR = 1.88 respectively). ARID5B rs10994982 alleles, genotypes and haplotypes are not associated with jSLE (p > 0.05). The ARID5B rs10821936 and rs10994982 genotypes showed non-significant associations with LN, proliferative versus non proliferative and biopsy grades (p > 0.05).

CONCLUSION

ARID5B rs10821936 SNP may be a susceptibility risk factor for juvenile SLE in the studied cohort of Egyptian children.

Notes on Three Decades of Methodology Workshops

Since 1987, a group of behavior geneticists have been teaching an annual methodology workshop on how to use state-of-the-art methods to analyze genetically informative data. In the early years, the focus was on analyzing twin and family data, using information of their known genetic relatedness to infer the role of genetic and environmental factors on phenotypic variation. With the rapid evolution of genotyping and sequencing technology and availability of measured genetic data, new methods to detect genetic variants associated with human traits were developed and became the focus of workshop teaching in alternate years. Over the years, many of the methodological advances in the field of statistical genetics have been direct outgrowths of the workshop, as evidence by the software and methodological publications authored by workshop faculty. We provide data and demographics of workshop attendees and evaluate the impact of the methodology workshops on scientific output in the field by evaluating the number of papers applying specific statistical genetic methodologies authored by individuals who have attended workshops.

LDscaff: LD-based scaffolding of de novo genome assemblies

Background

Genome assembly is fundamental for de novo genome analysis. Hybrid assembly, utilizing various sequencing technologies increases both contiguity and accuracy. While such approaches require extra costly sequencing efforts, the information provided millions of existed whole-genome sequencing data have not been fully utilized to resolve the task of scaffolding. Genetic recombination patterns in population data indicate non-random association among alleles at different loci, can provide physical distance signals to guide scaffolding.

Results

In this paper, we propose LDscaff for draft genome assembly incorporating linkage disequilibrium information in population data. We evaluated the performance of our method with both simulated data and real data. We simulated scaffolds by splitting the pig reference genome and reassembled them. Gaps between scaffolds were introduced ranging from 0 to 100 KB. The genome misassembly rate is 2.43% when there is no gap. Then we implemented our method to refine the Giant Panda genome and the donkey genome, which are purely assembled by NGS data. After LDscaff treatment, the resulting Panda assembly has scaffold N50 of 3.6 MB, 2.5 times larger than the original N50 (1.3 MB). The re-assembled donkey assembly has an improved N50 length of 32.1 MB from 23.8 MB.

Conclusions

Our method effectively improves the assemblies with existed re-sequencing data, and is an potential alternative to the existing assemblers required for the collection of new data.

Update on the genetics of primary open-angle glaucoma

Affecting nearly 80 million individuals, glaucoma is the number one cause of irreversible blindness in the world. This ocular disease describes a set of optic neuropathies of which primary open angle glaucoma (POAG) is the most common. POAG is associated with progressive visual field deterioration resulting from damage to the optic nerve and loss of retinal ganglion cells. Risk factors for POAG include elevated intraocular pressure, aging, African and Hispanic ancestry, and a positive family history of POAG. Multiple genes have been found to contribute to POAG. Much of POAG genetics and pathology has yet to be explained. Recent genome-wide association studies have identified a large number of novel loci associated with POAG and its endophenotypes. Genomic and proteomic profiling of biofluids has contributed to our knowledge of differential gene expression in POAG. Functional studies both in cell culture and animal models have confirmed the effects of variants and differential gene expression on ocular physiology while in silico analyses have increased our understanding of disease risk and progression so that we might better diagnose and treat this complex genetic illness.

Genomic data provide new insights on the demographic history and the extent of recent material transfers in Norway spruce

Primeval forests are today exceedingly rare in Europe, and transfer of forest reproductive material for afforestation and improvement has been very common, especially over the last two centuries. This can be a serious impediment when inferring past population movements in response to past climate changes such as the last glacial maximum (LGM), some 18,000 years ago. In the present study, we genotyped 1,672 individuals from three Picea species (P. abies, P. obovata, and P. omorika) at 400K SNPs using exome capture to infer the past demographic history of Norway spruce (P. abies) and estimate the amount of recent introduction used to establish the Norway spruce breeding program in southern Sweden. Most of these trees belong to P. abies and originate from the base populations of the Swedish breeding program. Others originate from populations across the natural ranges of the three species. Of the 1,499 individuals stemming from the breeding program, a large proportion corresponds to recent introductions from mainland Europe. The split of P. omorika occurred 23 million years ago (mya), while the divergence between P. obovata and P. abies began 17.6 mya. Demographic inferences retrieved the same main clusters within P. abies than previous studies, that is, a vast northern domain ranging from Norway to central Russia, where the species is progressively replaced by Siberian spruce (P. obovata) and two smaller domains, an Alpine domain and a Carpathian one, but also revealed further subdivision and gene flow among clusters. The three main domains divergence was ancient (15 mya), and all three went through a bottleneck corresponding to the LGM. Approximately 17% of P. abies Nordic domain migrated from P. obovata ~103K years ago, when both species had much larger effective population sizes. Our analysis of genomewide polymorphism data thus revealed the complex demographic history of Picea genus in Western Europe and highlighted the importance of material transfer in Swedish breeding program.

A comparative integrated gene-based linkage and locus ordering by linkage disequilibrium map for the Pacific white shrimp, Litopenaeus vannamei

The Pacific whiteleg shrimp, Litopenaeus vannamei, is the most farmed aquaculture species worldwide with global production exceeding 3 million tonnes annually. Litopenaeus vannamei has been the focus of many selective breeding programs aiming to improve growth and disease resistance. However, these have been based primarily on phenotypic measurements and omit potential gains by integrating genetic selection into existing breeding programs. Such integration of genetic information has been hindered by the limited available genomic resources, background genetic parameters and knowledge on the genetic architecture of commercial traits for L. vannamei. This study describes the development of a comprehensive set of genomic gene-based resources including the identification and validation of 234,452 putative single nucleotide polymorphisms in-silico, of which 8,967 high value SNPs were incorporated into a commercially available Illumina Infinium ShrimpLD-24 v1.0 genotyping array. A framework genetic linkage map was constructed and combined with locus ordering by disequilibrium methodology to generate an integrated genetic map containing 4,817 SNPs, which spanned a total of 4552.5 cM and covered an estimated 98.12% of the genome. These gene-based genomic resources will not only be valuable for identifying regions underlying important L. vannamei traits, but also as a foundational resource in comparative and genome assembly activities.

CYP2A6 Polymorphisms Associate with Outcomes of S-1 Plus Oxaliplatin Chemotherapy in Chinese Gastric Cancer Patients

Gastric carcinoma is a heterogeneous malignant disease involving genetic factors. To identify predictive markers for gastric cancer treatment in Chinese patients, we evaluated the association between polymorphisms of the gene encoding cytochrome P450 2A6 (CYP2A6) and outcomes of S-1 plus oxaliplatin (SOX) chemotherapy treatment. Clinical data on 60 consecutive gastric cancer patients receiving SOX regimen were collected prospectively. We sequenced all exons of CYP2A6 and a total of 22 different polymorphisms were detected in the present study. Comprehensive analyses of these genetic polymorphisms were performed to determine their association with both safety and efficacy of SOX regimen. Our results showed that polymorphisms of CYP2A6 were associated with the safety and efficacy of SOX treatment. Among them, missense mutations CYP2A6 rs60823196 and rs138978736 could be possible risk factors (P < 0.05) for severe diarrhea induced by SOX, whereas CYP2A6 rs138978736 could be a conceivable predictor for overall survival of patients treated with SOX adjuvant chemotherapy. Further large-scale randomized prospective studies are warranted to confirm these findings.

Estimating Disequilibrium Coefficients

Gametic phase disequilibrium is the nonrandom association of alleles within gametes. Linkage disequilibrium (LD) describes the special case of deviation from independence between alleles at two linked genetic loci. Estimation of allelic LD requires knowledge of haplotypes. Genotype-based LD measures dispense with the haplotype estimation step and avoid bias in LD estimation. In this chapter, the most important measures for allelic and genotypic LD are introduced. The use of software packages for LD estimation is illustrated.

UGT1A1*6, UGT1A7*3 and UGT1A9*1b polymorphisms are predictive markers for severe toxicity in patients with metastatic gastrointestinal cancer treated with irinotecan-based regimens

Irinotecan-induced severe neutropenia and diarrhea, which remain unpredictable, has restrained the dose and clinical efficiency of irinotecan administration. In the present study, a total of 70 irinotecan-treated patients with histologically confirmed metastatic gastrointestinal cancer were enrolled. Despite genotyping well-reported alleles, direct sequencing was specifically adopted to avoid ethnic heterogeneity and to identify novel variations. The promoter (-1000 bp) and exon 1 regions of UDP glucuronosyltransferase family 1 member A complex locus (UGT1A1) gene family members UGT1A1, UGT1A7 and UGT1A9 were sequenced, and comprehensive analysis of their genetic polymorphisms was performed to determine the association between inherited genetic variations and irinotecan-induced toxicity. A total of 23 different genetic variants were detected in the present study, including 2 novel polymorphisms. The results of the present study revealed that UGT1A1*6 and UGT1A7*3 are risk factors for irinotecan-induced severe neutropenia, and UGT1A9*1b is associated with severe diarrhea. These results may provide biomarkers for the selection of the optimal chemotherapy for Chinese patients with metastatic gastrointestinal cancer.

Effective Population Size and Signatures of Selection Using Bovine 50K SNP Chips in Korean Native Cattle (Hanwoo)

Inferring the effective population size and the pattern of selection signatures is of interest both from an evolutionary perspective and to improve models for mapping of quantitative trait genes. We used DNA samples of 61 sires and 486 progeny of the Hanwoo, genotyped by the Illumina Bovine SNP50 BeadChip, to analyze the genetic structure. Our study showed a persistent decline in effective population size throughout the period considered, but suggested a marked decline at one distinctive time point (100th generation) and two sharp decline intervals (50th-25th generation and 25th-10th generation). This pattern can be explained by Hanwoo formation and the modern breeding program. Our results revealed 95 regions exhibiting the footprint of recent positive selection at a threshold level of 0.01. We found an overlap of the 11 core regions presenting top P-values and those that had previously been identified as harboring quantitative trait loci from other breeds. The information generated from this study can be used to better understand the mechanism of selection in Hanwoo breeding, and provide important implications for the design and application of association studies in the Hanwoo population.

The Impact of Genetic Relationship and Linkage Disequilibrium on Genomic Selection

Genomic selection is a promising research area due to its practical application in breeding. In this study, impact of realized genetic relationship and linkage disequilibrium (LD) on marker density and training population size required was investigated and their impact on practical application was further discussed. This study is based on experimental data of two populations derived from the same two founder lines (B73, Mo17). Two populations were genotyped with different marker sets at different density: IBM Syn4 and IBM Syn10. A high-density marker set in Syn10 was imputed into the Syn4 population with low marker density. Seven different prediction scenarios were carried out with a random regression best linear unbiased prediction (RR-BLUP) model. The result showed that the closer the real genetic relationship between training and validation population, the fewer markers were required to reach a good prediction accuracy. Taken the short-term cost for consideration, relationship information is more valuable than LD information. Meanwhile, the result indicated that accuracies based on high LD between QTL and markers were more stable over generations, thus LD information would provide more robust prediction capacity in practical applications.

Systematic Validation of RNF213 Coding Variants in Japanese Patients With Moyamoya Disease

Background

A founder variant of RNF213, p.R4810K (c.14429G>A, rs112735431), was recently identified as a major genetic risk factor for moyamoya disease (MMD) in Japan. Although the association of p.R4810K was reported to be highly significant and reproducible, the disease susceptibility of other RNF213 variants remains largely unknown. In the present study, we systematically evaluated the coding variants detected in Japanese patients and controls for associations with MMD.

Methods and Results

To detect variants of RNF213, all coding exons were sequenced in 27 Japanese MMD patients without p.R4810K. We also validated all previously reported variants in our case–control samples and tested for associations in combination with previous Japanese study cohorts, including the 1000 Genomes Project data set, as population-based controls. Forty-six missense variants other than p.R4810K were identified among 370 combined patients and 279 combined controls in Japan. Sixteen of 46 variants were polymorphisms with minor allele frequency >1%, and, after conditioning on the p.R4810K genotype, were not associated with MMD. We conducted a variable threshold test using Combined Annotation-Dependent Depletion on the remaining 30 rare variants (minor allele frequency <1%), and the results showed that the frequency of potentially functional variants was significantly higher in patients than in controls (permutation, minimum P=0.045).

Conclusions

Not only p.4810K but also other functional missense variants of RNF213 conferred susceptibility to MMD. Our analysis also revealed that ≈20% of Japanese MMD patients did not harbor susceptibility variants of RNF213, indicating the presence of other susceptibility genes for MMD.

Genetic mapping of a new race specific resistance allele effective to Puccinia hordei at the Rph9/Rph12 locus on chromosome 5HL in barley

BACKGROUND

Barley is an important cereal crop cultivated for malt and ruminant feed and in certain regions it is used for human consumption. It is vulnerable to numerous foliar diseases including barley leaf rust caused by the pathogen Puccinia hordei.

RESULTS

A temporarily designated resistance locus RphCantala (RphC) identified in the Australian Hordeum vulgare L. cultivar 'Cantala' displayed an intermediate to low infection type (";12 = N") against the P. hordei pathotype 253P- (virulent on Rph1, Rph2, Rph4, Rph6, Rph8 and RphQ). Phenotypic assessment of a 'CI 9214' (susceptible) x 'Stirling' (RphC) (CI 9214/Stirling) doubled haploid (DH) population at the seedling stage using P. hordei pathotype 253P-, confirmed that RphC was monogenically inherited. Marker-trait association analysis of RphC in the CI 9214/Stirling DH population using 4,500 DArT-seq markers identified a highly significant (-log10Pvalue > 17) single peak on the long arm of chromosome 5H (5HL). Further tests of allelism determined that RphC was genetically independent of Rph3, Rph7, Rph11, Rph13 and Rph14, and was an allele of Rph12 (Rph9.z), which also maps to 5HL.

CONCLUSION

Multipathotype tests and subsequent pedigree analysis determined that 14 related Australian barley varieties (including 'Stirling' and 'Cantala') carry RphC and that the likely source of this resistance is via a Czechoslovakian landrace LV-Kvasice-NA-Morave transferred through common ancestral cultivars 'Hanna' and 'Abed Binder'. RphC is an allele of Rph12 (Rph9.z) and is therefore designated Rph9.am. Bioinformatic analysis using sequence arrays from DArT-seq markers in linkage disequilibrium with Rph9.am identified possible candidates for further gene cloning efforts and marker development at the Rph9/Rph12/Rph9.am locus.

Role of estrogen related receptor beta (ESRRB) in DFN35B hearing impairment and dental decay

BACKGROUND

Congenital forms of hearing impairment can be caused by mutations in the estrogen related receptor beta (ESRRB) gene. Our initial linkage studies suggested the ESRRB locus is linked to high caries experience in humans.

METHODS

We tested for association between the ESRRB locus and dental caries in 1,731 subjects, if ESRRB was expressed in whole saliva, if ESRRB was associated with the microhardness of the dental enamel, and if ESRRB was expressed during enamel development of mice.

RESULTS

Two families with recessive ESRRB mutations and DFNB35 hearing impairment showed more extensive dental destruction by caries. Expression levels of ESRRB in whole saliva samples showed differences depending on sex and dental caries experience.

CONCLUSIONS

The common etiology of dental caries and hearing impairment provides a venue to assist in the identification of individuals at risk to either condition and provides options for the development of new caries prevention strategies, if the associated ESRRB genetic variants are correlated with efficacy.

Accuracy of predicting genomic breeding values for residual feed intake in Angus and Charolais beef cattle

In beef cattle, phenotypic data that are difficult and/or costly to measure, such as feed efficiency, and DNA marker genotypes are usually available on a small number of animals of different breeds or populations. To achieve a maximal accuracy of genomic prediction using the phenotype and genotype data, strategies for forming a training population to predict genomic breeding values (GEBV) of the selection candidates need to be evaluated. In this study, we examined the accuracy of predicting GEBV for residual feed intake (RFI) based on 522 Angus and 395 Charolais steers genotyped on SNP with the Illumina Bovine SNP50 Beadchip for 3 training population forming strategies: within breed, across breed, and by pooling data from the 2 breeds (i.e., combined). Two other scenarios with the training and validation data split by birth year and by sire family within a breed were also investigated to assess the impact of genetic relationships on the accuracy of genomic prediction. Three statistical methods including the best linear unbiased prediction with the relationship matrix defined based on the pedigree (PBLUP), based on the SNP genotypes (GBLUP), and a Bayesian method (BayesB) were used to predict the GEBV. The results showed that the accuracy of the GEBV prediction was the highest when the prediction was within breed and when the validation population had greater genetic relationships with the training population, with a maximum of 0.58 for Angus and 0.64 for Charolais. The within-breed prediction accuracies dropped to 0.29 and 0.38, respectively, when the validation populations had a minimal pedigree link with the training population. When the training population of a different breed was used to predict the GEBV of the validation population, that is, across-breed genomic prediction, the accuracies were further reduced to 0.10 to 0.22, depending on the prediction method used. Pooling data from the 2 breeds to form the training population resulted in accuracies increased to 0.31 and 0.43, respectively, for the Angus and Charolais validation populations. The results suggested that the genetic relationship of selection candidates with the training population has a greater impact on the accuracy of GEBV using the Illumina Bovine SNP50 Beadchip. Pooling data from different breeds to form the training population will improve the accuracy of across breed genomic prediction for RFI in beef cattle.

Association between TGFB3 and nonsyndromic cleft lip with or without cleft palate in a Chilean population

OBJECTIVE

To assess the possible association between TGFB3 allele variants and nonsyndromic cleft lip with or without cleft palate in a Chilean population.

DESIGN

In our study we used a case-parents trios design. The sample consisted of 150 unrelated trios ascertained through probands affected with nonsyndromic cleft lip with or without cleft palate. Three TGFB3 polymorphisms were analyzed (rs2268626, rs2268625, and rs3917201). An allele/haplotype transmission disequilibrium test was used to evaluate the possible genotype-phenotype association.

RESULTS

An overtransmission from parents to affected progeny was observed for the A allele of rs3917201 (p = .03) and for the rs2268625-rs3917201 A-A haplotype (p = .022). A defect of transmission of rs2268625-rs3917201 G-G haplotype (p = .022) was observed also.

CONCLUSIONS

Allelic and haplotypic associations implicate a possible role of TGFB3 in nonsyndromic cleft lip with or without cleft palate in the Chilean population. Additional studies are needed in order to elucidate the possible mechanisms that can explain the role of TGFB3 genetic variants in the condition.

Ethnic heterogeneity of IRF6 AP-2a binding site promoter SNP association with nonsyndromic cleft lip and palate

OBJECTIVE

The goal of this study was to confirm the reported association between a noncoding SNP (rs642961) in IRF6 and nonsyndromic cleft lip and palate.

DESIGN, SETTING, AND PARTICIPANTS

Two SNPs in IRF6 (rs2235371 and rs64296) were genotyped in Hispanic and non-Hispanic white multiplex (122) and simplex (308) nonsyndromic cleft lip and palate families. Linkage and family-based association analyses were performed on the individual SNPs as well as the 2-SNP haplotype.

RESULTS

Only modest evidence was found for an association with rs642961 and the 2-SNP haplotype. In contrast, strong evidence was found for the association with rs2235371; this was most evident in the non-Hispanic white simplex families.

CONCLUSIONS

Although it is confirmed that variation in IRF6 is associated with nonsyndromic cleft lip and palate, the results do not support the reported association with SNP rs64296. Importantly, the association varies between ethnic groups. This finding underscores the need for evaluating additional variations in IRF6 across multiple populations to better determine its role in nonsyndromic cleft lip and palate.

Genomic dissection of variation in clutch size and egg mass in a wild great tit (Parus major) population

Clutch size and egg mass are life history traits that have been extensively studied in wild bird populations, as life history theory predicts a negative trade-off between them, either at the phenotypic or at the genetic level. Here, we analyse the genomic architecture of these heritable traits in a wild great tit (Parus major) population, using three marker-based approaches - chromosome partitioning, quantitative trait locus (QTL) mapping and a genome-wide association study (GWAS). The variance explained by each great tit chromosome scales with predicted chromosome size, no location in the genome contains genome-wide significant QTL, and no individual SNPs are associated with a large proportion of phenotypic variation, all of which may suggest that variation in both traits is due to many loci of small effect, located across the genome. There is no evidence that any regions of the genome contribute significantly to both traits, which combined with a small, nonsignificant negative genetic covariance between the traits, suggests the absence of genetic constraints on the independent evolution of these traits. Our findings support the hypothesis that variation in life history traits in natural populations is likely to be determined by many loci of small effect spread throughout the genome, which are subject to continued input of variation by mutation and migration, although we cannot exclude the possibility of an additional input of major effect genes influencing either trait.

Population-specific differences in gene conversion patterns between human SUZ12 and SUZ12P are indicative of the dynamic nature of interparalog gene conversion

Nonallelic homologous gene conversion (NAHGC) resulting from interparalog recombination without crossover represents an important influence on the evolution of duplicated sequences in the human genome. In 17q11.2, different paralogous sequences mediate large NF1 deletions by nonallelic homologous recombination with crossover (NAHR). Among these paralogs are SUZ12 and its pseudogene SUZ12P which harbour the breakpoints of type-2 (1.2-Mb) NF1 deletions. Such deletions are caused predominantly by mitotic NAHR since somatic mosaicism with normal cells is evident in most patients. Investigating whether SUZ12 and SUZ12P have also been involved in NAHGC, we observed gene conversion tracts between these paralogs in both Africans (AFR) and Europeans (EUR). Since germline type-2 NF1 deletions resulting from meiotic NAHR are very rare, the vast majority of the gene conversion tracts in SUZ12 and SUZ12P are likely to have resulted from mitotic recombination during premeiotic cell divisions of germ cells. A higher number of gene conversion tracts were noted within SUZ12 and SUZ12P in AFR as compared to EUR. Further, the distinctive signature of NAHGC (a high number of SNPs per paralog and a high number of shared SNPs between paralogs), a characteristic of many actively recombining paralogs, was observed in both SUZ12 and SUZ12P but only in AFR and not in EUR. A novel polymorphic 2.3-kb deletion in SUZ12P was identified which exhibited a high allele frequency in EUR. We postulate that this interparalog structural difference, together with low allelic recombination rates, could have caused a reduction in NAHGC between SUZ12 and SUZ12P during human evolution.

Genetic association analyses of nitric oxide synthase genes and neural tube defects vary by phenotype

Neural tube defects (NTDs) are caused by improper neural tube closure during the early stages of embryonic development. NTDs are hypothesized to have a complex genetic origin and numerous candidate genes have been proposed. The nitric oxide synthase 3 (NOS3) G594T polymorphism has been implicated in risk for spina bifida, and interactions between that single nucleotide polymorphism (SNP) and the methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism have also been observed. To evaluate other genetic variation in the NO pathway in the development of NTDs, we examined all three NOS genes: NOS1, NOS2, and NOS3. Using 3109 Caucasian samples in 745 families, we evaluated association in the overall dataset and within specific phenotypic subsets. Haplotype tagging SNPs in the NOS genes were tested for genetic association with NTD subtypes, both for main effects as well as for the presence of interactions with the MTHFR C677T polymorphism. Nominal main effect associations were found with all subtypes, across all three NOS genes, and interactions were observed between SNPs in all three NOS genes and MTHFR C677T. Unlike the previous report, the most significant associations in our dataset were with cranial subtypes and the AG genotype of rs4795067 in NOS2 (p = 0.0014) and the interaction between the rs9658490 G allele in NOS1 and MTHFR 677TT genotype (p = 0.0014). Our data extend the previous findings by implicating a role for all three NOS genes, independently and through interactions with MTHFR, in risk not only for spina bifida, but all NTD subtypes.

Genetic variants ofLRRK2 in Taiwanese Parkinson's disease

Genetic variants of leucine-rich repeat kinase 2 (LRRK2) were reported to alter the risk for Parkinson's disease (PD). However, the genetic spectrum of LRRK2 variants has not been clearly disclosed yet in Taiwanese population. Herein, we sequenced LRRK2 coding region in 70 Taiwanese early onset PD patients (age at onset ≤ 50), and found six amino acid-changing single nucleotide polymorphisms (SNPs, N551K, R1398H, R1628P, S1647T, G2385R and M2397T), one reported (R1441H) and 2 novel missense (R767H and S885N) mutations. We examined the frequency of identified LRRK2 variants by genotyping 573 Taiwanese patients with PD and 503 age-matched control subjects. The results showed that PD patients demonstrated a higher frequency of G2385R A allele (4.6%) than control subjects (2.1%; odds ratio = 2.27, 95% confidence interval: 1.38-3.88, P = 0.0017). Fewer PD patients (27.7%) carried the 1647T-2397T haplotype as compared with the control subjects (33.0%; odds ratio = 0.80, 95% confidence interval: 0.65-0.97, P = 0.0215). However, the frequency of 1647T-2385R-2397T haplotype (4.3%) in PD patients was still higher than in control subjects (1.9%, odds ratio: 2.15, 95% confidence interval: 1.27-3.78, P = 0.0058). While no additional subject was found to carry R767H and R1441H, one more patient was observed to carry the S885N variant. Our results indicate a robust risk association regarding G2385R and a new possible protective haplotype (1647T-2397T). Gene-environmental interaction and a larger cohort study are warranted to validate our findings. Additionally, two new missense mutations (R767H and S885N) regarding LRRK2 in PD patients were identified. Functional studies are needed to elucidate the effects of these LRRK2 variants on protein function.

A high-density SNP genetic linkage map for the silver-lipped pearl oyster, Pinctada maxima: a valuable resource for gene localisation and marker-assisted selection

Background

The silver-lipped pearl oyster, Pinctada maxima, is an important tropical aquaculture species extensively farmed for the highly sought "South Sea" pearls. Traditional breeding programs have been initiated for this species in order to select for improved pearl quality, but many economic traits under selection are complex, polygenic and confounded with environmental factors, limiting the accuracy of selection. The incorporation of a marker-assisted selection (MAS) breeding approach would greatly benefit pearl breeding programs by allowing the direct selection of genes responsible for pearl quality. However, before MAS can be incorporated, substantial genomic resources such as genetic linkage maps need to be generated. The construction of a high-density genetic linkage map for P. maxima is not only essential for unravelling the genomic architecture of complex pearl quality traits, but also provides indispensable information on the genome structure of pearl oysters.

Results

A total of 1,189 informative genome-wide single nucleotide polymorphisms (SNPs) were incorporated into linkage map construction. The final linkage map consisted of 887 SNPs in 14 linkage groups, spans a total genetic distance of 831.7 centimorgans (cM), and covers an estimated 96% of the P. maxima genome. Assessment of sex-specific recombination across all linkage groups revealed limited overall heterochiasmy between the sexes (i.e. 1.15:1 F/M map length ratio). However, there were pronounced localised differences throughout the linkage groups, whereby male recombination was suppressed near the centromeres compared to female recombination, but inflated towards telomeric regions. Mean values of LD for adjacent SNP pairs suggest that a higher density of markers will be required for powerful genome-wide association studies. Finally, numerous nacre biomineralization genes were localised providing novel positional information for these genes.

Conclusions

This high-density SNP genetic map is the first comprehensive linkage map for any pearl oyster species. It provides an essential genomic tool facilitating studies investigating the genomic architecture of complex trait variation and identifying quantitative trait loci for economically important traits useful in genetic selection programs within the P. maxima pearling industry. Furthermore, this map provides a foundation for further research aiming to improve our understanding of the dynamic process of biomineralization, and pearl oyster evolution and synteny.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-810) contains supplementary material, which is available to authorized users.

Folate metabolism genes, dietary folate and response to antidepressant medications in late-life depression

OBJECTIVE

The primary aims of this study were to (i) determine whether folate metabolism genetic polymorphisms predict age of onset and occurrence of late life depression; and (ii) determine whether folate metabolism genetic polymorphisms predict response to antidepressant medications in late-life depression.

METHODS

This study used the Conte Center for the Neuroscience of Depression and the Neurocognitive Outcomes of Depression in the Elderly Study database, which includes individuals aged ≥60. The folate nutrition assessment was determined by the Block Food Frequency Questionnaire. Genotype was evaluated for 15 single nucleotide polymorphisms from 10 folate metabolism genes. Logistic regression models were used to examine genetic polymorphisms and folate estimates with association with depression age of onset and remission status.

RESULTS

There were 304 Caucasians in the database, 106 of these were not depressed and 198 had a diagnosis of depression. There were no significant differences between remitters and non-remitters in age, sex or estimated folate intakes. There were no folate estimates or folate metabolism gene single nucleotide polymorphisms that significantly predicted age of onset of depression or occurrence of depression. Methionine synthase reductase (MTRR) A66G (rs1801394) was significantly associated with remission status (p = 0.0077) such that those with the AA genotype were 3.2 times as likely as those with the GG genotype to be in remission (p = 0.0020). Methylenetetrahydrofolate reductase A1298C (rs1801131) achieved a borderline significance for association with remission status (p = 0.0313).

CONCLUSION

The major finding from this study is that the MTRR A66G genotype predicts response to selective serotonin reuptake inhibitor antidepressants in late life depression.

Killer cell immunoglobulin-like receptor (KIR) genes and their HLA-C ligands in a Ugandan population

Killer cell immunoglobulin-like receptor (KIR) genes are expressed by natural killer cells and encoded by a family of genes exhibiting considerable haplotypic and allelic variation. HLA-C molecules, the dominant ligands for KIR, are present in all individuals and are discriminated by two KIR epitopes, C1 and C2. We studied the frequencies of KIR genes and HLA-C1 and C2 groups in a large cohort (n = 492) from Kampala, Uganda, East Africa and compared our findings with published data from other populations in sub-Saharan Africa (SSA) and several European populations. We find considerably more KIR diversity and weaker linkage disequilibrium in SSA compared to the European populations and describe several novel KIR genotypes. C1 and C2 frequencies were similar to other SSA populations with a higher frequency of the C2 epitope (54.9 %) compared to Europe (average 39.7 %). Analysis of this large cohort from Uganda in the context of other African populations reveals variations in KIR and HLA-C1 and C2 that are consistent with migrations within Africa and potential selection pressures on these genes. Our results will help understand how KIR/HLA-C interactions contribute to resistance to pathogens and reproductive success.

Epigenetic regulation of COL15A1 in smooth muscle cell replicative aging and atherosclerosis

Smooth muscle cell (SMC) proliferation is a hallmark of vascular injury and disease. Global hypomethylation occurs during SMC proliferation in culture and in vivo during neointimal formation. Regardless of the programmed or stochastic nature of hypomethylation, identifying these changes is important in understanding vascular disease, as maintenance of a cells' epigenetic profile is essential for maintaining cellular phenotype. Global hypomethylation of proliferating aortic SMCs and concomitant decrease of DNMT1 expression were identified in culture during passage. An epigenome screen identified regions of the genome that were hypomethylated during proliferation and a region containing Collagen, type XV, alpha 1 (COL15A1) was selected by 'genomic convergence' for characterization. COL15A1 transcript and protein levels increased with passage-dependent decreases in DNA methylation and the transcript was sensitive to treatment with 5-Aza-2'-deoxycytidine, suggesting DNA methylation-mediated gene expression. Phenotypically, knockdown of COL15A1 increased SMC migration and decreased proliferation and Col15a1 expression was induced in an atherosclerotic lesion and localized to the atherosclerotic cap. A sequence variant in COL15A1 that is significantly associated with atherosclerosis (rs4142986, P = 0.017, OR = 1.434) was methylated and methylation of the risk allele correlated with decreased gene expression and increased atherosclerosis in human aorta. In summary, hypomethylation of COL15A1 occurs during SMC proliferation and the consequent increased gene expression may impact SMC phenotype and atherosclerosis formation. Hypomethylated genes, such as COL15A1, provide evidence for concomitant epigenetic regulation and genetic susceptibility, and define a class of causal targets that sit at the intersection of genetic and epigenetic predisposition in the etiology of complex disease.

Novel genomic approaches unravel genetic architecture of complex traits in apple

Background

Understanding the genetic architecture of quantitative traits is important for developing genome-based crop improvement methods. Genome-wide association study (GWAS) is a powerful technique for mining novel functional variants. Using a family-based design involving 1,200 apple (Malus × domestica Borkh.) seedlings genotyped for an 8K SNP array, we report the first systematic evaluation of the relative contributions of different genomic regions to various traits related to eating quality and susceptibility to some physiological disorders. Single-SNP analyses models that accounted for population structure, or not, were compared with models fitting all markers simultaneously. The patterns of linkage disequilibrium (LD) were also investigated.

Results

A high degree of LD even at longer distances between markers was observed, and the patterns of LD decay were similar across successive generations. Genomic regions were identified, some of which coincided with known candidate genes, with significant effects on various traits. Phenotypic variation explained by the loci identified through a whole-genome scan ranged from 3% to 25% across different traits, while fitting all markers simultaneously generally provided heritability estimates close to those from pedigree-based analysis. Results from ‘Q+K’ and ‘K’ models were very similar, suggesting that the SNP-based kinship matrix captures most of the underlying population structure. Correlations between allele substitution effects obtained from single-marker and all-marker analyses were about 0.90 for all traits. Use of SNP-derived realized relationships in linear mixed models provided a better goodness-of-fit than pedigree-based expected relationships. Genomic regions with probable pleiotropic effects were supported by the corresponding higher linkage group (LG) level estimated genetic correlations.

Conclusions

The accuracy of artificial selection in plants species can be increased by using more precise marker-derived estimates of realized coefficients of relationships. All-marker analyses that indirectly account for population- and pedigree structure will be a credible alternative to single-SNP analyses in GWAS. This study revealed large differences in the genetic architecture of apple fruit traits, and the marker-trait associations identified here will help develop genome-based breeding methods for apple cultivar development.

Genome-wide linkage analysis of congenital heart defects using MOD score analysis identifies two novel loci

Background

Congenital heart defects (CHD) is the most common cause of death from a congenital structure abnormality in newborns and is often associated with fetal loss. There are many types of CHD. Human genetic studies have identified genes that are responsible for the inheritance of a particular type of CHD and for some types of CHD previously thought to be sporadic. However, occasionally different members of the same family might have anatomically distinct defects — for instance, one member with atrial septal defect, one with tetralogy of Fallot, and one with ventricular septal defect. Our objective is to identify susceptibility loci for CHD in families affected by distinct defects. The occurrence of these apparently discordant clinical phenotypes within one family might hint at a genetic framework common to most types of CHD.

Results

We performed a genome-wide linkage analysis using MOD score analysis in families with diverse CHD. Significant linkage was obtained in two regions, at chromosome 15 (15q26.3, P_empirical = 0.0004) and at chromosome 18 (18q21.2, P_empirical = 0.0005).

Conclusions

In these two novel regions four candidate genes are located: SELS, SNRPA1, and PCSK6 on 15q26.3, and TCF4 on 18q21.2. The new loci reported here have not previously been described in connection with CHD. Although further studies in other cohorts are needed to confirm these findings, the results presented here together with recent insight into how the heart normally develops will improve the understanding of CHD.

Role of TRAV locus in low caries experience

Caries is the most common chronic, multifactorial disease in the world today; and little is still known about the genetic factors influencing susceptibility. Our previous genome-wide linkage scan has identified five loci related to caries susceptibility: 5q13.3, 13q31.1, 14q11.2, 14q 24.3, and Xq27. In the present study, we fine mapped the 14q11.2 locus to identify genetic contributors to caries susceptibility. Four hundred seventy-seven subjects from 72 pedigrees with similar cultural and behavioral habits and limited access to dental care living in the Philippines were studied. An additional 387 DNA samples from unrelated individuals were used to determine allele frequencies. For replication purposes, a total of 1,446 independent subjects from four different populations were analyzed based on their caries experience (low versus high). Forty-eight markers in 14q11.2 were genotyped using TaqMan chemistry. Transmission disequilibrium test was used to detect over transmission of alleles in the Filipino families, and Chi-square, Fisher's exact and logistic regression were used to test for association between low caries experience and variant alleles in the replication data sets. We finally assessed the mRNA expression of TRAV4 in the saliva of 143 study subjects. In the Filipino families, statistically significant associations were found between low caries experience and markers in TRAV4. We were able to replicate these results in the populations studied that were characteristically from underserved areas. Direct sequencing of 22 subjects carrying the associated alleles detects one missense mutation (Y30R) that is predicted to be probably damaging. Finally, we observed higher expression in children and teenagers with low caries experience, correlating with specific alleles in TRAV4. Our results suggest that TRAV4 may have a role in protecting against caries.

Interpreting joint SNP analysis results: when are two distinct signals really two distinct signals?

In genetic association studies, much effort has focused on moving beyond the initial single-nucleotide polymorphism (SNP)-by-SNP analysis. One approach is to reanalyze a chromosomal region where an association has been detected, jointly analyzing the SNP thought to best represent that association with each additional SNP in the region. Such joint analyses may help identify additional, statistically independent association signals. However, it is possible for a single genetic effect to produce joint SNP results that would typically be interpreted as two distinct effects (e.g., both SNPs are significant in the joint model). We present a general approach that can (1) identify conditions under which a single variant could produce a given joint SNP result, and (2) use these conditions to identify variants from a list of known SNPs (e.g., 1000 Genomes) as candidates that could produce the observed signal. We apply this method to our previously reported joint result for smoking involving rs16969968 and rs588765 in CHRNA5. We demonstrate that it is theoretically possible for a joint SNP result suggestive of two independent signals to be produced by a single causal variant. Furthermore, this variant need not be highly correlated with the two tested SNPs or have a large odds ratio. Our method aids in interpretation of joint SNP results by identifying new candidate variants for biological causation that would be missed by traditional approaches. Also, it can connect association findings that may seem disparate due to lack of high correlations among the associated SNPs.

Association mapping of the high-grade myopia MYP3 locus reveals novel candidates UHRF1BP1L, PTPRR, and PPFIA2

PURPOSE

Myopia, or nearsightedness, is a common ocular genetic disease for which over 20 candidate genomic loci have been identified. The high-grade myopia locus, MYP3, has been reported on chromosome 12q21-23 by four independent linkage studies.

METHODS

We performed a genetic association study of the MYP3 locus in a family-based high-grade myopia cohort (n = 82) by genotyping 768 single-nucleotide polymorphisms (SNPs) within the linkage region. Qualitative testing for high-grade myopia (sphere ≤ -5 D affected, > -0.5 D unaffected) and quantitative testing on the average dioptric sphere were performed.

RESULTS

Several genetic markers were nominally significantly associated with high-grade myopia in qualitative testing, including rs3803036, a missense mutation in PTPRR (P = 9.1 × 10(-4)) and rs4764971, an intronic SNP in UHRF1BP1L (P = 6.1 × 10(-4)). Quantitative testing determined statistically significant SNPs rs4764971, also found by qualitative testing (P = 3.1 × 10(-6)); rs7134216, in the 3' untranslated region (UTR) of DEPDC4 (P = 5.4 × 10(-7)); and rs17306116, an intronic SNP within PPFIA2 (P < 9 × 10(-4)). Independently conducted whole genome expression array analyses identified protein tyrosine phosphatase genes PTPRR and PPFIA2, which are in the same gene family, as differentially expressed in normal rapidly growing fetal relative to normal adult ocular tissue (confirmed by RT-qPCR).

CONCLUSIONS

In an independent high-grade myopia cohort, an intronic SNP in UHRF1BP1L, rs4764971, was validated for quantitative association, and SNPs within PTPRR (quantitative) and PPFIA2 (qualitative and quantitative) approached significance. Three genes identified by our association study and supported by ocular expression and/or replication, UHRF1BP1L, PTPRR, and PPFIA2, are novel candidates for myopic development within the MYP3 locus that should be further studied.

Fine-mapping of 5q12.1-13.3 unveils new genetic contributors to caries

Caries is a multifactorial disease and little is still known about the host genetic factors influencing susceptibility. Our previous genome-wide linkage scan has identified the interval 5q12.1-5q13.3 as linked to low caries susceptibility in Filipino families. Here we fine-mapped this region in order to identify genetic contributors to caries susceptibility. Four hundred and seventy-seven subjects from 72 pedigrees with similar cultural and behavioral habits and limited access to dental care living in the Philippines were studied. DMFT scores and genotype data of 75 single-nucleotide polymorphisms were evaluated in the Filipino families with the Family-Based Association Test. For replication purposes, a total 1,467 independent subjects from five different populations were analyzed in a case-control format. In the Filipino cohort, statistically significant and borderline associations were found between low caries experience and four genes spanning 13 million base pairs (PART1, ZSWIM6, CCNB1, and BTF3). We were able to replicate these results in some of the populations studied. We detected PART1 and BTF3 expression in whole saliva, and the expression of BTF3 was associated with caries experience. Our results suggest BTF3 may have a functional role in protecting against caries.

Changes in uric acid levels following bariatric surgery are not associated with SLC2A9 variants in the Swedish Obese Subjects Study

Context and Objective

Obesity and SLC2A9 genotype are strong determinants of uric acid levels. However, data on SLC2A9 variants and weight loss induced changes in uric acid levels are missing. We examined whether the changes in uric acid levels two- and ten-years after weight loss induced by bariatric surgery were associated with SLC2A9 single nucleotide polymorphisms (SNPs) in the Swedish Obese Subjects study.

Methods

SNPs (N = 14) identified by genome-wide association studies and exonic SNPs in the SLC2A9 gene locus were genotyped. Cross-sectional associations were tested before (N = 1806), two (N = 1664) and ten years (N = 1201) after bariatric surgery. Changes in uric acid were compared between baseline and Year 2 (N = 1660) and years 2 and 10 (N = 1172). A multiple testing corrected threshold of P = 0.007 was used for statistical significance.

Results

Overall, 11 of the 14 tested SLC2A9 SNPs were significantly associated with cross-sectional uric acid levels at all three time points, with rs13113918 showing the strongest association at each time point (R² = 3.7−5.2%, 3.9×10⁻²²≤p≤7.7×10⁻¹¹). One SNP (rs737267) showed a significant association (R² = 0.60%, P = 0.002) with change in uric acid levels from baseline to Year 2, as common allele homozygotes (C/C, N = 957) showed a larger decrease in uric acid (−61.4 µmol/L) compared to minor allele carriers (A/X: −51.7 µmol/L, N = 702). No SNPs were associated with changes in uric acid from years 2 to 10.

Conclusions

SNPs in the SLC2A9 locus contribute significantly to uric acid levels in obese individuals, and the associations persist even after considerable weight loss due to bariatric surgery. However, we found little evidence for an interaction between genotype and weight change on the response of uric acid to bariatric surgery over ten years. Thus, the fluctuations in uric acid levels among the surgery group appear to be driven by the weight losses and gains, independent of SLC2A9 genotypes.

Enamel formation genes influence enamel microhardness before and after cariogenic challenge

There is evidence for a genetic component in caries susceptibility, and studies in humans have suggested that variation in enamel formation genes may contribute to caries. For the present study, we used DNA samples collected from 1,831 individuals from various population data sets. Single nucleotide polymorphism markers were genotyped in selected genes (ameloblastin, amelogenin, enamelin, tuftelin, and tuftelin interacting protein 11) that influence enamel formation. Allele and genotype frequencies were compared between groups with distinct caries experience. Associations with caries experience can be detected but they are not necessarily replicated in all population groups and the most expressive results was for a marker in AMELX (p=0.0007). To help interpret these results, we evaluated if enamel microhardness changes under simulated cariogenic challenges are associated with genetic variations in these same genes. After creating an artificial caries lesion, associations could be seen between genetic variation in TUFT1 (p=0.006) and TUIP11 (p=0.0006) with enamel microhardness. Our results suggest that the influence of genetic variation of enamel formation genes may influence the dynamic interactions between the enamel surface and the oral cavity.

Interaction between IRF6 and TGFA genes contribute to the risk of nonsyndromic cleft lip/palate

Previous evidence from tooth agenesis studies suggested IRF6 and TGFA interact. Since tooth agenesis is commonly found in individuals with cleft lip/palate (CL/P), we used four large cohorts to evaluate if IRF6 and TGFA interaction contributes to CL/P. Markers within and flanking IRF6 and TGFA genes were tested using Taqman or SYBR green chemistries for case-control analyses in 1,000 Brazilian individuals. We looked for evidence of gene-gene interaction between IRF6 and TGFA by testing if markers associated with CL/P were overtransmitted together in the case-control Brazilian dataset and in the additional family datasets. Genotypes for an additional 142 case-parent trios from South America drawn from the Latin American Collaborative Study of Congenital Malformations (ECLAMC), 154 cases from Latvia, and 8,717 individuals from several cohorts were available for replication of tests for interaction. Tgfa and Irf6 expression at critical stages during palatogenesis was analyzed in wild type and Irf6 knockout mice. Markers in and near IRF6 and TGFA were associated with CL/P in the Brazilian cohort (p<10⁻⁶). IRF6 was also associated with cleft palate (CP) with impaction of permanent teeth (p<10⁻⁶). Statistical evidence of interaction between IRF6 and TGFA was found in all data sets (p = 0.013 for Brazilians; p = 0.046 for ECLAMC; p = 10⁻⁶ for Latvians, and p = 0.003 for the 8,717 individuals). Tgfa was not expressed in the palatal tissues of Irf6 knockout mice. IRF6 and TGFA contribute to subsets of CL/P with specific dental anomalies. Moreover, this potential IRF6-TGFA interaction may account for as much as 1% to 10% of CL/P cases. The Irf6-knockout model further supports the evidence of IRF6-TGFA interaction found in humans.

Assessing footprints of selection in commercial Atlantic salmon populations using microsatellite data

Relatively large rates of response to traits of economic importance have been observed in different selection experiments in salmon. Several QTL have been mapped in the salmon genome, explaining unprecedented levels of phenotypic variation. Owing to the relatively large selection intensity, individual loci may be indirectly selected, leaving molecular footprints of selection, together with increased inbreeding, as its likely relatives will share the selected loci. We used population differentiation and levels of linkage disequilibrium in chromosomes known to be harbouring QTL for body weight, infectious pancreatic necrosis resistance and infectious salmon anaemia resistance to assess the recent selection history at the genomic level in Atlantic salmon. The results clearly suggest that the marker SSA0343BSFU on chromosome 3 (body weight QTL) showed strong evidence of directional selection. It is intriguing that this marker is physically mapped to a region near the coding sequence of DVL2 , making it an ideal candidate gene to explain the rapid evolutionary response of this chromosome to selection for growth in Salmo salar. Weak evidence of diversifying selection was observed in the QTL associated with infectious pancreatic necrosis and infectious salmon anaemia resistance. Overall, this study showed that artificial selection has produced important changes in the Atlantic salmon genome, validating QTL in commercial salmon populations used for production purposes according to the recent selection history.

Studies of TBX4 and chromosome 17q23.1q23.2: an uncommon cause of nonsyndromic clubfoot

Clubfoot is a common birth defect characterized by inward posturing and rigid downward displacement of one or both feet. The etiology of syndromic forms of clubfoot is varied and the causes of isolated clubfoot are not well understood. A microduplication of 2.2 Mb on chromosome 17q23.1q23.2 which includes T-box 4 (TBX4), a hindlimb-specific gene, and 16 other genes was recently identified in 3 of 66 families reported as nonsyndromic clubfoot, but additional non-foot malformations place them in the syndromic clubfoot category. Our study assesses whether variation in or around TBX4 contributes to nonsyndromic clubfoot. To determine whether this microduplication was a common cause of nonsyndromic clubfoot, 605 probands (from 148 multiplex and 457 simplex families) with nonsyndromic clubfoot were evaluated by copy number and oligonucleotide array CGH testing modalities. Only one multiplex family (0.68%) that had 16 with clubfoot and 9 with other foot anomalies, had a 350 kb microduplication, which included the complete duplication of TBX4 and NACA2 and partial duplication of BRIP1. The microduplication was transmitted in an autosomal dominant pattern and all with the microduplication had a range of phenotypes from short wide feet and toes to bilateral clubfoot. Minimal evidence was found for an association between TBX4 and clubfoot and no pathogenic sequence variants were identified in the two known TBX4 hindlimb enhancer elements. Altogether, these results demonstrate that variation in and around the TBX4 gene and the 17q23.1q23.2 microduplication are not a frequent cause of this common orthopedic birth defect and narrows the 17q23.1q23.2 nonsyndromic clubfoot-associated region.

Refining the genetic portrait of Portuguese Roma through X-chromosomal markers

Due to differences in transmission between X-chromosomal and autosomal DNA, the comparison of data derived from both markers allows deeper insight into the forces that shape the patterns of genetic diversity in populations. In this study, we applied this comparative approach to a sample of Portuguese Roma (Gypsies) by analyzing 43 X-chromosomal markers and 53 autosomal markers. Portuguese individuals of non-Gypsy ancestry were also studied. Compared with the host population, reduced levels of diversity on the X chromosome and autosomes were detected in Gypsies; this result was in line with known patterns of genetic diversity typical of Roma groups. As a consequence of the complex demographic past of the Roma, during which admixture and genetic drift played major roles, the amount of linkage disequilibrium (LD) on the X chromosome in Gypsies was considerably higher than that observed in non-Gypsies. When the pattern of differentiation on the X chromosome was compared with that of autosomes, there was evidence for asymmetries in female and male effective population sizes during the admixture between Roma and non-Roma. This result supplements previous data provided by mtDNA and the Y chromosome, underlining the importance of using combined information from the X chromosome and autosomes to dissect patterns of genetic diversity. Following the out-of-India dispersion, the Roma acquired a complex genetic pattern that was influenced by drift and introgression with surrounding populations, with important contributions from both males and females. We provide evidence that a sex-biased admixture with Europeans is probably associated with the founding of the Portuguese Gypsies.

Genomic selection for fruit quality traits in apple (Malus×domestica Borkh.)

The genome sequence of apple (Malus×domestica Borkh.) was published more than a year ago, which helped develop an 8K SNP chip to assist in implementing genomic selection (GS). In apple breeding programmes, GS can be used to obtain genomic breeding values (GEBV) for choosing next-generation parents or selections for further testing as potential commercial cultivars at a very early stage. Thus GS has the potential to accelerate breeding efficiency significantly because of decreased generation interval or increased selection intensity. We evaluated the accuracy of GS in a population of 1120 seedlings generated from a factorial mating design of four females and two male parents. All seedlings were genotyped using an Illumina Infinium chip comprising 8,000 single nucleotide polymorphisms (SNPs), and were phenotyped for various fruit quality traits. Random-regression best liner unbiased prediction (RR-BLUP) and the Bayesian LASSO method were used to obtain GEBV, and compared using a cross-validation approach for their accuracy to predict unobserved BLUP-BV. Accuracies were very similar for both methods, varying from 0.70 to 0.90 for various fruit quality traits. The selection response per unit time using GS compared with the traditional BLUP-based selection were very high (>100%) especially for low-heritability traits. Genome-wide average estimated linkage disequilibrium (LD) between adjacent SNPs was 0.32, with a relatively slow decay of LD in the long range (r² = 0.33 and 0.19 at 100 kb and 1,000 kb respectively), contributing to the higher accuracy of GS. Distribution of estimated SNP effects revealed involvement of large effect genes with likely pleiotropic effects. These results demonstrated that genomic selection is a credible alternative to conventional selection for fruit quality traits.