Comparison of Copy Number of HSF Genes in Two Buffalo Genomes

The copy number variation (CNV) is the number of copies of a particular gene in the genotype of an individual. Recent evidences show that the CNVs can vary in frequency and occurrence between breeds. These variations reportedly allowed different breeds to adapt to different environments. As copy number variations follow Mendelian pattern of inheritance, identification and distribution of these variants between populations can be used to infer the evolutionary history of the species. In this study, we have examined the absolute copy number of four Heat shock factor genes viz. HSF-1, 2, 4, and 5 in two different breeds of buffalo species using real-time PCR. Here, we report that the absolute copy number of HSF2 varies between the two breeds. In contrast no significant difference was observed in the copy number for HSF-1, 4, and 5 between the two breeds. Our results provide evidence for the presence of breed specific differences in HSF2 genomic copy number. This seems to be the first step in delineating the genetic factors underlying environmental adaptation between the two breeds. Nevertheless, a more detailed study is needed to characterize the functional consequence of this variation.     

Study of polymorphism in human rRNA gene clusters. Population and familial analysis of basic types of RFLP

The frequency and average amount of copy number per genome were defined for standard and a number of new variants of BamHI 5'-NTS RFLP from populations of Moscow, Riga and individuals with Down syndrome. It was demonstrated that the populations studied differ neither in population frequency nor in the average amount of copy number of the variants. New variants were detected in the EcoRI 3'-NTS RFLP system and their amplification, as well as discordance among MZ twins. Possible target for methylation in the HindII site of 3' end of 28S rRNA gene was revealed. Analysis of data obtained demonstrated inefficiency of using the RFLP systems in systematic mapping of NOR-chromosomes. Our data also suggested a possible role of amplification of one copy repeated unit rRNA genes in their evolution.     

The dynamics of transposable elements in structured populations

We analyzed the dynamics of transposable elements (TEs) according to Wright's island and continent-island models, assuming that selection tends to counter the deleterious effects of TEs. We showed that migration between host populations has no impact on either the existence or the stability of the TE copy number equilibrium points obtained in the absence of migration. However, if the migration rate is slower than the transposition rate or if selection is weak, then the TE copy numbers in all the populations can be expected to slowly become homogeneous, whereas a heterogeneous TE copy number distribution between populations is maintained if TEs are mobilized in some populations. The mean TE copy number is highly sensitive to the population size, but as a result of migration between populations, it decreases as the sum of the population sizes increases and tends to reach the same value in these populations. We have demonstrated the existence of repulsion between TE insertion sites, which is established by selection and amplified by drift. This repulsion is reduced as much as the migration rate is higher than the recombination rate between the TE insertion sites. Migration and demographic history are therefore strong forces in determining the dynamics of TEs within the genomes and the populations of a species.     

Mitochondrial DNA Length Variation and Heteroplasmy in Populations of White Sturgeon (Acipenser Transmontanus)

Southern blot analysis was used to quantify the extent of mtDNA D-loop length variation in two populations of white sturgeon, Acipenser transmontanus. Over 42% of individuals were heteroplasmic for up to six different mtDNA length variants attributable to varying copy numbers of an 82-bp repeat sequence. Chi-square analyses revealed that the frequencies of length genotypes and the incidence of heteroplasmy were significantly different between Fraser and Columbia River sturgeon populations but not between restriction site haplotypes. Heteroplasmic fish have, on average, higher copy number than homoplasmic fish. Forty-five of 101 homoplasmic individuals carry only a single copy of the repeat, while none of the 73 heteroplasmic fish has the single repeat as the predominant variant. On the basis of differences in frequency distributions of copy number within and between fish, we suggest that (1) heteroplasmy is maintained by high recurrent mutation of multiple copy genomes, favoring increased copy number and (2) the mutation pressure toward higher copy number heteroplasmy is partially offset by selection to reduced genome size and segregation to the homoplasmic condition.     

Sequence Variation within the KIV-2 Copy Number Polymorphism of the Human LPA Gene in African,Asian, and European Populations

Amazingly little sequence variation is reported for the kringle IV 2 copy number variation (KIV 2 CNV) in the human LPA gene. Apart from whole genome sequencing projects, this region has only been analyzed in some detail in samples of European populations. We have performed a systematic resequencing study of the exonic and flanking intron regions within the KIV 2 CNV in 90 alleles from Asian, European, and four different African populations. Alleles have been separated according to their CNV length by pulsed field gel electrophoresis prior to unbiased specific PCR amplification of the target regions. These amplicons covered all KIV 2 copies of an individual allele simultaneously. In addition, cloned amplicons from genomic DNA of an African individual were sequenced. Our data suggest that sequence variation in this genomic region may be higher than previously appreciated. Detection probability of variants appeared to depend on the KIV 2 copy number of the analyzed DNA and on the proportion of copies carrying the variant. Asians had a high frequency of so-called KIV 2 type B and type C (together 70% of alleles), which differ by three or two synonymous substitutions respectively from the reference type A. This is most likely explained by the strong bottleneck suggested to have occurred when modern humans migrated to East Asia. A higher frequency of variable sites was detected in the Africans. In particular, two previously unreported splice site variants were found. One was associated with non-detectable Lp(a). The other was observed at high population frequencies (10% to 40%). Like the KIV 2 type B and C variants, this latter variant was also found in a high proportion of KIV 2 repeats in the affected alleles and in alleles differing in copy numbers. Our findings may have implications for the interpretation of SNP analyses in other repetitive loci of the human genome.     

A mathematical model for plasmid replication and distribution in microbial populations

A mathematical model and a computer program for its implementation have been developed to predict the distribution of plasmid copy numbers in the individual cells of a microbial population. The kinetics of accumulation of plasmid-free cells. the copy number distribution within the population and the mean copy number can all be calculated using the computer program. The model has been shown to accurately predict these parameters for recombinant plasmids in yeast populations.     

Population Analysis of Large Copy Number Variants and Hotspots of Human Genetic Disease

Copy number variants (CNVs) contribute to human genetic and phenotypic diversity. However, the distribution of larger CNVs in the general population remains largely unexplored. We identify large variants in ~2500 individuals by using Illumina SNP data, with an emphasis on “hotspots” prone to recurrent mutations. We find variants larger than 500 kb in 5%–10% of individuals and variants greater than 1 Mb in 1%–2%. In contrast to previous studies, we find limited evidence for stratification of CNVs in geographically distinct human populations. Importantly, our sample size permits a robust distinction between truly rare and polymorphic but low-frequency copy number variation. We find that a significant fraction of individual CNVs larger than 100 kb are rare and that both gene density and size are strongly anticorrelated with allele frequency. Thus, although large CNVs commonly exist in normal individuals, which suggests that size alone can not be used as a predictor of pathogenicity, such variation is generally deleterious. Considering these observations, we combine our data with published CNVs from more than 12,000 individuals contrasting control and neurological disease collections. This analysis identifies known disease loci and highlights additional CNVs (e.g., 3q29, 16p12, and 15q25.2) for further investigation. This study provides one of the first analyses of large, rare (0.1%–1%) CNVs in the general population, with insights relevant to future analyses of genetic disease.     

Parasite load and MHC diversity in undisturbed and agriculturally modified habitats of the ornate dragon lizard

Major histocompatibility complex (MHC) gene polymorphism is thought to be driven by host–parasite co‐evolution, but the evidence for an association between the selective pressure from parasites and the number of MHC alleles segregating in a population is scarce and inconsistent. Here, we characterized MHC class I polymorphism in a lizard whose habitat preferences (rock outcrops) lead to the formation of well‐defined and stable populations. We investigated the association between the load of ticks, which were used as a proxy for the load of pathogens they transmit, and MHC class I polymorphism across populations in two types of habitat: undisturbed reserves and agricultural land. We hypothesized that the association would be positive across undisturbed reserve populations, but across fragmented agricultural land populations, the relationship would be distorted by the loss of MHC variation due to drift. After controlling for habitat, MHC diversity was not associated with tick number, and the habitats did not differ in this respect. Neither did we detect a difference between habitats in the relationship between MHC and neutral diversity, which was positive across all populations. However, there was extensive variation in the number of MHC alleles per individual, and we found that tick number was positively associated with the average number of alleles carried by lizards across reserve populations, but not across populations from disturbed agricultural land. Our results thus indicate that local differences in selection from parasites may contribute to MHC copy number variation within species, but habitat degradation can distort this relationship.     

Success and failure in replication of genotype–phenotype associations: How does replication help in understanding the genetic basis of phenotypic variation in outbred populations?

Recent developments in sequencing technologies have facilitated genomewide mapping of phenotypic variation in natural populations. Such mapping efforts face a number of challenges potentially leading to low reproducibility. However, reproducible research forms the basis of scientific progress. We here discuss the options for replication and the reasons for potential nonreproducibility. We then review the evidence for reproducible quantitative trait loci (QTL) with a focus on natural animal populations. Existing case studies of replication fall into three categories: (i) traits that have been mapped to major effect loci (including chromosomal inversion and supergenes) by independent research teams; (ii) QTL fine‐mapped in discovery populations; and (iii) attempts to replicate QTL across multiple populations. Major effect loci, in particular those associated with inversions, have been successfully replicated in several cases within and across populations. Beyond such major effect variants, replication has been more successful within than across populations, suggesting that QTL discovered in natural populations may often be population‐specific. This suggests that biological causes (differences in linkage patterns, allele frequencies or context‐dependencies of QTL) contribute to nonreproducibility. Evidence from other fields, notably animal breeding and QTL mapping in humans, suggests that a significant fraction of QTL is indeed reproducible in direction and magnitude at least within populations. However, there is also a large number of QTL that cannot be easily reproduced. We put forward that more studies should explicitly address the causes and context‐dependencies of QTL signals, in particular to disentangle linkage differences, allele frequency differences and gene‐by‐environment interactions as biological causes of nonreproducibility of QTL, especially between populations.     

Population genetics models of transposable elements

The control of transposable element copy number is of considerable theoretical and empirical interest. Under simple models, copy numbers may increase without limit. Mechanisms that can prevent such an increase include those in which the effect of selection increases with copy number, those in which the rate of transposition decreases with copy number, and those where unlimited increase in copy number is prevented by the consequences of functional heterogeneity in the transposable element family. Finite population sizes may attenuate the power of natural selection to act on transposable element copy number in a number of ways that may be of particular importance in laboratory populations. First, a small host population size will create occasional periods in which the variance between individuals in copy number is diminished, and with it the power of natural selection, even when the expected variance is Poisson. Second, small population sizes will produce high-frequency transposable element sites, systematically reducing the variance in copy number. The consequences will be particularly profound when the selective damage of transposable elements follows from their heterozygosity, as when ectopic exchange limits copy number. This revised version was published online in August 2006 with corrections to the Cover Date.     

Copy Number Variations in Tilapia Genomes

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R ² > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.     

Segregation of Naturally Occurring Mitochondrial DNA Variants in a Mini-Pig Model

The maternally inherited mitochondrial genome (mtDNA) is present in multimeric form within cells and harbors sequence variants (heteroplasmy). While a single mtDNA variant at high load can cause disease, naturally occurring variants likely persist at low levels across generations of healthy populations. To determine how naturally occurring variants are segregated and transmitted, we generated a mini-pig model, which originates from the same maternal ancestor. Following next-generation sequencing, we identified a series of low-level mtDNA variants in blood samples from the female founder and her daughters. Four variants, ranging from 3% to 20%, were selected for validation by high-resolution melting analysis in 12 tissues from 31 animals across three generations. All four variants were maintained in the offspring, but variant load fluctuated significantly across the generations in several tissues, with sex-specific differences in heart and liver. Moreover, variant load was persistently reduced in high-respiratory organs (heart, brain, diaphragm, and muscle), which correlated significantly with higher mtDNA copy number. However, oocytes showed increased heterogeneity in variant load, which correlated with increased mtDNA copy number during in vitro maturation. Altogether, these outcomes show that naturally occurring mtDNA variants segregate and are maintained in a tissue-specific manner across generations. This segregation likely involves the maintenance of selective mtDNA variants during organogenesis, which can be differentially regulated in oocytes and preimplantation embryos during maturation.     

Partitioning of genetic variation across the genome using multimarker methods in a wild bird population

The underlying basis of genetic variation in quantitative traits, in terms of the number of causal variants and the size of their effects, is largely unknown in natural populations. The expectation is that complex quantitative trait variation is attributable to many, possibly interacting, causal variants, whose effects may depend upon the sex, age and the environment in which they are expressed. A recently developed methodology in animal breeding derives a value of relatedness among individuals from high‐density genomic marker data, to estimate additive genetic variance within livestock populations. Here, we adapt and test the effectiveness of these methods to partition genetic variation for complex traits across genomic regions within ecological study populations where individuals have varying degrees of relatedness. We then apply this approach for the first time to a natural population and demonstrate that genetic variation in wing length in the great tit (Parus major) reflects contributions from multiple genomic regions. We show that a polygenic additive mode of gene action best describes the patterns observed, and we find no evidence of dosage compensation for the sex chromosome. Our results suggest that most of the genomic regions that influence wing length have the same effects in both sexes. We found a limited amount of genetic variance in males that is attributed to regions that have no effects in females, which could facilitate the sexual dimorphism observed for this trait. Although this exploratory work focuses on one complex trait, the methodology is generally applicable to any trait for any laboratory or wild population, paving the way for investigating sex‐, age‐ and environment‐specific genetic effects and thus the underlying genetic architecture of phenotype in biological study systems.     

A coalescent model of ancestry for a rare allele

In disequilibrium mapping from data on a rare allele, interest may focus on the ancestry of a random sample of current descendants of a mutation. The mutation is assumed to have been introduced into the population as a single copy a known time ago and to have reached a given copy number within the population. Theory has been developed to describe the ancestral distribution under arbitrary patterns of population expansion. Further results permit convenient realization of the ancestry for a random sample of copies of a rare allele within populations of constant size or within populations growing or shrinking at constant exponential rate. In this article, we present an efficient approximate method for realizing coalescence times under more general patterns of population growth. We also apply diagnostics, checking the age of the mutation. In the course of the derivation, some additional insight is gained into the dynamics of the descendants of the mutation.     

Demography and the Age of Rare Variants

Large whole-genome sequencing projects have provided access to much rare variation in human populations, which is highly informative about population structure and recent demography. Here, we show how the age of rare variants can be estimated from patterns of haplotype sharing and how these ages can be related to historical relationships between populations. We investigate the distribution of the age of variants occurring exactly twice ( variants) in a worldwide sample sequenced by the 1000 Genomes Project, revealing enormous variation across populations. The median age of haplotypes carrying variants is 50 to 160 generations across populations within Europe or Asia, and 170 to 320 generations within Africa. Haplotypes shared between continents are much older with median ages for haplotypes shared between Europe and Asia ranging from 320 to 670 generations. The distribution of the ages of haplotypes is informative about their demography, revealing recent bottlenecks, ancient splits, and more modern connections between populations. We see the effect of selection in the observation that functional variants are significantly younger than nonfunctional variants of the same frequency. This approach is relatively insensitive to mutation rate and complements other nonparametric methods for demographic inference.     

Rates of recombination in the ribosomal DNA of apomictically propagated Daphnia obtusa lines

Ribosomal (r)DNA undergoes concerted evolution, the mechanisms of which are unequal crossing over and gene conversion. Despite the fundamental importance of these mechanisms to the evolution of rDNA, their rates have been estimated only in a few model species. We estimated recombination rate in rDNA by quantifying the relative frequency of intraindividual length variants in an expansion segment of the 18S rRNA gene of the cladoceran crustacean, Daphnia obtusa, in four apomictically propagated lines. We also used quantitative PCR to estimate rDNA copy number. The apomictic lines were sampled every 5 generations for 90 generations, and we considered each significant change in the frequency distribution of length variants between time intervals to be the result of a recombination event. Using this method, we calculated the recombination rate for this region to be 0.02-0.06 events/generation on the basis of three different estimates of rDNA copy number. In addition, we observed substantial changes in rDNA copy number within and between lines. Estimates of haploid copy number varied from 53 to 233, with a mean of 150. We also measured the relative frequency of length variants in 30 lines at generations 5, 50, and 90. Although length variant frequencies changed significantly within and between lines, the overall average frequency of each length variant did not change significantly between the three generations sampled, suggesting that there is little or no bias in the direction of change due to recombination.     

A Genomic Portrait of Haplotype Diversity and Signatures of Selection in Indigenous Southern African Populations

We report a study of genome-wide, dense SNP (∼900K) and copy number polymorphism data of indigenous southern Africans. We demonstrate the genetic contribution to southern and eastern African populations, which involved admixture between indigenous San, Niger-Congo-speaking and populations of Eurasian ancestry. This finding illustrates the need to account for stratification in genome-wide association studies, and that admixture mapping would likely be a successful approach in these populations. We developed a strategy to detect the signature of selection prior to and following putative admixture events. Several genomic regions show an unusual excess of Niger-Kordofanian, and unusual deficiency of both San and Eurasian ancestry, which were considered the footprints of selection after population admixture. Several SNPs with strong allele frequency differences were observed predominantly between the admixed indigenous southern African populations, and their ancestral Eurasian populations. Interestingly, many candidate genes, which were identified within the genomic regions showing signals for selection, were associated with southern African-specific high-risk, mostly communicable diseases, such as malaria, influenza, tuberculosis, and human immunodeficiency virus/AIDs. This observation suggests a potentially important role that these genes might have played in adapting to the environment. Additionally, our analyses of haplotype structure, linkage disequilibrium, recombination, copy number variation and genome-wide admixture highlight, and support the unique position of San relative to both African and non-African populations. This study contributes to a better understanding of population ancestry and selection in south-eastern African populations; and the data and results obtained will support research into the genetic contributions to infectious as well as non-communicable diseases in the region.