Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Taxonomía y caracterización del arbusto forrajero Atriplex halimus L.

El objetivo de esta Tesis Doctoral ha sido la revisión de la taxonomía infraespecífica de la especie Atriplex halimus L., la delimitación del área geográfica que ocupan sus taxones y su caracterización morfológica. En base a las citas bibliográficas y material de herbario consultado se ha trazado un mapa de distribución de la especie, que delimita su distribución como planta espontánea a la Cuenca Mediterránea, desierto de Siria, Macizo de Hoggar y Sahara Occidental. El resto debe considerarse poblaciones naturalizadas o cultivadas. La taxonomía actual, basada en relaciones filogenéticas, nos indica que lo adecuado es separar las poblaciones de A. halimus en dos clados que se caracterizan por tener diferente nivel de ploidía. Para definir el área de distribución de los clados diploide y tetraploide se ha determinado la ploidía de 19 poblaciones en el espacio comprendido entre las dos áreas de distribución conocidas. Y como resultado el área de distribución de ambos clados puede definirse trazando una recta que une el estrecho de Gibraltar con Estambul, las poblaciones al sur de esta línea son tetraploides a las que habría que añadir las poblaciones del este y sur de la isla de Cerdeña que son también tetraploides. Se han determinado cuales son los parámetros morfológicos que permiten caracterizar ambos clados, para lo cual se han estudiado 52 poblaciones que comprenden toda el área de distribución de la especie. Se han evaluado caracteres cualitativos y cuantitativos de valvas fructíferas, plántulas cultivadas en cámara y arbustos de 2 a 4 años situados en la misma parcela. El único elemento fiable para determinar la pertenencia al clado diploide o tetraploide es la medición del nivel de ploidía. Ninguno de los caracteres morfológicos estudiados es plenamente fiable para determinar si los individuos de una población pertenecen al clado diploide o tetraploide a pesar de realizarse las mediciones sobre poblaciones cultivadas en igualdad de condiciones de suelo y clima, y en condiciones muy favorables para la expresión de caracteres cualitativos y cuantitativos. Los caracteres cualitativos observados que guardan una mayor correlación con el nivel de ploidía son: la base de las hojas de las plántulas, el porte del arbusto, la ramificación, el ángulo de inserción de las ramas y la longitud de la inflorescencia. Respecto a los datos cuantitativos los más relevantes han sido, la longitud de los cotiledones y la longitud del limbo de la 1ª, 3ª o 5ª hoja de la plántula o la razón entre este valor y la anchura del limbo. ABSTRACT The aim of this Thesis has been in the revision of the infraspecies taxa of Atriplex halimus L., the delimitation of their natural growing area and their morphological characterization. Using references to relevant literature and consultation of herbariums, it has been possible to map the distribution of the species as spontaneous plant in the Mediterranean Basin, Desert of Syria, Massif of Hoggar and Western Sahara. The remaining populations must be considered as naturalized or cultivated. The current taxonomy, based on phylogenetic relationships, indicates that A.halimus should be divided in two clades. These clades are characterized by their different ploidy level. Samples of nineteen populations collected from localizations between the defined distribution area of both clades have been characterized for their ploidy level. As a result, it can be seen that the line of clade separation can be trazed from Istambul tothe Strait of Gibraltar, where the populations north to this line are classified as diploid and those south of the line are clsassified as tetraploid. The only exception of this rule is in Sardinia, hwre populations to the West and South of the Island are also classified as tetraploid. The morphological parameters, which include the characterization of both clades, have been defined, and the study of fifty-two populations from the entire distribution area according to these parameters have been analysed. The qualitative and quantitative characteristics of the fruit valves have been evaluated, with seedlings grown in a culture chamber, and 2-4 years old shrubs grown in the same orchard. None of the morphological characters under study are fully reliable to determine whether individuals in a population belong to clade diploid or tetraploid. This is despite measurements being made on populations grown under the same conditions of soil and climate, and in very favorable conditions for the expression of qualitative and quantitative traits. The qualitative characteristics that show a higher correlation with the ploidy level are: the base of the leaves of the seedlings, the bearing bush, the branch arrangement, the angle of insertion of the branches and the inflorescence length. With regards to the quantitative data the most relevant were the length of the cotyledons and the length of the blade of the 1st, 3rd or 5th leaves of the seedling, and the ratio of any of this values and the leaf blade width. The only reliable way to determine the diploid or tetraploid clade membership is the measurement of ploidy level.

Type 2 diabetes: Evidence for linkage on chromosome 20 in 716 Finnish affected sib pairs

We are conducting a genome scan at an average resolution of 10 centimorgans (cM) for type 2 diabetes susceptibility genes in 716 affected sib pairs from 477 Finnish families. To date, our best evidence for linkage is on chromosome 20 with potentially separable peaks located on both the long and short arms. The unweighted multipoint maximum logarithm of odds score (MLS) was 3.08 on 20p (location, x̂ = 19.5 cM) under an additive model, whereas the weighted MLS was 2.06 on 20q (x̂ = 57 cM, recurrence risk, λ̂s = 1.25, P = 0.009). Weighted logarithm of odds scores of 2.00 (x̂ = 69.5 cM, P = 0.010) and 1.92 (x̂ = 18.5 cM, P = 0.013) were also observed. Ordered subset analyses based on sibships with extreme mean values of diabetes-related quantitative traits yielded sets of families who contributed disproportionately to the peaks. Two-hour glucose levels in offspring of diabetic individuals gave a MLS of 2.12 (P = 0.0018) at 9.5 cM. Evidence from this and other studies suggests at least two diabetes-susceptibility genes on chromosome 20. We have also screened the gene for maturity-onset diabetes of the young 1, hepatic nuclear factor 4-a (HNF-4α) in 64 affected sibships with evidence for high chromosomal sharing at its location on chromosome 20q. We found no evidence that sequence changes in this gene accounted for the linkage results we observed.

The estrogen receptor locus is associated with a major gene influencing litter size in pigs.

Identification of individual major genes affecting quantitative traits in livestock species has been limited to date. By using a candidate gene approach and a divergent breed cross involving the Chinese Meishan pig, we have shown that a specific allele of the estrogen receptor (ER) locus is associated with increased litter size. Female pigs from synthetic lines with a 50% Meishan background that were homozygous for this beneficial allele produced 2.3 more pigs in first parities and 1.5 more pigs averaged over all parities than females from the same synthetic lines and homozygous for the undesirable allele. This beneficial ER allele was also found in pigs with Large White breed ancestory. Analysis of females with Large White breed background showed an advantage for females homozygous for the beneficial allele as compared to females homozygous for the other allele of more than 1 total pig born. Analyses of growth performance test records detected no significant unfavorable associations of the beneficial allele with growth and developmental traits. Mapping of the ER gene demonstrated that the closest known genes or markers were 3 centimorgans from ER. To our knowledge, one of these, superoxide dismutase gene (SOD2), was mapped for the first time in the pig. Analysis of ER and these linked markers indicated that ER is the best predictor of litter size differences. Introgression of the beneficial allele into commercial pig breeding lines, in which the allele was not present, and marker-assisted selection for the beneficial allele in lines with Meishan and Large White background have begun.

Mapeamento associativo para múltiplos ambientes e múltiplos locos, visando tolerância à seca em milho

A seca é um dos estresses abióticos mais importantes na cultura do milho, o qual ocasiona reduções significativas na produção de grãos. A arquitetura genética da tolerância à seca é complexa, fazendo-se necessária a melhor compreensão desse caráter. Estudos envolvendo mapeamento associativo são úteis por explorarem a variação genética de caracteres quantitativos e, adicionalmente, levam em conta informações acerca de genótipos, ambientes e interações genótipo por ambiente (G × E). Ao considerar efeitos de G × E em modelos de mapeamento associativo há possibilidade de identificar regiões no genoma associadas à condições e ambientes específicos. Este trabalho teve como objetivo detectar associações relacionadas à tolerância à seca em milho por meio de um modelo de mapeamento associativo para múltiplos ambientes e múltiplos locos, o qual permitiu distinguir associações com efeitos ambiente-específico daquelas com efeitos principais e de interação associação por ambiente (QEI). O painel associativo foi composto por 190 linhagens, classificadas de acordo com os grupos heteróticos quanto ao tipo de grão. Marcadores SNPs (∼500k) foram utilizados para a genotipagem do painel associativo. Duas linhagens (L228-3 e L3) foram usadas como testadores comuns e os híbridos obtidos foram avaliados em duas localidades (Janaúba-MG e Teresina-PI), dois anos agrícolas (2010 e 2011), sob duas condições de tratamento (irrigado e não irrigado). Ao total, consideraram-se seis caracteres: peso de grãos, intervalo de florescimento, florescimento feminino e masculino, altura de planta e de espiga. Consideraram-se dois grupos de mapeamento, agrupados de acordo com os testadores utilizados. SNPs foram úteis para testar associações ao longo do genoma do milho e investigar o relacionamento genético entre indivíduos. O modelo de mapeamento associativo, com inclusão de informações sobre interação G × E, detectou o total de 179 associações, e o maior número de associações foram relacionadas aos caracteres de florescimento. A maioria das associações (168) apresentaram QEI significativo, sendo que o tamanho e a magnitude desses efeitos distinguiram-se de acordo com o ambiente em avaliação. Apenas o caráter florescimento feminino não apresentou associações com efeitos estáveis ao longo dos ambientes em estudo. A detecção de algumas associações em posições próximas do genoma evidenciam possíveis efeitos de pleiotropia. Algumas associações foram co-localizadas em regiões do genoma do milho relacionadas à tolerância à seca, sendo que algumas dessas associações estavam envolvidas a fatores pertencentes à vias metabólicas de interesse. O presente estudo forneceu informações úteis para a compreensão da base genética da tolerância à seca em milho sob os ambientes específicos em avaliação.

Análise da estrutura populacional e do desequilíbrio de ligação de um painel de acessos de sorgo: uma abordagem usando teoria da coalescência

A estrutura populacional e o desequilíbrio de ligação são dois processos fundamentais para estudos evolutivos e de mapeamento associativo. Tradicionalmente, ambos têm sido investigados por meio de métodos clássicos comumente utilizados. Tais métodos certamente forneceram grandes avanços no entendimento dos processos evolutivos das espécies. No entanto, em geral, nenhum deles utiliza uma visão genealógica de forma a considerar eventos genéticos ocorridos no passado, dificultando a compreensão dos padrões de variação observados no presente. Uma abordagem que possibilita a investigação retrospectiva com base no atual polimorfismo observado é a teoria da coalescência. Assim, o objetivo deste trabalho foi analisar, com base na teoria da coalescência, a estrutura populacional e o desequilíbrio de ligação de um painel mundial de acessos de sorgo (Sorghum bicolor). Para tanto, análises de mutação, migração com fluxo gênico e recombinação foram realizadas para cinco regiões genômicas relacionadas à altura de plantas e maturidade (Dw1, Dw2, Dw4, Ma1 e Ma3) e sete populações previamente selecionadas. Em geral, elevado fluxo gênico médio (Μ = m/μ = 41,78 − 52,07) foi observado entre as populações considerando cada região genômica e todas elas simultaneamente. Os padrões sugeriram intenso intercâmbio de acessos e história evolutiva específica para cada região genômica, mostrando a importância da análise individual dos locos. A quantidade média de migrantes por geração (Μ) não foi simétrica entre pares recíprocos de populações, de acordo com a análise individual e simultânea das regiões. Isso sugere que a forma pela qual as populações se relacionaram e continuam interagindo evolutivamente não é igual, mostrando que os métodos clássicos utilizados para investigar estrutura populacional podem ser insatisfatórios. Baixas taxas médias de recombinação (ρL = 2Ner = 0,030 − 0,246) foram observadas utilizando o modelo de recombinação constante ao longo da região. Baixas e altas taxas médias de recombinação (ρr = 2Ner = 0,060 − 3,395) foram estimadas utilizando o modelo de recombinação variável ao longo da região. Os métodos tradicional (r2) e via coalescência (E[r2 rhomap]) utilizados para a estimação do desequilíbrio de ligação mostraram resultados próximos para algumas regiões genômicas e populações. No entanto, o r2 sugeriu padrões descontínuos de desequilíbrio em várias ocasiões, dificultando o entendimento e a caracterização de possíveis blocos de associação. O método via coalescência (E[r2 rhomap]) forneceu resultados que pareceram ter sido mais consistentes, podendo ser uma estratégia eventualmente importante para um refinamento dos padrões não-aleatórios de associação. Os resultados aqui encontrados sugerem que o mapeamento genético a partir de um único pool gênico pode ser insuficiente para detectar associações causais importantes para características quantitativas em sorgo.

Deriva genética de caracteres quantitativos em milho

A obtenção de genótipos superiores no melhoramento de plantas depende da existência de variabilidade genética. A existência de coleções de germoplasma representativas e a utilização de um tamanho adequado de amostra são fundamentais para a preservação das frequências alélicas e genotípicas, diminuindo a perda de variabilidade genética e postergando o aparecimento dos efeitos da deriva genética. Assim, teve-se como objetivo avaliar os efeitos da deriva genética em caracteres quantitativos em subpopulações de milho. Este estudo foi realizado a partir das populações originais BR-105 e BR-106, das quais 10 subpopulações foram obtidas em cada um dos cinco ciclos sucessivos de amostragem com tamanho efetivo reduzido, totalizando 50 subpopulações para cada população original, as quais foram posteriormente autofecundadas, gerando um nível a mais de endogamia. Os tratamentos foram constituídos de 10 amostras da população original sem autofecundação, 10 amostras com autofecundação, 50 subpopulações obtidas da população original e 50 subpopulações autofecundadas, totalizando 120 tratamentos para cada população, avaliados separadamente. Utilizou-se o delineamento em blocos casualizados no esquema de parcelas subdivididas em faixas hierárquico, em quatro ambientes com duas repetições por ambiente. Os caracteres avaliados foram produção de grãos (PG), prolificidade (PROL), comprimento e diâmetro de espigas (CE e DE), número de fileiras por espiga (NFE), número de grãos por fileira (NGF), altura de planta e espiga (AP e AE), florescimento masculino e feminino (FM e FF) e número de ramificações do pendão (NRP). Foram estimados os efeitos da deriva genética entre as médias das subpopulações nos dois níveis de endogamia e os efeitos da depressão por endogamia nas subpopulações dentro dos ciclos. Posteriormente, realizaram-se análises de regressão linear para as subpopulações nos dois níveis de endogamia, separadamente, e em conjunto. Foi verificada uma grande variação nas médias das subpopulações ao longo dos ciclos, indicando que a deriva genética causou diferenciação entre as mesmas e que estas se diferenciaram das populações originais. Detectaram-se efeitos significativos da deriva genética nas populações não autofecundadas para todos os caracteres avaliados, em maior número para PG, já que este caráter é mais sensível à deriva genética por possuir maior grau de dominância que os demais. Houve diminuição no número de estimativas de deriva significativas para as populações autofecundadas, incluindo mudanças na magnitude e no sinal das mesmas em relação às populações não autofecundadas. Para as estimativas de depressão por endogamia, os caracteres PG, NGF, FM e FF apresentaram maior quantidade de estimativas significativas que os demais. Para a maioria dos caracteres, a regressão linear explicou a maior parte da variação encontrada com o aumento dos coeficientes de endogamia. As populações BR-105 e BR-106, por terem estruturas genéticas distintas, apresentaram performances diferentes quanto aos efeitos da deriva genética. Enfim, como a deriva genética interfere na integridade genética das populações, torna-se importante considerar seus efeitos na coleta e manutenção dos bancos de germoplasma e nas populações utilizadas no melhoramento genético de plantas.

A genome scan for eye color in 502 twin families: Most variation is due to a QTL on chromosome 15q

We have rated eye color on a 3-point scale (1=blue/grey, 2=hazel/green, 3=brown) in 502 twin families and carried out a 5-10 cM genome scan (400-757 markers). We analyzed eye color as a threshold trait and performed multipoint sib pair linkage analysis using variance components analysis in Mx. A lod of 19.2 was found at the marker D15S1002, less than 1 cM from OCA2, which has been previously implicated in eye color variation. We estimate that 74% of variance in eye color liability is due to this QTL and a further 18% due to polygenic effects. However, a large shoulder on this peak suggests that other loci affecting eye color may be telomeric of OCA2 and inflating the QTL estimate. No other peaks reached genome-wide significance, although lods >2 were seen on 5p and 14q and lods >1 were additionally seen on chromosomes 2, 3, 6, 7, 8, 9, 17 and 18. Most of these secondary peaks were reduced or eliminated when we repeated the scan as a two locus analysis with the 15q linkage included, although this does not necessarily exclude them as false positives. We also estimated the interaction between the 15q QTL and the other marker locus but there was only minor evidence for additive x additive epistasis. Elaborating the analysis to the full two-locus model including non-additive main effects and interactions did not strengthen the evidence for epistasis. We conclude that most variation in eye color in Europeans is due to polymorphism in OCA2 but that there may be modifiers at several other loci.

Copulas in QTL mapping

The standard variance components method for mapping quantitative trait loci is derived on the assumption of normality. Unsurprisingly, statistical tests based on this method do not perform so well if this assumption is not satisfied. We use the statistical concept of copulas to relax the assumption of normality and derive a test that can perform well under any distribution of the continuous trait. In particular, we discuss bivariate normal copulas in the context of sib-pair studies. Our approach is illustrated by a linkage analysis of lipoprotein(a) levels, whose distribution is highly skewed. We demonstrate that the asymptotic critical levels of the test can still be calculated using the interval mapping approach. The new method can be extended to more general pedigrees and multivariate phenotypes in a similar way as the original variance components method.

Effects of dopamine receptor D4 variation on alcohol and tobacco use and on novelty seeking: Multivariate linkage and association analysis

The dopamine D4 receptor gene contains a polymorphic sequence consisting of a variable number of 48-base-pair (bp) repeats, and there have been a number of reports that this polymorphism is associated with variation in novelty seeking or in substance abuse and addictive behaviors. In this study we have assessed the linkage and association of DRD4 genotype with novelty seeking, alcohol use, and smoking in a sample of 377 dizygotic twin pairs and 15 single twins recruited from the Australian Twin Registry (ATR). We found no evidence of linkage or association of the DRD4 locus with any of the phenotypes. We made use of repeated measures for some phenotypes to increase power by multivariate genetic analysis, but allelic effects were still non-significant. Specifically, it has been suggested that the DRD4 7-repeat allele is associated with increased novelty seeking in males but we found no evidence for this, despite considerable power to do so. We conclude that DRD4 variation does not have an effect on use of alcohol and the problems that arise from it, on smoking, or on novelty seeking behavior. (C) 2003 Wiley-Liss, Inc.

Segregation distortion for seed testa color in mungbean (Vigna radiata L. Wilcek)

Genetic segregation experiments with plant species are commonly used for understanding the inheritance of traits. A basic assumption in these experiments is that each gamete developed from megasporogenesis has an equal chance of fusing with a gamete developed from microsporogenesis, and every zygote formed has an equal chance of survival. If gametic and/or zygotic selection occurs whereby certain gametes or zygotic combinations have a reduced chance of survival, progeny distributions are skewed and are said to exhibit segregation distortion. In this study, inheritance data are presented for the trait seed testa color segregating in large populations (more than 200 individuals) derived from closely related mungbean (Vigna radiata L. Wilcek) taxa. Segregation ratios suggested complex inheritance, including dominant and recessive epistasis. However, this genetic model was rejected in favor of a single-gene model based on evidence of segregation distortion provided by molecular marker data. The segregation distortion occurred after each generation of self-pollination from F-1 thru F-7 resulting in F-7 phenotypic frequencies of 151:56 instead of the expected 103.5:103.5. This study highlights the value of molecular markers for understanding the inheritance of a simply inherited trait influenced by segregation distortion.

Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within- family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.

GENETIC DIVERSITY AMONG SORGHUM LANDRACES OF SOUTHWESTERN HIGHLANDS OF UGANDA

Sorghum ( Sorghum bicolor L. Moench) is an economic and staple crop in sub-Saharan Africa. The genetic diversity in its germplasm is an invaluable aid for its crop improvement. The objective of this study was to assess the existing genetic diversity among sorghum landraces in the southwestern highlands of Uganda. A total of 47 sorghum landraces, collected from southwestern highlands of Uganda, were characterised using 12 qualitative and 13 quantitative traits. The study was conducted at Kachwekano Research Farm in Kabale District, at an altitude of 2,223 m above sea level, during growing season of December 2014 to August 2015. Panicle shape and compactness were the most varied qualitative traits. Grain yield (1.23 to 11.31 t ha-1) and plant height (144.7 to 351.6 cm) were among quantitative traits that showed high variability. Days to 50% flowering (115 to 130 days) showed the least variability. Results of UPGMA cluster analysis generated a dendrogram with three clusters. Panicle weight, leaf width, stem girth, exertion length, peduncle length, panicle shape and compactness, glume colour and threshability were major traits responsible for the observed clustering (P<0.001). Principal Component Analysis revealed the largest variation contributors.

Natural selection and quantitative genetics of life-history traits in Western women: A twin study

Whether contemporary human populations are still evolving as a result of natural selection has been hotly debated. For natural selection to cause evolutionary change in a trait, variation in the trait must be correlated with fitness and be genetically heritable and there must be no genetic constraints to evolution. These conditions have rarely been tested in human populations. In this study, data from a large twin cohort were used to assess whether selection Will cause a change among women in contemporary Western population for three life-history traits: age at menarche, age at first reproduction, and age at menopause. We control for temporal variation in fecundity (the baby boom phenomenon) and differences between women in educational background and religious affiliation. University-educated women have 35% lower fitness than those with less than seven years education, and Roman Catholic women have about 20% higher fitness than those of other religions. Although these differences were significant, education and religion only accounted for 2% and 1% of variance in fitness, respectively. Using structural equation modeling, we reveal significant genetic influences for all three life-history traits, with heritability estimates of 0.50, 0.23, and 0.45, respectively. However, strong genetic covariation with reproductive fitness could only be demonstrated for age at first reproduction, with much weaker covariation for age at menopause and no significant covariation for age at menarche. Selection may, therefore, lead to the evolution of earlier age at first reproduction in this population. We also estimate substantial heritable variation in fitness itself, with approximately 39% of the variance attributable to additive genetic effects, the remainder consisting of unique environmental effects and small effects from education and religion. We discuss mechanisms that could be maintaining such a high heritability for fitness. Most likely is that selection is now acting on different traits from which it did in pre-industrial human populations.

Quantitative trait loci for agronomic and seed quality traits in an F-2 and F-4 : 6 soybean population

Molecular breeding is becoming more practical as better technology emerges. The use of molecular markers in plant breeding for indirect selection of important traits can favorably impact breeding efficiency. The purpose of this research is to identify quantitative trait loci (QTL) on molecular linkage groups (MLG) which are associated with seed protein concentration, seed oil concentration, seed size, plant height, lodging, and maturity, in a population from a cross between the soybean cultivars 'Essex' and 'Williams.' DNA was extracted from F-2 generation soybean leaves and amplified via polymerase chain reaction (PCR) using simple sequence repeat (SSR) markers. Markers that were polymorphic between the parents were analyzed against phenotypic trait data from the F-2 and F-4:6 generation. For the F-2 population, significant additive QTL were Satt540 (MLG M, maturity, r(2)=0.11; height, r(2)=0.04, seed size, r(2)=0.061, Satt373 (MLG L, seed size, r(2)=0.04; height, r(2)=0.14), Satt50 (MLG A1, maturity r(2)=0.07), Satt14 (MLG D2, oil, r(2)=0.05), and Satt251 (protein r(2)=0.03, oil, r(2)=0.04). Significant dominant QTL for the F-2 population were Satt540 (MLG M, height, r(2)=0.04; seed size, r(2)=0.06) and Satt14 (MLG D2, oil, r(2)=0.05). In the F-4:6 generation significant additive QTL were Satt239 (MLG I, height, r(2)=0.02 at Knoxville, TN and r(2)=0.03 at Springfield, TN), Satt14 (MLG D2, seed size, r(2)=0.14 at Knoxville, TN), Satt373 (MLG L, protein, r(2)=0.04 at Knoxville, TN) and Satt251 (MLG B I, lodging r(2)=0.04 at Springfield, TN). Averaged over both environments in the F-4:6 generation, significant additive QTL were identified as Satt251 (MLG B 1, protein, r(2)=0.03), and Satt239 (MLG 1, height, r(2)=0.03). The results found in this study indicate that selections based solely on these QTL would produce limited gains (based on low r(2) values). Few QTL were detected to be stable across environments. Further research to identify stable QTL over environments is needed to make marker-assisted approaches more widely adopted by soybean breeders.