Relationships Between Wright's F ST and F IS Statistics in a Context of Wahlund Effect.

Waples (2015) has suggested a formula for the Wahlund effect in a case of unequal contribution of samples from genetically different populations that relates Wright's inbreeding coefficient, F IS , and normalized variance in allele frequencies between populations, F ST . I generalize this relationship to a case of multiple alleles and multiple populations not assuming Hardy-Weinberg ratios prior to mixing. This can help to evaluate the impact of a Wahlund effect on heterozygote deficiency relative to other factors such as null alleles, nonrandom mating, or selection. It is suggested that Wahlund effect cannot be an important factor of deviations from Hardy-Weinberg proportions in natural populations in the majority of instances, but it can have a substantial contribution to heterozygote deficiency in a population that has low genetic diversity compared to that among immigrants or in mixed samples that contain comparable fractions of individuals from genetically different populations.

Separating the post-Glacial coancestry of European and Asian Y chromosomes within haplogroup R1a.

Human Y-chromosome haplogroup structure is largely circumscribed by continental boundaries. One notable exception to this general pattern is the young haplogroup R1a that exhibits post-Glacial coalescent times and relates the paternal ancestry of more than 10% of men in a wide geographic area extending from South Asia to Central East Europe and South Siberia. Its origin and dispersal patterns are poorly understood as no marker has yet been described that would distinguish European R1a chromosomes from Asian. Here we present frequency and haplotype diversity estimates for more than 2000 R1a chromosomes assessed for several newly discovered SNP markers that introduce the onset of informative R1a subdivisions by geography. Marker M434 has a low frequency and a late origin in West Asia bearing witness to recent gene flow over the Arabian Sea. Conversely, marker M458 has a significant frequency in Europe, exceeding 30% in its core area in Eastern Europe and comprising up to 70% of all M17 chromosomes present there. The diversity and frequency profiles of M458 suggest its origin during the early Holocene and a subsequent expansion likely related to a number of prehistoric cultural developments in the region. Its primary frequency and diversity distribution correlates well with some of the major Central and East European river basins where settled farming was established before its spread further eastward. Importantly, the virtual absence of M458 chromosomes outside Europe speaks against substantial patrilineal gene flow from East Europe to Asia, including to India, at least since the mid-Holocene.

Decreased rate of evolution in Y chromosome STR loci of increased size of the repeat unit.

BACKGROUND: Polymorphic Y chromosome short tandem repeats (STRs) have been widely used in population genetic and evolutionary studies. Compared to di-, tri-, and tetranucleotide repeats, STRs with longer repeat units occur more rarely and are far less commonly used. PRINCIPAL FINDINGS: In order to study the evolutionary dynamics of STRs according to repeat unit size, we analysed variation at 24 Y chromosome repeat loci: 1 tri-, 14 tetra-, 7 penta-, and 2 hexanucleotide loci. According to our results, penta- and hexanucleotide repeats have approximately two times lower repeat variance and diversity than tri- and tetranucleotide repeats, indicating that their mutation rate is about half of that of tri- and tetranucleotide repeats. Thus, STR markers with longer repeat units are more robust in distinguishing Y chromosome haplogroups and, in some cases, phylogenetic splits within established haplogroups. CONCLUSIONS: Our findings suggest that Y chromosome STRs of increased repeat unit size have a lower rate of evolution, which has significant relevance in population genetic and evolutionary studies.

Extended Y chromosome haplotypes resolve multiple and unique lineages of the Jewish priesthood.

It has been known for over a decade that a majority of men who self report as members of the Jewish priesthood (Cohanim) carry a characteristic Y chromosome haplotype termed the Cohen Modal Haplotype (CMH). The CMH has since been used to trace putative Jewish ancestral origins of various populations. However, the limited number of binary and STR Y chromosome markers used previously did not provide the phylogenetic resolution needed to infer the number of independent paternal lineages that are encompassed within the Cohanim or their coalescence times. Accordingly, we have genotyped 75 binary markers and 12 Y-STRs in a sample of 215 Cohanim from diverse Jewish communities, 1,575 Jewish men from across the range of the Jewish Diaspora, and 2,099 non-Jewish men from the Near East, Europe, Central Asia, and India. While Cohanim from diverse backgrounds carry a total of 21 Y chromosome haplogroups, 5 haplogroups account for 79.5% of Cohanim Y chromosomes. The most frequent Cohanim lineage (46.1%) is marked by the recently reported P58 T->C mutation, which is prevalent in the Near East. Based on genotypes at 12 Y-STRs, we identify an extended CMH on the J-P58* background that predominates in both Ashkenazi and non-Ashkenazi Cohanim and is remarkably absent in non-Jews. The estimated divergence time of this lineage based on 17 STRs is 3,190 +/- 1,090 years. Notably, the second most frequent Cohanim lineage (J-M410*, 14.4%) contains an extended modal haplotype that is also limited to Ashkenazi and non-Ashkenazi Cohanim and is estimated to be 4.2 +/- 1.3 ky old. These results support the hypothesis of a common origin of the CMH in the Near East well before the dispersion of the Jewish people into separate communities, and indicate that the majority of contemporary Jewish priests descend from a limited number of paternal lineages.

Developing STR databases on structured populations: the native South Siberian population versus the Russian population.

Developing a forensic DNA database on a population that consists of local ethnic groups separated by physical and cultural barriers is questionable as it can be genetically subdivided. On the other side, small sizes of ethnic groups, especially in alpine regions where they are sub-structured further into small villages, prevent collecting a large sample from each ethnic group. For such situations, we suggest to obtain both a total population database on allele frequencies across ethnic groups and a list of theta-values between the groups and the total data. We have genotyped 558 individuals from the native population of South Siberia, consisting of nine ethnic groups, at 17 autosomal STR loci of the kit packages AmpFlSTR SGM Plus i, Cyrillic AmpFlSTR Profiler Plus. The groups differentiate from each other with average theta-values of around 1.1%, and some reach up to three to four percent at certain loci. There exists between-village differentiation as well. Therefore, a database for the population of South Siberia is composed of data on allele frequencies in the pool of ethnic groups and data on theta-values that indicate variation in allele frequencies across the groups. Comparison to additional data on northeastern Asia (the Chukchi and Koryak) shows that differentiation in allele frequencies among small groups that are separated by large geographic distance can be even greater. In contrast, populations of Russians that live in large cities of the European part of Russia are homogeneous in allele frequencies, despite large geographic distance between them, and thus can be described by a database on allele frequencies alone, without any specific information on theta-values.

An STR database on the Volga-Ural population.

This work develops a reference STR database on the Volga-Ural population, Russia, comprised of 640 individuals that were sampled from eight ethnic groups (Finno-Ugric Mari, Mordva-Moksha, Mordva-Erzja, Komi-Permjak, and Udmurt, and Turkic-speaking Bashkir, Tatar-Mishary, and Chuvash) and typed with 10 autosomal STR markers: TH01, CSF1P0, FGA, vWA, D3S1358, TPOX, D16S539, D8S1179, D13S317, FES. The groups differentiate in allele frequencies, and therefore we computed theta-values between allele frequencies in each ethnic groups and those in the database as a measure of their differentiation. Nevertheless, the Volga-Ural ethnic groups form a relatively compact cluster that greatly deviate from the Romanic Moldovans and the Turkic Yakuts, taken for comparison, and are closer to the Slavic Russians, Belarusians, and Ukrainians, although significantly differ from those as well.

A reference data base on STR allele frequencies in the Belarus population developed from paternity cases.

The study considers data on 15 STR-loci from paternity cases across Belarus districts from 2004 to the beginning of 2006 (set #1, 580 individuals) and since that to the beginning of 2007 (set #2, 530 individuals); the data majorly consist of ethnic Belarusians. No significant difference was found between the sets, as well as between the country districts in which the cases occurred. The allele frequencies obtained are very similar to those based on population survey at common loci. Therefore, a data base can be constructed of data from wide survey on paternity cases. Pooling the sets together provides a reference data base on the Belarus population. Additionally, we compared the allele profiles to those in other Slavic groups from the former USSR: Russians from Moscow and eastern Ukrainians from Kharkov city.

Dissecting the molecular architecture and origin of Bayash Romani patrilineages: genetic influences from South-Asia and the Balkans.

The Bayash are a branch of Romanian speaking Roma living dispersedly in Central, Eastern, and Southeastern Europe. To better understand the molecular architecture and origin of the Croatian Bayash paternal gene pool, 151 Bayash Y chromosomes were analyzed for 16 SNPs and 17 STRs and compared with European Romani and non-Romani majority populations from Europe, Turkey, and South Asia. Two main layers of Bayash paternal gene pool were identified: ancestral (Indian) and recent (European). The reduced diversity and expansion signals of H1a patrilineages imply descent from closely related paternal ancestors who could have settled in the Indian subcontinent, possibly as early as between the eighth and tenth centuries AD. The recent layer of the Bayash paternal pool is dominated by a specific subset of E1b1b1a lineages that are not found in the Balkan majority populations. At least two private mutational events occurred in the Bayash during their migrations from the southern Balkans toward Romania. Additional admixture, evident in the low frequencies of typical European haplogroups, J2, R1a, I1, R1b1b2, G, and I2a, took place primarily during the early Bayash settlement in the Balkans and the Romani bondage in Romania. Our results indicate two phenomena in the Bayash and analyzed Roma: a significant preservation of ancestral H1a haplotypes as a result of considerable, but variable level of endogamy and isolation and differential distribution of less frequent, but typical European lineages due to different patterns of the early demographic history in Europe marked by differential admixture and genetic drift.

Y-chromosome diversity characterizes the Gulf of Oman.

Arabia has served as a strategic crossroads for human disseminations, providing a natural connection between the distant populations of China and India in the east to the western civilizations along the Mediterranean. To explore this region's critical role in the migratory episodes leaving Africa to Eurasia and back, high-resolution Y-chromosome analysis of males from the United Arab Emirates (164), Qatar (72) and Yemen (62) was performed. The role of the Levant in the Neolithic dispersal of the E3b1-M35 sublineages is supported by the data, and the distribution and STR-based analyses of J1-M267 representatives points to their spread from the north, most likely during the Neolithic. With the exception of Yemen, southern Arabia, South Iran and South Pakistan display high diversity in their Y-haplogroup substructure possibly a result of gene flow along the coastal crescent-shaped corridor of the Gulf of Oman facilitating human dispersals. Elevated rates of consanguinity may have had an impact in Yemen and Qatar, which experience significant heterozygote deficiencies at various hypervariable autosomal STR loci.

Combining p-values in large-scale genomics experiments.

In large-scale genomics experiments involving thousands of statistical tests, such as association scans and microarray expression experiments, a key question is: Which of the L tests represent true associations (TAs)? The traditional way to control false findings is via individual adjustments. In the presence of multiple TAs, p-value combination methods offer certain advantages. Both Fisher's and Lancaster's combination methods use an inverse gamma transformation. We identify the relation of the shape parameter of that distribution to the implicit threshold value; p-values below that threshold are favored by the inverse gamma method (GM). We explore this feature to improve power over Fisher's method when L is large and the number of TAs is moderate. However, the improvement in power provided by combination methods is at the expense of a weaker claim made upon rejection of the null hypothesis - that there are some TAs among the L tests. Thus, GM remains a global test. To allow a stronger claim about a subset of p-values that is smaller than L, we investigate two methods with an explicit truncation: the rank truncated product method (RTP) that combines the first K-ordered p-values, and the truncated product method (TPM) that combines p-values that are smaller than a specified threshold. We conclude that TPM allows claims to be made about subsets of p-values, while the claim of the RTP is, like GM, more appropriately about all L tests. GM gives somewhat higher power than TPM, RTP, Fisher, and Simes methods across a range of simulations.