Haplotype analysis of the mitochondrial DNA d‐loop region reveals the maternal origin and historical dynamics among the indigenous goat populations in east and west of the Democratic Republic of Congo

This study aimed at assessing haplotype diversity and population dynamics of three Congolese indigenous goat populations that included Kasai goat (KG), small goat (SG), and dwarf goat (DG) of the Democratic Republic of Congo (DRC). The 1169 bp d‐loop region of mitochondrial DNA (mtDNA) was sequenced for 339 Congolese indigenous goats. The total length of sequences was used to generate the haplotypes and evaluate their diversities, whereas the hypervariable region (HVI, 453 bp) was analyzed to define the maternal variation and the demographic dynamic. A total of 568 segregating sites that generated 192 haplotypes were observed from the entire d‐loop region (1169 bp d‐loop). Phylogenetic analyses using reference haplotypes from the six globally defined goat mtDNA haplogroups showed that all the three Congolese indigenous goat populations studied clustered into the dominant haplogroup A, as revealed by the neighbor‐joining (NJ) tree and median‐joining (MJ) network. Nine haplotypes were shared between the studied goats and goat populations from Pakistan (1 haplotype), Kenya, Ethiopia and Algeria (1 haplotype), Zimbabwe (1 haplotype), Cameroon (3 haplotypes), and Mozambique (3 haplotypes). The population pairwise analysis (F_ST ) indicated a weak differentiation between the Congolese indigenous goat populations. Negative and significant (p‐value <.05) values for Fu''s Fs (−20.418) and Tajima''s (−2.189) tests showed the expansion in the history of the three Congolese indigenous goat populations. These results suggest a weak differentiation and a single maternal origin for the studied goats. This information will contribute to the improvement of the management strategies and long‐term conservation of indigenous goats in DRC.     

Genome‐wide selective detection of Mile red‐bone goat using next‐generation sequencing technology

The ecotype population of goats (Capra hircus) was created by long‐term artificial selection and natural adaptation. Mile red‐bone goat is an indigenous breed with visible red bones, and its special bone structure has received extensive attention. This study aimed to identify genetic variants and candidate genes associated with specific bone phenotypes using next‐generation sequencing technology (NGS). The results revealed that 31,828,206 single nucleotide polymorphisms (SNPs) were obtained from 72 goats (20 Mile red‐bone goats and 52 common goats) by NGS. A total of 100 candidate genes were identified on the basis top 1% window interaction from nucleotide diversity (π), π ratio (π _A/π _B), and pairwise fixation index (F _ST). Exactly 77 known signaling pathways were enriched. Specifically, three coding genes (NMNAT2, LOC102172983, and PNLIP) were annotated in the vitamin metabolism signaling pathways, and NCF2 was annotated to the osteoclast (OC) differentiation pathway. Furthermore, 5862 reliable copy number variations (CNVs) were obtained, and 14 and 24 genes were annotated with the top 1‰ CNV based on F _ST (>0.490) and V _ST (>0.527), respectively. Several pathways related to bone development and metabolism of exogenous substances in vivo, including calcium signaling pathway, OC differentiation, and glycerophospholipid metabolism, were annotated. Specifically, six genes from 19 candidate CNVs, which were obtained by interaction of the top 1‰ CNVs with F _ST and V _ST, were annotated to mucin‐type O‐glycan biosynthesis and metabolic pathways. Briefly, the results implied that pseudopurpurin and specific genetic variants work together to contribute to the red‐bone color and specific bone structure of Mile red‐bone goat. This study is helpful to understanding the genetic basis of the unique bone phenotype of Mile red‐bone goats.     

The role of anthropogenic dispersal in shaping the distribution and genetic composition of a widespread North American tree species

Dispersal and colonization are among the most important ecological processes for species persistence as they allow species to track changing environmental conditions. During the last glacial maximum (LGM), many cold‐intolerant Northern Hemisphere plants retreated to southern glacial refugia. During subsequent warming periods, these species expanded their ranges northward. Interestingly, some tree species with limited seed dispersal migrated considerable distances after the LGM ~19,000 years before present (YBP). It has been hypothesized that indigenous peoples may have dispersed valued species, in some cases beyond the southern limits of the Laurentide Ice Sheet. To investigate this question, we employed a molecular genetics approach on a widespread North American understory tree species whose fruit was valued by indigenous peoples. Twenty putative anthropogenic (near pre‐Columbian habitations) and 62 wild populations of Asimina triloba (pawpaw), which produces the largest edible fruit of any North American tree, were genetically assayed with nine microsatellite loci. Putative anthropogenic populations were characterized by reduced genetic diversity and greater excess heterozygosity relative to wild populations. Anthropogenic populations in regions that were glaciated during the LGM had profiles consistent with founder effects and reduced gene flow, and shared rare alleles with wild populations hundreds of kilometers away (mean = 723 km). Some of the most compelling evidence for human‐mediated dispersal is that putative anthropogenic and wild populations sharing rare alleles were separated by significantly greater distances (mean = 695 km) than wild populations sharing rare alleles (mean = 607 km; p = .014). Collectively, the genetic data suggest that long‐distance dispersal played an important role in the distribution of pawpaw and is consistent with the hypothesized role of indigenous peoples.     

The assembly of caprine Y chromosome sequence reveals a unique paternal phylogenetic pattern and improves our understanding of the origin of domestic goat

The mammalian Y chromosome offers a unique perspective on the male reproduction and paternal evolutionary histories. However, further understanding of the Y chromosome biology for most mammals is hindered by the lack of a Y chromosome assembly. This study presents an integrated in silico strategy for identifying and assembling the goat Y‐linked scaffolds using existing data. A total of 11.5 Mb Y‐linked sequences were clustered into 33 scaffolds, and 187 protein‐coding genes were annotated. We also identified high abundance of repetitive elements. A 5.84 Mb subset was further ordered into an assembly with the evidence from the goat radiation hybrid map (RH map). The existing whole‐genome resequencing data of 96 goats (worldwide distribution) were utilized to exploit the paternal relationships among bezoars and domestic goats. Goat paternal lineages were clearly divided into two clades (Y1 and Y2), predating the goat domestication. Demographic history analyses indicated that maternal lineages experienced a bottleneck effect around 2,000 YBP (years before present), after which goats belonging to the A haplogroup spread worldwide from the Near East. As opposed to this, paternal lineages experienced a population decline around the 10,000 YBP. The evidence from the Y chromosome suggests that male goats were not affected by the A haplogroup worldwide transmission, which implies sexually unbalanced contribution to the goat trade and population expansion in post‐Neolithic period.     

Plasticity and artificial selection for developmental mode in a poecilogonous sea slug

The contribution of phenotypically plastic traits to evolution depends on the degree of environmental influence on the target of selection (the phenotype) as well as the underlying genetic structure of the trait and plastic response. Likewise, maternal effects can help or hinder evolution through affects to the response to selection. The sacoglossan sea slug Alderia willowi exhibits intraspecific variation for developmental mode (= poecilogony) that is environmentally modulated with populations producing more yolk‐feeding (lecithotrophic) larvae during the summer, and more planktonic‐feeding (planktotrophic) larvae in the winter. I found significant family‐level variation in the reaction norms between 17 maternal families of A. willowi when reared in a split‐brood design in low (16 ppt) versus high (32 ppt) salinity, conditions which mimic seasonal variation in salinity of natural populations. I documented a significant response to selection for lecithotrophic larvae in high and low salinity. The slope of the reaction norm was maintained following one generation of selection for lecithotrophy. When the maternal environment was controlled in the laboratory, I found significant maternal effects, which reduced the response to selection. These results suggest there is standing genetic variation for egg‐mass type in A. willowi, but the ability of selection to act on that variation may depend on the environment in which the phenotype is expressed in preceding generations.     

Phylogeography and genetic diversity of the widespread katydid Ducetia japonica (Thunberg, 1815) across China

Habitat fragmentation can lower migration rates and genetic connectivity among remaining populations of native species. Ducetia japonica is one of the most widespread katydids in China, but little is known about its genetic structure and phylogeographic distribution. We combined the five‐prime region of cytochrome c oxidase subunit I (COI‐5P), 11 newly developed microsatellite loci coupled with an ecological niche model (ENM) to examine the genetic diversity and population structure of D. japonica in China and beyond to Laos and Singapore. Both Bayesian inference (BI) and haplotype network methods revealed six mitochondrial COI‐5P lineages. The distribution of COI‐5P haplotypes may not demonstrate significant phylogeographic structure (N _ST > G _ST, p > .05). The STRUCTURE analysis based on microsatellite data also revealed six genetic clusters, but discordant with those obtained from COI‐5P haplotypes. For both COI‐5P and microsatellite data, Mantel tests revealed a significant positive correlation between geographic and genetic distances in mainland China. Bayesian skyline plot (BSP) analyses indicated that the population size of D. japonica''s three major mitochondrial COI‐5P lineages were seemingly not affected by last glacial maximum (LGM, 0.015–0.025 Mya). The ecological niche models showed that the current distribution of D. japonica was similar to the species’ distribution during the LGM period and only slightly extended in northern China. Further phylogeographic studies based on more extensive sampling are needed to identify specific locations of glacial refugia in northern China.     

Phylogeography of sugar kelp: Northern ice‐age refugia in the Gulf of Alaska

Many Northeast (NE) Pacific fishes and invertebrates survived Pleistocene glaciations in northern refugia, but the extent that kelps survived in northern areas is uncertain. Here, we test the hypothesis that populations of sugar kelp (Saccharina latissima) persisted in the Gulf of Alaska during ice‐age maxima when the western margin of the Cordilleran ice sheet covered coastal areas around the NE Pacific Ocean. We estimated genetic diversities within and phylogeographical relationships among 14 populations along 2,800 km in the NE Pacific and Bering Sea with partial sequences of mitochondrial DNA 5′‐cytochrome oxidase subunit I (COI, bp = 624, n = 543), chloroplast DNA ribulose‐1,5‐bisphosphate carboxylase large subunit‐3′ (rbcL, bp = 735, n = 514), and 11 microsatellite loci. Concatenated sequences of rbcL and COI showed moderate levels of within‐population genetic diversity (mean h = 0.200) but substantial differences among populations (Φ_ST = 0.834, p < .0001). Microsatellites showed moderate levels of heterozygosity within populations (mean H _E = 0.391). Kelps in the same organellar lineage tended to cluster together, regardless of geographic origins, as indicated in a principal coordinate analysis (PCoA) of microsatellite genotypes. The PCoA also showed evidence of nuclear hybridizations between co‐occurring organellar lineages. Individual admixture plots with population clusters of K = 2, 6, and 9 showed increasing complexity with considerable historical admixture between some clusters. A time‐calibrated phylogeny placed divergences between rbcL‐COI lineages at 1.4 million years at most. The time frames of mutation in the rbcL‐COI lineages and microsatellite population clusters differed among locations. The existence of ancient lineages in the Gulf of Alaska, moderate levels of genetic diversity, and the absence of departures from neutrality are consistent with northern refugia during multiple Croll‐Milankovitch climate cycles in the Pleistocene Epoch.     

Genetic Variations of 13 Indigenous Chinese Goat Breeds Based on Cytochrome b Gene Sequences

Phylogenetic relationships among and genetic variability within 13 Chinese indigenous goat breeds and Boer goat were analyzed using cytochrome b gene sequences. There were 44 variable sites found in a 642 bp sequence, and 46 Cyt b haplotypes were subsequently defined. The phylogeny analysis of haplotypes in combination with goat Cyt b sequences from GenBank shows that Chinese goats are obviously separated from wild goats and might come from Capra aegagrus. Further analysis indicated that indigenous Chinese goats might descend from at least two lineages; most of the individuals analyzed could be classified into lineage A as defined by Luikart, but five other goats were of uncertain lineage. The Tibet plateau is a possible place of origin for Chinese goats. The neighbor-joining tree based on pairwise differences among populations shows that most Tibetan goats, except the Middle Tibet type, cluster closely with North China goats, and then with South China goats. This result confirms that differences in genetic structure exist among goats in different geographic locations. Nucleotide diversity varied among populations. Tibet and North China goats had higher genetic diversity than South China goats. The fixation index (F _st=87.72%) suggested that most of the total genetic variation was due to variation within populations. In addition, the results indicate that Cyt b gene sequence information alone might not be enough for phylogeny analysis among breeds within species, as shown by fewer polymorphic sites and lower bootstrap values on the neighbor-joining tree.     

Genetic differentiation and connectivity of morphological types of the broadcast‐spawning coral Galaxea fascicularis in the Nansei Islands,Japan

Population connectivity resulting from larval dispersal is essential for the maintenance or recovery of populations in marine ecosystems, including coral reefs. Studies of species diversity and genetic connectivity within species are essential for the conservation of corals and coral reef ecosystems. We analyzed mitochondrial DNA sequence types and microsatellite genotypes of the broadcast‐spawning coral, Galaxea fascicularis, from four regions in the subtropical Nansei Islands in the northwestern Pacific Ocean. Two types (soft and hard types) of nematocyst morphology are known in G. fascicularis and are significantly correlated with the length of a mitochondrial DNA noncoding sequence (soft type: mt‐L; hard type: mt‐S type). Using microsatellites, significant genetic differentiation was detected between the mitochondrial DNA sequence types in all regions. We also found a third genetic cluster (mt‐L+), and this unexpected type may be a cryptic species of Galaxea. High clonal diversity was detected in both mt‐L and mt‐S types. Significant genetic differentiation, which was found among regions within a given type (F _ST = 0.009–0.024, all Ps ≤ 0.005 in mt‐L; 0.009–0.032, all Ps ≤ 0.01 in mt‐S), may result from the shorter larval development than in other broadcast‐spawning corals, such as the genus Acropora. Nevertheless, intraspecific genetic diversity and connectivity have been maintained, and with both sexual and asexual reproduction, this species appears to have a potential for the recovery of populations after disturbance.     

RAD sequencing of common whelk,Buccinum undatum,reveals fine‐scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F _ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation‐by‐distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome‐wide association studies.     

How rapidly do self‐compatible populations evolve selfing? Mating system estimation within recently evolved self‐compatible populations of Azorean Tolpis succulenta (Asteraceae)

Genome‐wide genotyping and Bayesian inference method (BORICE) were employed to estimate outcrossing rates and paternity in two small plant populations of Tolpis succulenta (Asteraceae) on Graciosa island in the Azores. These two known extant populations of T. succulenta on Graciosa have recently evolved self‐compatibility. Despite the expectation that selfing would occur at an appreciable rate (self‐incompatible populations of the same species show low but nonzero selfing), high outcrossing was found in progeny arrays from maternal plants in both populations. This is inconsistent with an immediate transition to high selfing following the breakdown of a genetic incompatibility system. This finding is surprising given the small population sizes and the recent colonization of an island from self‐incompatible colonists of T. succulenta from another island in the Azores, and a potential paucity of pollinators, all factors selecting for selfing through reproductive assurance. The self‐compatible lineage(s) likely have high inbreeding depression (ID) that effectively halts the evolution of increased selfing, but this remains to be determined. Like their progeny, all maternal plants in both populations are fully outbred, which is consistent with but not proof of high ID. High multiple paternity was found in both populations, which may be due in part to the abundant pollinators observed during the flowering season.     

Genetic diversity and sex‐biased dispersal in the brown spotted pitviper (Protobothrops mucrosquamatus): Evidence from microsatellite markers

Dispersal plays a vital role in the geographical distribution, population genetic structure, quantity dynamics, and evolution of a species. Sex‐biased dispersal is common among vertebrates and many studies have documented a tendency toward male‐biased dispersal in mammals and female‐biased dispersal in birds. However, dispersal patterns in reptiles remain poorly understood. In this study, we explored the genetic diversity and dispersal patterns of the widely distributed Asian pitviper Protobothrops mucrosquamatus. In total, 16 polymorphic microsatellite loci were screened in 150 snakes (48 males, 44 females, 58 samples without sex information) covering most of their distribution. Microsatellite analysis revealed high genetic diversity in P. mucrosquamatus. Bayesian clustering of population assignment identified two major clusters for all populations, somewhat inconsistent with the mitochondrial DNA phylogeny of P. mucrosquamatus reported in previous research. Analyses based on 92 sex‐determined and 37 samples of P. mucrosquamatus from three small sites in Sichuan, China (Mingshan, Yibin, and Zizhong) consistently suggested female‐biased dispersal in P. mucrosquamatus, which is the first example of this pattern in snakes. The female‐biased dispersal patterns in P. mucrosquamatus may be explained by local resource competition.     

Genetic studies of various Prosopis species (Leguminosae,Section Algarobia) co‐occurring in oases of the Atacama Desert (northern Chile)

In the Atacama Desert from northern Chile (19–24°S), Prosopis (Leguminosae) individuals are restricted to oases that are unevenly distributed and isolated from each other by large stretches of barren landscape constituting an interesting study model as the degree of connectivity between natural populations depends on their dispersal capacity and the barriers imposed by the landscape. Our goal was to assess the genetic diversity and the degree of differentiation among groups of Prosopis individuals of different species from Section Algarobia and putative hybrids (hereafter populations) co‐occurring in these isolated oases from the Atacama Desert and determine whether genetic patterns are associated with dispersal barriers. Thirteen populations were sampled from oases located on three hydrographic basins (Pampa del Tamarugal, Rio Loa, and Salar de Atacama; northern, central, and southern basins, respectively). Individuals genotyped by eight SSRs show high levels of genetic diversity (H _O = 0.61, A _r = 3.5) and low but significant genetic differentiation among populations (F _ST = 0.128, F _ST‐ENA = 0.129, D _JOST = 0.238). The AMOVA indicates that most of the variation occurs within individuals (79%) and from the variance among individuals (21%); almost, the same variation can be found between basins and between populations within basins. Differentiation and structure results were not associated with the basins, retrieving up to four genetic clusters and certain admixture in the central populations. Pairwise differentiation comparisons among populations showed inconsistencies considering their distribution throughout the basins. Genetic and geographic distances were significantly correlated at global and within the basins considered (p < .02), but low correlation indices were obtained (r < .37). These results are discussed in relation to the fragmented landscape, considering both natural and non‐natural (humans) dispersal agents that may be moving Prosopis in the Atacama Desert.     

High genetic and epigenetic variation of transposable elements: Potential drivers to rapid adaptive evolution for the noxious invasive weed Mikania micrantha

Why invasive species can rapidly adapt to novel environments is a puzzling question known as the genetic paradox of invasive species. This paradox is explainable in terms of transposable elements (TEs) activity, which are theorized to be powerful mutational forces to create genetic variation. Mikania micrantha, a noxious invasive weed, in this sense provides an excellent opportunity to test the explanation. The genetic and epigenetic variation of 21 invasive populations of M. micrantha in southern China have been examined by using transposon display (TD) and transposon methylation display (TMD) techniques to survey 12 TE superfamilies. Our results showed that M. micrantha populations maintained an almost equally high level of TE‐based genetic and epigenetic variation and they have been differentiated into subpopulations genetically and epigenetically. A similar positive spatial genetic and epigenetic structure pattern was observed within 300 m. Six and seven TE superfamilies presented significant genetic and epigenetic isolation by distance (IBD) pattern. In total, 59 genetic and 86 epigenetic adaptive TE loci were identified. Of them, 51 genetic and 44 epigenetic loci were found to correlate with 25 environmental variables (including precipitation, temperature, vegetation coverage, and soil metals). Twenty‐five transposon‐inserted genes were sequenced and homology‐based annotated, which are found to be involved in a variety of molecular and cellular functions. Our research consolidates the importance of TE‐associated genetic and epigenetic variation in the rapid adaptation and invasion of M. micrantha.     

Life in the desert: The impact of geographic and environmental gradients on genetic diversity and population structure of Ivesia webberi

For range‐restricted species with disjunct populations, it is critical to characterize population genetic structure, gene flow, and factors that influence functional connectivity among populations in order to design effective conservation programs. In this study, we genotyped 314 individuals from 16 extant populations of Ivesia webberi, a United States federally threatened Great Basin Desert using six microsatellite loci. We assessed the effects of Euclidean distance, landscape features, and ecological dissimilarity on the pairwise genetic distance of the sampled populations, while also testing for a potential relationship between I. webberi genetic diversity and diversity in the vegetative communities. The results show low levels of genetic diversity overall (H _e = 0.200–0.441; H _o = 0.192–0.605) and high genetic differentiation among populations. Genetic diversity was structured along a geographic gradient, congruent with patterns of isolation by distance. Populations near the species’ range core have relatively high genetic diversity, supporting in part a central‐marginal pattern, while also showing some evidence for a metapopulation dynamic. Peripheral populations have lower genetic diversity, significantly higher genetic distances, and higher relatedness. Genotype cluster admixture results suggest a complex dispersal pattern among populations with dispersal direction and distance varying on the landscape. Pairwise genetic distance strongly correlates with elevation, actual evapotranspiration, and summer seasonal precipitation, indicating a role for isolation by environment, which the observed phenological mismatches among the populations also support. The significant correlation between pairwise genetic distance and floristic dissimilarity in the germinated soil seed bank suggests that annual regeneration in the plant communities contribute to the maintenance of genetic diversity in I. webberi.     

Mitochondrial DNA variation reveals maternal origins and demographic dynamics of Ethiopian indigenous goats

The Horn of Africa forms one of the two main historical entry points of domestics into the continent and Ethiopia is particularly important in this regard. Through the analysis of mitochondrial DNA (mtDNA) d‐loop region in 309 individuals from 13 populations, we reveal the maternal genetic variation and demographic dynamics of Ethiopian indigenous goats. A total of 174 variable sites that generated 231 haplotypes were observed. They defined two haplogroups that were present in all the 13 study populations. Reference haplotypes from the six globally defined goat mtDNA haplogroups show the two haplogroups present in Ethiopia to be A and G, the former being the most predominant. Although both haplogroups are characterized by an increase in effective population sizes (N_e) predating domestication, they also have experienced a decline in N_e at different time periods, suggesting different demographic histories. We observed seven haplotypes, six were directly linked to the central haplotypes of the two haplogroups and one was central to haplogroup G. The seven haplotypes were common between Ethiopia, Kenya, Egypt, and Saudi Arabia populations, suggesting common maternal history and the introduction of goats into East Africa via Egypt and the Arabian Peninsula, respectively. While providing new mtDNA data from a historically important region, our results suggest extensive intermixing of goats mediated by human socio‐cultural and economic interactions. These have led to the coexistence of the two haplogroups in different geographic regions in Ethiopia resulting in a large caprine genetic diversity that can be exploited for genetic improvement.     

Patterns of genomic and allochronic strain divergence in the fall armyworm,Spodoptera frugiperda (J.E. Smith)

Speciation is the process through which reproductive isolation develops between distinct populations. Because this process takes time, speciation studies often necessarily examine populations within a species that are at various stages of divergence. The fall armyworm, Spodoptera frugiperda (J.E. Smith), is comprised of two strains (R = Rice & C = Corn) that serve as a novel system to explore population divergence in sympatry. Here, we use ddRADSeq data to show that fall armyworm strains in the field are largely genetically distinct, but some interstrain hybridization occurs. Although we detected F1 hybrids of both R‐ and C‐strain maternal origin, only hybrids with R‐strain mtDNA were found to contribute to subsequent generations, possibly indicating a unidirectional barrier to gene flow. Although these strains have been previously defined as “host plant‐associated,” we recovered an equal proportion of R‐ and C‐strain moths in fields dominated by C‐strain host plants. As an alternative to host‐associated divergence, we tested the hypothesis that differences in nightly activity patterns could account for reproductive isolation by genotyping temporally collected moths. Our data indicates that strains exhibit a significant shift in the timing of their nightly activities in the field. This divergence in phenology creates a prezygotic reproductive barrier that likely maintains the genetic isolation between strains. Thus, we conclude that it may be ecologically inaccurate to refer to the C‐ and R‐ strain as “host‐associated” and they should more appropriately be considered “allochronic strains.”     

Genetic signature of the natural gene pool of Tilia cordata Mill. in Lithuania: Compound evolutionary and anthropogenic effects

Tilia cordata Mill. is a valuable tree species enriching the ecological values of the coniferous‐dominated boreal forests in Europe. Following the historical decline, spreading of Tilia sp. is challenged by the elevated inbreeding and habitat fragmentation. We studied the geographical distribution of genetic diversity of Tilia cordata populations in Lithuania. We used 14 genomic microsatellite markers to genotype 543 individuals from 23 wild‐growing populations. We found that Tilia cordata retained high levels of genetic diversity (population F _IS = 0–0.15, H _o = 0.53–0.69, H _e = 0.56–0.75). AMOVA, Bayesian clustering, and Monmonier''s barrier detection indicate weak but significant differentiation among the populations (F _ST = 0.037***) into geographically interpretable clusters of (a) western Lithuania with high genetic heterogeneity but low genetic diversity, bottleneck effects, (b) relatively higher genetic diversity of Tilia cordata on rich and most soils of midland lowland, and (c) the most differentiated populations on poor soils of the coolest northeastern highland possessing the highest rare allele frequency but elevated inbreeding and bottleneck effects. Weak genetic differentiation among the Tilia cordata populations in Lithuania implies common ancestry, absence of strong adaptive gradients, and effective genetic exchange possible mediated via the riparian networks. A hypothesis on riparian networks as gene flow mediators in Tilia cordata was raised based on results of this study.     

Investigation of the genetic diversity of domestic Capra hircus breeds reared within an early goat domestication area in Iran

Background

Iran is an area of particular interest for investigating goat diversity. Archaeological remains indicate early goat domestication (about 10 000 years ago) in the Iranian Zagros Mountains as well as in the high Euphrates valley and southeastern Anatolia. In addition, mitochondrial DNA data of domestic goats and wild ancestors (C. aegagrusor bezoar) suggest a pre-domestication management of wild populations in southern Zagros and central Iranian Plateau. In this study genetic diversity was assessed in seven Iranian native goat breeds, namely Markhoz, Najdi, Taleshi, Khalkhali, Naini, native Abadeh and Turki-Ghashghaei. A total of 317 animals were characterized using 14 microsatellite loci. Two Pakistani goat populations, Pahari and Teddy, were genotyped for comparison.

Results

Iranian goats possess a remarkable genetic diversity (average expected heterozygosity of 0.671 across loci, 10.7 alleles per locus) mainly accounted for by the within-breed component (G_ST = 5.9%). Positive and highly significant F_IS values in the Naini, Turki-Ghashghaei, Abadeh and Markhoz breeds indicate some level of inbreeding in these populations. Multivariate analyses cluster Iranian goats into northern, central and western groups, with the western breeds relatively distinct from the others. Pakistani breeds show some relationship with Iranian populations, even if their position is not consistent across analyses. Gene flow was higher within regions (west, north, central) compared to between regions but particularly low between the western and the other two regions, probably due to the isolating topography of the Zagros mountain range. The Turki-Ghashghaei, Najdi and Abadeh breeds are reared in geographic areas where mtDNA provided evidence of early domestication. These breeds are highly variable, located on basal short branches in the neighbor-joining tree, close to the origin of the principal component analysis plot and, although highly admixed, they are quite distinct from those reared on the western side of the Zagros mountain range.

Conclusions

These observations call for further investigation of the nuclear DNA diversity of these breeds within a much wider geographic context to confirm or re-discuss the current hypothesis (based on maternal lineage data) of an almost exclusive contribution of the eastern Anatolian bezoar to the domestic goat gene pool.     

Georeferenced phylogenetic analysis of a global collection of wild and cultivated Citrullus species

The geographical origin of watermelon (Citrullus lanatus) remains debated. While a first hypothesis suggests the center of origin to be West Africa, where the endemic sister species C. mucosospermus thrives, a second hypothesis suggests northeastern Africa where the white‐fleshed Sudanese Kordophan melon is cultivated. In this study, we infer biogeographical and haplotype genealogy for C. lanatus, C. mucosospermus, C. amarus, and C. colocynthis using noncoding cpDNA sequences (trnT‐trnL and ndhF‐rpl32 regions) from a global collection of 135 accessions. In total, we identified 38 haplotypes in C. lanatus, C. mucosospermus, C. amarus, and C. colocynthis; of these, 21 were found in Africa and 17 appear endemic to the continent. The least diverse species was C. mucosospermus (5 haplotypes) and the most diverse was C. colocynthis (16 haplotypes). Some haplotypes of C. mucosospermus were nearly exclusive to West Africa, and C. lanatus and C. mucosospermus shared haplotypes that were distinct from those of both C. amarus and C. colocynthis. The results support previous findings that revealed C. mucosospermus to be the closest relative to C. lanatus (including subsp. cordophanus). West Africa, as a center of endemism of C. mucosospermus, is an area of interest in the search of the origin of C. lanatus. This calls for further historical and phylogeographical investigations and wider collection of samples in West and northeastern Africa.