Nucleotide diversity patterns of three divergent soybean populations: evidences for population‐dependent linkage disequilibrium and taxonomic status of Glycine gracilis

The level of linkage disequilibrium (LD) is a major factor to determine DNA polymorphism pattern of a population and to construct high‐resolution maps useful in localizing and gene cloning of complicated traits. Here, we investigated LD level of three soybean populations with different genetic backgrounds and taxonomic status of G. gracilis by comparing the DNA polymorphism patterns of four high‐diversity single‐copy nuclear genes. A total of 152, 22, and 77 accessions of G. soja, G. gracilis, and G. max were observed. The results indicated that G. max retained only 75.3 (π) and 39% (θ) of the nucleotide polymorphism found in G. soja. Four gene loci evolved according to neutrality in both G. max and G. gracilis populations, and three gene loci evolved according to neutrality in G. soja population by Tajima's and Fu and Li's test. However, one gene locus deviated from neutrality by Fu and Li's test in the G. soja population. Further, medial level of LD (average r² = 0.2426) was found in intragene in G. max and G. gracilis populations, but unexpected low level of LD (r² ≤ 0.0539) was found in G. soja population. Significant genetic differentiation was detected between G. max and G. soja populations and also between G. max and G. gracilis populations; however, nonsignificant genetic differentiation was found between G. gracilis and G. soja populations. The results suggest that LD level depends on genetic background of soybean population, and implicit that G. gracilis should be regarded as the variant of G. soja, not as an independent species.     

Genomic conflict drives patterns of X‐linked population structure in Drosophila neotestacea

Intragenomic conflict has the potential to cause widespread changes in patterns of genetic diversity and genome evolution. In this study, we investigate the consequences of sex‐ratio (SR) drive on the population genetic patterns of the X‐chromosome in Drosophila neotestacea. An SR X‐chromosome prevents the maturation of Y‐bearing sperm during male spermatogenesis and thus is transmitted to ~100% of the offspring, nearly all of which are daughters. Selection on the rest of the genome to suppress SR can be strong, and the resulting conflict over the offspring sex ratio can result in the accumulation of multiple loci on the X‐chromosome that are necessary for the expression of drive. We surveyed variation at 12 random X‐linked microsatellites across 16 populations of D. neotestacea that range in SR frequency from 0% to 30%. First, every locus was differentiated between SR and wild‐type chromosomes, and this drives genetic structure at the X‐chromosome. Once the association with SR is accounted for, the patterns of differentiation among populations are similar to the autosomes. Second, within wild‐type chromosomes, the relative heterozygosity is reduced in populations with an increased prevalence of drive, and the heterozygosity of SR chromosomes is higher than expected based on its prevalence. The combination of the relatively high prevalence of SR drive and the structuring of polymorphism between the SR and wild‐type chromosomes suggests that genetic conflict because of SR drive has had significant consequences on the patterns of X‐linked polymorphism and thus also probably affects the tempo of X‐chromosome evolution in D. neotestacea.     

Population structure and linkage disequilibrium in a large collection of Fusarium oxysporum strains analysed through iPBS markers

In the current study, 160 pathogenic strains of Fusarium oxysporum collected from tomato, eggplant and pepper were studied. Eighteen inter‐primer binding site (iPBS)‐retrotransposon primers were used, and these primers generated 205 scorable polymorphic bands. The number of polymorphic bands per primer varied between 9 and 19, with a mean of 11 bands per primer. The highest polymorphism information content (PIC) value was determined as 0.27, and the lowest was 0.05. The unweighted pair‐group method with arithmetic averages (UPGMA) dendrogram including a heat map revealed that the 160 pathogenic strains of F. oxysporum were divided into two main clusters. The first cluster mainly included F. oxysporum f. sp. capsici (FOC) and F. oxysporum f. sp. melongenae (FOMG) isolates. The second cluster mainly comprised F. oxysporum f. sp. lycopersici (FOL) and F. oxysporum f. sp. radicis lycopersici (FORL) isolates. The highest percentage of loci in significant linkage disequilibrium (LD) was detected for FOL, whereas the lowest level of LD was found for FOC, and 95.2%, 99.4%, 99.1% and 97.4% of the relative kinship estimates were less than 0.4 for FOL, FOMG, FORL and FOC, respectively. LD differences were detected among formae speciales, and LD was higher in FOL as compare to FOC species. The findings of this study confirm that iPBS‐retrotransposon markers are highly polymorphic at the intraspecific level in Fusarium spp.     

Reconstructing the population history of Nicaragua by means of mtDNA,Y‐chromosome STRs,and autosomal STR markers

Before the arrival of the Spaniards in Nicaragua, diverse Native American groups inhabited the territory. In colonial times, Native Nicaraguan populations interacted with Europeans and slaves from Africa. To ascertain the extent of this genetic admixture and provide genetic evidence about the origin of the Nicaraguan ancestors, we analyzed the mitochondrial control region (HVSI and HVSII), 17 Y chromosome STRs, and 15 autosomal STRs in 165 Mestizo individuals from Nicaragua. To carry out interpopulation comparisons, HVSI sequences from 29 American populations were compiled from the literature. The results reveal a close relationship between Oto‐manguean, Uto‐Aztecan, Mayan groups from Mexico, and a Chibchan group to Nicaraguan lineages. The Native American contribution to present‐day Nicaraguan Mestizos accounts for most of the maternal lineages, whereas the majority of Nicaraguan Y chromosome haplogroups can be traced back to a West Eurasian origin. Pairwise Fst distances based on Y‐STRs between Nicaragua and European, African and Native American populations show that Nicaragua is much closer to Europeans than the other populations. Additionally, admixture proportions based on autosomal STRs indicate a predominantly Spanish contribution. Our study reveals that the Nicaraguan Mestizo population harbors a high proportion of European male and Native American female substrate. Finally, the amount of African ancestry is also interesting, probably because of the contribution of Spanish conquerors with NorthAfrican genetic traces or that of West African slaves. Am J Phys Anthropol, 2010. © 2010 Wiley‐Liss, Inc.     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

Validation of whole genome linkage‐linkage disequilibrium and association results,and identification of markers to predict genetic merit for twinning

A previous genome‐wide search with a moderate density 10K marker set identified many marker associations with twinning rate, either through single‐marker analysis or combined linkage‐linkage disequilibrium (LLD; haplotype) analysis. The objective of the current study was to validate putative marker associations using an independent set of phenotypic data. Holstein bulls (n = 921) from 100 paternal half‐sib families were genotyped. Twinning rate predicted transmitting abilities were calculated using calving records from 1994 to 1998 (Data I) and 1999 to 2006 (Data II), and the underlying liability scores from threshold model analysis were used as the trait in marker association analyses. The previous analysis used 201 bulls with daughter records in Data I. In the current analysis, this was increased to 434, providing a revised estimate of effect and significance. Bulls with daughter records in Data II totaled 851, and analysis of this data provided the validation of results from analysis of Data I. Single nucleotide polymorphisms (SNPs) were selected to validate previously significant single‐marker associations and LLD results. Bulls were genotyped for a total of 306 markers. Nine of 13 LLD regions located on chromosomes 1, 2, 3, 6, 9, 22, 23(2) and 26 were validated, showing significant results for both Data I and II. Association analysis revealed 55 of 174 markers validated, equating to a single‐marker validation rate of 31%. Stepwise backward elimination and cross‐validation analyses identified 18 SNPs for use in a final reduced marker panel explaining 34% of the genetic variation, and to allow prediction of genetic merit for twinning rate.     

Population structure,seasonal genotypic differentiation,and clonal diversity of weedy dandelions in three Boston area populations (Taraxacum sp.)

Weedy dandelions have a worldwide distribution and thrive in urban environments despite a lack of sexual reproduction throughout most of its range. North American dandelions, introduced from Eurasia, are believed to be primarily, if not exclusively, apomictic triploids. In some European populations, apomicts co‐occur with diploid sexual individuals and hybridizations can create genetically unique apomicts, which may subsequently disperse and establish new populations globally. Using six nuclear microsatellite markers and a cpDNA intergenic spacer, we investigate the impact of this unusual natural history on population structure and diversity in three urban Boston area dandelion populations. Our results show high levels of genetic diversity within populations, spatial population structure, and seasonal genotypic differentiation in flowering times. We find evidence that sexual reproduction and recombination, presumably in Europe, and extensive gene flow drive these patterns of diversity and create the appearance of panmixia despite the lack of evidence for local sexual reproduction.     

Autosomal genetic diversity in non‐breed horses from eastern Eurasia provides insights into historical population movements

Many events in the history of eastern Eurasia, including the process of domestication itself, the initial spread of domestic horses and subsequent movements, are believed to have affected the genetic structure of domestic horse populations in this area. We investigated levels of within‐ and between‐population genetic diversity in ‘non‐breed horses’ (working horses sampled in remote areas) from 17 locations in Asia and parts of Eastern Europe, using 26 autosomal microsatellite loci. Non‐breed horses have not been subject to the same intensity of artificial selection and closed breeding as have most breed animals and are thus expected to better reflect the population history of domestic horses. Despite geographic distances of between 300 and 7000 km between sampling locations, pairwise F _ST was very low (range: <0.001 to ?0.033), suggesting historically high levels of gene flow. Our analyses of non‐breed horses revealed a pattern of isolation by distance and a significant decline in genetic diversity (expected heterozygosity and allelic richness) from east to west, consistent with a westward expansion of horses out of East Asia. Although the timing of this putative expansion is unclear, our results highlight the benefit of studying animals that do not belong to particular breeds when investigating aspects of a population's history.     

The non‐hierarchical,non‐uniformly branching topology of a leuconoid sponge aquiferous system revealed by 3D reconstruction and morphometrics using corrosion casting and X‐ray microtomography

Hammel, J.U., Filatov, M.V., Herzen, J, Beckmann, F., Kaandorp, J.A. and Nickel M. 2011. The non‐hierarchical, non‐uniformly branching topology of a leuconoid sponge aquiferous system revealed by 3D reconstruction and morphometrics using corrosion casting and X‐ray microtomography. —Acta Zoologica (Stockholm) 00 :1–12. As sessile filter feeders, sponges rely on a highly efficient fluid transport system. Their physiology depends on efficient water exchange, which is performed by the aquiferous system. This prominent poriferan anatomical character represents a dense network of incurrent and excurrent canals on which we lack detailed 3D models. To overcome this, we investigated the complex leucon‐type architecture in the demosponge Tethya wilhelma using corrosion casting, microtomography, and 3D reconstructions. Our integrative qualitative and quantitative approach allowed us to create, for the first time, high‐resolution 3D representations of entire canal systems which were used for detailed geometric and morphometric measurements. Canal diameters lack distinct size classes, and bifurcations are non‐uniformly ramified. A relatively high number of bifurcations show previously unknown and atypical cross‐sectional area ratios. Scaling properties and topological patterns of the canals indicate a more complex overall architecture than previously assumed. As a consequence, it might be more convenient to group canals into functional units rather than hierarchical clusters. Our data qualify the leucon canal system architecture of T. wilhelma as a highly efficient fluid transport system adapted toward minimal flow resistance. Our results and approach are relevant for a better understanding of sponge biology and cultivation techniques.     

Inference of sex‐specific expansion patterns in human populations from Y‐chromosome polymorphism

Studying the current distribution of genetic diversity in humans has important implications for our understanding of the history of our species. We analyzed a set of linked STR and SNP loci from the paternally inherited Y chromosome to infer the past demography of 55 African and Eurasian populations, using both the parametric and nonparametric coalescent‐based methods implemented in the BEAST application. We inferred expansion events in most sedentary farmer populations, while we found constant effective population sizes for both nomadic hunter‐gatherers and seminomadic herders. Our results differed, on several aspects, from previous results on mtDNA and autosomal markers. First, we found more recent expansion patterns in Eurasia than in Africa. This discrepancy, substantially stronger than the ones found with the other kind of markers, may result from a lower effective population size for men, which might have made male‐transmitted markers more sensitive to the out‐of‐Africa bottleneck. Second, we found expansion signals only for sedentary farmers but not for nomadic herders in Central Asia, while these signals were found for both kind of populations in this area when using mtDNA or autosomal markers. Expansion signals in this area may result from spatial expansion processes and may have been erased for the Y chromosome among the herders because of restricted male gene flow. Am J Phys Anthropol 157:217–225, 2015. © 2015 Wiley Periodicals, Inc.     

speed‐ne: Software to simulate and estimate genetic effective population size (Ne) from linkage disequilibrium observed in single samples

The genetic effective population size, N_e, can be estimated from the average gametic disequilibrium () between pairs of loci, but such estimates require evaluation of assumptions and currently have few methods to estimate confidence intervals. speed‐ne is a suite of matlab computer code functions to estimate from with a graphical user interface and a rich set of outputs that aid in understanding data patterns and comparing multiple estimators. speed‐ne includes functions to either generate or input simulated genotype data to facilitate comparative studies of estimators under various population genetic scenarios. speed‐ne was validated with data simulated under both time‐forward and time‐backward coalescent models of genetic drift. Three classes of estimators were compared with simulated data to examine several general questions: what are the impacts of microsatellite null alleles on , how should missing data be treated, and does disequilibrium contributed by reduced recombination among some loci in a sample impact . Estimators differed greatly in precision in the scenarios examined, and a widely employed estimator exhibited the largest variances among replicate data sets. speed‐ne implements several jackknife approaches to estimate confidence intervals, and simulated data showed that jackknifing over loci and jackknifing over individuals provided ~95% confidence interval coverage for some estimators and should be useful for empirical studies. speed‐ne provides an open‐source extensible tool for estimation of from empirical genotype data and to conduct simulations of both microsatellite and single nucleotide polymorphism (SNP) data types to develop expectations and to compare estimators.     

Population structure in the highly fragmented range of Tor douronensis (Cyprinidae) in Sarawak, Malaysia revealed by microsatellite DNA markers

1. Understanding the structure of threatened populations, particularly those that exist in degraded or fragmented habitats is crucial for their effective management and conservation. Recently developed methods of individual-based analysis of genetic data provide an unprecedented opportunity to understand the relationships amongst fragmented populations.
2. In the present study, population structure of an important cyprinid species ( Tor douronensis ), which is indigenous to Sarawak, Malaysia, is investigated as part of an ongoing conservation effort to restore threatened wild populations of the species. The population structure inferred using data from seven autosomal microsatellite loci was generally consistent with geography and habitat fragmentation.
3. The results indicate that there are two well-defined clusters of T. douronensis in Sarawak, namely the 'northeastern' and the 'southwestern' clusters. In addition, a further subdivision was observed in each of the clusters distributed between river systems. Low levels of gene flow were also observed and migrants between habitat fragments were identified, possibly resulting from human-mediated translocations.
4. Implications of the findings for management and conservation of T. douronensis are discussed.     

Genome‐wide variation in nucleotides and retrotransposons in alpine populations of Arabis alpina (Brassicaceae)

Advances in high‐throughput sequencing have promoted the collection of reference genomes and genome‐wide diversity. However, the assessment of genomic variation among populations has hitherto mainly been surveyed through single‐nucleotide polymorphisms (SNPs) and largely ignored the often major fraction of genomes represented by transposable elements (TEs). Despite accumulating evidence supporting the evolutionary significance of TEs, comprehensive surveys remain scarce. Here, we sequenced the full genomes of 304 individuals of Arabis alpina sampled from four nearby natural populations to genotype SNPs as well as polymorphic long terminal repeat retrotransposons (polymorphic TEs; i.e., presence/absence of TE insertions at specific loci). We identified 291,396 SNPs and 20,548 polymorphic TEs, comparing their contributions to genomic diversity and divergence across populations. Few SNPs were shared among populations and overall showed high population‐specific variation, whereas most polymorphic TEs segregated among populations. The genomic context of these two classes of variants further highlighted candidate adaptive loci having a putative impact on functional genes. In particular, 4.96% of the SNPs were identified as nonsynonymous or affecting start/stop codons. In contrast, 43% of the polymorphic TEs were present next to Arabis genes enriched in functional categories related to the regulation of reproduction and responses to biotic as well as abiotic stresses. This unprecedented data set, mapping variation gained from SNPs and complementary polymorphic TEs within and among populations, will serve as a rich resource for addressing microevolutionary processes shaping genome variation.     

Hidden genetic variance contributes to increase the short‐term adaptive potential of selfing populations

Standing genetic variation is considered a major contributor to the adaptive potential of species. The low heritable genetic variation observed in self‐fertilizing populations has led to the hypothesis that species with this mating system would be less likely to adapt. However, a non‐negligible amount of cryptic genetic variation for polygenic traits, accumulated through negative linkage disequilibrium, could prove to be an important source of standing variation in self‐fertilizing species. To test this hypothesis, we simulated populations under stabilizing selection subjected to an environmental change. We demonstrate that, when the mutation rate is high (but realistic), selfing populations are better able to store genetic variance than outcrossing populations through genetic associations, notably due to the reduced effective recombination rate associated with predominant selfing. Following an environmental shift, this diversity can be partially remobilized, which increases the additive variance and adaptive potential of predominantly (but not completely) selfing populations. In such conditions, despite initially lower observed genetic variance, selfing populations adapt as readily as outcrossing ones within a few generations. For low mutation rates, purifying selection impedes the storage of diversity through genetic associations, in which case, as previously predicted, the lower genetic variance of selfing populations results in lower adaptability compared to their outcrossing counterparts. The population size and the mutation rate are the main parameters to consider, as they are the best predictors of the amount of stored diversity in selfing populations. Our results and their impact on our knowledge of adaptation under high selfing rates are discussed.     

Analysis of genetic structure in soil populations of Rhizobium leguminosarum recovered from the USA and the UK

Although many studies have shown that animal-associated bacterial species exhibit linkage disequilibrium at chromosomal loci, recent studies indicate that both animal-associated and soil-borne bacterial species can display a nonclonal genetic structure in which alleles at chromosomal loci are in linkage equilibrium. To examine the situation in soil-borne species further, we compared genetic structure in two soil populations of Rhizobium leguminosarum bv. trifolii and two populations of R. leguminosarum bv. viciae from two sites in Oregon, with genetic structure in R. leguminosarum bv. viciae populations recovered from peas grown at a site in Washington, USA, and at a site in Norfolk, UK. A total of 234 chromosomal types (ET) were identified among 682 strains analysed for allelic variation at 13 enzyme-encoding chromosomal loci by multilocus enzyme electrophoresis (MLEE). Chi-square tests for heterogeneity of allele frequencies showed that the populations were not genetically uniform. A comparison of the genetic diversity within combined and individual populations confirmed that the Washington population was the primary cause of genetic differentiation between the populations. Each individual population exhibited linkage disequilibrium, with the magnitude of the disequilibrium being greatest in the Washington population and least in the UK population of R. leguminosarum bv. viciae. Linkage disequilibrium in the UK population was created between two clusters of 9 and 23 ETs, which, individually, were in linkage equilibrium. Strong linkage disequilibrium between the two major clusters of 8 and 12 ETs in the Washington population was caused by the low genetic diversity of the ETs within each cluster relative to the inter-cluster genetic distance. Because neither the magnitude of genetic diversity nor of linkage disequilibrium increased as hierarchical combinations of the six local populations were analysed, we conclude that the populations have not been isolated from each other for sufficient time, nor have they been exposed to enough selective pressure to develop unique multilocus genetic structure.     

Comparison of SNPs and microsatellites for assessing the genetic structure of chicken populations

Many studies in human genetics compare informativeness of single‐nucleotide polymorphisms (SNPs) and microsatellites (single sequence repeats; SSR) in genome scans, but it is difficult to transfer the results directly to livestock because of different population structures. The aim of this study was to determine the number of SNPs needed to obtain the same differentiation power as with a given standard set of microsatellites. Eight chicken breeds were genotyped for 29 SSRs and 9216 SNPs. After filtering, only 2931 SNPs remained. The differentiation power was evaluated using two methods: partitioning of the Euclidean distance matrix based on a principal component analysis (PCA) and a Bayesian model‐based clustering approach. Generally, with PCA‐based partitioning, 70 SNPs provide a comparable resolution to 29 SSRs. In model‐based clustering, the similarity coefficient showed significantly higher values between repeated runs for SNPs compared to SSRs. For the membership coefficients, reflecting the proportion to which a fraction segment of the genome belongs to the ith cluster, the highest values were obtained for 29 SSRs and 100 SNPs respectively. With a low number of loci (29 SSRs or ≤100 SNPs), neither marker types could detect the admixture in the Gödöllö Nhx population. Using more than 250 SNPs allowed a more detailed insight into the genetic architecture. Thus, the admixed population could be detected. It is concluded that breed differentiation studies will substantially gain power even with moderate numbers of SNPs.     

Spatio‐temporal genetic structure and the effects of long‐term fishing in two partially sympatric offshore demersal fishes

Environmental gradients have been shown to disrupt gene flow in marine species, yet their influence in structuring populations at depth remains poorly understood. The Cape hakes (Merluccius paradoxus and M. capensis) are demersal species co‐occurring in the Benguela Current system, where decades of intense fishing resulted in severely depleted stocks in the past. Previous studies identified conflicting mtDNA genetic substructuring patterns and thus contrasting evolutionary trajectories for both species. Using 10 microsatellite loci, the control region of mtDNA and employing a seascape genetics approach, we investigated genetic connectivity and the impact of prolonged exploitation in the two species, which are characterized by different patterns of fishing pressure. Three consecutive years were sampled covering the entire distribution (N = 2100 fishes). Despite large estimated population sizes, both species exhibited low levels of contemporary genetic diversity (0.581 < H_E < 0.692), implying that fishing has had a significant impact on their genetic composition and evolutionary trajectories. Further, for M. paradoxus, significant temporal, but not spatial, divergence points to the presence of genetic chaotic patchiness. In contrast, M. capensis exhibited a clear latitudinal cline in genetic differentiation between Namibia and South Africa (F_ST = 0.063, P < 0.05), with low (0.2% per generation) estimates of contemporary gene flow. Seascape analyses reveal an association with bathymetry and upwelling events, suggesting that adaptation to local environmental conditions may drive genetic differentiation in M. capensis. Importantly, our results highlight the need for temporal sampling in disentangling the complex factors that impact population divergence in marine fishes.