Replicated analysis of the genetic architecture of quantitative traits in two wild great tit populations

Currently, there is much debate on the genetic architecture of quantitative traits in wild populations. Is trait variation influenced by many genes of small effect or by a few genes of major effect? Where is additive genetic variation located in the genome? Do the same loci cause similar phenotypic variation in different populations? Great tits (Parus major) have been studied extensively in long‐term studies across Europe and consequently are considered an ecological ‘model organism’. Recently, genomic resources have been developed for the great tit, including a custom SNP chip and genetic linkage map. In this study, we used a suite of approaches to investigate the genetic architecture of eight quantitative traits in two long‐term study populations of great tits—one in the Netherlands and the other in the United Kingdom. Overall, we found little evidence for the presence of genes of large effects in either population. Instead, traits appeared to be influenced by many genes of small effect, with conservative estimates of the number of contributing loci ranging from 31 to 310. Despite concordance between population‐specific heritabilities, we found no evidence for the presence of loci having similar effects in both populations. While population‐specific genetic architectures are possible, an undetected shared architecture cannot be rejected because of limited power to map loci of small and moderate effects. This study is one of few examples of genetic architecture analysis in replicated wild populations and highlights some of the challenges and limitations researchers will face when attempting similar molecular quantitative genetic studies in free‐living populations.     

Recombination rates across porcine autosomes inferred from high-density linkage maps

Studies of the variation in recombination rate across the genome provide a better understanding of evolutionary genomics and are also an important step towards mapping and dissecting complex traits in domestic animals. With the recent completion of the porcine genome sequence and the availability of a high‐density porcine single nucleotide polymorphism (SNP) array, it is now possible to construct a high‐density porcine linkage map and estimate recombination rate across the genome. A total of 416 animals were genotyped with the Porcine SNP60BeadChip, and high‐density chromosome linkage maps were constructed using CRI‐MAP, assuming the physical order of the Sscrofa10 assembly. The total linkage map length was 2018.79 cM, using 658 meioses and 14 503 SNPs. The estimated average recombination rate across the porcine autosomes was 0.86 cM/Mb. However, a large variation in recombination rate was observed among chromosomes. The estimated average recombination rates (cM/Mb) per chromosome ranged from 0.48 in SSC1 to 1.48 in SSC10, displaying a significant negative correlation with the chromosome sizes. In addition, the analysis of the variation in the recombination rates taking 1‐Mb sliding windows has allowed us to demonstrate the variation in recombination rates within chromosomes. In general, a larger recombination rate was observed in the extremes than in the centre of the chromosome. Finally, the ratio between female and male recombination rates was also inferred, obtaining a value of 1.38, with the heterogametic sex having the least recombination.     

A microsatellite linkage map for Drosophila montana shows large variation in recombination rates,and a courtship song trait maps to an area of low recombination

Current advances in genetic analysis are opening up our knowledge of the genetics of species differences, but challenges remain, particularly for out‐bred natural populations. We constructed a microsatellite‐based linkage map for two out‐bred lines of Drosophila montana derived from divergent populations by taking advantage of the Drosophila virilis genome and available cytological maps of both species. Although the placement of markers was quite consistent with cytological predictions, the map indicated large heterogeneity in recombination rates along chromosomes. We also performed a quantitative trait locus (QTL) analysis on a courtship song character (carrier frequency), which differs between populations and is subject to strong sexual selection. Linkage mapping yielded two significant QTLs, which explained 3% and 14% of the variation in carrier frequency, respectively. Interestingly, as in other recent studies of traits which can influence speciation, the strongest QTL mapped to a genomic region partly covered by an inversion polymorphism.     

A linkage map of the porcine genome from a large-scale White Duroc × Erhualian resource population and evaluation of factors affecting recombination rates

A porcine genome linkage map composed of 194 microsatellite markers was constructed with a large-scale White Duroc × Erhualian resource population. The marker order on this linkage map was consistent with the USDA-MARC reference map except for two markers on SSC3, two markers on SSC13 and two markers on SSCX. The length of the sex-averaged map (2344.9 cM) was nearly the same as that of the USDA-MARC and NIAI map. Highly significant heterogeneity in recombination rates between sexes was observed. Except for SSC1 and SSC13, the female autosomes had higher average recombination rates than the male autosomes. Moreover, recombination rates in the pseudoautosomal region were greater in males than in females. These observations are consistent with those of previous reports. The recombination rates on each paternal and maternal chromosome of F₂ animals were calculated. Recombination rates were not significantly affected by the age (in days) or parity of the F₁ animals. However, recombination rates on paternal chromosomes were affected by the mating season of the F₁ animals. This could represent an effect of environmental temperature on spermatogenesis.     

Spatial genetic analyses reveal strong genetic structure in two populations of the outcrossing tree fern Alsophila firma (Cyatheaceae)

The development of spatial genetic structure (SGS) in seed plants has been linked to several biological attributes of species, such as breeding system and life form. However, little is known about SGS in ferns, which together with lycopods are unique among land plants in having two free‐living life stages. We combined spatial aggregation statistics and spatial genetic autocorrelation analyses using five plastid microsatellites and one nuclear gene to investigate SGS in two populations of the outcrossing tree fern Alsophila firma (Cyatheaceae). We assessed how the observed patterns compare with those estimated for other ferns and seed plants. Populations of A. firma exhibited strong SGS, spatial clustering of individuals, substantial clonal diversity and no inbreeding. SGS in ferns appears to be higher than in most seed plants analysed to date. Contrary to our expectations, an outcrossing breeding system, wind dispersal and an arborescent life form did not translate into weak or no SGS. In ferns, SGS is probably being affected by the life cycle with two free‐living life stages. The reproductive biology of ferns appears to be more complex than previously thought. This implies that SGS in ferns is affected by some factors that cannot be inferred from the study of flowering plants. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 177 , 439–449.     

Construction of a high-density genetic linkage map and QTL analysis of morphological traits in an F1 Malusdomestica × Malus baccata hybrid

Apple is considered the most commonly grown fruit crop in temperate regions that brings great economic profits to fruit growers. Dwarfing rootstocks have been extensively used in apple breeding as well as commercial orchards, but the molecular and genetic basis of scion dwarfing and other morphological traits induced by them is still unclear. At present, we report a genetic map of Malusdomestica × Malus baccata with high density. The F₁ population was sequenced by a specific length amplified fragment (SLAF). In the genetic map, 5064 SLAF markers spanning 17 linkage groups (LG) were included. Dwarf-related and other phenotypic traits of the scion were evaluated over a 3-year growth period. Based on quantitative trait loci (QTL) evaluation of plant height and trunk diameter, two QTL clusters were found on LG 11, which exhibited remarkable influences on dwarfing of the scion. In this analysis, QTL DW2, which was previously reported as a locus that controls dwarfing, was confirmed. Moreover, three novel QTLs for total flower number and branching flower number were detected on LG2 and LG4, exhibited the phenotypic variation that has been explained by QTL ranging from 8.80% to 34.80%. The findings of the present study are helpful to find scion dwarfing and other phenotypes induced by rootstock in the apple.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12298-021-01069-0.     

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.     

High-density physical maps reveal that the dominant male-sterile gene Ms3 is located in a genomic region of low recombination in wheat and is not amenable to map-based cloning

In wheat it is essential to know whether a gene is located in a high or low recombination region of the genome before initiating a map-based cloning approach. The objective of this study was to explore the potential feasibility of map-based cloning of the dominant male-sterile gene Ms3 of wheat. High-density physical maps of the short arms of the group-5 chromosomes (5AS, 5BS, and 5DS) of Triticum aestivum L. were constructed by mapping 40 DNA markers on a set of 17 homozygous deletion lines. One hundred RFLP loci were mapped: 35 on 5AS, 37 on 5BS, and 28 on 5DS. A consensus physical map was colinearly aligned with a consensus genetic map of the group-5 short arms. Sixteen of the 17 markers in the consensus genetic map encompass a genetic distance of 25 cM and correspond to the distal region (FL 0.56–0.97) of the consensus physical map. Two rice probes, RG463 and RG901, previously identified to be linked to markers CDO344 and CDO749 (group-5 short arm of wheat), respectively, in the genetic map of rice chromosome 12, map between FL 0.56 and 0.63 in the consensus map. Thus at least a part of the group-5 short arm is homoeologous to a region of chromosome 12 of rice. The genetic map of chromosome arm 5AS was constructed using a population of 139 BC₁ plants derived from a cross between the euploid wheat ”Chris” carrying a dominant male-sterile gene Ms3 and a disomic substitution line in which chromosome 5A of T. aestivum cv Chinese Spring was substituted by chromosome 5A from Triticum turgidum ssp. dicoccoides. The map has a genetic length of 53.4 cM with 11 DNA markers. The initial map showed that the gene Ms3 cosegregated with three markers, WG341, BCD1130 and CDO677. High-resolution mapping using an additional 509 BC₁ plants indicated that the marker WG341 was closely linked to Ms3 at a genetic distance of 0.8 cM. The Ms3 was mapped physically in the region spanning 40% of the arm length from the centromere of 5AS. Therefore, map-based cloning of the Ms3 is not feasible, although WG341 can be used as a useful tag for the Ms3 gene for breeding purposes. Received: 12 December 2000 / Accepted: 26 January 2001     

Linkage disequilibrium patterns and genetic structure of Amerindian and non‐Amerindian Brazilian populations revealed by long‐range X‐STR markers

The extent of X‐chromosome linkage disequilibrium (LD) was studied in a southern Brazilian population, and in a pool of samples from Amerindian populations. For this purpose, 11 microsatellites, located mostly in a Xq region comprising ～86 Mb was investigated. The lower Amerindian gene diversity associated with significant differences between the populations studied indicated population structure as the main cause for the higher LD values in the Amerindian pool. On the other hand, the LD levels of the non‐Amerindian Brazilian sample, although less extensive than that of the Amerindians, were probably determined by admixture events. Our results indicated that different demographic histories have significant effects on LD levels of human populations, and provide a first approach to the X‐chromosome ancestry of Amerindian and non‐Amerindian Brazilian populations, being valuable for future studies involving mapping and population genetic studies. Am J Phys Anthropol 2009. © 2009 Wiley‐Liss, Inc.     

Gene variation and genetic differentiation among populations of the solitary mud dauber wasp Trypoxylon (Trypargilum) albitarse Fabricius 1804 (Hymenoptera, Crabronidae)

Trypoxylon is a genus of solitary crabronid wasps whose populationgenetics is poorly known. The purpose of the present study was to investigate thegenetic variation and differentiation among five populations of Trypoxylonalbitarse, a species widely distributed throughout the Neotropics, withrecords from Panama to northern Argentina. Eight species-specific microsatellite lociwere used for genotyping 96 adult wasps (one female per nest) sampled at five sitesin Brazil. The analysis of allelic richness and private alleles indicated highgenetic diversity in the populations sampled. Pairwise comparisons using theF_st and D_est indices revealed significant differentiation for all, but one pair ofpopulations. F_st, D_est, AMOVA and assignment test values pointed to inter-population differentiation.Additionally, the analysis of population structure using Bayesian and PCA methodscharacterized two alternative genetic groups. The Mantel test indicated nocorrelation between genetic and geographic distances. Despite evidence ofconsiderable dispersal capacity for T. albitarse, the data indicatelow to moderate population structuring in this species.     

RFLP genetic linkage maps from four F(2.3) populations and a joinmap of Gossypium hirsutum L

An RFLP genetic linkage joinmap was constructed from four different mapping populations of cotton (Gossypium hirsutum L.). Genetic maps from two of the four populations have been previously reported. The third genetic map was constructed from 199 bulk-sampled plots of an F_2.3 (HQ95–6×’MD51ne’) population. The map comprises 83 loci mapped to 24 linkage groups with an average distance between markers of 10.0 centiMorgan (cM), covering 830.1 cM or approximately 18% of the genome. The fourth genetic map was developed from 155 bulk-sampled plots of an F_2.3 (119– 5 sub-okra×’MD51ne’) population. This map comprises 56 loci mapped to 16 linkage groups with an average distance between markers of 9.3 cM, covering 520.4 cM or approximately 11% of the cotton genome. A core of 104 cDNA probes was shared between populations, yielding 111 RFLP loci. The constructed genetic linkage joinmap from the above four populations comprises 284 loci mapped to 47 linkage groups with the average distance between markers of 5.3 cM, covering 1,502.6 cM or approximately 31% of the total recombinational length of the cotton genome. The linkage groups contained from 2 to 54 loci each and ranged in distance from 1.0 to 142.6 cM. The joinmap provided further knowledge of competitive chromosome arrangement, parental relationships, gene order, and increased the potential to map genes for the improvement of the cotton crop. This is the first genetic linkage joinmap assembled in G. hirsutum with a core of RFLP markers assayed on different genetic backgrounds of cotton populations (Acala, Delta, and Texas plain). Research is ongoing for the identification of quantitative trait loci for agronomic, physiological and fiber quality traits on these maps, and the identification of RFLP loci lineage for G. hirsutum from its diploid progenitors (the A and D genomes). Received: 23 February 2001 / Accepted: 8 June 2001     

Vital rates of two small populations of brown bears in Canada and range‐wide relationship between population size and trend

Identifying mechanisms of population change is fundamental for conserving small and declining populations and determining effective management strategies. Few studies, however, have measured the demographic components of population change for small populations of mammals (<50 individuals). We estimated vital rates and trends in two adjacent but genetically distinct, threatened brown bear (Ursus arctos) populations in British Columbia, Canada, following the cessation of hunting. One population had approximately 45 resident bears but had some genetic and geographic connectivity to neighboring populations, while the other population had <25 individuals and was isolated. We estimated population‐specific vital rates by monitoring survival and reproduction of telemetered female bears and their dependent offspring from 2005 to 2018. In the larger, connected population, independent female survival was 1.00 (95% CI: 0.96–1.00) and the survival of cubs in their first year was 0.85 (95% CI: 0.62–0.95). In the smaller, isolated population, independent female survival was 0.81 (95% CI: 0.64–0.93) and first‐year cub survival was 0.33 (95% CI: 0.11–0.67). Reproductive rates did not differ between populations. The large differences in age‐specific survival estimates resulted in a projected population increase in the larger population (λ = 1.09; 95% CI: 1.04–1.13) and population decrease in the smaller population (λ = 0.84; 95% CI: 0.72–0.95). Low female survival in the smaller population was the result of both continued human‐caused mortality and an unusually high rate of natural mortality. Low cub survival may have been due to inbreeding and the loss of genetic diversity common in small populations, or to limited resources. In a systematic literature review, we compared our population trend estimates with those reported for other small populations (<300 individuals) of brown bears. Results suggest that once brown bear populations become small and isolated, populations rarely increase and, even with intensive management, recovery remains challenging.     

Whole‐genome sequencing of two North American Drosophila melanogaster populations reveals genetic differentiation and positive selection

The prevailing demographic model for Drosophila melanogaster suggests that the colonization of North America occurred very recently from a subset of European flies that rapidly expanded across the continent. This model implies a sudden population growth and range expansion consistent with very low or no population subdivision. As flies adapt to new environments, local adaptation events may be expected. To describe demographic and selective events during North American colonization, we have generated a data set of 35 individual whole‐genome sequences from inbred lines of D. melanogaster from a west coast US population (Winters, California, USA) and compared them with a public genome data set from Raleigh (Raleigh, North Carolina, USA). We analysed nuclear and mitochondrial genomes and described levels of variation and divergence within and between these two North American D. melanogaster populations. Both populations exhibit negative values of Tajima's D across the genome, a common signature of demographic expansion. We also detected a low but significant level of genome‐wide differentiation between the two populations, as well as multiple allele surfing events, which can be the result of gene drift in local subpopulations on the edge of an expansion wave. In contrast to this genome‐wide pattern, we uncovered a 50‐kilobase segment in chromosome arm 3L that showed all the hallmarks of a soft selective sweep in both populations. A comparison of allele frequencies within this divergent region among six populations from three continents allowed us to cluster these populations in two differentiated groups, providing evidence for the action of natural selection on a global scale.     

Genetic experiments with model populations: Fallacies in genetic analysis performed from samples of recombinants selected for two markers

Abstract: Samples selected for two markers (one from each parent) from the progeny of a genetic cross are altered both in the genotype and in the frequency distribution as compared to the original progeny population. The consequences of these alterations were analyzed in selected samples obtained from model progeny populations of hypothetical genetic crosses. In complete progeny populations, distribution of genotypes (pattern of genotypes), sequence of markers, and frequency distribution of individual genotypes exhibit genetically intelligent relationship to each other only under a sole condition: at the correct sequence of markers. In selected samples, the singularity of the interrelationship is relaxed. Consequences of the relaxation in the genetic analysis are multiple, the most striking of which is an insinuation of circularity for the linear gene map.     

Two colonisation stages generate two different patterns of genetic diversity within native and invasive ranges of Ulex europaeus

Genetic diversity and the way a species is introduced influence the capacity of populations of invasive species to persist in, and adapt to, their new environment. The diversity of introduced populations affects their evolutionary potential, which is particularly important for species that have invaded a wide range of habitats and climates, such as European gorse, Ulex europaeus. This species originated in the Iberian peninsula and colonised Europe in the Neolithic; over the course of the past two centuries it was introduced to, and has become invasive in, other continents. We characterised neutral genetic diversity and its structure in the native range and in invaded regions. By coupling these results with historical data, we have identified the way in which gorse populations were introduced and the consequences of introduction history on genetic diversity. Our study is based on the genotyping of individuals from 18 populations at six microsatellite loci. As U. europaeus is an allohexaploid species, we used recently developed tools that take into account genotypic ambiguity. Our results show that genetic diversity in gorse is very high and mainly contained within populations. We confirm that colonisation occurred in two stages. During the first stage, gorse spread out naturally from Spain towards northern Europe, losing some genetic diversity. During the second stage, gorse was introduced by humans into different regions of the world, from northern Europe. These introductions resulted in the loss of rare alleles but did not significantly reduce genetic diversity and thus the evolutionary potential of this invasive species.     

Using a reference population yardstick to calibrate and compare genetic diversity reported in different studies: an example from the brown bear

In species with large geographic ranges, genetic diversity of different populations may be well studied, but differences in loci and sample sizes can make the results of different studies difficult to compare. Yet, such comparisons are important for assessing the status of populations of conservation concern. We propose a simple approach of using a single well-studied reference population as a ‘yardstick'' to calibrate results of different studies to the same scale, enabling comparisons. We use a well-studied large carnivore, the brown bear (Ursus arctos), as a case study to demonstrate the approach. As a reference population, we genotyped 513 brown bears from Slovenia using 20 polymorphic microsatellite loci. We used this data set to calibrate and compare heterozygosity and allelic richness for 30 brown bear populations from 10 different studies across the global distribution of the species. The simplicity of the reference population approach makes it useful for other species, enabling comparisons of genetic diversity estimates between previously incompatible studies and improving our understanding of how genetic diversity is distributed throughout a species range.     

A genetic linkage map of Pacific white shrimp (<Emphasis Type="Italic">Litopenaeus vannamei</Emphasis>): sex-linked microsatellite markers and high recombination rates

Pacific white shrimp (Litopenaeus vannamei) is the leading species farmed in the Western Hemisphere and an economically important aquaculture species in China. In this project, a genetic linkage map was constructed using amplified fragment length polymorphism (AFLP) and microsatellite markers. One hundred and eight select AFLP primer combinations and 30 polymorphic microsatellite markers produced 2071 markers that were polymorphic in either of the parents and segregated in the progeny. Of these segregating markers, 319 were mapped to 45 linkage groups of the female framework map, covering a total of 4134.4 cM; and 267 markers were assigned to 45 linkage groups of the male map, covering a total of 3220.9 cM. High recombination rates were found in both parental maps. A sex-linked microsatellite marker was mapped on the female map with 6.6 cM to sex and a LOD of 17.8, two other microsatellite markers were also linked with both 8.6 cM to sex and LOD score of 14.3 and 16.4. The genetic maps presented here will serve as a basis for the construction of a high-resolution genetic map, quantitative trait loci (QTLs) detection, marker-assisted selection (MAS) and comparative genome mapping.     

Relevant genetic differentiation among Brazilian populations of Anastrepha fraterculus (Diptera,Tephritidae)

We used a population genetic approach to detect the presence of genetic diversity among six populations of Anastrepha fraterculus across Brazil. To this aim, we used Simple Sequence Repeat (SSR) markers, which may capture the presence of differentiative processes across the genome in distinct populations. Spatial analyses of molecular variance were used to identify groups of populations that are both genetically and geographically homogeneous while also being maximally differentiated from each other. The spatial analysis of genetic diversity indicates that the levels of diversity among the six populations vary significantly on an eco-geographical basis. Particularly, altitude seems to represent a differentiating adaptation, as the main genetic differentiation is detected between the two populations present at higher altitudes and the other four populations at sea level. The data, together with the outcomes from different cluster analyses, identify a genetic diversity pattern that overlaps with the distribution of the known morphotypes in the Brazilian area.     

Assessing performance of single‐sample molecular genetic methods to estimate effective population size: empirical evidence from the endangered Gochu Asturcelta pig breed

Estimating effective population size (N_e) using linkage disequilibrium (LD) information (N_e(LD)) has the operational advantage of using a single sample. However, N_e(LD) estimates assume discrete generations and its performance are constrained by demographic issues. However, such concerns have received little empirical attention so far. The pedigree of the endangered Gochu Asturcelta pig breed includes individuals classified into discrete filial generations and individuals with generations overlap. Up to 780 individuals were typed with a set of 17 microsatellites. Performance of N_e(LD) was compared with N_e estimates obtained using genealogical information, molecular coancestry (N_e(M)) and a temporal (two‐sample) method (N_e(JR)). Molecular‐based estimates of N_e exceeded those obtained using pedigree data. Estimates of N_e(LD) for filial generations F3 and F4 (17.0 and 17.3, respectively) were lower and steadier than those obtained using yearly or biannual samplings. N_e(LD) estimated for samples including generations overlap could only be compared with those obtained for the discrete filial generations when sampling span approached a generation interval and demographic correction for bias was applied. Single‐sample N_e(M) estimates were lower than their N_e(LD) counterparts. N_e(M) estimates are likely to partially reflect the number of founders rather than population size. In any case, estimates of LD and molecular coancestry tend to covary and, therefore, N_e(M) and N_e(LD) can hardly be considered independent. Demographically adjusted estimates of N_e(JR) and N_e(LD) took comparable values when: (1) the two samples used for the former were separated by one equivalent to discrete generations in the pedigree and (2) sampling span used for the latter approached a generation interval. Overall, the empirical evidence given in this study suggested that the advantage of using single‐sample methods to obtain molecular‐based estimates of N_e is not clear in operational terms. Estimates of N_e obtained using methods based in molecular information should be interpreted with caution.