Expected early genetic gain from selection for milk yield in dairy cattle

Summary A matrix program to predict short term genetic gain from single trait selection for milk yield was developed. Rate of genetic gain was calculated as the annual change in the mean breeding value of all producing females. Several parameters sets representing various selection policies were used to examine situations pertinent to dairy populations of the United States. Approach to the asymptotic rates of genetic gain within the model varied with the choice of parameters, but even with consistent selection policies, predicted total genetic gain in the first 10 years was only half of the expected from classical theory. Considerable year to year variation in the rate of gain occurred. Early gains were more dependent on female selection decisions than gains during the steady state. In a two-phase model, the approach to the linear rate of gain in the second phase was accelerated by starting with an ongoing improvement program, but considerable delays still existed. Selection for sex- limited traits such as milk yield, which require pedigree selection and a waiting time for progeny test results reached asymptotic rates more slowly than previously assumed.     

A universal mechanism generating clusters of differentiated loci during divergence‐with‐migration

Genome‐wide patterns of genetic divergence reveal mechanisms of adaptation under gene flow. Empirical data show that divergence is mostly concentrated in narrow genomic regions. This pattern may arise because differentiated loci protect nearby mutations from gene flow, but recent theory suggests this mechanism is insufficient to explain the emergence of concentrated differentiation during biologically realistic timescales. Critically, earlier theory neglects an inevitable consequence of genetic drift: stochastic loss of local genomic divergence. Here, we demonstrate that the rate of stochastic loss of weak local differentiation increases with recombination distance to a strongly diverged locus and, above a critical recombination distance, local loss is faster than local “gain” of new differentiation. Under high migration and weak selection, this critical recombination distance is much smaller than the total recombination distance of the genomic region under selection. Consequently, divergence between populations increases by net gain of new differentiation within the critical recombination distance, resulting in tightly linked clusters of divergence. The mechanism responsible is the balance between stochastic loss and gain of weak local differentiation, a mechanism acting universally throughout the genome. Our results will help to explain empirical observations and lead to novel predictions regarding changes in genomic architectures during adaptive divergence.     

Current versus historical population sizes in vertebrate species with high gene flow: a comparison based on mitochondrial DNA lineages and inbreeding theory for neutral mutations

Using inbreeding theory as applied to neutral alleles inherited maternally, we generate expected probability distributions of times to identity by descent for random pairs of mitochondrial genotypes within a population or within an entire species characterized by high gene flow. For comparisons with these expectations, empirical distributions of times to most recent common ancestry were calculated (by conventional mtDNA clock calibrations) from mtDNA haplotype distances observed within each of three vertebrate species--American eels, hardhead catfish, and redwinged blackbirds. These species were chosen for analysis because census population size in each is currently large and because both genetic and life-history data are consistent with the postulate that historical gene flow within these species has been high. The observed molecular distances among mtDNA lineages were two to three orders of magnitude lower than predicted from census sizes of breeding females, suggesting that rate of mtDNA evolution is decelerated in these species and/or that long-term effective population size is vastly smaller than present-day population size. Several considerations point to the latter possibility as most likely. The genetic structure of any species is greatly influenced by historical demography; even for species that are currently abundant, mtDNA gene lineages appear to have been channeled through fairly small numbers of ancestors.      

GENES WITH SOCIAL EFFECTS ARE EXPECTED TO HARBOR MORE SEQUENCE VARIATION WITHIN AND BETWEEN SPECIES

The equilibrium sequence diversity of genes within a population and the rate of sequence divergence between populations or species depends on a variety of factors, including expression pattern, mutation rate, nature of selection, random drift, and mating system. Here, we extend population genetic theory developed for maternal-effect genes to predict the equilibrium polymorphism within species and sequence divergence among species for genes with social effects on fitness. We show how the fitness effects of genes, mating system, and genetic system affect predicted gene polymorphism. We find that, because genes with indirect social effects on fitness effectively experience weaker selection, they are expected to harbor higher levels of polymorphism relative to genes with direct fitness effects. The relative increase in polymorphism is proportional to the inverse of the genetic relatedness between individuals expressing the gene and their social partners that experience the fitness effects of the gene. We find a similar pattern of more rapid divergence between populations or species for genes with indirect social effects relative to genes with direct effects. We focus our discussion on the social insects, organisms with diverse indirect genetic effects, mating and genetic systems, and we suggest specific examples for testing our predictions with emerging sociogenomic tools.     

Exome sequencing and targeted gene panels: a simulated comparison of diagnostic yield using data from 158 patients with rare diseases

Next-generation sequencing (NGS) has altered clinical genetic testing by widening the access to molecular diagnosis of genetically determined rare diseases. However, physicians may face difficulties selecting the best diagnostic approach. Our goal is to estimate the rate of possible molecular diagnoses missed by different targeted gene panels using data from a cohort of patients with rare genetic diseases diagnosed with exome sequencing (ES). For this purpose, we simulated a comparison between different targeted gene panels and ES: the list of genes harboring clinically relevant variants from 158 patients was used to estimate the theoretical rate of diagnoses missed by NGS panels from 53 different NGS panels from eight different laboratories. Panels presented a mean rate of missed diagnoses of 64% (range 14%-100%) compared to ES, representing an average predicted sensitivity of 36%. Metabolic abnormalities represented the group with highest mean of missed diagnoses (86%), while seizure represented the group with lowest mean (46%). Focused gene panels are restricted in covering select sets of genes implicated in specific diseases and they may miss molecular diagnoses of rare diseases compared to ES. However, their role in genetic diagnosis remains important especially for well-known genetic diseases with established genetic locus heterogeneity.     

The role of recombination,niche‐specific gene pools and flexible genomes in the ecological speciation of bacteria

Bacteria diversify into genetic clusters analogous to those observed in sexual eukaryotes, but the definition of bacterial species is an ongoing problem. Recent work has focused on adaptation to distinct ecological niches as the main driver of clustering, but there remains debate about the role of recombination in that process. One view is that homologous recombination occurs too rarely for gene flow to constrain divergent selection. Another view is that homologous recombination is frequent enough in many bacterial populations that barriers to gene flow are needed to permit divergence. Niche‐specific gene pools have been proposed as a general mechanism to limit gene flow. We use theoretical models to evaluate additional hypotheses that evolving genetic architecture, specifically the effect sizes of genes and gene gain and loss, can limit gene flow between diverging populations. Our model predicts that (a) in the presence of gene flow and recombination, ecological divergence is concentrated in few loci of large effect and (b) high rates of gene flow plus recombination promote gene loss and favor the evolution of niche‐specific genes. The results show that changing genetic architecture and gene loss can facilitate ecological divergence, even without niche‐specific gene pools. We discuss these results in the context of recent studies of sympatric divergence in microbes.     

Overcoming redundancy: an RNAi enhancer screen for morphogenesis genes in Caenorhabditis elegans

Morphogenesis is an important component of animal development. Genetic redundancy has been proposed to be common among morphogenesis genes, posing a challenge to the genetic dissection of morphogenesis mechanisms. Genetic redundancy is more generally a challenge in biology, as large proportions of the genes in diverse organisms have no apparent loss of function phenotypes. Here, we present a screen designed to uncover redundant and partially redundant genes that function in an example of morphogenesis, gastrulation in Caenorhabditis elegans. We performed an RNA interference (RNAi) enhancer screen in a gastrulation-sensitized double-mutant background, targeting genes likely to be expressed in gastrulating cells or their neighbors. Secondary screening identified 16 new genes whose functions contribute to normal gastrulation in a nonsensitized background. We observed that for most new genes found, the closest known homologs were multiple other C. elegans genes, suggesting that some may have derived from rounds of recent gene duplication events. We predict that such genes are more likely than single copy genes to comprise redundant or partially redundant gene families. We explored this prediction for one gene that we identified and confirmed that this gene and five close relatives, which encode predicted substrate recognition subunits (SRSs) for a CUL-2 ubiquitin ligase, do indeed function partially redundantly with each other in gastrulation. Our results implicate new genes in C. elegans gastrulation, and they show that an RNAi-based enhancer screen in C. elegans can be used as an efficient means to identify important but redundant or partially redundant developmental genes.     

High gene flow promotes the genetic homogeneity of arable weed populations at the landscape level

To determine whether weed populations growing in neighbouring fields were genetically isolated, we investigated the genetic contamination of Alopecurus myosuroides populations in organic fields by populations in conventional fields. Herbicide resistance was used as a marker for gene flow to organic populations, which are not under herbicide selective pressure. Organic fields contained on average 74.5% herbicide-resistant plants (80.1% in conventional fields). Identical resistance alleles were found in neighbouring organic and conventional fields. AFLP analysis revealed that populations from organic and conventional fields had similar neutral genetic diversity. Massive pollen flow from conventional fields is likely chiefly responsible for the genetic makeup of A. myosuroides populations from organic fields. Using a demo-genetic model, we propose that demographic collapses of populations due to effective weed control enhance gene flow towards these populations. Fields with a low weed density could act as ‘genetic sinks’ that would facilitate the diffusion of genes from neighbouring, dense weed populations. Populations of allogamous wind-pollinated weed species like A. myosuroides occurring in neighbouring fields are therefore clearly not independent units. Adaptive evolution and management of such weeds should thus be considered on a scale broader than the field.     

Population genetic structure and history of a generalist parasite infecting multiple sympatric host species

Host specificity is predicted to shape patterns of parasite gene flow between host species; specialist parasites should have low gene flow between host species, while generalists are predicted to have high gene flow between species. However, even for generalist parasites external forces, including ecological differences between host species may sometimes intervene to limit gene flow and create genetic structure. To investigate the potential for cryptic parasite genetic structure to arise under such circumstances, we examined the population genetic structure and history of the generalist nematode, Trichostrongylus axei, infecting six sympatric wild ungulate species in North America. Using genotypes for 186 T. axei larvae at two mitochondrial genes, cox1 and nad4, we found that T. axei was completely panmictic across host species, with 0% of genetic variation structured between host species and 97% within individual hosts. In addition, T. axei showed no evidence of recent genetic bottlenecks, had high nucleotide diversities (above 2%), and an effective population size estimated to be in the tens of millions. Our result that T. axei maintains high rates of gene flow between multiple sympatric host species adds to a growing body of information on trichostrongylid population genetic structure in different ecological contexts. Furthermore, the high rates of gene flow, coupled with high levels of genetic diversity and large effective population size which we observed in T. axei, point to a potentially broad capacity for rapid evolutionary change in this parasite.     

Loss of MHC and neutral variation in Peary caribou: genetic drift is not mitigated by balancing selection or exacerbated by MHC allele distributions

Theory and empirical results suggest that the rate of loss of variation at Mhc and neutral microsatellite loci may differ because selection influences Mhc genes, and because a high proportion of rare alleles at Mhc loci may result in high rates of loss via drift. Most published studies compare Mhc and microsatellite variation in various contemporary populations to infer the effects of population size on genetic variation, even though different populations are likely to have different demographic histories that may also affect contemporary genetic variation. We directly compared loss of variation at Mhc and microsatellite loci in Peary caribou by comparing historical and contemporary samples. We observed that similar proportions of genetic variation were lost over time at each type of marker despite strong evidence for selection at Mhc genes. These results suggest that microsatellites can be used to estimate genome-wide levels of variation, but also that adaptive potential is likely to be lost following population bottlenecks. However, gene conversion and recombination at Mhc loci may act to increase variation following bottlenecks.     

THE GENETIC STRUCTURE OF A GYNODIOECIOUS PLANT: NUCLEAR AND CYTOPLASMIC GENES

Sex expression in gynodioecious plants is often determined by an interaction between biparentally and maternally inherited genes. Their relative rates of gene flow should be considered when modeling the evolution of the sex ratio in structured populations. In order to understand patterns of gene flow in Silene vulgaris, a gynodioecious plant, genetic structure was estimated from biparentally inherited genetic markers (allozymes) and a maternally inherited marker (chloroplast DNA) using Wright's F_st. Based on data from 16 local populations, chloroplast DNA showed considerably more genetic structure than did allozymes (F_st values of 0.62 and 0.22, respectively). This suggests that the rate of gene flow is about three times greater for nuclear genes.     

Joint effect of selection, linkage and partial inbreeding on additive genetic variance in an infinite population

Under the inifinitesimal model of gene effects, selection reduces the additive genetic variance by inducing negative linkage disequilibrium among selected genes. If the selected genes are linked, the decay of linkage disequilibrium is delayed, and the reduction of additive genetic variance is enhanced. Inbreeding in an infinite population also alters the additive genetic variance through the generation of positive association among genes within a locus. In the present study, the joint effect of selection, linkage and partial inbreeding (partial selfing or partial full-sib mating) on the additive genetic variance was modeled. The recurrence relations of the additive genetic variance between successive generations and the prediction equation of the asymptotic additive genetic variance were derived. Numerical computation showed that although partially inbred populations initially maintain larger genetic variances, the accumulated effect of selection overrides the effect of inbreeding. Stochastic simulation was carried out to check the precision of prediction, showing that the obtained equations give a satisfactory prediction during initial generations. However, the predicted values always overestimate the simulated values, especially in later generations. Based on these results, possible extensions and perspectives of the assumed model were discussed.     

Relationships between migration rates and landscape resistance assessed using individual-based simulations

Linking landscape effects on gene flow to processes such as dispersal and mating is essential to provide a conceptual foundation for landscape genetics. It is particularly important to determine how classical population genetic models relate to recent individual-based landscape genetic models when assessing individual movement and its influence on population genetic structure. We used classical Wright-Fisher models and spatially explicit, individual-based, landscape genetic models to simulate gene flow via dispersal and mating in a series of landscapes representing two patches of habitat separated by a barrier. We developed a mathematical formula that predicts the relationship between barrier strength (i.e., permeability) and the migration rate (m) across the barrier, thereby linking spatially explicit landscape genetics to classical population genetics theory. We then assessed the reliability of the function by obtaining population genetics parameters (m, F(ST) ) using simulations for both spatially explicit and Wright-Fisher simulation models for a range of gene flow rates. Next, we show that relaxing some of the assumptions of the Wright-Fisher model can substantially change population substructure (i.e., F(ST) ). For example, isolation by distance among individuals on each side of a barrier maintains an F(ST) of ～0.20 regardless of migration rate across the barrier, whereas panmixia on each side of the barrier results in an F(ST) that changes with m as predicted by classical population genetics theory. We suggest that individual-based, spatially explicit modelling provides a general framework to investigate how interactions between movement and landscape resistance drive population genetic patterns and connectivity across complex landscapes.     

Prediction of response with overlapping generations accounting for multistage selection

Summary A generalization of Hill's equations predicting response to selection is developed that accounts for multiple stage selection in either or both sexes. The method accounts for the flow of genes for animals selected at later stages. This allows for the use of genetic gains from later stages, which explains the reduction in variance due to previous selection. Genetic gains from different selection differentials in each reproductive pathway are incorporated into the equations. The asymptotic response to a single cycle of selection is shown to agree with classical selection theory.The method is applied to a dairy progeny testing scheme representative of an artificial insemination organization in the USA. Two models were compared: (1) the first model accounted for two-stage selection of males, the first stage being based on pedigree information and the second stage on both pedigree and progeny test information; and (2) the second model assumed single-stage male selection. Selection was based on milk volume, milk fat, and milk protein yields. The predicted asymptotic rates for a single cycle of selection were overestimated by 6% and the cumulative response to continuous selection over 20 years was overestimated by 8% by assuming singlestage male selection.Journal Paper No. J14146 of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa; Project No. 1053     

Measuring Gene Flow among Populations Having High Levels of Genetic Fragmentation

We present an analysis of the genetic structures of 22 species of salamanders, with regard to levels of gene flow among populations. We estimate the gene flow parameter, Nm (the product of the effective population number and rate of migration among populations) using two alternative methods described by Wright and Slatkin. For most species, these two methods give approximately congruent estimates of Nm; when estimates differ, the method of Wright produces values slightly larger than those derived by the method of Slatkin. We analyze these results in light of independently derived historical inferences of the fragmentation of populations. This analysis suggests that the Nm values calculated from protein polymorphisms may contain information more relevant to historical patterns of gene exchange than to the current population dynamics; moderately large values of Nm may be calculated for species containing populations known to be no longer exchanging genes. Application of a method for estimating the maximum possible rate of gene exchange among populations indicates that, for most species studied here, gene flow among populations probably is no greater than the mutation rate. We suggest that most plethodontid species cannot be viewed as units whose cohesion is maintained by continuing gene exchange. Furthermore, we suggest that phenotypic uniformity among populations is not easily explained by hypotheses of continual stabilizing selection and propose that future work concentrate upon clarification of the genetic and epigenetic factors conferring self-maintenance or autopoietic properties on living systems.     

An Approach to Population and Evolutionary Genetic Theory for Genes in Mitochondria and Chloroplasts, and Some Results

We developed population genetic theory for organelle genes, using an infinite alleles model appropriate for molecular genetic data, and considering the effects of mutation and random drift on the frequencies of selectively neutral alleles. The effects of maternal inheritance and vegetative segregation of organelle genes are dealt with by defining new effective gene numbers, and substituting these for 2N(e) in classical theory of nuclear genes for diploid organisms. We define three different effective gene numbers. The most general is N(lambda), defined as a function of population size, number of organelle genomes per cell, and proportions of genes contributed by male and female gametes to the zygote. In many organisms, vegetative segregation of organelle genomes and intracellular random drift of organelle gene frequencies combine to produce a predominance of homoplasmic cells within individuals in the population. Then, the effective number of organelle genes is N(eo), a simple function of the numbers of males and females and of the maternal and paternal contributions to the zygote. Finally, when the paternal contribution is very small, N( eo) is closely approximated by the number of females, N( f). Then if the sex ratio is 1, the mean time to fixation or loss of new mutations is approximately two times longer for nuclear genes than for organelle genes, and gene diversity is approximately four times greater. The difference between nuclear and organelle genes disappears or is reversed in animals in which males have large harems. The differences between nuclear and organelle gene behavior caused by maternal inheritance and vegetative segregation are generally small and may be overshadowed by differences in mutation rates to neutral alleles. For monoecious organisms, the effective number of organelle genes is approximately equal to the total population size N. We also show that a population can be effectively subdivided for organelle genes at migration rates which result in panmixis for nuclear genes, especially if males migrate more than females.     

Widespread genetic linkage of mating signals and preferences in the Hawaiian cricket Laupala

The evolution of novel sexual communication systems is integral to the process of speciation, as it discourages gene flow between incipient species. Physical linkage between genes underlying male-female communication (i.e. sexual signals and preferences for them) facilitates both rapid and coordinated divergence of sexual communication systems between populations and reduces recombination in the face of occasional hybridization between diverging populations. Despite these ramifications of the genetic architecture of sexual communication for sexual selection and speciation, few studies have examined this relationship empirically. Previous studies of the closely related Hawaiian crickets Laupala paranigra and Laupala kohalensis have indirectly suggested that many of the genes underlying the difference in pulse rate of male song are physically linked with genes underlying the difference in female preference for pulse rate. Using marker-assisted introgression, we moved 'slow pulse rate' alleles from L. paranigra at five known quantitative trait loci (QTL) underlying male pulse rate into the 'fast pulse rate' genetic background of L. kohalensis and assessed the effect of these loci on female preference. An astounding four out of five song QTL predicted the preferences of female fourth-generation backcrosses, providing direct evidence for the extensive genetic linkage of song and preference in one of the fastest diversifying genera currently known.     

Genetic population structure and contemporary dispersal patterns of a recent European invader, the Chinese mitten crab, Eriocheir sinensis

Genetic studies of recently established populations are challenging because the assumption of equilibrium underlying many analyses is likely to be violated. Using microsatellites, we investigated determinants of genetic structure and migration among invasive European-Chinese mitten crab populations, applying a combination of traditional population genetic analyses and nonequilibrium Bayesian methods. Consistent with their recent history, invasive populations showed much lower levels of genetic diversity than a native Chinese population, indicative of recent bottlenecks. Population differentiation was generally low but significant and especially pronounced among recently established populations. Significant differentiation among cohorts from the same geographical location (River Thames) suggests the low effective population size and associated strong genetic drift that would be anticipated from a very recent colonization. An isolation-by-distance pattern appears to be driven by an underlying correlation between geographical distance and population age, suggesting that cumulative homogenizing gene flow reduces founder bottleneck-associated genetic differentiation between longer-established populations. This hypothesis was supported by a coalescent analysis, which supported a drift + gene flow model as more likely than a model excluding gene flow. Furthermore, admixture analysis identified several recent migrants between the UK and Continental European population clusters. Admixture proportions were significantly predicted by the volume of shipping between sites, indicating that human-mediated transport remains a significant factor for dispersal of mitten crabs after the initial establishment of populations. Our study highlights the value of nonequilibrium methods for the study of invasive species, and also the importance of evaluating nonequilibrium explanations for isolation by distance patterns.     

Prospects for gene therapy

Experiments in animal models and human cells in vitro suggest that gene transfer using retroviral vectors may be useful to treat genetic diseases and to gain information that may improve treatment of other common diseases such as cancer. The approach to treatment of genetic diseases by inserting genes into bone marrow cells and experimental models, and a novel application of gene transfer technology to cancer research are discussed herein.     

Annotation of the domestic dog genome sequence: finding the missing genes

There are over 350 genetically distinct breeds of domestic dog that present considerable variation in morphology, physiology, and disease susceptibility. The genome sequence of the domestic dog was assembled and released in 2005, providing an estimated 20,000 protein-coding genes that are a great asset to the scientific community that uses the dog system as a genetic biomedical model and for comparative and evolutionary studies. Although the canine gene set had been predicted using a combination of ab initio methods, homology studies, motif analysis, and similarity-based programs, it still requires a deep annotation of noncoding genes, alternative splicing, pseudogenes, regulatory regions, and gain and loss events. Such analyses could benefit from new sequencing technologies (RNA-Seq) to better exploit the advantages of the canine genetic system in tracking disease genes. Here, we review the catalog of canine protein-coding genes and the search for missing genes, and we propose rationales for an accurate identification of noncoding genes though next-generation sequencing.