A comparison of genetic variation in two sibling species pairs of haplodiploid insects

Sibling species pairs of sweat bees (Halictus confusus andH. tumulorum) and pine sawflies (Neodiprion pratti andN. maurus) were surveyed for genetic variability using enzyme electrophoresis. Levels of heterozygosity were found to be within the ranges earlier recorded for Hymenoptera. Expected heterozygosities were not significantly higher in the sawflies than in the sweat bees. Estimates of genetic identity between the sibling species were not lower than those generally found for diplodiploid insect species: no evidence was found for an increased rate of evolution in these haplodiploids. Genetic identity data among populations ofH. confusus and betweenHalictus species were within the range expected for conspecific populations and sibling species, respectively. InNeodiprion all genetic distances were low but the two populations ofN. pratti had similar genetic distances as each did toN. maurus, indicating the necessity for further systematic studies of the genus. The research reported here was supported by NSERC operating grants to the junior author. Collections of samples were made through NSERC funding previously available to Drs. R. E. Owen and G. Knerer, L. R. was supported by an Internationaliserringsstipendium from the University of Copenhagen, and L. P. by an NSERC University Research Fellowship.     

Post-glacial partitioning of mitochondrial genetic variation in the field vole

Genetic markers are often used to examine population history. There is considerable debate about the behaviour of molecular clock rates around the population-species transition. Nevertheless, appropriate calibration is critical to any inference regarding the absolute timing and scale of demographic changes. Here, we use a mitochondrial cytochrome b gene genealogy, based entirely on modern sequences and calibrated from recent geophysical events, to date the post-glacial expansion of the Eurasian field vole (Microtus agrestis), a widespread temperate mammal species. The phylogeographic structure reflects the subsequent expansion of populations that went through bottlenecks at the time of the Younger Dryas (ca 12,000 years BP) rather than the Last Glacial Maximum (LGM, ca 24,000 years BP), which is usually seen as the time when present-day patterns were determined. The nucleotide substitution rate that was estimated here, ca 4 × 10(-7) substitutions/site/year, remains extremely high throughout the relevant time frame. Calibration with similarly high population-based substitution rates, rather than long-term rates derived from species divergence times, will show that post-LGM climatic events generated current phylogeographic structure in many other organisms from temperate latitudes.     

Age-related variation in genetic control of height growth in Douglas-fir

Summary The development of genetic variances in height growth of Douglas-fir over a 53-year period is analyzed and found to fall into three periods. In the juvenile period, variances in environmental error increase logarithmically, genetic variance within populations exists at moderate levels, and variance among populations is low but increasing. In the early reproductive period, the response to environmental sources of error variance is restricted, genetic variance within populations disappears, and populational differences strongly emerge but do not increase as expected. In the later period, environmental error again increases rapidly, but genetic variance within populations does not reappear and population differences are maintained at about the same level as established in the early reproductive period. The change between the juvenile and early reproductive periods is perhaps associated with the onset of ecological dominance and significant allocations of energy to reproduction.This paper is published with the approval of the Director of Research, North Carolina Agricultural Experiment Station, as No. 3361 of the Journal Series. The computing services in this project were supported by NIH Grant GM-11 546, held by the Institute of Statistics, North Carolina State University at Raleigh.     

Composite statistics for QTL mapping with moderately discordant sibling pairs

Extreme discordant sibling-pair (EDSP) designs have been shown in theory to be very powerful for mapping quantitative-trait loci (QTLs) in humans. However, their practical applicability has been somewhat limited by the need to phenotype very large populations to find enough pairs that are extremely discordant. In this paper, we demonstrate that there is also substantial power in pairs that are only moderately discordant, and that designs using moderately discordant pairs can yield a more practical balance between phenotyping and genotyping efforts. The power we demonstrate for moderately discordant pairs stems from a new statistical result. Statistical analysis in discordant-pair studies is generally done by testing for reduced identity by descent (IBD) sharing in the pairs. By contrast, the most commonly-used statistical methods for more standard QTL mapping are Haseman-Elston regression and variance-components analysis. Both of these use statistics that are functions of the trait values given IBD information for the pedigree. We show that IBD sharing statistics and "trait value given IBD" statistics contribute complementary rather than redundant information, and thus that statistics of the two types can be combined to form more powerful tests of linkage. We propose a simple composite statistic, and test it with simulation studies. The simulation results show that our composite statistic increases power only minimally for extremely discordant pairs. However, it boosts the power of moderately discordant pairs substantially and makes them a very practical alternative. Our composite statistic is straightforward to calculate with existing software; we give a practical example of its use by applying it to a Genetic Analysis Workshop (GAW) data set.     

Genome scan of Diabrotica virgifera virgifera for genetic variation associated with crop rotation tolerance

Abstract:  Crop rotation has been a valuable technique for control of Diabrotica virgifera virgifera for almost a century. However, during the last two decades, crop rotation has ceased to be effective in an expanding area of the US corn belt. This failure appears to be due to a change in the insect's oviposition behaviour, which, in all probability, has an underlying genetic basis. A preliminary genome scan using 253 amplified fragment-length polymorphism (AFLP) markers sought to identify genetic variation associated with the circumvention of crop rotation. Samples of D. v. virgifera from east-central Illinois, where crop rotation is ineffective, were compared with samples from Iowa at locations that the behavioural variant has yet to reach. A single AFLP marker showed signs of having been influenced by selection for the circumvention of crop rotation. However, this marker was not diagnostic. The lack of markers strongly associated with the trait may be due to an insufficient density of marker coverage throughout the genome. A weak but significant general heterogeneity was observed between the Illinois and Iowa samples at microsatellite loci and AFLP markers. This has not been detected in previous population genetic studies of D. v. virgifera and may indicate a reduction in gene flow between variant and wild-type beetles.     

Mapping genetic loci that determine leukocyte telomere length in a large sample of unselected female sibling pairs

Telomeres play a central role in cellular senescence and cancer pathobiology and are associated with age-related diseases such as atherosclerosis and dementia. Telomere length varies between individuals of the same age, is influenced by DNA-damaging factors such as oxidative stress, and is heritable. We performed a quantitative-trait linkage analysis using an approximate 10-cM genomewide map for mean leukocyte terminal-restriction fragment (TRF) lengths measured by Southern blotting, in 2,050 unselected women aged 18-80 years, comprising 1,025 complete dizygotic twin pairs. Heritability of mean batch-adjusted TRF was 36% (95% confidence interval [CI] 18%-48%), with a large common environmental effect of 49% (95% CI 40%-58%). Significant linkage was observed on chromosome 14 (LOD 3.9) at 14q23.2, and suggestive linkage at 10q26.13 (LOD 2.4) and 3p26.1 (LOD 2.7). This is the first report of loci, mapped in a sample of healthy individuals, that influence mean telomere variation in humans.     

Inter- and intraspecific genetic and morphological variation in a sibling pair of carabid species

Background

Pogonus littoralis and Pogonus chalceus are very close related species with quite different ecological preferences within salt marshes. We study the evolutionary processes in and between these presumably young species. Therefore, we compare the variation in ecologically relevant characters and the genetic variation within one of the species (intraspecific differentiation) with the variation of the two types of characters between the two species (interspecific variation). Data are compared between two independent sets of populations, one set at a small geographical scale (the ecologically diverse Guérande area in France) and the other set at a Atlantic-Mediterranean scale.

Results

Body and relative wing size and IDH1 allozyme data show that the intraspecific variation in P. chalceus is high and in the same range as the interspecific variation (P. chalceus versus P. littoralis). Based on neutral markers (other allozymes and mitochondrial DNA) on the other hand, the intraspecific variation in P. chalceus is much lower in comparison to the interspecific variation.

Conclusion

The different ecotypes in the highly polytypic species P. chalceus are as highly differentiated in ecological characters as true species, but are not recognised as such by screening neutral DNA polymorphisms. This can be interpreted as a case of ongoing speciation driven by natural selection adapting each ecotype to its respective ecological niche. The same ecological process can be recognised in the differentiation between the two sister species, where en plus reproductive isolation between the two gene pools occurred, allowing independent drift and mutation accumulation in neutral genetic characters.     

A powerful and robust new linkage statistic for discordant sibling pairs

Previously, Szatkiewicz and colleagues evaluated the performance of a wide variety of statistics for quantitative-trait-locus linkage, using discordant sibling pairs. They found that the most powerful statistics, in general, were a score statistic and a "composite statistic." However, whereas these two statistics have equal power under ideal conditions, each has limitations that reduce its power in certain circumstances. The score statistic depends on estimates of trait parameters and can lose a lot of power if those estimates are incorrect. The composite statistic is not sensitive to trait-parameter estimates but does depend on arbitrary weights that must be chosen on the basis of the ascertainment scheme. In this report, we elucidate the algebraic relationship between the score and composite statistics and then use that relationship to suggest a new statistic that combines the best properties of both. We call our new statistic the "robust discordant pair" (RDP) statistic. We report simulation studies to show that the RDP statistic does, indeed, have all of the strengths and none of the weaknesses of the score and composite statistics.     

Linkage analysis of extremely discordant and concordant sibling pairs identifies quantitative trait loci influencing variation in human menopausal age

Age at natural menopause may be used as parameter for evaluating the rate of ovarian aging. Environmental factors determine only a small part of the large variation in menopausal age. Studies have shown that genetic factors are likely to be involved in variation in menopausal age. To identify quantitative-trait loci for this trait, we performed a genomewide linkage study with age at natural menopause as a continuous quantitative phenotype in Dutch sister pairs, through use of a selective sampling scheme. A total of 165 families were ascertained using extreme selected sampling and were genotyped for 417 markers. Data were analyzed by Haseman-Elston regression and by an adjusted variance-components analysis. Subgroup analyses for early and late menopausal age were conducted by Haseman-Elston regression. In the adjusted variance-components analysis, 12 chromosomes had a LOD score of > or =1.0. Two chromosomal regions showed suggestive linkage: 9q21.3 (LOD score 2.6) and Xp21.3 (LOD score 3.1). Haseman-Elston regression showed rather similar locations of the peaks but yielded lower LOD scores. A permutation test to obtain empirical P values resulted in a significant peak on the X chromosome. To our knowledge, this is the first study to attempt to identify loci responsible for variability in menopausal age and in which several chromosomal regions were identified with suggestive and significant linkage. Although the finding of the region on the X chromosome comes as no surprise, because of its widespread involvement in premature ovarian failure, the definition of which particular gene is involved is of great interest. The region on chromosome 9 deserves further consideration. Both findings require independent confirmation.     

Adaptation from standing genetic variation

Populations adapt to novel environments in two distinct ways: selection on pre-existing genetic variation and selection on new mutations. These alternative sources of beneficial alleles can result in different evolutionary dynamics and distinct genetic outcomes. Compared with new mutations, adaptation from standing genetic variation is likely to lead to faster evolution, the fixation of more alleles of small effect and the spread of more recessive alleles. There is potential to distinguish between adaptation from standing variation and that from new mutations by differences in the genomic signature of selection. Here we review these approaches and possible examples of adaptation from standing variation in natural populations. Understanding how the source of genetic variation affects adaptation will be integral for predicting how populations will respond to changing environments.     

Recent advances in human quantitative-trait-locus mapping: comparison of methods for discordant sibling pairs

Extreme discordant sibling pairs (EDSPs) are theoretically powerful for the mapping of quantitative-trait loci (QTLs) in humans. EDSPs have not been used much in practice, however, because of the need to screen very large populations to find enough pairs that are extreme and discordant. Given appropriate statistical methods, another alternative is to use moderately discordant sibling pairs (MDSPs)--pairs that are discordant but not at the far extremes of the distribution. Such pairs can be powerful yet far easier to collect than extreme discordant pairs. Recent work on statistical methods for QTL mapping in humans has included a number of methods that, though not developed specifically for discordant pairs, may well be powerful for MDSPs and possibly even EDSPs. In the present article, we survey the new statistics and discuss their applicability to discordant pairs. We then use simulation to study the type I error and the power of various statistics for EDSPs and for MDSPs. We conclude that the best statistic(s) for discordant pairs (moderate or extreme) is (are) to be found among the new statistics. We suggest that the new statistics are appropriate for many other designs as well-and that, in fact, they open the way for the exploration of entirely novel designs.     

Recent advances in human quantitative-trait-locus mapping: comparison of methods for selected sibling pairs

During the past few years, there has been a great deal of new work on methods for mapping quantitative-trait loci by use of sibling pairs and sibships. There are several new methods based on linear regression, as well as several more that are based on score statistics. In theory, most of the new methods should be relatively robust to violations of distributional assumptions and to selected sampling, but, in practice, there has been little evaluation of how the methods perform on selected samples. We survey most of the new regression-based statistics and score statistics and propose a few minor variations on the score statistics. We use simulation to evaluate the type I error and the power of all of the statistics, considering (a) population samples of sibling pairs and (b) sibling pairs ascertained on the basis of at least one sibling with a trait value in the top 10% of the distribution. Most of the statistics have correct type I error for selected samples. The statistics proposed by Xu et al. and by Sham and Purcell are generally the most powerful, along with one of our score statistic variants. Even among the methods that are most powerful for "nice" data, some are more robust than others to non-Gaussian trait models and/or misspecified trait parameters.     

Linkage analysis of extremely discordant and concordant sibling pairs identifies quantitative-trait loci that influence variation in the human personality trait neuroticism

Several theoretical studies have suggested that large samples of randomly ascertained siblings can be used to ascertain phenotypically extreme individuals and thereby increase power to detect genetic linkage in complex traits. Here, we report a genetic linkage scan using extremely discordant and concordant sibling pairs, selected from 34,580 sibling pairs in the southwest of England who completed a personality questionnaire. We performed a genomewide scan for quantitative-trait loci (QTLs) that influence variation in the personality trait of neuroticism, or emotional stability, and we established genomewide empirical significance thresholds by simulation. The maximum pointwise P values, expressed as the negative logarithm (base 10), were found on 1q (3.95), 4q (3.84), 7p (3.90), 12q (4.74), and 13q (3.81). These five loci met or exceeded the 5% genomewide significance threshold of 3.8 (negative logarithm of the P value). QTLs on chromosomes 1, 12, and 13 are likely to be female specific. One locus, on chromosome 1, is syntenic with that reported from QTL mapping of rodent emotionality, an animal model of neuroticism, suggesting that some animal and human QTLs influencing emotional stability may be homologous.     

The probability of parallel genetic evolution from standing genetic variation

Parallel evolution is often assumed to result from repeated adaptation to novel, yet ecologically similar, environments. Here, we develop and analyse a mathematical model that predicts the probability of parallel genetic evolution from standing genetic variation as a function of the strength of phenotypic selection and constraints imposed by genetic architecture. Our results show that the probability of parallel genetic evolution increases with the strength of natural selection and effective population size and is particularly likely to occur for genes with large phenotypic effects. Building on these results, we develop a Bayesian framework for estimating the strength of parallel phenotypic selection from genetic data. Using extensive individual‐based simulations, we show that our estimator is robust across a wide range of genetic and evolutionary scenarios and provides a useful tool for rigorously testing the hypothesis that parallel genetic evolution is the result of adaptive evolution. An important result that emerges from our analyses is that existing studies of parallel genetic evolution frequently rely on data that is insufficient for distinguishing between adaptive evolution and neutral evolution driven by random genetic drift. Overcoming this challenge will require sampling more populations and the inclusion of larger numbers of loci.     

Body height of mummified pharaohs supports historical suggestions of sibling marriages

Body height is an important factor in reconstructing health conditions and it serves as an indicator of socio‐economic status. Researchers rely on ancient data to analyze evolutionary aspects of human health and its interrelation with environmental influences. This study presents body height estimates from all periods of ancient Egyptian history and compares the general population with the existing mummies of the members of royal families. A sample of 259 adult Egyptian mummies originating from various collections and published sources with body lengths (long bone measures or/and overall measurements, CT data) were analyzed, and royal mummies were scored with respect to the level of consanguinity. Male royals were taller than males in the general ancient Egyptian population, while female royals were shorter than females in the general population. The body height variation of the royals is significantly reduced when compared with a pool of non‐royal mummies. This provides evidence for inbreeding resulting from consanguineous marriages. However, there appears to be no correlation between the level of inbreeding and individual body height. The random sample of general population does not show signs of inbreeding. Due to the present lack of larger, technically and ethically challenging genetic studies, the selected non‐invasive approach of body height is the most reliable indicator of sibling marriages of pharaohs based on direct physical evidence. Am J Phys Anthropol 157:519–525, 2015. © 2015 Wiley Periodicals, Inc.     

DNA methylation profiles in sibling pairs discordant for intrauterine exposure to maternal gestational diabetes

Intrauterine exposure to hyperglycemia is reported to confer increased metabolic risk in later life, supporting the ‘developmental origins of health and disease’ hypothesis. Epigenetic alterations are suggested as one of the possible underlying mechanisms. In this study, we compared pairwise DNA methylation differences between siblings whose intrauterine exposure to maternal gestational diabetes (GDM) were discordant. Methylation of peripheral blood DNA of 18 sibling pairs was measured using Infinium HumanMethylation450 BeadChip assays. Of the 465,447 CpG sites analyzed, 12 showed differential methylation (false discovery rate <0.15), including markers within genes associated with monogenic diabetes (HNF4A) or obesity (RREB1). The overall methylation at HNF4A showed inverse correlations with mRNA expression levels, though non significant. In a gene set enrichment analysis, metabolism and signal transduction pathways were enriched. In conclusion, we found DNA methylation markers associated with intrauterine exposure to maternal GDM, including those within genes previously implicated in diabetes or obesity.     

Fine-scale spatial partitioning of genetic variation and evolutionary contestability in the invasive estuarine mussel Xenostrobus securis

Initial studies of Australian populations of the widely invasive mussel Xenostrobus securis raised the possibility of non-random spatial partitioning within estuaries of mitochondrial DNA (mtDNA) clades that have dispersed intercontinentally and those that have not. A fine-scale phylogeographic investigation was made to examine this possibility using cytochrome c oxidase DNA sequences and data from five microsatellite loci developed here. X. securis in the study region was found to comprise multiple robustly supported mtDNA clades, many of which were genetically very distinct. Frequently, the clades comprised numerous haplotypes that have narrow ranges and may have evolved in situ within drainages or nearby. Microsatellite data reveal significant intra- and inter-drainage differentiation but at lower levels than mtDNA. Variation in mtDNA appears to reflect not only recolonization after infrequent local extinction within drainages but also the low probability of migratory haplotypes successfully establishing during evolutionary contests, with resident haplotypes having selective or demographic advantages derived from local occupation. The patterns of mtDNA haplotype distribution suggest that central New South Wales is the source of the internationally invasive lineages of X. securis.     

The impact of transmission-ratio distortion on allele sharing in affected sibling pairs

There have been recent reports of transmission-ratio distortion (TRD) or segregation distortion in families not selected for genetic disease. If TRD exists but is ignored, linkage studies searching for disease genes in affected relatives may be misinterpreted. We show that the identical-by-descent sharing patterns for affected sib pairs are strongly affected by TRD and, further, that the estimated statistical significance of a sib-pair linkage study may be extremely biased. However, we also show that, if TRD is suspected during the planning stage of a study, the planned sample size of the study needs to be increased by only a small amount to maintain the desired power.