Estimates of genetic parameters for growth and wood properties in Eucalyptus pellita F. Muell. to support tree breeding in Vietnam

Eucalyptus pellita demonstrated good growth and wood quality traits in this study, with young plantation grown timber being suitable for both solid and pulp wood products. All traits examined were under moderate levels of genetic control with little genotype by environment interaction when grown on two contrasting sites in Vietnam. Eucalyptus pellita currently has a significant role in reforestation in the tropics. Research to support expanded of use of this species is needed: particularly, research to better understand the genetic control of key traits will facilitate the development of genetically improved planting stock. This study aimed to provide estimates of the heritability of diameter at breast height over bark, wood basic density, Kraft pulp yield, modulus of elasticity and microfibril angle, and the genetic correlations among these traits, and understand the importance of genotype by environment interactions in Vietnam. Data for diameter and wood properties were collected from two 10-year-old, open-pollinated progeny trials of E. pellita in Vietnam that evaluated 104 families from six native range and three orchard sources. Wood properties were estimated from wood samples using near-infrared (NIR) spectroscopy. Data were analysed using mixed linear models to estimate genetic parameters (heritability, proportion of variance between seed sources and genetic correlations). Variation among the nine sources was small compared to additive variance. Narrow-sense heritability and genetic correlation estimates indicated that simultaneous improvements in most traits could be achieved from selection among and within families as the genetic correlations among traits were either favourable or close to zero. Type B genetic correlations approached one for all traits suggesting that genotype by environment interactions were of little importance. These results support a breeding strategy utilizing a single breeding population advanced by selecting the best individuals across all seed sources. Both growth and wood properties have been evaluated. Multi-trait selection for growth and wood property traits will lead to more productive populations of E. pellita both with improved productivity and improved timber and pulp properties.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

Gene-environment interactions in Autism Spectrum Disorders (ASD)

Spectrum Disorder (ASD), is a heterogeneous neurodevelopmental disorder with na estimated global prevalence rate of 17:10000, and a male to female ratio of 4:1. Patients with ASD presente language and communication difficulties and stereotyped behaviours. Comorbidity with other disorders, such as Intelectual Disability, Fragile-X syndrome (FXS) epilepsy and tuberous sclerosis frequently occurs. ASD presents amultifactorial etiopathology, and genetic factos alone are not suficiente to explain how the syndrome arises, with recente studies establishing ASD heritability at approximately 50%. Pre-, peri- and post-natal exposure to toxic environmental factos has been implicated in the development of ASD. Involvement of epigenetic regulatory mechanisms has been suggested, supported by the occurrence of autistic symptoms in patients with disorders aris ing from epigenetic mutations, such as FXS. A polygenic and epistatic model is a strong hypothesis to explain ASD. The main goal of this project is to identify specific exposure patterns to environmental toxicants in children diagnosed with ASD and integrate the results with genetic and epigenetic data.

Neuropsychiatric disease (NPD) clustering in families with Autism Spectrum Disorder (ASD): genetic, epigenetic and environmental issues

Autism Spectrum Disorder (ASD) is a neurodevelopmental disorder characterized by deficits in social communication/interaction and by unusual repetitive and restricted behaviors and interests. ASD often co-occurs in the same families with other neuropsychiatric diseases (NPD), such as intellectual disability, schizophrenia, epilepsy, depression and attention deficit hyperactivity disorder. Genetic factors have an important role in ASD etiology. Multiple copy number variants (CNVs) and single nucleotide variants (SNVs) in candidate genes have been associated with an increased risk to develop ASD. Nevertheless, recent heritability estimates and the high genotypic and phenotypic heterogeneity characteristic of ASD indicate a role of environmental and epigenetic factors, such as long noncoding RNA (lncRNA) and microRNA (miRNA), as modulators of genetic expression and further clinical presentation. Both miRNA and lncRNA are functional RNA molecules that are transcribed from DNA but not translated into proteins, instead they act as powerful regulators of gene expression. While miRNA are small noncoding RNAs with 22-25 nucleotides in length that act at the post-transcriptional level of gene expression, the lncRNA are bigger molecules (>200 nucleotides in length) that are capped, spliced, and polyadenylated, similar to messenger RNA. Although few lncRNA were well characterized until date, there is a great evidence that they are implicated in several levels of gene expression (transcription/post-transcription/post-translation, organization of protein complexes, cell– cell signaling as well as recombination) as shown in figure 1.

Genetic variation in growth and wood-quality traits of Corymbia citriodora subsp. variegata across three sites in south-east Queensland, Australia

Ten growth or wood-quality traits were assessed in three nearby Corymbia citriodora subsp. variegata (CCV) open-pollinated family-within-provenance trials (18 provenances represented by a total of 374 families) to provide information for the development of a breeding program targeting both pulp and solid-wood products. Growth traits (diameter at breast high over bark [DBH], height and conical volume) were assessed at 3 and 7 years of age. Wood-quality traits (density [DEN], Kraft pulp yield [KPY], modulus of elasticity [MoE] and microfibril angle [MfA]) were predicted using near-infrared spectroscopy on wood samples collected from these trials when aged between 10 and 12 years. The high average KPY, DEN and MoE, and low average MfA observed indicates CCV is very suitable for both pulp and timber products. All traits were under moderate to strong genetic control. In across- trials analyses, high (>0.4) heritability estimates were observed for height, DEN, MoE and MfA, while moderate heritability estimates (0.24 to 0.34) were observed for DBH, volume and KPY. Most traits showed very low levels of genotype × site interaction. Estimated age–age genetic correlations for growth traits were strong at both the family (0.97) and provenance (0.99) levels. Relationships among traits (additive genetic correlation estimates) were favourable, with strong and positive estimates between growth traits (0.84 to 0.98), moderate and positive values between growth and wood-quality traits (0.32 to 0.68), moderate and positive between KPY and MoE (0.64), and high and positive between DEN and MoE (0.82). However, negative (but favourable) correlations were detected between MfA and all other evaluated traits (−0.31 to −0.96). The genetic correlation between the same trait expressed on two different sites, at family level, ranged from 0.24 to 0.42 for growth traits, and from 0.29 to 0.53 for wood traits. Therefore simultaneous genetic improvement of growth and wood property traits in CCV for the target environment in south-east Queensland should be possible, given the moderate to high estimates of heritability and favourable correlations amongst all traits studied, unless genotype × site interactions are greater than was evident. © 2016 NISC (Pty) Ltd

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

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

GENETIC EVALUATION OF SEED TRAITS FROM INTRASPECIFIC CROSSING OF GENETICALLY DISTINCT WATERMELON VARIETIES

Citrullus lanatus (Thunb.) Matsumura and Nakai (Cucurbitaceae) is an important cucurbit crop worldwide. Global production of watermelon is about 90 million metric tonnes per annum, making it among the top five most consumed fresh fruits. The objective of this study was to evaluate seed variability in different segregating populations, and determine heritability of traits of watermelon. Interspecific crosses were made between two cultivars of C. lanatus (Bebu and Wlêwlê Small Seeds (WSS) were performed at Research Station of Nangui Abrogoua University in Abidjan, Côte d’Ivoire. There was wide variability between parental, F1, BC1 (first generation of back-crossing) and F2 seeds. Seeds of all hybrid populations were intermediate versus those of the parents. Also, crossing did not affect F1 and F2 seed characters, but affected those of BC1 because of maternal effects. Thus, back-crossing on Bebu cultivar produced seeds which looked like those of Bebu; while back-crossing on WSS cultivar produced seeds similar to those of WSS. Principal Component Analysis (PCA) and individuals repartitioning revealed that Bebu and WSS cultivars were genetically distinct and showed three main groups: two groups from each parental line and one from a recombinant line (hybrids). F2 population had a wide individual’s dispersion, and contained seeds of all other populations. High heritability was observed for all evaluated characters.

Low-frequency and common genetic variation in ischemic stroke: The METASTROKE collaboration

Erratum in: Low-frequency and common genetic variation in ischemic stroke: The METASTROKE collaboration. [Neurology. 2016]

Estimates of genetic parameters for visual scores and their correlation with production and reproductive traits in Brahman cattle

The objective of this study was to evaluate the association of visual scores of body structure, precocity and muscularity with production (body weight at 18 months and average daily gain) and reproductive (scrotal circumference) traits in Brahman cattle in order to determine the possible use of these scores as selection criteria to improve carcass quality. Covariance components were estimated by the restricted maximum likelihood method using an animal model that included contemporary group as fixed effect. A total of 1,116 observations of body structure, precocity and muscularity were used. Heritability was 0.39, 043 and 0.40 for body structure, precocity and muscularity, respectively. The genetic correlations were 0.79 between body structure and precocity, 0.87 between body structure and muscularity, and 0.91 between precocity and muscularity. The genetic correlations between visual scores and body weight at 18 months were positive (0.77, 0.57 and 0.59 for body structure, precocity and muscularity, respectively). Similar genetic correlations were observed between average daily gain and visual scores (0.60, 0.57 and 0.48, respectively), whereas the genetic correlations between scrotal circumference and these scores were low (0.13, 0.02, and 0.13). The results indicate that visual scores can be used as selection criteria in Brahman breeding programs. Favorable correlated responses should be seen in average daily gain and body weight at 18 months. However, no correlated response is expected for scrotal circumference.

In two waters: contemporary evolution of lagoonal and marine white seabream (Diplodus sargus) populations

Brackish water ecosystems are often exposed to wide variations in environmental variables, including temperature and salinity, which may cause strong selective pressures on organisms modifying the genetic patterns of species. The aim of this work was to test whether there is a ‘divergence-with-gene flow’ in coastal lagoon populations of white seabream (Diplodus sargus) (Ria Formosa, S Portugal and Mar Menor, SE Spain) respect to four marine populations, by using partial sequences of cyt b mitochondrial gene and information from nine microsatellite loci. Genetic diversity was highest in both coastal lagoons (Mar Menor and Ria Formosa) considering mitochondrial and nuclear markers. Although some of FST population pairwise comparisons were not significant, analyses of molecular variance (AMOVAs) detected differences between groups (coastal lagoon and marine) close to significance. Also, only two haplotypes (Cytb-17 and Cytb-18) were detected in both coastal lagoon sampling sites and these localities (Mar Menor and Ria Formosa) showed the highest number of singletons, some of them with a high number of mutations, as has been already described for other Mar Menor populations (Pomatochistus marmoratus and Holothuria polii). Also, several tests detected significant positive and balancing selection considering mtDNA and microsatellite data. These data support the hypothesis of selection as one of the drivers of the genetic differences found between coastal lagoon and marine populations. The life strategy adopted by Diplodus sargus in coastal lagoons allows it to decrease its mortality rate and improve the heritability of its genes. Also, the increase time spent in coastal lagoons with different temperatures and salinities favours the fitness selection and the maintenance of exclusive haplotypes and genotypes in coastal lagoon inhabitants favouring the ‘divergence-with-gene-flow’.

Modelos linear e não linear em análises genéticas para sobrevivência de crias de ovinos da raça Santa Inês.

Resumo: Registros de sobrevivência do nascimento ao desmame de 3846 crias de ovinos da raça Santa Inês foram analisados por modelos de reprodutor linear e não linear (modelo de limiar), para estimar componentes de variância e herdabilidade. Os modelos usados para sobrevivência, analisada como característica da cria, incluíram os efeitos fixos de sexo, da combinação tipo de nascimento-criação da cria e da idade da ovelha ao parto, efeito da covariável peso da cria ao nascer e efeitos aleatórios de reprodutor, da classe rebanho-ano-estação e do resíduo. Componentes de variância para o modelo linear foram estimados pelo método da máxima verossimilhança restrita (REML) e para o modelo não linear por uma aproximação da máxima verossimilhança marginal (MML), pelo programa CMMAT2. O coeficiente de herdabilidade (h2) estimado pelo modelo de limiar foi de 0,29, e pelo modelo linear, 0,14. A correlação de ordem de Spearman entre as capacidades de transmissão dos reprodutores, com base nos dois modelos foi de 0,96. As estimativas de h2 obtidas indicam a possibilidade de se obter, por seleção, ganho genético para sobrevivência. [Linear and nonlinear models in genetic analyses of lamb survival in the Santa Inês hair sheep breed]. Abstract: Records of 3,846 lambs survival from birth to weaning of Santa Inês hair sheep breed, were analyzed by linear and non linear sire models (threshold model) to estimate variance components and heritability (h2). The models that were used to analyze survival, considered in this study as a lamb trait, included the fixed effects of sex of the lamb, combination of type of birth-rearing of lamb, and age of ewe, birth weight of lamb as covariate, and random effects of sire, herd-year-season and residual. Variance components were obtained using restricted maximum likelihood (REML), in linear model and marginal maximum likelihood in threshold model through CMMAT2 program. Estimate of heritability (h2) obtained by threshold model was 0.29 and by linear model was 0.14. Rank correlation of Spearman, between sire solutions based on the two models was 0.96. The obtained estimates in this study indicate that it is possible to acquire genetic gain to survival by selection.

A genetic perspective on rapid evolution in cane toads (Rhinella marina)

The process of biological invasion exposes a species to novel pressures, in terms of both the environments it encounters and the evolutionary consequences of range expansion. Several invaders have been shown to exhibit rapid evolutionary changes in response to those pressures, thus providing robust opportunities to clarify the processes at work during rapid phenotypic transitions. The accelerating pace of invasion of cane toads (Rhinella marina) in tropical Australia during its 80-year history has been well characterized at the phenotypic level, including common-garden experiments that demonstrate heritability of several dispersal-relevant traits. Individuals from the invasion front (and their progeny) show distinctive changes in morphology, physiology and behaviour that, in combination, result in far more rapid dispersal than is true of conspecifics from long-colonized areas. The extensive body of work on cane toad ecology enables us to place into context studies of the genetic basis of these traits. Our analyses of differential gene expression from toads from both ends of this invasion-history transect reveal substantial upregulation of many genes, notably those involved in metabolism and cellular repair. Clearly, then, the dramatically rapid phenotypic evolution of cane toads in Australia has been accompanied by substantial shifts in gene expression, suggesting that this system is well suited to investigating the genetic underpinnings of invasiveness.

Tripartite motif-containing 55 identified as functional candidate for spontaneous cardiac hypertrophy in the rat locus cardiac mass 22

BACKGROUND: Left ventricular (LV) hypertrophy is a risk factor for cardiovascular death, but the genetic factors determining LV size and predisposition to hypertrophy are not well understood. We have previously linked the quantitative trait locus cardiac mass 22 (Cm22) on chromosome 2 with cardiac hypertrophy independent of blood pressure in the spontaneously hypertensive rat. From an original cross of spontaneously hypertensive rat with F344 rats, we derived a normotensive polygenic model of spontaneous cardiac hypertrophy, the hypertrophic heart rat (HHR) and its control strain, the normal heart rat (NHR).

METHODS AND RESULTS: To identify the genes and molecular mechanisms underlying spontaneous LV hypertrophy we sequenced the HHR genome with special focus on quantitative trait locus Cm22. For correlative analyses of function, we measured global RNA transcripts in LV of neonatal HHR and NHR and 198 neonatal rats of an HHR × NHR F2 crossbred population. Only one gene within locus Cm22 was differentially expressed in the parental generation: tripartite motif-containing 55 (Trim55), with mRNA downregulation in HHR (P < 0.05) and reduced protein expression. Trim55 mRNA levels were negatively correlated with LV mass in the F2 cross (r = -0.16, P = 0.025). In exon nine of Trim55 in HHR, we found one missense mutation that functionally alters protein structure. This mutation was strongly associated with Trim55 mRNA expression in F2 rats (F = 10.35, P < 0.0001). Similarly, in humans, we found reduced Trim55 expression in hearts of subjects with idiopathic dilated cardiomyopathy.

CONCLUSION: Our study suggests that the Trim55 gene, located in Cm22, is a novel candidate gene for polygenic LV hypertrophy independent of blood pressure.

Multiscale evaluation of thermal dependence in the glucocorticoid response of vertebrates

Environmental temperature has profound effects on animal physiology, ecology, and evolution. Glucocorticoid (GC) hormones, through effects on phenotypic performance and life history, provide fundamental vertebrate physiological adaptations to environmental variation, yet we lack a comprehensive understanding of how temperature influences GC regulation in vertebrates. Using field studies and metaand comparative phylogenetic analyses, we investigated how acute change and broadscale variation in temperature correlated with baseline and stress-induced GC levels. Glucocorticoid levels were found to be temperature and taxon dependent, but generally, vertebrates exhibited strong positive correlations with acute changes in temperature. Furthermore, reptile baseline, bird baseline, and capture stressinduced GC levels to some extent covaried with broadscale environmental temperature. Thus, vertebrate GC function appears clearly thermally influenced. However, we caution that lack of detailed knowledge of thermal plasticity, heritability, and the basis for strong phylogenetic signal in GC responses limits our current understanding of the role of GC hormones in species’ responses to current and future climate variation.

Behavioral and molecular genetics of reading-related AM and FM detection thresholds

Auditory detection thresholds for certain frequencies of both amplitude modulated (AM) and frequency modulated (FM) dynamic auditory stimuli are associated with reading in typically developing and dyslexic readers. We present the first behavioral and molecular genetic characterization of these two auditory traits. Two extant extended family datasets were given reading tasks and psychoacoustic tasks to determine FM 2 Hz and AM 20 Hz sensitivity thresholds. Univariate heritabilities were significant for both AM (h2 = 0.20) and FM (h2 = 0.29). Bayesian posterior probability of linkage (PPL) analysis found loci for AM (12q, PPL = 81 %) and FM (10p, PPL = 32 %; 20q, PPL = 65 %). Bivariate heritability analyses revealed that FM is genetically correlated with reading, while AM was not. Bivariate PPL analysis indicates that FM loci (10p, 20q) are not also associated with reading.