POOR PERFORMANCE OF BOOTSTRAP CONFIDENCE INTERVALS FOR THE LOCATION OF A QUANTITATIVE TRAIT LOUCS

The aim of many genetic studies is to locate the genomic regions (called quantitative trait loci, QTLs) that contribute to variation in a quantitative trait (such as body weight). Confidence intervals for the locations of QTLs are particularly important for the design of further experiments to identify the gene or genes responsible for the effect. Likelihood support intervals are the most widely used method to obtain confidence intervals for QTL location, but the non-parametric bootstrap has also been recommended. Through extensive computer simulation, we show that bootstrap confidence intervals are poorly behaved and so should not be used in this context. The profile likelihood (or LOD curve) for QTL location has a tendency to peak at genetic markers, and so the distribution of the maximum likelihood estimate (MLE) of QTL location has the unusual feature of point masses at genetic markers; this contributes to the poor behavior of the bootstrap. Likelihood support intervals and approximate Bayes credible intervals, on the other hand, are shown to behave appropriately

Applied Plasma Research

Contains reports on two research projects.National Science Foundation (Grant GK-37979X)U. S. Army - Research Office - Durham (Contract DAHC04-72-C-0044

Fast and accurate relatedness estimation from high-throughput sequencing data in the presence of inbreeding.

BACKGROUND: The estimation of relatedness between pairs of possibly inbred individuals from high-throughput sequencing (HTS) data has previously not been possible for samples where we cannot obtain reliable genotype calls, as in the case of low-coverage data. RESULTS: We introduce ngsRelateV2, a major revision of ngsRelateV1, a program that originally allowed for estimation of relatedness from HTS data among non-inbred individuals only. The new revised version takes into account the possibility of individuals being inbred by estimating the 9 condensed Jacquard coefficients along with various other relatedness statistics. The program is threaded and scales linearly with the number of cores allocated to the process. CONCLUSION: The program is available as an open source C/C++ program under the GPL license and hosted at https://github.com/ANGSD/ngsRelate. To facilitate easy analysis, the program is able to work directly on the most commonly used container formats for raw sequence (BAM/CRAM) and summary data (VCF/BCF)

Distribution of the

OBJECTIVE: To study whether reported, but inconsistent, associations between the FMR1 CGG repeat lengths in the intermediate, high normal, or low normal range differentiate women diagnosed with diminished ovarian reserve (DOR) from population controls and whether associations vary by race/ethnic group. DESIGN: Case-control study. SETTING: Academic and private fertility clinics. PATIENT(S): DOR cases (n = 129; 95 Whites, 22 Asian, 12 other) from five U.S. fertility clinics were clinically diagnosed, with regular menses and no fragile X syndrome family history. Normal fertility controls (n = 803; 386 Whites, 219 African-Americans, 102 Japanese, 96 Chinese) from the United States-based SWAN Study had one or more menstrual period in the 3 months pre-enrollment, one or more pregnancy, no history of infertility or hormone therapy, and menopause ≥46 years. Previously, the SWAN Chinese and Japanese groups had similar FMR1 CGG repeat lengths, thus they were combined. INTERVENTION(S): None. MAIN OUTCOME MEASURE(S): FMR1 CGG repeat lengths. RESULT(S): Median CGG repeats were nearly identical by case/control group. DOR cases had fewer CGG repeats in the shorter FMR1 allele than controls among Whites, but this was not significant among Asians. White cases had fewer CGG repeats in the shorter allele than Asian cases. No significant differences were found in the high normal/intermediate range between cases and controls or by race/ethnic group within cases in the longer allele. CONCLUSION(S): This study refutes prior reports of an association between DOR and high normal/intermediate repeats and confirms an association between DOR and low normal repeats in Whites

Identifying variants that contribute to linkage for dichotomous and quantitative traits in extended pedigrees

Compared to genome-wide association analysis, linkage analysis is less influenced by allelic heterogeneity. The use of linkage information in large families should provide a great opportunity to identify less frequent variants. We perform a linkage scan for both dichotomous and quantitative traits in eight extended families. For the dichotomous trait, we identified one linkage region on chromosome 4q. For quantitative traits, we identified two regions on chromosomes 4q and 6p for Q1 and one region on chromosome 6q for Q2. To identify variants that contribute to these linkage signals, we performed standard association analysis in genomic regions of interest. We also screened less frequent variants in the linkage region based on the risk ratio and phenotypic distribution among carriers. Two rare variants at VEGFC and one common variant on chromosome 4q conferred the greatest risk for the dichotomous trait. We identified two rare variants on chromosomes 4q (VEGFC) and 6p (VEGFA) that explain 12.4% of the total phenotypic variance of trait Q1. We also identified four variants (including one at VNN3) on chromosome 6q that are able to drop the linkage LOD from 3.7 to 1.0. These results suggest that the use of classical linkage and association methods in large families can provide a useful approach to identifying variants that are responsible for diseases and complex traits in families

Polygenic transcriptome risk scores (PTRS) can improve portability of polygenic risk scores across ancestries

Background: Polygenic risk scores (PRS) are valuable to translate the results of genome-wide association studies (GWAS) into clinical practice. To date, most GWAS have been based on individuals of European-ancestry leading to poor performance in populations of non-European ancestry. Results: We introduce the polygenic transcriptome risk score (PTRS), which is based on predicted transcript levels (rather than SNPs), and explore the portability of PTRS across populations using UK Biobank data. Conclusions: We show that PTRS has a significantly higher portability (Wilcoxon p=0.013) in the African-descent samples where the loss of performance is most acute with better performance than PRS when used in combination

Applied Plasma Research

Contains reports on two research projects.National Science Foundation (Grant GK-28282X1)National Science Foundation (Grant GK-33843)U. S. Army - Research Office - Durham (Contract DAHC04-72-C-0044

Quantitative trait loci for bone traits segregating independently of those for growth in an F-2 broiler X layer cross

An F broiler-layer cross was phenotyped for 18 skeletal traits at 6, 7 and 9 weeks of age and genotyped with 120 microsatellite markers. Interval mapping identified 61 suggestive and significant QTL on 16 of the 25 linkage groups for 16 traits. Thirty-six additional QTL were identified when the assumption that QTL were fixed in the grandparent lines was relaxed. QTL with large effects on the lengths of the tarsometatarsus, tibia and femur, and the weights of the tibia and femur were identified on GGA4 between 217 and 249 cM. Six QTL for skeletal traits were identified that did not co-locate with genome wide significant QTL for body weight and two body weight QTL did not coincide with skeletal trait QTL. Significant evidence of imprinting was found in ten of the QTL and QTL x sex interactions were identified for 22 traits. Six alleles from the broiler line for weight- and size-related skeletal QTL were positive. Negative alleles for bone quality traits such as tibial dyschondroplasia, leg bowing and tibia twisting generally originated from the layer line suggesting that the allele inherited from the broiler is more protective than the allele originating from the layer

To what extent can headteachers be held to account in the practice of social justice leadership?

Internationally, leadership for social justice is gaining prominence as a global travelling theme. This article draws from the Scottish contribution to the International School Leadership Development Network (ISLDN) social justice strand and presents a case study of a relatively small education system similar in size to that of New Zealand, to explore one system's policy expectations and the practice realities of headteachers (principals) seeking to address issues around social justice. Scottish policy rhetoric places responsibility with headteachers to ensure socially just practices within their schools. However, those headteachers are working in schools located within unjust local, national and international contexts. The article explores briefly the emerging theoretical analyses of social justice and leadership. It then identifies the policy expectations, including those within the revised professional standards for headteachers in Scotland. The main focus is on the headteachers' perspectives of factors that help and hinder their practice of leadership for social justice. Macro systems-level data is used to contextualize equity and outcomes issues that headteachers are working to address. In the analysis of the dislocation between policy and reality, the article asks, 'to what extent can headteachers be held to account in the practice of social justice leadership?

Deep-coverage whole genome sequences and blood lipids among 16,324 individuals.

Large-scale deep-coverage whole-genome sequencing (WGS) is now feasible and offers potential advantages for locus discovery. We perform WGS in 16,324 participants from four ancestries at mean depth >29X and analyze genotypes with four quantitative traits-plasma total cholesterol, low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol, and triglycerides. Common variant association yields known loci except for few variants previously poorly imputed. Rare coding variant association yields known Mendelian dyslipidemia genes but rare non-coding variant association detects no signals. A high 2M-SNP LDL-C polygenic score (top 5th percentile) confers similar effect size to a monogenic mutation (~30 mg/dl higher for each); however, among those with severe hypercholesterolemia, 23% have a high polygenic score and only 2% carry a monogenic mutation. At these sample sizes and for these phenotypes, the incremental value of WGS for discovery is limited but WGS permits simultaneous assessment of monogenic and polygenic models to severe hypercholesterolemia