A new species of Sturisoma Swainson, 1838 (Loricariidae: Loricariinae), from the Madeira River basin, with a discussion of historical biogeography of western Amazonas and Paraguay River basins

A new loricariin species of Sturisoma is described from the Cautário, Guaporé, Mamoré, Machado and Soteiro rivers, Madeira River basin, in Bolivia and Brazil. The new species is distinguished from its congeners by the presence of a middorsal longitudinal, thin dark brown stripe on the caudal peduncle, extending from two or three plates posterior to the dorsal-fin base, reaching the origin of the caudal fin, or one or two plates anterior to the origin of the caudal fin; small squarish anteriormost abdominal plates; and a middorsal longitudinal dark-brown stripe from first predorsal plate to near the dorsal-fin origin. Furthermore, the new species is diagnosed from congeners by plate morphology, counts on the median series, coalescent plates and ventrolateral thoracic plates, in addition to measurements related to body and head structures. An analysis of genetic distances using cytochrome C oxidase subunit 1 gene marker of the mitochondrial genome between the new species and several congeners is presented, in addition to a likelihood analysis to illustrate the position of the new taxon within Sturisoma. An identification key for species of the genus currently recorded at the upper Amazonas River basin is provided.

Genetic diversity in 10 populations of domestic Turkeys by using microsatellites markers

The domestic turkey is a native breed in danger of extinction due to the introduction of new breeds specializing in meat production and yield. Turkeys have lost some prominence in urban areas, and only certain breeds of turkeys are preserved in rural areas. Wild and domestic turkeys are different; rural or indigenous turkeys, with black plumage, were domesticated from Mexican turkeys and have been reproduced throughout Latin America. Some of them were taken to Europe in the 16th century and later arrived in North America, where they crossed with another wild species, from which the bronze turkey emerged: the ancestor of all commercial turkeys. The objective of the present work was to evaluate the genetic diversity in 10 populations of domestic turkeys worldwide by using breeds from Europe: Spain and Italy; America: Mexico, United States and Brazil; and the Near East: Iran and Egypt. A total of 522 blood samples of both sexes were collected from domestic turkey populations. Thirty-four microsatellites were used to obtain genetic parameters, and genetic diversity was evaluated. All microsatellites used were polymorphic, and a total of 427 alleles were detected across the 34 markers investigated. In this study, a mean number of 13.44 alleles was found. The four most diverse breeds were from the Andalusia, Mexico, United States, and wild populations, which had the highest mean heterozygosity expected (0.619, 0.612, 0.650, and 0.773) and heterozygosity observed (0.422, 0.521, 0.429, and 0.627), respectively. The MNT348 marker deviated from the HWE in all populations. Our study has shown that the populations close to the species origin are more diverse than those resulting from posterior expansions. Mexican birds were the most diverse, followed by the Spanish populations because Spain imported a large number of turkeys coming from America. Such information can be complementary to other genotypic data required to validate the evolutionary relationships among turkey populations.

Comparative Genetic Analysis of Durum Wheat Landraces and Cultivars Widespread in Tunisia

The durum wheat (Triticum turgidum L. ssp. durum Desf.) landraces constitute a useful natural germplasm to increase the genetic diversity in the modern durum cultivars. The Tunisian durum germplasm constitutes 28 accessions conserved in Genebank of Tunisia, which are still unexplored. In this study, a comparative genetic analysis was performed to investigate the relationships between the Tunisian durum lines and the modern cultivars and detect divergent loci involved in breeding history. The genetic diversity analyses carried out using nine morphological descriptors and the 25K single-nucleotide polymorphism (SNP) array allowed us to distinguish two groups of Tunisian landraces and one of durum cultivars. The analysis of molecular variance and diversity indices confirmed the genetic variability among the groups. A total of 529 SNP loci were divergent between Tunisian durum landraces and modern cultivars. Candidate genes related to plant and spike architecture, including FLOWERING LOCUS T (FT-B1), zinc finger CONSTANS, and AP2/EREBPs transcription factors, were identified. In addition, divergent genes involved in grain composition and biotic stress nucleotide-binding site and leucine-reach repeats proteins and disease resistance proteins (NBS-LRR and RPM) were found, suggesting that the Tunisian durum germplasm may represent an important source of favorable alleles to be used in future durum breeding programs for developing well-adapted and resilient cultivars.

Comparative Analysis of Eight Mitogenomes of Bark Beetles and Their Phylogenetic Implications

Simple Summary

Many bark beetles are destructive pests in coniferous forests and cause extensive ecological and economic losses worldwide. Comparative studies of the structural characteristics of mitogenomes and phylogenetic relationships of bark beetles can improve our understanding of mitogenome evolution. In this study, we sequenced eight mitogenomes of bark beetles. Our results show that the use of start and stop codons, the abundance of amino acids, and the relative frequency of codon use are conserved among the eight bark beetles. Different regions of tRNA exhibit different degrees of conservatism. Together with the analysis of evolutionary rates and genetic distance among bark beetle species, our results reveal phylogenetic relationships among bark beetles of the subfamily Scolytinae.

Abstract

Many bark beetles of the subfamily Scolytinae are the most economically important insect pests of coniferous forests worldwide. In this study, we sequenced the mitochondrial genomes of eight bark beetle species, including Dendroctonus micans, Orthotomicus erosus, Polygraphus poligraphus, Dryocoetes hectographus, Ips nitidus, Ips typographus, Ips subelongatus, and Ips hauseri, to examine their structural characteristics and determine their phylogenetic relationships. We also used previously published mitochondrial genome sequence data from other Scolytinae species to identify and localize the eight species studied within the bark beetle phylogeny. Their gene arrangement matched the presumed ancestral pattern of these bark beetles. Start and stop codon usage, amino acid abundance, and the relative codon usage frequencies were conserved among bark beetles. Genetic distances between species ranged from 0.037 to 0.418, and evolutionary rates of protein-coding genes ranged from 0.07 for COI to 0.69 for ND2. Our results shed light on the phylogenetic relationships and taxonomic status of several bark beetles in the subfamily Scolytinae and highlight the need for further sequencing analyses and taxonomic revisions in additional bark beetle species.

Proposed Update to the Taxonomy of Pestiviruses: Eight Additional Species within the Genus Pestivirus, Family Flaviviridae

Pestiviruses are plus-stranded RNA viruses belonging to the family Flaviviridae. They comprise several important pathogens like classical swine fever virus and bovine viral diarrhea virus that induce economically important animal diseases. In 2017, the last update of pestivirus taxonomy resulted in demarcation of 11 species designated Pestivirus A through Pestivirus K. Since then, multiple new pestiviruses have been reported including pathogens associated with disease in pigs or small ruminants. In addition, pestivirus sequences have been found during metagenomics analysis of different non-ungulate hosts (bats, rodents, whale, and pangolin), but the consequences of this pestivirus diversity for animal health still need to be established. To provide a systematic classification of the newly discovered viruses, we analyzed the genetic relationship based on complete coding sequences (cds) and deduced polyprotein sequences and calculated pairwise distances that allow species demarcation. In addition, phylogenetic analysis was performed based on a highly conserved region within the non-structural protein NS5B. Taking into account the genetic relationships observed together with available information about antigenic properties, host origin, and characteristics of disease, we propose to expand the number of pestivirus species to 19 by adding eight additional species designated Pestivirus L through Pestivirus S.

New species of the genus Inversidens Haas, 1911 (Unionoida, Unionidae, Gonideinae) from Jiangxi Province, China

We diagnose and describe a new freshwater mussel species of the genus Inversidens, I.rentianensis sp. nov. from Jiangxi Province, China based on morphological characters and molecular data. This paper includes a morphological description and photograph of the holotype, and partial sequences of mitochondrial COI as DNA barcode data.

Two New Species of the Genus Eucorydia (Blattodea: Corydiidae) from the Nansei Islands in Southwest Japan

Two new species of the cockroach genus Eucorydia Hebard, 1929 from the Nansei Islands in Southwest Japan were compared to two closely related congeners, Eucorydia yasumatsui Asahina, 1971 and Eucorydia dasytoides (Walker, 1868). Eucorydia donanensis Yanagisawa, Sakamaki, and Shimano sp. nov. from Yonaguni-jima Island was characterized by an overall length of 12.5-14.5 mm in males. The dorsal side of the male abdomen was entirely dark purple and there was an obscure orange band running down the middle of the tegmen. Eucorydia tokaraensis Yanagisawa, Sakamaki, and Shimano sp. nov. was characterized by an overall length of 12.0-13.0 mm in males and a distinct orange band running down the middle of the tegmen. Eucorydia yasumatsui, E. donanensis, E. tokaraensis and the zonata population of E. dasytoides were divided into four lineages in a maximum-likelihood tree generated from a dataset concatenated from five (two nuclear, 28S rRNA, histone H3, and three mitochondrial, COII, 12S rRNA, 16S rRNA) genes. We recognized the three Japanese lineages E. yasumatsui, E. donanensis, and E. tokaraensis as distinct species, which were also supported by the pairwise genetic distances (5.4-7.8%, K2P) of the COI sequences. Morphometric analysis was performed on the genitalia. A principal component analysis plot revealed that the sizes of the genitalia in the three Japanese species were similar to each other and smaller than that of the zonata population of E. dasytoides. The analysis also revealed that the three Japanese species were distinguished from each other by combinations of the sizes of L3 and L7 sclerites and the shape of R2 sclerite, with some overlapping exceptions.

Convergent Evolution of Himalayan Marmot with Some High-Altitude Animals through ND3 Protein

The Himalayan marmot (Marmota himalayana) mainly lives on the Qinghai-Tibet Plateau and it adopts multiple strategies to adapt to high-altitude environments. According to the principle of convergent evolution as expressed in genes and traits, the Himalayan marmot might display similar changes to other local species at the molecular level. In this study, we obtained high-quality sequences of the CYTB gene, CYTB protein, ND3 gene, and ND3 protein of representative species (n = 20) from NCBI, and divided them into the marmot group (n = 11), the plateau group (n = 8), and the Himalayan marmot (n = 1). To explore whether plateau species have convergent evolution on the microscale level, we built a phylogenetic tree, calculated genetic distance, and analyzed the conservation and space structure of Himalayan marmot ND3 protein. The marmot group and Himalayan marmots were in the same branch of the phylogenetic tree for the CYTB gene and CYTB protein, and mean genetic distance was 0.106 and 0.055, respectively, which was significantly lower than the plateau group. However, the plateau group and the Himalayan marmot were in the same branch of the phylogenetic tree, and the genetic distance was only 10% of the marmot group for the ND3 protein, except Marmota flaviventris. In addition, some sites of the ND3 amino acid sequence of Himalayan marmots were conserved from the plateau group, but not the marmot group. This could lead to different structures and functional diversifications. These findings indicate that Himalayan marmots have adapted to the plateau environment partly through convergent evolution of the ND3 protein with other plateau animals, however, this protein is not the only strategy to adapt to high altitudes, as there may have other methods to adapt to this environment.

Breeding system of diploid sexuals within the Ranunculus auricomus complex and its role in a geographical parthenogenesis scenario

The larger distribution area of asexuals compared with their sexual relatives in geographical parthenogenesis (GP) scenarios has been widely attributed to the advantages of uniparental reproduction and polyploidy. However, potential disadvantages of sexuals due to their breeding system have received little attention so far. Here, we study the breeding system of five narrowly distributed sexual lineages of Ranunculus notabilis s.l. (R. auricomus complex) and its effects on outcrossing, inbreeding, female fitness, and heterozygosity. We performed selfing and intra- and interlineage crossings by bagging 481 flowers (59 garden individuals) followed by germination experiments. We compared seed set and germination rates, and related them to genetic distance and genome-wide heterozygosity (thousands of RADseq loci). Selfings (2.5%) unveiled a significantly lower seed set compared with intra- (69.0%) and interlineage crossings (69.5%). Seed set of intra- (65%) compared to interpopulation crossings (78%) was significantly lower. In contrast, all treatments showed comparable germination rates (32%-43%). Generalized linear regressions between seed set and genetic distance revealed positive relationships in general and between lineages, and a negative one within lineages. Seed set was the main decisive factor for female fitness. Germination rates were not related to genetic distance at any level, but were positively associated with heterozygosity in interlineage crossings. Experiments confirmed full crossability and predominant outcrossing among sexual R. notabilis s.l. lineages. However, up to 5% (outliers 15%-31%) of seeds were formed by selfing, probably due to semi-self-compatibility in a multi-locus gametophytic SI system. Less seed set in intrapopulation crossings, and higher seed set and germination rates from crossings of genetically more distant and heterozygous lineages (interlineage) indicate negative inbreeding and positive outbreeding effects. In GP scenarios, sexual species with small and/or isolated populations can suffer from decreased female fitness due to their breeding system. This factor, among others, probably limits range expansion of sexuals.

genodive version 3.0: Easy-to-use software for the analysis of genetic data of diploids and polyploids

genodive version 3.0 is a user-friendly program for the analysis of population genetic data. This version presents a major update from the previous version and now offers a wide spectrum of different types of analyses. genodive has an intuitive graphical user interface that allows direct manipulation of the data through transformation, imputation of missing data, and exclusion and inclusion of individuals, population and/or loci. Furthermore, genodive seamlessly supports 15 different file formats for importing or exporting data from or to other programs. One major feature of genodive is that it supports both diploid and polyploid data, up to octaploidy (2n = 8x) for some analyses, but up to hexadecaploidy (2n = 16x) for other analyses. The different types of analyses offered by genodive include multiple statistics for estimating population differentiation (φ_ST , F_ST , F'_ST , G_ST , G'_ST , G''_ST , D_est , R_ST , ρ), analysis of molecular variance-based K-means clustering, Hardy-Weinberg equilibrium, hybrid index, population assignment, clone assignment, Mantel test, Spatial Autocorrelation, 23 ways of calculating genetic distances, and both principal components and principal coordinates analyses. A unique feature of genodive is that it can also open data sets with nongenetic variables, for example environmental data or geographical coordinates that can be included in the analysis. In addition, genodive makes it possible to run several external programs (lfmm, structure, instruct and vegan) directly from its own user interface, avoiding the need for data reformatting and use of the command line. genodive is available for computers running Mac OS X 10.7 or higher and can be downloaded freely from: http://www.patrickmeirmans.com/software.

HLA diversity in Brazil

Brazil is the fifth largest country in the world in area and the fifth most populous. The Brazilian voluntary Bone Marrow Donor Registry is the third largest in terms of number of donors in the world, being a valuable source of HLA genetics to characterize the donor population of Brazil as well. The genetic background of the Brazilian population is quite heterogeneous, resulting from 5 centuries of admixture among Native Americans, Europeans and Africans, making the Brazilian population unique in terms of genetic ancestry. The unique characteristics of populations in different Brazilian regions make them an exciting focus for genetic diversity studies. Studies on HLA genetic diversity of Brazilian populations have been conducted since the late 1980s and, in this review, we highlight the main findings from studies carried out in Brazil based on classical HLA. In addition, we calculated the genetic distance from the molecular data of the studies included in this review in order to have a broader view of the HLA diversity in Brazilian populations. We emphasize that characterization of HLA diversity is not only important for transplantation programs, but can shed a light on ancestry, history and other demographic patterns with or without association with autoimmune disease.

Dispersal of Amur tiger from spatial distribution and genetics within the eastern Changbai mountain of China

Population dispersal and migration often indicate an expanded habitat and reduced inbreeding probability, and to some extend reflects improvement in the condition of the population. The Amur tiger population in the northern region of the Changbai mountain in China mostly distributes along the Sino-Russian border, next to the population in southwest Primorye in Russia. The successful dispersal westward and transboundary movement are crucial for the persistence of the Amur tiger in this area. This study explored the spatial dispersal of the population, transboundary migration, and the genetic condition of the Amur tiger population within the northern Changbai mountain in China, using occurrence data and fecal samples. Our results from 2003 to 2016 showed that the Amur tiger population in this area was spreading westward at a speed of 12.83 ± 4.41 km every three years. Genetic diversity of the Amur tiger populations in southwest Primorye was slightly different than the population in our study area, and the potential individual migration rate between these two populations was shown to be about 13.04%. Furthermore, the relationships between genetic distances and spatial distances indicated the existence of serious limitations to the dispersal of the Amur tiger in China. This study provided important information about spatial dispersal, transboundary migration, and the genetic diversity of Amur tigers in China, showed the urgent need for Amur tiger habitat restoration, and suggested some important conservation measures, such as corridor construction to eliminate dispersal barriers and joint international conservation to promote trans-boundary movement.

Coalescing molecular evolution and DNA barcoding

The DNA barcoding concept (Woese et al. ; Hebert et al. ) has considerably boosted taxonomy research by facilitating the identification of specimens and discovery of new species. Used alone or in combination with DNA metabarcoding on environmental samples (Taberlet et al. ), the approach is becoming a standard for basic and applied research in ecology, evolution and conservation across taxa, communities and ecosystems (Scheffers et al. ; Kress et al. ). However, DNA barcoding suffers from several shortcomings that still remain overlooked, especially when it comes to species delineation (Collins & Cruickshank ). In this issue of Molecular Ecology, Barley & Thomson () demonstrate that the choice of models of sequence evolution has substantial impacts on inferred genetic distances, with a propensity of the widely used Kimura 2-parameter model to lead to underestimated species richness. While DNA barcoding has been and will continue to be a powerful tool for specimen identification and preliminary taxonomic sorting, this work calls for a systematic assessment of substitution models fit on barcoding data used for species delineation and reopens the debate on the limitation of this approach.

SYMPATRIC CONVERGENCE AND ENVIRONMENTAL CORRELATION BETWEEN TWO LAND-SNAIL SPECIES

In the southern Appalachian region of North America, the phylogenetically convergent shells of the polygyrid snails Triodopsinae Neohelix major (Binney) and Polygyrinae Mesodon normalis (Pilsbry) are even more convergent in size and shape in sympatry (7 sites) than in allopatry (23 and 10 sites). Environmental correlations account for 34% and 30% of size and shape variations in N. major (larger, taller, and more loosely coiled at northern, high-altitude, sheltered sites), but for only 14% and 9% in M. normalis (larger, flatter, and more loosely coiled at south-facing, exposed sites). The statistical significance of the sympatric convergence dropped out when these correlations were removed. This phenomenon helps account for the many cases in eastern North America of nearly identical land-snail shells in sympatry and questions the importance of competitive character displacement in the evolution of land-snail shell morphology. This apparently nonmimetic case of sympatric convergence provides an unusually precise and well-delimited, naturally replicated experiment in evolutionary morphology, which is analyzed for controlling factors in a follow-up paper.

Karyotype characteristics, larval morphology and chromosomal polymorphism peculiarities of Glyptotendipes salinus Michailova, 1983 (Diptera, Chironomidae) from Tambukan Lake, Central Caucasus

Data on the karyotype characteristics, larval morphology and features of chromosomal polymorphism of a population of Glyptotendipes salinus Michailova, 1983 (Diptera, Chironomidae) from Tambukan Lake (on the northern macroslope of the central Caucasus) are presented. It was found that diagnostic larval characters of Glyptotendipes salinus from Caucasus in general are similar to those described in previous studies, but with some significant differences. By some morphological characteristics Caucasian larvae appeared to be closer to Glyptotendipes barbipes than to ones provided for European larvae of Glyptotendipes salinus by Contreras-Lichtenberg (1999). Obtained morphological data make possible to conclude that Caucasian population of Glyptotendipes salinus can be a markedly diverged population of the species, probably even subspecies. In the Caucasian population 12 banding sequences were found: two in arms A, B, C, E, and G, and one in arms D and F. Eight of these are already known for this species, and four, salA2, salB2, salEX, and salG3, are described for the first time. Genetic distances between all the previously studied populations of Glyptotendipes salinus were measured using Nei criteria (1972). The population of the central Caucasus occupies a distinct position on the dendrogram compared with populations from Altai and Kazakhstan. All the obtained morphological and cytogenetic data can indicate the plausible relative isolation and complexity of the Caucasus from the viewpoint of microevolution. More researches are required in other parts of Caucasus and other geographically distant regions for more specific allegations.

Assessing the performance of DNA barcoding using posterior predictive simulations

Accurate estimates of biodiversity are required for research in a broad array of biological subdisciplines including ecology, evolution, systematics, conservation and biodiversity science. The use of statistical models and genetic data, particularly DNA barcoding, has been suggested as an important tool for remedying the large gaps in our current understanding of biodiversity. However, the reliability of biodiversity estimates obtained using these approaches depends on how well the statistical models that are used describe the evolutionary process underlying the genetic data. In this study, we utilize data from the Barcode of Life Database and posterior predictive simulations to assess the performance of DNA barcoding under commonly used substitution models. We demonstrate that the success of DNA barcoding varies widely across DNA substitution models and that model choice has a substantial impact on the number of operational taxonomic units identified (changing results by ~4-31%). Additionally, we demonstrate that the widely followed practice of a priori assuming the Kimura 2-parameter model for DNA barcoding is statistically unjustified and should be avoided. Using both data-based and inference-based test statistics, we detect variation in model performance across taxonomic groups, clustering algorithms, genetic divergence thresholds and substitution models. Taken together, these results illustrate the importance of considering both model selection and model adequacy in studies quantifying biodiversity.

Molecular-based estimate of species number, phylogenetic relationships and divergence times for the genus Stenotaenia (Chilopoda, Geophilomorpha) in the Italian region

Stenotaenia is one of the largest and most widespread genera of geophilid centipedes in the Western Palearctic, with a very uniform morphology and about fifteen species provisionally recognized. For a better understanding of Stenotaenia species-level taxonomy, we have explored the possibility of using molecular data. As a preliminary assay, we sampled twelve populations, mainly from the Italian region, and analyzed partial sequences of the two genes COI and 28S. We employed a DNA-barcoding approach, complemented by a phylogenetic analysis coupled with divergence time estimation. Assuming a barcoding gap of 10-16% K2P pairwise distances, we found evidence for the presence of at least six Stenotaenia species in the Italian region, which started diverging about 50 million years ago, only partially matching with previously recognized species. We found that small-sized oligopodous species belong to a single clade that originated about 33 million years ago, and obtained some preliminary evidence of the related genus Tuoba being nested within Stenotaenia.

Mantel test in population genetics

The comparison of genetic divergence or genetic distances, estimated by pairwise FST and related statistics, with geographical distances by Mantel test is one of the most popular approaches to evaluate spatial processes driving population structure. There have been, however, recent criticisms and discussions on the statistical performance of the Mantel test. Simultaneously, alternative frameworks for data analyses are being proposed. Here, we review the Mantel test and its variations, including Mantel correlograms and partial correlations and regressions. For illustrative purposes, we studied spatial genetic divergence among 25 populations of Dipteryx alata ("Baru"), a tree species endemic to the Cerrado, the Brazilian savannas, based on 8 microsatellite loci. We also applied alternative methods to analyze spatial patterns in this dataset, especially a multivariate generalization of Spatial Eigenfunction Analysis based on redundancy analysis. The different approaches resulted in similar estimates of the magnitude of spatial structure in the genetic data. Furthermore, the results were expected based on previous knowledge of the ecological and evolutionary processes underlying genetic variation in this species. Our review shows that a careful application and interpretation of Mantel tests, especially Mantel correlograms, can overcome some potential statistical problems and provide a simple and useful tool for multivariate analysis of spatial patterns of genetic divergence.

Mantel test in population genetics

The comparison of genetic divergence or genetic distances, estimated by pairwise FST and related statistics, with geographical distances by Mantel test is one of the most popular approaches to evaluate spatial processes driving population structure. There have been, however, recent criticisms and discussions on the statistical performance of the Mantel test. Simultaneously, alternative frameworks for data analyses are being proposed. Here, we review the Mantel test and its variations, including Mantel correlograms and partial correlations and regressions. For illustrative purposes, we studied spatial genetic divergence among 25 populations of Dipteryx alata ("Baru"), a tree species endemic to the Cerrado, the Brazilian savannas, based on 8 microsatellite loci. We also applied alternative methods to analyze spatial patterns in this dataset, especially a multivariate generalization of Spatial Eigenfunction Analysis based on redundancy analysis. The different approaches resulted in similar estimates of the magnitude of spatial structure in the genetic data. Furthermore, the results were expected based on previous knowledge of the ecological and evolutionary processes underlying genetic variation in this species. Our review shows that a careful application and interpretation of Mantel tests, especially Mantel correlograms, can overcome some potential statistical problems and provide a simple and useful tool for multivariate analysis of spatial patterns of genetic divergence.

Mitochondrial D-loop analysis reveals low diversity in Mangalica pigs and their relationship to historical specimens

The genetic relationship between 195 Mangalica and 79 non-Mangalica pigs was studied using mitochondrial D-loop SNP genotyping. Altogether, 35 polymorphic sites and 27 haplotypes were identified. Of the haplotypes, eight and 16 are Mangalica and non-Mangalica specific, respectively, while three contain both Mangalica and non-Mangalica individuals. Genetic distance values and phylogenetic analysis indicate that Mangalica individuals are very closely related, and five haplotypes represent approximately 92% of the Mangalica pigs involved in the study, thus determining the major maternal lineages. In contrast to previous microsatellite studies, individuals of Mangalica could not be distinguished as three separate breeds using mtDNA genotyping. Comparing modern and archaeological mtDNA sequences revealed that present day Mangalica is related to pigs that lived in the Carpathian basin where postulated ancestors of Mangalica also lived. This is the first DNA-based genetic evidence to support the described breeding history of Mangalica.

Unraveling reticulate evolution in North American Dryopteris (Dryopteridaceae)

BACKGROUND

The thirteen species of Dryopteris in North America have long been suspected of having undergone a complicated history of reticulate evolution via allopolyploid hybridization. Various explanations for the origins of the allopolyploid taxa have been suggested, and though most lines of evidence have supported the so-called "semicristata" hypothesis, contention over the group's history has continued in several recent, conflicting studies.

RESULTS

Sequence data from nine plastid and two nuclear markers were collected from 73 accessions representing 35 species of Dryopteris. Sequences from each of the allopolyploids are most closely related to their progenitor species as predicted by the "semicristata" hypothesis. Allotetraploid D. campyloptera appears to be derived from a hybrid between diploid D. expansa and D. intermedia; D. celsa, from diploid D. ludoviciana x D. goldiana; and D. carthusiana and D. cristata, from diploid "D. semicristata" x D. intermedia and D. ludoviciana, respectively. Allohexaploid D. clintoniana appears to be derived from D. cristata x D.goldiana. The earliest estimated dates of formation of the allopolyploids, based on divergence time analyses, were within the last 6 Ma. We found no evidence for recurrent formation of any of the allopolyploids. The sexual allopolyploid taxa are derived from crosses between parents that show intermediate levels of genetic divergence relative to all pairs of potential progenitors. In addition, the four allotetraploids are transgressive with respect to geographic range relative to one or both of their parents (their ranges extend beyond those of the parents), suggesting that ecological advantages in novel habitats or regions may promote long-term regional coexistence of the hybrid taxa with their progenitors.

CONCLUSIONS

This study provides the first thorough evaluation of the North American complex of woodferns using extensive sampling of taxa and genetic markers. Phylogenies produced from each of three datasets (one plastid and two nuclear) support the "semicristata" hypothesis, including the existence of a missing diploid progenitor, and allow us to reject all competing hypotheses. This study demonstrates the value of using multiple, biparentally inherited markers to evaluate reticulate complexes, assess the frequency of recurrent polyploidization, and determine the relative importance of introgression vs. hybridization in shaping the histories of such groups.