Empirical evaluation of genetic clustering methods using multilocus genotypes from 20 chicken breeds

Clumppling: cluster matching and permutation program with integer linear programming

MOTIVATION

In the mixed-membership unsupervised clustering analyses commonly used in population genetics, multiple replicate data analyses can differ in their clustering solutions. Combinatorial algorithms assist in aligning clustering outputs from multiple replicates so that clustering solutions can be interpreted and combined across replicates. Although several algorithms have been introduced, challenges exist in achieving optimal alignments and performing alignments in reasonable computation time.

RESULTS

We present Clumppling, a method for aligning replicate solutions in mixed-membership unsupervised clustering. The method uses integer linear programming for finding optimal alignments, embedding the cluster alignment problem in standard combinatorial optimization frameworks. In example analyses, we find that it achieves solutions with preferred values of a desired objective function relative to those achieved by Pong and that it proceeds with less computation time than Clumpak. It is also the first method to permit alignments across replicates with multiple arbitrary values of the number of clusters K.

AVAILABILITY AND IMPLEMENTATION

Clumppling is available at https://github.com/PopGenClustering/Clumppling.

SNP data reveals the complex and diverse evolutionary history of the blue-ringed octopus genus (Octopodidae: Hapalochlaena) in the Asia-Pacific

The blue-ringed octopus species complex (Hapalochlaena spp.), known to occur from Southern Australia to Japan, currently contains four formally described species (Hapalochlaena maculosa, Hapalochlaena fasciata, Hapalochlaena lunulata and Hapalochlaena nierstraszi). These species are distinguished based on morphological characters (iridescent blue rings and/or lines) along with reproductive strategies. However, the observation of greater morphological diversity than previously captured by the current taxonomic framework indicates that a revision is required. To examine species boundaries within the genus we used mitochondrial (12S rRNA, 16S rRNA, cytochrome c oxidase subunit 1 [COI], cytochrome c oxidase subunit 3 [COIII] and cytochrome b [Cytb]) and genome-wide SNP data (DaRT seq) from specimens collected across its geographic range including variations in depth from 3m to >100m. This investigation indicates substantially greater species diversity present within the genus Hapalochlaena than is currently described. We identified 10,346 SNPs across all locations, which when analysed support a minimum of 11 distinct clades. Bayesian phylogenetic analysis of the mitochondrial COI gene on a more limited sample set dates the diversification of the genus to ∼30mya and corroborates eight of the lineages indicated by the SNP analyses. Furthermore, we demonstrate that the diagnostic lined patterning of H. fasciata found in North Pacific waters and NSW, Australia is polyphyletic and therefore likely the result of convergent evolution. Several "deep water" (> 100m) lineages were also identified in this study with genetic convergence likely to be driven by external selective pressures. Examination of morphological traits, currently being undertaken in a parallel morphological study, is required to describe additional species within the complex.

Crimp: An efficient tool for summarizing multiple clusterings in population structure analysis and beyond

When a data set is repeatedly clustered using unsupervised techniques, the resulting clusterings, even if highly similar, may list their clusters in different orders. This so-called 'label-switching' phenomenon obscures meaningful differences between clusterings, complicating their comparison and summary. The problem often arises in the context of population structure analysis based on multilocus genotype data. In this field, a variety of popular tools apply model-based clustering, assigning individuals to a prespecified number of ancestral populations. Since such methods often involve stochastic components, it is a common practice to perform multiple replicate analyses based on the same input data and parameter settings. Available postprocessing tools allow to mitigate label switching, but leave room for improvements, in particular, regarding large input data sets. In this work, I present Crimp, a lightweight command-line tool, which offers a relatively fast and scalable heuristic to align clusters across replicate clusterings consisting of the same number of clusters. For small problem sizes, an exact algorithm can be used as an alternative. Additional features include row-specific weights, input and output files similar to those of CLUMPP (Jakobsson & Rosenberg, 2007) and the evaluation of a given solution in terms of CLUMPP as well as its own objective functions. Benchmark analyses show that Crimp, especially when applied to larger data sets, tends to outperform alternative tools considering runtime requirements and various quality measures. While primarily targeting population structure analysis, Crimp can be used as a generic tool to correct multiple clusterings for label switching. This facilitates their comparison and allows to generate an averaged clustering. Crimp's computational efficiency makes it even applicable to relatively large data sets while offering competitive solution quality.

A Dirichlet model of alignment cost in mixed-membership unsupervised clustering

Mixed-membership unsupervised clustering is widely used to extract informative patterns from data in many application areas. For a shared data set, the stochasticity and unsupervised nature of clustering algorithms can cause difficulties in comparing clustering results produced by different algorithms, or even multiple runs of the same algorithm, as outcomes can differ owing to permutation of the cluster labels or genuine differences in clustering results. Here, with a focus on inference of individual genetic ancestry in population-genetic studies, we study the cost of misalignment of mixed-membership unsupervised clustering replicates under a theoretical model of cluster memberships. Using Dirichlet distributions to model membership coefficient vectors, we provide theoretical results quantifying the alignment cost as a function of the Dirichlet parameters and the Hamming permutation difference between replicates. For fixed Dirichlet parameters, the alignment cost is seen to increase with the Hamming distance between permutations. Data sets with low variance across individuals of membership coefficients for specific clusters generally produce high misalignment costs-so that a single optimal permutation has far lower cost than suboptimal permutations. Higher variability in data, as represented by greater variance of membership coefficients, generally results in alignment costs that are similar between the optimal permutation and suboptimal permutations. We demonstrate the application of the theoretical results to data simulated under the Dirichlet model, as well as to membership estimates from inference of human-genetic ancestry. The results can contribute to improving cluster alignment algorithms that seek to find optimal permutations of replicates.

Genetic Diversity and Population Structure of Spirobolus bungii as Revealed by Mitochondrial DNA Sequences

Soil macrofauna, such as Spirobolus bungii, are an important component of ecosystems. However, systematic studies of the genetic diversity, population genetic structure, and the potential factors affecting the genetic differentiation of S. bungii are lacking. We performed a population genetic study of 166 individuals from the mountains to the south of the Yangtze River, north of the Yangtze River in Nanjing city, and near Tianjin city, in order to investigate the correlations between geographical distance and genetic diversity. A total of 1182 bp of COX2 and Cytb gene sequences of mitochondrial DNA, and 700 bp of the 18S rRNA gene sequence were analyzed. There were two haplotypes and one variable site in the 18S rRNA gene, and 28 haplotypes and 78 variable sites in the COX2 and Cytb genes. In this study, the 18S rRNA gene was used for species identification, and mtDNA (concatenated sequences with Cytb and COX2) was used for population genetic analysis. Structure cluster analysis indicated that the genetic structures of the different populations of S. bungii tended to be consistent at small geographical scales. Phylogenetic trees revealed that the haplotypes were clearly divided into three branches: the area south of the Yangtze River, the area to the north of the Yangtze River in Nanjing, and the area in Tianjin. Large geographical barriers and long geographical distance significantly blocked gene flow between populations of S. bungii. Our results provide a basic theoretical basis for subsequent studies of millipede taxonomy and population genetic evolution.

Genetic diversity and population structure of pigeonpea (Cajanus cajan [L.] Millspaugh) landraces grown in Benin revealed by Genotyping-By-Sequencing

Genetic diversity studies provide important details on target trait availability and its variability, for the success of breeding programs. In this study, GBS approach was used to reveal a new structuration of genetic diversity and population structure of pigeonpea in Benin. We used a total of 688 high-quality Single Nucleotide Polymorphism markers for a total of 44 pigeonpea genotypes. The distribution of SNP markers on the 11 chromosomes ranged from 14 on chromosome 5 to 133 on chromosome 2. The Polymorphism Information Content and gene diversity values were 0.30 and 0.34 respectively. The analysis of population structure revealed four clear subpopulations. The Weighted Neighbor Joining tree agreed with structure analyses by grouping the 44 genotypes into four clusters. The PCoA revealed that genotypes from subpopulations 1, 2 and 3 intermixed among themselves. The Analysis of Molecular Variance showed 7% of the total variation among genotypes while the rest of variation (93%) was within genotypes from subpopulations indicating a high gene exchange (Nm = 7.13) and low genetic differentiation (PhiPT = 0.07) between subpopulations. Subpopulation 2 presented the highest mean values of number of different alleles (Na = 1.57), number of loci with private alleles (Pa = 0.11) and the percentage of polymorphic loci (P = 57.12%). We discuss our findings and demonstrate how the genetic diversity and the population structure of this specie can be used through the Genome Wide Association Studies and Marker-Assisted Selection to enhance genetic gain in pigeonpea breeding programs in Benin.

Genetic Diversity and Structure of Quercus hondae, a Rare Evergreen Oak Species in Southwestern Japan

Conservation of rare species is essential for maintaining ecosystem function. Quercus hondae is a rare evergreen oak species (Cyclobalanopsis) endemic to Japan. This species is found in several locations in Southwestern Japan; small populations remain in the tutelary forests of the Japanese shrine. To evaluate the genetic diversity and phylogeographic structure of this rare species, 11 microsatellite loci and chloroplast DNA sequences are analyzed for 12 populations of Q. hondae and 8 populations of the more widespread congeneric species, Q. glauca. It is found that heterozygosity at both the population and species level is substantially lower in Q. hondae than in Q. glauca. Genetic differentiation among populations of Q. hondae was high, in contrast to Q. glauca, in which populations exhibit largely insignificant differentiation. STRUCTURE analysis shows that at K = 7, the clusters largely corresponded to major predefined populations. This study suggests that there is little gene flow among extant Q. hondae populations and that Q. hondae is genetically differentiated due to the greater effect of genetic drift in small populations. This pattern is in sharp contrast to that of a more common congeneric species, which will be an important consideration in the conservation of Q. hondae.

Gene Pool Boundaries for the Yosemite Toad (Anaxyrus canorus) Reveal Asymmetrical Migration Within Meadow Neighborhoods

The Yosemite toad (Anaxyrus [Bufo] canorus) is a federally threatened species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. The species is one of the first amphibians to undergo a large demographic collapse that was well-documented, and is reputed to remain in low abundance throughout its range. Recent phylogeographic work has demonstrated that Pleistocene toad lineages diverged and then admixed to differing extents across an elevational gradient. Although lineage divisions may have significant effects on evolutionary trajectories over large spatial and temporal scales, present-day population dynamics must be delineated in order to manage and conserve the species effectively. In this study, we used a double-digest RADseq dataset to address three primary questions: (1) Are single meadows or neighborhoods of nearby meadows most correlated with population boundaries? (2) Does asymmetrical migration occur among neighborhoods of nearby meadows? (3) What topographic or hydrological variables predict such asymmetrical migration in these meadow neighborhoods? Hierarchical STRUCTURE and AMOVA analyses suggested that populations are typically circumscribed by a single meadow, although 84% of meadows exist in neighborhoods of at least two meadows connected by low levels of migration, and over half (53%) of neighborhoods examined display strong asymmetrical migration. Meadow neighborhoods often contain one or more large and flat “hub” meadows that experience net immigration, surrounded by smaller and topographically rugged “satellite” meadows with net emigration. Hubs tend to contain more genetic diversity and could be prioritized for conservation and habitat management and as potential sources for reestablishment efforts.

Evaluación de la utilidad de marcadores microsatélites en la población avícola rustipollos

Los marcadores moleculares son una herramienta de gran utilidad para estudios de diversidad genética, que permite identificar poblaciones con características genéticas particulares, que soportan el establecimiento de programas de conservación y mejoramiento genético. El objetivo de este estudio fue evaluar el grado de información generada por un panel de 30 marcadores microsatélites en la población avícola Rustipollos. Se obtuvieron muestras de sangre de 50 individuos, la amplificación de fragmentos se realizó mediante PCR, utilizando 30 microsatélites recomendados por la FAO-ISAG para estudios de biodiversidad en gallinas. La estimación de los tamaños de los fragmentos se realizó en un secuenciador automático ABI Prism 377. Fueron determinados el número de alelos por locus y el Contenido de Información Polimórfica (PIC), mediante el programa Microsatellite-Toolkit. El número total de alelos reportados fue de 99 en los 30 marcadores microsatélites, con un valor medio de 3.3 ±1.06 alelos por locus. La determinación del PIC registró un promedio de 0.46, con un rango de 0.18 a 0.76 en los marcadores MCW016 y ADL278, respectivamente. El 43% de los marcadores empleados resultaron altamente informativos para la población evaluada. En general, los marcadores microsatélites demostraton ser útiles para estudios genéticos en la población avícola Rustipollos.

Genetic Diversity and Population Structure of Metaphire vulgaris Based on the Mitochondrial COI Gene and Microsatellites

The earthworm species Metaphire vulgaris (a member of the Clitellata class) is widely distributed across China, and has important ecological functions and medicinal value. However, investigations into its genetic diversity and differentiation are scarce. Consequently, we evaluated the genetic diversity of five populations of M. vulgaris (GM, HD, NYYZ, QDDY, and QDY) in Yancheng, China via the mitochondrial COI gene and the novel microsatellites developed there. A total of nine haplotypes were obtained by sequencing the mitochondrial COI gene, among which NYYZ and QDDY populations had the greatest number of haplotypes (nh = 5). Further, the nucleotide diversity ranged from 0.00437 to 0.1243. The neighbor-joining trees and the TCS network of haplotypes indicated that earthworm populations within close geographical range were not genetically isolated at these small scale distances. Results of the identification of microsatellite molecular markers revealed that the allele number in 12 microsatellite loci ranged from 4 to 13. The observed heterozygosity ranged from 0.151 to 0.644, whereas the expected heterozygosity ranged from 0.213 to 0.847. The polymorphism data content of most sites was >0.5, which indicated that the designed sites had high polymorphism. Structural analysis results indicated that GM, HD, and NYYZ had similar genetic structures across the five populations. The Nei's genetic distance between HD and NYYZ populations was the smallest (D s = 0.0624), whereas that between HD and QDY populations was the largest (D s = 0.2364). The UPGMA tree showed that HD were initially grouped with NYYZ, followed by GM, and then with QDDY. Furthermore, cross-species amplification tests were conducted for Metaphire guillelmi, which indicated that the presented markers were usable for this species. This study comprised a preliminary study on the genetic diversity of M. vulgaris, which provides basic data for future investigations into this species.

Temporal dynamics of QTL effects on vegetative growth in Arabidopsis thaliana

We assessed early vegetative growth in a population of 382 accessions of Arabidopsis thaliana using automated non-invasive high-throughput phenotyping. All accessions were imaged daily from 7 d to 18 d after sowing in three independent experiments and genotyped using the Affymetrix 250k SNP array. Projected leaf area (PLA) was derived from image analysis and used to calculate relative growth rates (RGRs). In addition, initial seed size was determined. The generated datasets were used jointly for a genome-wide association study that identified 238 marker-trait associations (MTAs) individually explaining up to 8% of the total phenotypic variation. Co-localization of MTAs occurred at 33 genomic positions. At 21 of these positions, sequential co-localization of MTAs for 2-9 consecutive days was observed. The detected MTAs for PLA and RGR could be grouped according to their temporal expression patterns, emphasizing that temporal variation of MTA action can be observed even during the vegetative growth phase, a period of continuous formation and enlargement of seemingly similar rosette leaves. This indicates that causal genes may be differentially expressed in successive periods. Analyses of the temporal dynamics of biological processes are needed to gain important insight into the molecular mechanisms of growth-controlling processes in plants.

Addressing the Linnean shortfall in a cryptic species complex

Biodiversity is being lost at an alarming rate, but the rate of this loss is likely to be underestimated as a result of a deficit in taxonomic knowledge (i.e. the Linnean shortfall). This knowledge gap is more extensive for morphologically indistinct taxa. The advancement of molecular techniques and delimitation methods has facilitated the identification of such cryptic species, but a majority of these taxa remain undescribed. To investigate the effects of taxonomic uncertainty on understanding of biodiversity, we applied the general lineage concept of species to an amphipod species complex, the Gammaruslacustris lineage that occupies springs of the northern Chihuahuan Desert, which is emerging in contemporary times. We investigated species boundaries using a validation-based approach and examined genetic structure of the lineage using a suite of microsatellite markers to identify independently evolving metapopulations. Our results show that each spring contains a genetically distinct population that is geographically isolated from other springs, suggesting evolutionary independence and status as separate species. Additionally, we observed subtle interspecific morphological variation among the putative species. We used multiple lines of evidence to formally describe four new species (Gammarus langi sp. nov., G. percalacustris sp. nov., G. colei sp. nov. and G. malpaisensis sp. nov.) endemic to the northern Chihuahuan Desert. Cryptic speciation is likely to be high in other aquatic taxa within these ecosystems, and across arid landscapes throughout North America and elsewhere, suggesting that the magnitude of the Linnean shortfall is currently underestimated in desert springs worldwide.

Genetic diversity and population structure of Camellia huana (Theaceae), a limestone species with narrow geographic range, based on chloroplast DNA sequence and microsatellite markers

Camellia huana is an endangered species with a narrow distribution in limestone hills of northern Guangxi and southern Guizhou provinces, China. We used one chloroplast DNA (cpDNA) fragment and 12 pairs of microsatellite (simple sequence repeat; SSR) markers to assess the genetic diversity and structure of 12 C. huana populations. A total of 99 alleles were detected for 12 polymorphic loci, and eight haplotypes and nine polymorphic sites were detected within 5200 bp of cpDNA. C. huana populations showed a low level of genetic diversity (n = 8, Hd = 0.759, Pi = 0.00042 for cpDNA, N _A  = 3.931, H _E  = 0.466 for SSRs), but high genetic differentiation between populations (F _ST  = 0.2159 for SSRs, F _ST  = 0.9318 for cpDNA). This can be attributed to the narrow distribution and limestone habitat of C. huana. STRUCTURE analysis divided natural C. huana populations into two groups, consistent with their geographical distribution. Thus, we suggest that five natural C. huana populations should be split into two units to be managed effectively.

Characterization of 25 new microsatellite markers for the fin whale (Balaenoptera physalus) and cross-species amplification in other cetaceans

Cetaceans are large mammals widely distributed on Earth. The fin whale, Balaenoptera physalus, is the second largest living animal. In the 20th century, commercial whaling reduced its global population by 70%, and in the Mediterranean Sea not only was their overall population depleted but the migration between the Mediterranean Sea and the Atlantic Ocean was reduced. Previous genetic studies identified isolation between these two regions, with a limited gene-flow between these adjacent populations based on nuclear and mitochondrial markers. However, only limited information exists for the Mediterranean population as genetic diversity and abundance trends are still unknown. In this study, 39 highly polymorphic microsatellite markers were tested, including 25 markers developed de novo together with 14 markers previously published. An average allelic diversity of 8.3 alleles per locus was reported, ranging from 3 to 15 alleles per locus, for B. physalus. Expected heterozygosity was variable among loci and ranged from 0.34 to 0.91. Only two markers in the new set were significantly deviant from the Hardy Weinberg equilibrium. Cross-species amplification was tested in four other cetacean species. A total of 27 markers were successfully amplified in the four species (Balaenoptera acutorostrata, Megaptera novaeangliae, Physeter macrocephalus and Globicephala melas). A multivariate analysis on the multilocus genotypes successfully discriminated the five species. This new set of microsatellite markers will not only provide a useful tool to identify and understand the genetic diversity and the evolution of the B. physalus population, but it will also be relevant for other cetacean species, and will allow further parentage analyses. Eventually, this new set of microsatellite markers will provide critical data that will shed light on important biological data within a conservation perspective.

Mining whole genome sequence data to efficiently attribute individuals to source populations

Whole genome sequence (WGS) data could transform our ability to attribute individuals to source populations. However, methods that efficiently mine these data are yet to be developed. We present a minimal multilocus distance (MMD) method which rapidly deals with these large data sets as well as methods for optimally selecting loci. This was applied on WGS data to determine the source of human campylobacteriosis, the geographical origin of diverse biological species including humans and proteomic data to classify breast cancer tumours. The MMD method provides a highly accurate attribution which is computationally efficient for extended genotypes. These methods are generic, easy to implement for WGS and proteomic data and have wide application.

Genetic diversity among early provitamin A quality protein maize inbred lines and the performance of derived hybrids under contrasting nitrogen environments

BACKGROUND

Information on population structure and genetic diversity of germplasm in a breeding programme is useful because it enhances judicious utilisation of genetic resources to achieve breeding objectives. Seventy early maturing provitamin A (PVA) quality protein maize (QPM) inbreds developed by the IITA- maize improvement programme were genotyped using 8171 DArTseq markers. Furthermore, 96 hybrids derived from 24 selected inbreds plus four checks were evaluated under low-N and optimal environments in Nigeria during 2016 and 2017. Genotypic and phenotypic data of inbreds and hybrids respectively, were analysed to (i) assess the level of genetic dissimilarities and population structure of the inbreds, and (ii) investigate the grain yield performance of derived hybrids under low-N, optimal and across environments.

RESULTS

Genetic diversity among the seventy inbreds was high varying from 0.042 to 0.500 with an average of 0.357. Sixty-six inbred lines with probabilities ≥0.70 were assigned to a single group. The population structure analysis, the UPGMA phylogeny, and the principal Coordinate Analysis (PCoA) of the DArTseq markers revealed a clear separation of five groups and each followed pedigree records. Clustered inbreds displayed common characteristics including high PVA levels, and drought and low-N tolerance. The top performing hybrid, TZEIORQ 40 × TZEIORQ 26 out-yielded the best hybrid control, TZEIOR 127 × TZEIOR 57 by 8, 3, and 9% under low-N, optimal, and across environments, respectively. High repeatability estimates were detected for grain yield under each and across environments. Similarly, high breeding efficiency of 71, 70 and 72% were computed under low-N, optimal, and across environments, respectively.

CONCLUSIONS

The UPGMA clustering, the structure analysis, and the PCoA consistently revealed five groups which largely followed pedigree information indicating the existence of genetically distinct groups in the inbred lines. High repeatability and breeding efficiency values estimated for grain yield of hybrids under low-N, optimal and across environments demonstrated that high productive hybrids could be developed using inbreds from the opposing clusters identified by the DArTseq markers. The 15 top performing hybrids identified, particularly TZEIORQ 40 × TZEIORQ 26 and TZEIORQ 29 × TZEIORQ 43 should be further evaluated for release and commercialization in SSA.

Low genetic diversity and limited genetic structure across the range of the critically endangered Mexican howler monkey (Alouatta palliata mexicana)

Genetic diversity provides populations with the possibility to persist in ever-changing environments, where selective regimes change over time. Therefore, the long-term survival of a population may be affected by its level of genetic diversity. The Mexican howler monkey (Alouatta palliata mexicana) is a critically endangered primate restricted to southeast Mexico. Here, we evaluate the genetic diversity and population structure of this subspecies based on 83 individuals from 31 groups sampled across the distribution range of the subspecies, using 29 microsatellite loci. Our results revealed extremely low genetic diversity (H_O  = 0.21, H_E  = 0.29) compared to studies of other A. palliata populations and to other Alouatta species. Principal component analysis, a Bayesian clustering method, and analyses of molecular variance did not detect strong signatures of genetic differentiation among geographic populations of this subspecies. Although we detect small but significant F_ST values between populations, they can be explained by a pattern of isolation by distance. These results and the presence of unique alleles in different populations highlight the importance of implementing conservation efforts in multiple populations across the distribution range of A. p. mexicana to preserve its already low genetic diversity. This is especially important given current levels of population isolation due to the extreme habitat fragmentation across the distribution range of this primate.

Genetic Diversity and Population Structure of Bermudagrass [Cynodon dactylon (L.) Pers.] along Latitudinal Gradients and the Relationship with Polyploidy Level

Understanding the population genetic pattern and process of gene flow requires a detailed knowledge of how landscape characteristics structure populations. Although Cynodon dactylon (L.) Pers. (common bermudagrass) is widely distributed in the world, information on its genetic pattern and population structure along latitudinal gradients is limited. We tried to estimate the genetic diversity and genetic structure of C. dactylon along a latitudinal gradient across China. Genetic diversity among different ploidy levels was also compared in the study. The material used consisted of 296 C. dactylon individuals sampled from 16 geographic sites from 22°35′ N to 36°18′ N. Genetic diversity was estimated using 153 expressed sequence tag-derived simple sequence repeat (EST-SSR) loci. Higher within-population genetic diversity appeared at low-latitude, as well as having positive correlation with temperature and precipitation. The genetic diversity increased with the ploidy level of C. dactylon, suggesting polyploidy creates higher genetic diversity. No isolation by distance and notable admixture structure existed among populations along latitudes. Both seed dispersal (or vegetative organs) and extrinsic pollen played important roles for gene flow in shaping the spatial admixture population structure of C. dactylon along latitudes. In addition, populations were separated into three clusters according to ploidy levels. C. dactylon has many such biological characters of perennial growth, wind-pollination, polyploidy, low genetic differentiation among populations, sexual and asexual reproduction leading to higher genetic diversity, which gives it strong adaptability with its genetic patterns being very complex across all the sampled latitudes. The findings of this study are related to landscape population evolution, polyploidy speciation, preservation, and use of bermudagrass breeding.

Assessment of morphological and genetic variability through genic microsatellite markers for essential oil in Sandalwood (Santalum album L.)

Sandalwood (Santalum album L; family Santalaceae) is a highly significant aromatic oil yielding tree. It is valued for two important traits, heartwood and essential oil obtained from the heartwood. This study was proposed to assess the morphological and genetic variability of sandalwood accessions. For this, genotypes were randomly selected (n = 177) from the 14 populations from three states in southern India. The total heartwood oil content and quality was estimated by UV method and GC-MS. Total 14 oil-specific genic SSR markers were procured to evaluate the genetic diversity among the sandalwood accessions. Total core size, heartwood content, and oil of S. album ranged from 4.4 to 19.1 cm; 0.0 to 17.3 cm; and 0.0 to 5.96% with covariance 27.61, 85.25, and 73.12% followed by mean 9.74, 3.77, and 2.71, respectively. Genetic diversity estimates were highly polymorphic in terms of Na 7.28, Ne 5.89, He 8.0 PIC 0.891, with little Ho, and F-0.922. AMOVA revealed that minimal genetic variation among populations and highest variation was found among individuals with Nm (58.4). The UPGMA reveals the cluster favored the grouping pattern by the PCA analysis. Structure and PCA analysis clustered the entire populations into two major groups with F _ST 0.046 in which population of Kerala and Karnataka were pure and Telangana accessions were found admixtures. No significant correlation (r ² = 0.23, P = 0.00) was observed between heartwood oil and genetic structures. A high degree of transferability of genic markers would facilitate the assessment of novel genotypes for future tree improvement and conservation of Sandalwood populations.

Antioxidant properties and structured biodiversity in a diverse set of wild cranberry clones

Wild germplasm with elevated antioxidants are a useful resource for using directly and in a breeding program. In a study with 136 wild clones and two cranberry cultivars, phenolic, flavonoid and anthocyanin contents varied 2.79, 2.70 and 17.46 times, respectively. The antioxidant activity ranged from 1.17 ± 0.01 to 2.53 ± 0.05 mg/g and varied 2.16 times. Seventy-five of wild clones and the cultivar Franklin were grouped into five distinct classes by molecular structure analysis using inter simple sequence repeat (ISSR), expressed sequence tag-simple sequence repeat (EST-SSR) and EST-polymerase chain reaction (PCR) markers. Grouping with DNA markers did not coincide with that of based on antioxidant properties. Present study indicates that genetic diversity analysis combined with antioxidant properties of wild germplasm play a significant role for conservation and in selecting diverse genotypes for future berry crop improvement.

Direct-to-Consumer Genetic Testing's Red Herring: "Genetic Ancestry" and Personalized Medicine

The growth in the direct-to-consumer genetic testing industry poses a number of challenges for healthcare practice, among a number of other areas of concern. Several companies providing this service send their customers reports including information variously referred to as genetic ethnicity, genetic heritage, biogeographic ancestry, and genetic ancestry. In this article, we argue that such information should not be used in healthcare consultations or to assess health risks. Far from representing a move toward personalized medicine, use of this information poses risks both to patients as individuals and to racialized ethnic groups because of the way it misrepresents human genetic diversity.

Genetic Differentiation among Commercial Lines of Laying-type Japanese Quail

Recently in Japan, approximately six million quails were primarily being reared for commercial egg production. It is believed that almost all commercial quails in the country became extinct during World War II, and that the present commercial gene pool was restored from the limited number of surviving birds. The present study evaluates the genetic diversity within and differentiation between 12 laying-type Japanese quail lines on the basis of 45 microsatellite genotypes. The mean number of alleles per locus and the expected heterozygosity within a quail line were 5.22–5.69 and 0.601–0.618, respectively. These results showed that laying-type quail lines in the present study exhibited a higher degree of genetic diversity than experimental quail lines in a previous study. Pairwise genetic differentiations (F_ST) between lines were significant but weak (F_ST=0.0028–0.0254; 57.6%), and no significant differentiations were found between the remainder. This was also confirmed by genetic clustering analyses, in which individuals did not form independent clusters consistent with their line origins. The results of the present study indicate relatively high genetic diversity within and no clear genetic differentiation between laying-type quail lines. Absence of genetic differentiation may reflect the breeding history of laying-type quails.

Genetic diversity of 21 experimental chicken lines with diverse origins and genetic backgrounds

The genetic characteristics and diversity of 21 experimental chicken lines registered with the National BioResource Project of Japan were examined using mitochondrial D-loop sequences and 54 microsatellite DNA markers. A total of 12 haplotypes were detected in the 500-bp mitochondrial DNA sequences of the hypervariable segment I for 349 individuals of 21 lines. The 12 haplotypes belonged to three (A, D, and E) haplogroups, out of the eight (A‒H) common haplogroups in domestic chickens and red junglefowls. The haplogroups A and D were widely represented in indigenous chickens in the Asian and Pacific regions, and the haplogroup E was the most prevalent in domestic chickens. Genetic clustering by discriminant analysis of principal components with microsatellite markers divided 681 individuals of 21 lines into three groups that consisted of Fayoumi-, European-, and Asian- derived lines. In each of the cladograms constructed with Nei's genetic distances based on allele frequencies and the membership coefficients provided by STRUCTURE and with the genetic distance based on the proportion of shared alleles, the genetic relationships coincided well with the breeding histories of the lines. Microsatellite markers showed remarkably lower genetic heterozygosities (less than 0.1 observed heterozygosity) for eight lines (GSP, GSN/1, YL, PNP, BM-C, WL-G, BL-E, and #413), which have been maintained as closed colonies for more than 40 years (except for #413), indicating their usefulness as experimental chicken lines in laboratory animal science research.

A tutorial on how not to over-interpret STRUCTURE and ADMIXTURE bar plots

Genetic clustering algorithms, implemented in programs such as STRUCTURE and ADMIXTURE, have been used extensively in the characterisation of individuals and populations based on genetic data. A successful example is the reconstruction of the genetic history of African Americans as a product of recent admixture between highly differentiated populations. Histories can also be reconstructed using the same procedure for groups that do not have admixture in their recent history, where recent genetic drift is strong or that deviate in other ways from the underlying inference model. Unfortunately, such histories can be misleading. We have implemented an approach, badMIXTURE, to assess the goodness of fit of the model using the ancestry “palettes” estimated by CHROMOPAINTER and apply it to both simulated data and real case studies. Combining these complementary analyses with additional methods that are designed to test specific hypotheses allows a richer and more robust analysis of recent demographic history.

Clustering methods such as STRUCTURE and ADMIXTURE are widely used in population genetic studies to investigate ancestry. Here, the authors provide a tutorial on how to interpret results of these analyses and a tool to test the goodness of fit of the model.

Possible mathematical definitions of the biological term “breed”

Abstract. For scientific discussions it is necessary to clearly define the terms used. Otherwise, scientific statements are open to interpretation and remain unclear which hampers scientific progress. A clear specification of scientific terms can be reached using tools such as mathematical definitions. So far, the biological term “breed” lacks a mathematical definition. Based on common biological literature, four mathematical definitions of the term “breed” for gonochoric species are proposed in this paper. These definitions were tested on the dog breed “whippet”. It was discussed if bully whippets and their offspring can be considered to be part of the whippet dog breed. Although all mathematical definitions proposed are consistent with common word-based definitions, the results of the bully whippet test differ. According to two definitions bully whippets may be considered part of the whippet breed. However, according to the other definitions bully whippets are not part of the whippet breed. Under certain conditions the introduced mathematical structure may also be used to define other biological terms like “variety” or “species”.

Distinguishing industrial meat from that of indigenous chickens with molecular markers

The aim of investigation was to evaluate a traceability system to detect industrial chicken meat among indigenous products, considering issues that could affect assignment accuracy. The dataset included 2 Italian indigenous meat breeds, namely Bionda Piemontese (2 ecotypes) and Bianca di Saluzzo, one broiler line, and 3 layer lines. Assignment tests were performed using a standard panel of 28 microsatellite loci. To evaluate effects of inbreeding and substructure on assignment accuracy, a simulated dataset was prepared. Broilers and layers belong to homogeneous populations and never enter the clusters of indigenous breeds. Ambiguity or misallocation are expected between the Bionda ecotypes and between the 2 indigenous breeds, but it is unlikely that niche products provided by Bionda and Bianca will compete with one another. Non-random mating reduces accuracy, but only populations having weak genetic differentiation are involved, namely those that are less interesting to discriminate. The dataset can be used as a reference population to distinguish commercial meat from indigenous meat with great accuracy. Misallocations increase as number of loci decreases, but only within or between the indigenous breeds. A subpanel of the most resolving 14 loci keeps sufficient informative content to provide accuracy and to correctly allocate additional test samples within the reference population. This analytical tool is economically sustainable as a method to detect fraud or mislabeling. Adoption of a monitoring system should increase the value of typical products because the additional burden of molecular analyses would improve commercial grade and perception of quality.

Genetic Diversity Analysis Reveals Genetic Differentiation and Strong Population Structure in Calotropis Plants

The genus Calotropis (Asclepiadaceae) is comprised of two species, C. gigantea and C. procera, which both show significant economic potential for use of their seed fibers in the textile industry, and of their bioactive compounds as new medicinal resources. The available wild-sourced germplasm contains limited genetic information that restricts further germplasm exploration for the purposes of domestication. We here developed twenty novel EST-SSR markers and applied them to assess genetic diversity, population structure and differentiation within Calotropis. The polymorphic information index of these markers ranged from 0.102 to 0.800; indicating that they are highly informative. Moderate genetic diversity was revealed in both species, with no difference between species in the amount of genetic diversity. Population structure analysis suggested five main genetic groups (K = 5) and relatively high genetic differentiation (F_ST = 0.528) between the two species. Mantel test analysis showed strong correlation between geographical and genetic distance in C. procera (r = 0.875, p = 0.020) while C. gigantea showed no such correlation (r = 0.390, p = 0.210). This study provides novel insights into the genetic diversity and population structure of Calotropis, which will promote further resource utilization and the development of genetic improvement strategies for Calotropis.

Origin of Pest Lineages of the Colorado Potato Beetle (Coleoptera: Chrysomelidae)

Colorado potato beetle (Leptinotarsa decemlineata Say [Coleoptera: Chrysomelidae]) is a pest of potato throughout the Northern Hemisphere, but little is known about the beetle's origins as a pest. We sampled the beetle from uncultivated Solanum host plants in Mexico, and from pest and non-pest populations in the United States and used mitochondrial DNA and nuclear loci to examine three hypotheses on the origin of the pest lineages: 1) the pest beetles originated from Mexican populations, 2) they descended from hybridization between previously divergent populations, or 3) they descended from populations that are native to the Plains states in the United States. Mitochondrial haplotypes of non-pest populations from Mexico and Arizona differed substantially from beetles collected from the southern plains and potato fields in the United States, indicating that beetles from Mexico and Arizona did not contribute to founding the pest lineages. Similar results were observed for AFLP and microsatellite data . In contrast, non-pest populations from the states of Colorado, Kansas, Nebraska, New Mexico, and Texas were genetically similar to U.S. pest populations, indicating that they contributed to the founding of the pest lineages. Most of the pest populations do not show a significant reduction in genetic diversity compared to the plains populations in the United States. We conclude that genetically heterogeneous beetle populations expanded onto potato from native Solanum hosts. This mode of host range expansion may have contributed to the abundant genetic diversity of contemporary populations, perhaps contributing to the rapid evolution of climate tolerance, host range, and insecticide resistance.

Trait Mapping Approaches Through Association Analysis in Plants

Previously, association mapping (AM) methodology was used to unravel genetic complications in animal science by measuring the complex traits for candidate and non-candidate genes. Nowadays, this statistical approach is widely used to clarify the complexity in plant breeding program-based genome-wide breeding strategies, marker development, and diversity analysis. This chapter is particularly focused on methodologies with limitations and provides an overview of AM models and software used up to now. Association or linkage disequilibrium mapping has become a very popular method for discovering candidate and non-candidate genes and confirmation of quantitative trait loci (QTL) on various parts of the genome and in marker-assisted selection for breeding. Previously, various QTL investigations were carried out for different plants exclusively by linkage mapping. To help to understand the basics of modern molecular genetic techniques, in this chapter we summarize previous studies done on different crops. AM offers high-resolution power when there is large genotypic diversity and low linkage disequilibrium (LD) for the germplasm being investigated. The benefits of AM, compared with traditional QTL mapping, include a relatively detailed mapping resolution and a far less time-consuming approach since no mapping populations need to be generated. The advancements in genotyping and computational techniques have encouraged the use of AM. AM provides a fascinating approach for genetic investigation of QTLs, due to its resolution and the possibility to study the various genomic areas at the same time without construction of mapping populations. In this chapter we also discuss the advantages and disadvantages of AM, especially in the dicotyledonous crops Fabaceae and Solanaceae, with various genome-size reproductive strategies (clonal vs. sexual), and statistical models. The main objective of this chapter is to highlight the uses of association genetics in major and minor crop species that have trouble being analyzed for dissection of complex traits by identification of the factor responsible for controlling the effect of trait. Graphical Abstract.

Low levels of hybridization in two species of African driver ants

Hybridization in ants can have consequences different from those observed in most other species, with many of the potential deleterious effects being mitigated due to haplodiploidy and eusociality. In some species where colonies are either headed by multiple queens or single queens that mate with many males, hybridization is associated with genetic caste determination, where hybrids develop into workers and purebred individuals develop into queens. A previous study suggested that hybridization occurs between two Dorylus army ant species with multiply mated queens. However, the extent and exact pattern of hybridization have remained unclear, and its possible effect on caste determination has not been investigated. In this study, we aimed to determine the extent and direction of hybridization by measuring how frequently hybrids occur in colonies of both species, and to investigate the possibility of genetic caste determination. We show that hybridization is bidirectional and occurs at equal rates in both species. Hybrid workers make up only 1-2% of the population, and successful interspecific matings represent approximately 2% of all matings in both species. This shows that, although interspecific matings that give rise to worker offspring occur regularly, they are much rarer than intraspecific mating. Finally, we find no evidence of an association between hybridization and genetic caste determination in this population. This means that genetic caste determination is not a necessary outcome of hybridization in ants, even in species where queens mate with multiple males.

The 22Rv1 prostate cancer cell line carries mixed genetic ancestry: Implications for prostate cancer health disparities research using pre-clinical models

BACKGROUND

Understanding how biological factors contribute to prostate cancer (PCa) health disparities requires mechanistic functional analysis of specific genes or pathways in pre-clinical cellular and animal models of this malignancy. The 22Rv1 human prostatic carcinoma cell line was originally derived from the parental CWR22R cell line. Although 22Rv1 has been well characterized and used in numerous mechanistic studies, no racial identifier has ever been disclosed for this cell line. In accordance with the need for racial diversity in cancer biospecimens and recent guidelines by the NIH on authentication of key biological resources, we sought to determine the ancestry of 22RV1 and authenticate previously reported racial identifications for four other PCa cell lines.

METHODS

We used 29 established Ancestry Informative Marker (AIM) single nucleotide polymorphisms (SNPs) to conduct DNA ancestry analysis and assign ancestral proportions to a panel of five PCa cell lines that included 22Rv1, PC3, DU145, MDA-PCa-2b, and RC-77T/E.

RESULTS

We found that 22Rv1 carries mixed genetic ancestry. The main ancestry proportions for this cell line were 0.41 West African (AFR) and 0.42 European (EUR). In addition, we verified the previously reported racial identifications for PC3 (0.73 EUR), DU145 (0.63 EUR), MDA-PCa-2b (0.73 AFR), and RC-77T/E (0.74 AFR) cell lines.

CONCLUSIONS

Considering the mortality disparities associated with PCa, which disproportionately affect African American men, there remains a burden on the scientific community to diversify the availability of biospecimens, including cell lines, for mechanistic studies on potential biological mediators of these disparities. This study is beneficial by identifying another PCa cell line that carries substantial AFR ancestry. This finding may also open the door to new perspectives on previously published studies using this cell line.

Genetic diversity of Afrikaner cattle in southern Africa

The Afrikaner is an indigenous South African breed of "Sanga" type beef cattle along with breeds such as the Drakensberger and Nguni. Six composite breeds have been developed from crosses with the Afrikaner. Additionally, Afrikaner has been the base from which exotic breeds were established in South Africa through backcrossing. The study examined genetic diversity of Afrikaner cattle by genotyping 1257 animals from 27 herds in different geographic areas of South Africa and Namibia using 11 microsatellite markers. Multiple-locus assignment, performed using the Bayesian clustering algorithm of STRUCTURE, revealed three underlying genotypic groups. These groups were not geographically localized. Across herds and markers, the proportion of unbiased heterozygosity ranged from 0.49 to 0.72 averaging 0.57; mean number of alleles per locus ranged from 3.18 to 7.09, averaging 4.81; and allelic richness ranged from 2.35 to 3.38, averaging 2.67. It is concluded that a low inbreeding level of 2.7% and a moderate to high degree of variation still persists within the Afrikaner cattle breed, despite the recent decline in numbers of animals.

Isolation and gene flow in a speciation continuum in newts

Because reproductive isolation often evolves gradually, differentiating lineages may retain the potential for genetic exchange for prolonged periods, providing an opportunity to quantify and to understand the fundamental role of gene flow during speciation. Here we delimit evolutionary lineages, reconstruct the phylogeny and infer gene flow in newts of the Lissotriton vulgaris species complex based on 74 nuclear markers sampled from 127 localities. We demonstrate that distinct lineages along the speciation continuum in newts exchange nontrivial amounts of genes, affecting their evolutionary trajectories. By integrating a wide array of methods, we delimit nine evolutionary lineages and show that two principal factors have driven their genetic differentiation: time since the last common ancestor determining levels of shared ancestral polymorphism, and shifts in geographic distributions determining the extent of secondary contact. Post-divergence gene flow, indicative of evolutionary non-independence, has been most extensive in Central Europe, while four southern European lineages have acquired the population-genetic hallmarks of independent species (L. graecus, L. kosswigi, L. lantzi, L. schmidtleri). We obtained strong statistical support for widespread mtDNA introgression following secondary contact, previously suggested by discordance between mtDNA phylogeny and morphology. Our study reveals long-term evolutionary persistence of evolutionary lineages that may periodically exchange genes with one another: although some of these lineages may become extinct or fuse, others will acquire complete reproductive isolation and will carry signatures of this complex history in their genomes.

Intercontinental and intracontinental biogeography of the eastern Asian - Eastern North American disjunct Panax (the ginseng genus, Araliaceae), emphasizing its diversification processes in eastern Asia

The intercontinental biogeography between eastern Asia and eastern North America has attracted much attention from evolutionary biologists. Further insights into understanding the evolution of the intercontinental disjunctions have been hampered by the lack of studies on the intracontinental biogeography in eastern Asia, a region with complex geology, geography, climates and habitats. Herein we studied the biogeographic history of the eastern Asian-eastern North American disjunct genus Panax with special emphasis on the investigation of its uneven diversification in Asia. This study reconstructs the diversification history of Panax and also emphasizes a large clade of Panax taxa, which has a wide distribution in eastern Asia, but was unresolved in previous studies. We examined the noncoding plastid DNA fragments of trnH-psbA, rps16, and psbM-trnD, the mitochondrial b/c intron of NAD1, and the nuclear ribosomal internal transcribed spacer (ITS) region of 356 samples from 47 populations. The results revealed the subtropical Northern Hemisphere origin (Asia or Asia and North America) of Panax in the Paleocene. Intercontinental disjunctions between eastern Asia and eastern North America formed twice in Panax, once estimated in early Eocene for the split of P. trifolius and another in mid-Miocene for the divergence of P. quinquefolius. Intercontinental diversifications in Panax showed temporal correlation with the increase of global temperature. The evolutionary radiation of the P. bipinnatifidus species complex occurred around the boundary of Oligocene and Miocene. Strong genetic structure among populations of the species complex was detected and the populations may be isolated by distance. The backbone network and the Bayesian clustering analysis revealed a major evolutionary radiation centered in the Hengduan Mountains of western China. Our results suggested that the evolutionary radiation of Panax was promoted by geographic barriers, including mountain ranges (Hengduan Mountains, Nanling Mountains and Wuyishan Mountains), oceans and altitudinal shifts, which further contribute to the knowledge of the uneven species diversification between eastern Asia and North America.

Genetic structure and gene flow of the flea Xenopsylla cheopis in Madagascar and Mayotte

Background

The flea Xenopsylla cheopis (Siphonaptera: Pulicidae) is a vector of plague. Despite this insect’s medical importance, especially in Madagascar where plague is endemic, little is known about the organization of its natural populations. We undertook population genetic analyses (i) to determine the spatial genetic structure of X. cheopis in Madagascar and (ii) to determine the potential risk of plague introduction in the neighboring island of Mayotte.

Results

We genotyped 205 fleas from 12 sites using nine microsatellite markers. Madagascan populations of X. cheopis differed, with the mean number of alleles per locus per population ranging from 1.78 to 4.44 and with moderate to high levels of genetic differentiation between populations. Three distinct genetic clusters were identified, with different geographical distributions but with some apparent gene flow between both islands and within Malagasy regions. The approximate Bayesian computation (ABC) used to test the predominant direction of flea dispersal implied a recent population introduction from Mayotte to Madagascar, which was estimated to have occurred between 1993 and 2012. The impact of this flea introduction in terms of plague transmission in Madagascar is unclear, but the low level of flea exchange between the two islands seems to keep Mayotte free of plague for now.

Conclusion

This study highlights the occurrence of genetic structure among populations of the flea vector of plague, X. cheopis, in Madagascar and suggests that a flea population from Mayotte has been introduced to Madagascar recently. As plague has not been reported in Mayotte, this introduction is unlikely to present a major concern for plague transmission. Nonetheless, evidence of connectivity among flea populations in the two islands indicates a possibility for dispersal by fleas in the opposite direction and thus a risk of plague introduction to Mayotte.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2290-6) contains supplementary material, which is available to authorized users.

Effective Dispersal of Caribbean Reef Fish is Smaller than Current Spacing Among Marine Protected Areas

The oceans are deteriorating at a fast pace. Conservation measures, such as Marine Protected Areas, are being implemented to relieve areas from local stressors and allow populations to restore to natural levels. Successful networks of MPAs operate if the space among MPAs is smaller than the dispersal capacity of the species under protection. We studied connectivity patterns across populations in a series of MPAs in the common yellowhead Jawfish, Opistognathus aurifrons. Using the power of genome-wide variation, we estimated that the maximum effective dispersal is 8.3 km. We found that MPAs exchange migrants likely via intermediate unprotected habitats through stepping stone dispersal. At scales >50 km such connectivity is decreased, particularly across the Mona Passage. The MPA network studied would be unable to maintain connectivity of these small benthic fishes if habitat in between them is extirpated. Our study highlights the power of SNPs to derive effective dispersal distance and the ability of SNPs to make inferences from single individuals. Given that overall reef fish diversity is driven by species with life histories similar to that of the yellowhead jawfish, managers face a challenge to develop strategies that allow connectivity and avoid isolation of populations and their possible extinction.

Multiple origins and a narrow genepool characterise the African tea germplasm: concordant patterns revealed by nuclear and plastid DNA markers

Despite the highly economic value of tea in Africa, its genetic and geographic origins remain largely unexplored. Here we address this by collecting 439 samples across 11 countries in Africa and Asia to investigate the origin and genepool composition of African tea based on 23 nuclear microsatellites loci (nSSRs) and three cpDNA intergenic spacer regions. Our results indicated that the African tea represents a potpourri originating from multiple introductions over time. The nSSR analysis revealed that the majority (79%) of tea accessions collected in Africa belong to Indian Assam tea which have likely originated from India and/or Sri Lanka. The patterns of nSSR variation also showed that Chinese Assam tea is genetically distinct from Indian Assam tea, and has rarely been used in African tea breeding efforts since only 4% of the African tea accessions possessed this genotype. We found a total of 22 cpDNA haplotypes, which grouped into three main geographic clades that were concordant with the distribution of microsatellite genotypes. Several private cpDNA haplotypes were identified in Chinese Assam tea in Southern Yunnan province of China. Therefore Chinese Assam tea will be important for the enrichment of African tea gene pools. Our results is a useful guide in future tea breeding programmes in Africa.

Molecular Diversity Analysis and Genetic Mapping of Pod Shatter Resistance Loci in Brassica carinata L

Seed lost due to easy pod dehiscence at maturity (pod shatter) is a major problem in several members of Brassicaceae family. We investigated the level of pod shatter resistance in Ethiopian mustard (Brassica carinata) and identified quantitative trait loci (QTL) for targeted introgression of this trait in Ethiopian mustard and its close relatives of the genus Brassica. A set of 83 accessions of B. carinata, collected from the Australian Grains Genebank, was evaluated for pod shatter resistance based on pod rupture energy (RE). In comparison to B. napus (RE = 2.16 mJ), B. carinata accessions had higher RE values (2.53 to 20.82 mJ). A genetic linkage map of an F₂ population from two contrasting B. carinata selections, BC73526 (shatter resistant with high RE) and BC73524 (shatter prone with low RE) comprising 300 individuals, was constructed using a set of 6,464 high quality DArTseq markers and subsequently used for QTL analysis. Genetic analysis of the F₂ and F_2:3 derived lines revealed five statistically significant QTL (LOD ≥ 3) that are linked with pod shatter resistance on chromosomes B1, B3, B8, and C5. Herein, we report for the first time, identification of genetic loci associated with pod shatter resistance in B. carinata. These characterized accessions would be useful in Brassica breeding programs for introgression of pod shatter resistance alleles in to elite breeding lines. Molecular markers would assist marker-assisted selection for tracing the introgression of resistant alleles. Our results suggest that the value of the germplasm collections can be harnessed through genetic and genomics tools.

Genome-Wide Association Study Identifies Candidate Loci Underlying Agronomic Traits in a Middle American Diversity Panel of Common Bean

Common bean ( L.) breeding programs aim to improve both agronomic and seed characteristics traits. However, the genetic architecture of the many traits that affect common bean production are not completely understood. Genome-wide association studies (GWAS) provide an experimental approach to identify genomic regions where important candidate genes are located. A panel of 280 modern bean genotypes from race Mesoamerica, referred to as the Middle American Diversity Panel (MDP), were grown in four US locations, and a GWAS using >150,000 single-nucleotide polymorphisms (SNPs) (minor allele frequency [MAF] ≥ 5%) was conducted for six agronomic traits. The degree of inter- and intrachromosomal linkage disequilibrium (LD) was estimated after accounting for population structure and relatedness. The LD varied between chromosomes for the entire MDP and among race Mesoamerica and Durango-Jalisco genotypes within the panel. The LD patterns reflected the breeding history of common bean. Genome-wide association studies led to the discovery of new and known genomic regions affecting the agronomic traits at the entire population, race, and location levels. We observed strong colocalized signals in a narrow genomic interval for three interrelated traits: growth habit, lodging, and canopy height. Overall, this study detected ∼30 candidate genes based on a priori and candidate gene search strategies centered on the 100-kb region surrounding a significant SNP. These results provide a framework from which further research can begin to understand the actual genes controlling important agronomic production traits in common bean.

Insights into the Genetic Relationships and Breeding Patterns of the African Tea Germplasm Based on nSSR Markers and cpDNA Sequences

Africa is one of the key centers of global tea production. Understanding the genetic diversity and relationships of cultivars of African tea is important for future targeted breeding efforts for new crop cultivars, specialty tea processing, and to guide germplasm conservation efforts. Despite the economic importance of tea in Africa, no research work has been done so far on its genetic diversity at a continental scale. Twenty-three nSSRs and three plastid DNA regions were used to investigate the genetic diversity, relationships, and breeding patterns of tea accessions collected from eight countries of Africa. A total of 280 African tea accessions generated 297 alleles with a mean of 12.91 alleles per locus and a genetic diversity (H S) estimate of 0.652. A STRUCTURE analysis suggested two main genetic groups of African tea accessions which corresponded well with the two tea types Camellia sinensis var. sinensis and C. sinensis var. assamica, respectively, as well as an admixed "mosaic" group whose individuals were defined as hybrids of F2 and BC generation with a high proportion of C. sinensis var. assamica being maternal parents. Accessions known to be C. sinensis var. assamica further separated into two groups representing the two major tea breeding centers corresponding to southern Africa (Tea Research Foundation of Central Africa, TRFCA), and East Africa (Tea Research Foundation of Kenya, TRFK). Tea accessions were shared among countries. African tea has relatively lower genetic diversity. C. sinensis var. assamica is the main tea type under cultivation and contributes more in tea breeding improvements in Africa. International germplasm exchange and movement among countries within Africa was confirmed. The clustering into two main breeding centers, TRFCA, and TRFK, suggested that some traits of C. sinensis var. assamica and their associated genes possibly underwent selection during geographic differentiation or local breeding preferences. This study represents the first step toward effective utilization of differently inherited molecular markers for exploring the breeding status of African tea. The findings here will be important for planning the exploration, utilization, and conservation of tea germplasm for future breeding efforts in Africa.

Genetic structure reveals a history of multiple independent origins followed by admixture in the allopolyploid weed Salsola ryanii

It has recently become clear that many invasive species have evolved in situ via hybridization or polyploidy from progenitors which themselves are introduced species. For species formed by hybridization or polyploidy, genetic diversity within the newly formed species is influenced by the number of independent evolutionary origins of the species. For recently formed species, an analysis of genetic structure can provide insight into the number of independent origin events involved in the formation of the species. For a putative invasive allopolyploid species, the number of origins involved in the species formation, the genetic diversity present within these origins, and the level of gene flow between independent origins determines the genetic composition of the neospecies. Here we analyze the genetic structure of the newly formed allopolyploid species, Salsola ryanii, a tumbleweed which evolved within the last 20–100 years in California. We utilize the genetic structure analysis to determine that this new species is the result of at least three independent allopolyplodization events followed by gene flow between the descendants of independent origins.

Genetic variability of two Italian indigenous chicken breeds inferred from microsatellite marker analysis

The objective of this study was to determine the genetic structure and variability of Bionda Piemontese and Bianca di Saluzzo (Piedmont, Northwest Italy) using an international set of microsatellite loci (AVIANDIV-FAO). Differences compared with commercial lines and other Italian breeds were verified to justify the implementation of conservation programmes. Flock contribution to genetic variability was assessed following the approach implemented in the MolKin software. Comparison was performed using the fixation index and the Reynolds genetic distance. The most likely number of different populations was estimated using the clustering procedure implemented in STRUCTURE. The molecular information suggests that management practices could have prevented random mating and produced inbreeding and heterogeneity across flocks. In this respect, Bionda and Bianca show substructuring and are more similar to British breeds than other continental European breeds. Bionda and Bianca fit into the European breeds provided with the highest number of alleles and expected heterozygosity. There is a clear distinction between the Piedmont breeds and the other populations. The Piedmont poultry differ from both commercial lines and other Italian breeds and retain a high level of genetic variability. As for other indigenous breeds, Bionda and Bianca could make an original contribution to the industry in the future. A collective planned approach to restoration is essential, because the flocks are managed with poor regulation. Enhancing connection between breeders with an efficient replacement interchange and mating plan is the right way of controlling inbreeding, preventing substructuring and increasing variability within the flocks.