A single origin and moderate bottleneck during domestication of soybean (Glycine max): implications from microsatellites and nucleotide sequences

Background and Aims

It is essential to illuminate the evolutionary history of crop domestication in order to understand further the origin and development of modern cultivation and agronomy; however, despite being one of the most important crops, the domestication origin and bottleneck of soybean (Glycine max) are poorly understood. In the present study, microsatellites and nucleotide sequences were employed to elucidate the domestication genetics of soybean.

Methods

The genomes of 79 landrace soybeans (endemic cultivated soybeans) and 231 wild soybeans (G. soja) that represented the species-wide distribution of wild soybean in East Asia were scanned with 56 microsatellites to identify the genetic structure and domestication origin of soybean. To understand better the domestication bottleneck, four nucleotide sequences were selected to simulate the domestication bottleneck.

Key Results

Model-based analysis revealed that most of the landrace genotypes were assigned to the inferred wild soybean cluster of south China, South Korea and Japan. Phylogeny for wild and landrace soybeans showed that all landrace soybeans formed a single cluster supporting a monophyletic origin of all the cultivars. The populations of the nearest branches which were basal to the cultivar lineage were wild soybeans from south China. The coalescent simulation detected a bottleneck severity of K′ = 2 during soybean domestication, which could be explained by a foundation population of 6000 individuals if domestication duration lasted 3000 years.

Conclusions

As a result of integrating geographic distribution with microsatellite genotype assignment and phylogeny between landrace and wild soybeans, a single origin of soybean in south China is proposed. The coalescent simulation revealed a moderate genetic bottleneck with an effective wild soybean population used for domestication estimated to be ≈2 % of the total number of ancestral wild soybeans. Wild soybeans in Asia, especially in south China contain tremendous genetic resources for cultivar improvement.     

Genetic variants in root architecture-related genes in a Glycine soja accession,a potential resource to improve cultivated soybean

Background

Root system architecture is important for water acquisition and nutrient acquisition for all crops. In soybean breeding programs, wild soybean alleles have been used successfully to enhance yield and seed composition traits, but have never been investigated to improve root system architecture. Therefore, in this study, high-density single-feature polymorphic markers and simple sequence repeats were used to map quantitative trait loci (QTLs) governing root system architecture in an inter-specific soybean mapping population developed from a cross between Glycine max and Glycine soja.

Results

Wild and cultivated soybean both contributed alleles towards significant additive large effect QTLs on chromosome 6 and 7 for a longer total root length and root distribution, respectively. Epistatic effect QTLs were also identified for taproot length, average diameter, and root distribution. These root traits will influence the water and nutrient uptake in soybean. Two cell division-related genes (D type cyclin and auxin efflux carrier protein) with insertion/deletion variations might contribute to the shorter root phenotypes observed in G. soja compared with cultivated soybean. Based on the location of the QTLs and sequence information from a second G. soja accession, three genes (slow anion channel associated 1 like, Auxin responsive NEDD8-activating complex and peroxidase), each with a non-synonymous single nucleotide polymorphism mutation were identified, which may also contribute to changes in root architecture in the cultivated soybean. In addition, Apoptosis inhibitor 5-like on chromosome 7 and slow anion channel associated 1-like on chromosome 15 had epistatic interactions for taproot length QTLs in soybean.

Conclusion

Rare alleles from a G. soja accession are expected to enhance our understanding of the genetic components involved in root architecture traits, and could be combined to improve root system and drought adaptation in soybean.

Electronic supplementary material

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

The Population Genetics of Cultivation: Domestication of a Traditional Chinese Medicine,Scrophularia ningpoensis Hemsl. (Scrophulariaceae)

Background

Domestic cultivation of medicinal plants is an important strategy for protecting these species from over harvesting. Some species of medicinal plants have been brought into cultivation for more than hundreds years. Concerns about severe loss of genetic diversity and sustainable cultivation can potentially limit future use of these valuable plants. Genetic studies with comprehensive sampling of multiple medicinal species by molecular markers will allow for assessment and management of these species. Here we examine the population genetic consequences of cultivation and domestication in Scrophularia ningpoensis Hemsl. We used chloroplast DNA and genomic AFLP markers to clarify not only the effects of domestication on genetic diversity, but also determine the geographic origins of cultivars and their genetic divergence from native populations. These results will allow both better management of cultivated populations, but also provide insights for crop improvement.

Results

Twenty-one cpDNA haplotypes of S. ningpoensis were identified. Wild populations contain all haplotypes, whereas only three haplotypes were found in cultivated populations with wild populations having twice the haplotype diversity of cultivated populations. Genetic differentiation between cultivated populations and wild populations was significant. Genomic AFLP markers revealed similar genetic diversity patterns. Furthermore, Structure analysis grouped all wild populations into two gene pools; two of which shared the same gene pool with cultivated S. ningpoensis. The result of Neighbor-Joining analysis was consistent with the structure analysis. In principal coordinate analysis, three cultivated populations from Zhejiang Province grouped together and were separated from other cultivated populations.

Conclusions

These results suggest that cultivated S. ningpoensis has experienced dramatic loss of genetic diversity under anthropogenic influence. We postulate that strong artificial selection for medicinal quality has resulted in genetic differentiation between cultivated and wild populations. Furthermore, it appears that wild populations in Jiangxi-Hunan area were involved in the origin of cultivated S. ningpoensis.     

Multi-locus analysis reveals a different pattern of genetic diversity for mitochondrial and nuclear DNA between wild and domestic pigs in East Asia

Background

A major reduction of genetic diversity in mtDNA occurred during the domestication of East Asian pigs. However, the extent to which genetic diversity has been lost in the nuclear genome is uncertain. To reveal levels and patterns of nucleotide diversity and to elucidate the genetic relationships and demographic history of domestic pigs and their ancestors, wild boars, we investigated 14 nuclear markers (including 8 functional genes, 2 pseudogenes and 4 intergenic regions) from 11 different chromosomes in East Asia-wide samples and pooled them with previously obtained mtDNA data for a combined analysis.

Principal Findings

The results indicated that domestic pigs and wild boars possess comparable levels of nucleotide diversity across the nuclear genome, which is inconsistent with patterns that have been found in mitochondrial genome.

Conclusions

This incongruence between the mtDNA and nuclear genomes is suggestive of a large-scale backcross between male wild boars and female domestic pigs in East Asia. Our data reveal the impacts of founder effects and backcross on the pig genome and help us better understand the complex demographic histories of East Asian pigs, which will be useful for future work on artificial selection.     

A comparison of rpoB and 16S rRNA as markers in pyrosequencing studies of bacterial diversity

Background

The 16S rRNA gene is the gold standard in molecular surveys of bacterial and archaeal diversity, but it has the disadvantages that it is often multiple-copy, has little resolution below the species level and cannot be readily interpreted in an evolutionary framework. We compared the 16S rRNA marker with the single-copy, protein-coding rpoB marker by amplifying and sequencing both from a single soil sample. Because the higher genetic resolution of the rpoB gene prohibits its use as a universal marker, we employed consensus-degenerate primers targeting the Proteobacteria.

Methodology/Principal Findings

Pyrosequencing can be problematic because of the poor resolution of homopolymer runs. As these erroneous runs disrupt the reading frame of protein-coding sequences, removal of sequences containing nonsense mutations was found to be a valuable filter in addition to flowgram-based denoising. Although both markers gave similar estimates of total diversity, the rpoB marker revealed more species, requiring an order of magnitude fewer reads to obtain 90% of the true diversity. The application of population genetic methods was demonstrated on a particularly abundant sequence cluster.

Conclusions/Significance

The rpoB marker can be a complement to the 16S rRNA marker for high throughput microbial diversity studies focusing on specific taxonomic groups. Additional error filtering is possible and tests for recombination or selection can be employed.     

Molecular characterization and genetic diversity analysis of soybean (Glycine max (L.) Merr.) germplasm accessions in India

Molecular characterization and genetic diversity among 82 soybean accessions was carried out by using 44 simple sequence repeat (SSR) markers. Of the 44 SSR markers used, 40 markers were found polymorphic among 82 soybean accessions. These 40 polymorphic markers produced a total of 119 alleles, of which five were unique alleles and four alleles were rare. The allele number for each SSR locus varied between two to four with an average of 2.97 alleles per marker. Polymorphic information content values of SSRs ranged from 0.101 to 0.742 with an average of 0.477. Jaccard’s similarity coefficient was employed to study the molecular diversity of 82 soybean accessions. The pairwise genetic similarity among 82 soybean accessions varied from 0.28 to 0.90. The dendrogram constructed based on genetic similarities among 82 soybean accessions identified three major clusters. The majority of genotypes including four improved cultivars were grouped in a single subcluster IIIa of cluster III, indicating high genetic resemblance among soybean germplasm collection in India.

Electronic supplementary material

The online version of this article (doi:10.1007/s12298-014-0266-y) contains supplementary material, which is available to authorized users.     

Amplified fragment length polymorphism (AFLP) in soybean: species diversity, inheritance, and near-isogenic line analysis

Amplified fragment length polymorphism (AFLP) analysis is a PCR-based technique capable of detecting more than 50 independent loci in a single PCR reaction. The objectives of the present study were to: (1) assess the extent of AFLP variation in cultivated (Gycine max L. Merr.) and wild soybean (G. soja Siebold & Zucc.), (2) determine genetic relationships among soybean accessions using AFLP data, and (3) evaluate the usefulness of AFLPs as genetic markers. Fifteen AFLP primer pairs detected a total of 759 AFLP fragments in a sample of 23 accessions of wild and cultivated soybean, with an average of 51 fragments produced per primer pair per accession. Two-hundred and seventy four fragments (36% of the total observed) were polymorphic, among which 127 (17%) were polymorphic in G. max and 237 (31%) were polymorphic in G. soja. F₂ segregation analysis of six AFLP fragments indicated that they segregate as stable Mendelian loci. The number of polymorphic loci detected per AFLP primer pair in a sample of 23 accessions ranged from 9 to 27. The AFLP phenotypic diversity values were greater in wild than in cultivated soybean. Cluster and principal component analyses using AFLP data clearly separated G. max and G. soja accessions. Within the G. max group, adapted soybean cultivars were tightly clustered, illustrating the relatively low genetic diversity present in cultivated soybean. AFLP analysis of four soybean near-isogenic lines (NILs) identified three AFLP markers putatively linked to a virus resistance gene from two sources. The capacity of AFLP analysis to detect thousands of independent genetic loci with minimal cost and time requirements makes them an ideal marker for a wide array of genetic investigations.     

Rapid loss of genetic variation in a founding population of Primula elatior (Primulaceae) after colonization

Background and Aims

Land-use changes and associated extinction/colonization dynamics can have a large impact on population genetic diversity of plant species. The aim of this study was to investigate genetic diversity in a founding population of the self-incompatible forest herb Primula elatior and to elucidate the processes that affect genetic diversity shortly after colonization.

Methods

AFLP markers were used to analyse genetic diversity across three age classes and spatial genetic structure within a founding population of P. elatior in a recently established stand in central Belgium. Parentage analyses were used to assess the amount of gene flow from outside the population and to investigate the contribution of mother plants to future generations.

Results

The genetic diversity of second and third generation plants was significantly reduced compared with that of first generation plants. Significant spatial genetic structure was observed. Parentage analyses showed that <20 % of the youngest individuals originated from parents outside the study population and that >50 % of first and second generation plants did not contribute to seedling recruitment.

Conclusions

These results suggest that a small effective population size and genetic drift can lead to rapid decline of genetic diversity of offspring in founding populations shortly after colonization. This multigenerational study also highlights that considerable amounts of gene flow seem to be required to counterbalance genetic drift and to sustain high levels of genetic diversity after colonization in recently established stands.Key words: AFLP, colonization, forest regeneration, genetic diversity, genetic drift, parentage analysis, spatial genetic structure     

Estimating genetic diversity across the neutral genome with the use of dense marker maps

Background

With the advent of high throughput DNA typing, dense marker maps have become available to investigate genetic diversity on specific regions of the genome. The aim of this paper was to compare two marker based estimates of the genetic diversity in specific genomic regions lying in between markers: IBD-based genetic diversity and heterozygosity.

Methods

A computer simulated population was set up with individuals containing a single 1-Morgan chromosome and 1665 SNP markers and from this one, an additional population was produced with a lower marker density i.e. 166 SNP markers. For each marker interval based on adjacent markers, the genetic diversity was estimated either by IBD probabilities or heterozygosity. Estimates were compared to each other and to the true genetic diversity. The latter was calculated for a marker in the middle of each marker interval that was not used to estimate genetic diversity.

Results

The simulated population had an average minor allele frequency of 0.28 and an LD (r²) of 0.26, comparable to those of real livestock populations. Genetic diversities estimated by IBD probabilities and by heterozygosity were positively correlated, and correlations with the true genetic diversity were quite similar for the simulated population with a high marker density, both for specific regions (r = 0.19-0.20) and large regions (r = 0.61-0.64) over the genome. For the population with a lower marker density, the correlation with the true genetic diversity turned out to be higher for the IBD-based genetic diversity.

Conclusions

Genetic diversities of ungenotyped regions of the genome (i.e. between markers) estimated by IBD-based methods and heterozygosity give similar results for the simulated population with a high marker density. However, for a population with a lower marker density, the IBD-based method gives a better prediction, since variation and recombination between markers are missed with heterozygosity.     

Diversity,expression and mRNA targeting abilities of Argonaute-targeting miRNAs among selected vascular plants

Background

Micro (mi)RNAs are important regulators of plant development. Across plant lineages, Dicer-like 1 (DCL1) proteins process long ds-like structures to produce micro (mi) RNA duplexes in a stepwise manner. These miRNAs are incorporated into Argonaute (AGO) proteins and influence expression of RNAs that have sequence complementarity with miRNAs. Expression levels of AGOs are greatly regulated by plants in order to minimize unwarranted perturbations using miRNAs to target mRNAs coding for AGOs. AGOs may also have high promoter specificity-sometimes expression of AGO can be limited to just a few cells in a plant. Viral pathogens utilize various means to counter antiviral roles of AGOs including hijacking the host encoded miRNAs to target AGOs. Two host encoded miRNAs namely miR168 and miR403 that target AGOs have been described in the model plant Arabidopsis and such a mechanism is thought to be well conserved across plants because AGO sequences are well conserved.

Results

We show that the interaction between AGO mRNAs and miRNAs is species-specific due to the diversity in sequences of two miRNAs that target AGOs, sequence diversity among corresponding target regions in AGO mRNAs and variable expression levels of these miRNAs among vascular plants. We used miRNA sequences from 68 plant species representing 31 plant families for this analysis. Sequences of miR168 and miR403 are not conserved among plant lineages, but surprisingly they differ drastically in their sequence diversity and expression levels even among closely related plants. Variation in miR168 expression among plants correlates well with secondary structures/length of loop sequences of their precursors.

Conclusions

Our data indicates a complex AGO targeting interaction among plant lineages due to miRNA sequence diversity and sequences of miRNA targeting regions among AGO mRNAs, thus leading to the assumption that the perturbations by viruses that use host miRNAs to target antiviral AGOs can only be species-specific. We also show that rapid evolution and likely loss of expression of miR168 isoforms in tobacco is related to the insertion of MITE-like transposons between miRNA and miRNA* sequences, a possible mechanism showing how miRNAs are lost in few plant lineages even though other close relatives have abundantly expressing miRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1049) contains supplementary material, which is available to authorized users.     

Functional significance of allelic variation at methuselah, an aging gene in Drosophila

Background

Longevity and age-specific patterns of mortality are complex traits that vary within and among taxa. Multiple candidate genes for aging have been identified in model systems by extended longevity mutant phenotypes, including the G-protein coupled receptor methuselah (mth) in D. melanogaster. These genes offer important insights into the mechanisms of lifespan determination and have been major targets of interest in the biology of aging. However, it is largely unknown whether these genes contribute to genetic variance for lifespan in natural populations, and consequently contribute to lifespan evolution.

Methodology/Principle Findings

For a gene to contribute to genetic variance for a particular trait, it must meet two criteria: natural allelic variation and functional differences among variants. Previous work showed that mth varies significantly among wild populations; here we assess the functional significance of wild-derived mth alleles on lifespan, fecundity and stress resistance using a quantitative complementation scheme. Our results demonstrate that mth alleles segregating in nature have a functional effect on all three traits.

Conclusions/Significance

These results suggest that allelic variation at mth contributes to observed differences in lifespan and correlated phenotypes in natural populations, and that evaluation of genetic diversity at candidate genes for aging can be a fruitful approach to identifying loci contributing to lifespan evolution.     

A Validated Methodology for Genetic Identification of Tuna Species (Genus Thunnus)

Background

Tuna species of the genus Thunnus, such as the bluefin tunas, are some of the most important and yet most endangered trade fish in the world. Identification of these species in traded forms, however, may be difficult depending on the presentation of the products, which may hamper conservation efforts on trade control. In this paper, we validated a genetic methodology that can fully distinguish between the eight Thunnus species from any kind of processed tissue.

Methodology

After testing several genetic markers, a complete discrimination of the eight tuna species was achieved using Forensically Informative Nucleotide Sequencing based primarily on the sequence variability of the hypervariable genetic marker mitochondrial DNA control region (mtDNA CR), followed, in some specific cases, by a second validation by a nuclear marker rDNA first internal transcribed spacer (ITS1). This methodology was able to distinguish all tuna species, including those belonging to the subgenus Neothunnus that are very closely related, and in consequence can not be differentiated with other genetic markers of lower variability. This methodology also took into consideration the presence of introgression that has been reported in past studies between T. thynnus, T. orientalis and T. alalunga. Finally, we applied the methodology to cross-check the species identity of 26 processed tuna samples.

Conclusions

Using the combination of two genetic markers, one mitochondrial and another nuclear, allows a full discrimination between all eight tuna species. Unexpectedly, the genetic marker traditionally used for DNA barcoding, cytochrome oxidase 1, could not differentiate all species, thus its use as a genetic marker for tuna species identification is questioned.     

Genetic Differentiation of the Mitochondrial Cytochrome Oxidase c Subunit I Gene in Genus Paramecium (Protista,Ciliophora)

Background

The mitochondrial cytochrome c oxidase subunit I (COI) gene is being used increasingly for evaluating inter- and intra-specific genetic diversity of ciliated protists. However, very few studies focus on assessing genetic divergence of the COI gene within individuals and how its presence might affect species identification and population structure analyses.

Methodology/Principal findings

We evaluated the genetic variation of the COI gene in five Paramecium species for a total of 147 clones derived from 21 individuals and 7 populations. We identified a total of 90 haplotypes with several individuals carrying more than one haplotype. Parsimony network and phylogenetic tree analyses revealed that intra-individual diversity had no effect in species identification and only a minor effect on population structure.

Conclusions

Our results suggest that the COI gene is a suitable marker for resolving inter- and intra-specific relationships of Paramecium spp.     

Development of DArT markers and assessment of diversity in Fusarium oxysporum f. sp. ciceris,wilt pathogen of chickpea (Cicer arietinum L.)

Background

Fusarium oxysporum f. sp. ciceris (Foc), the causal agent of Fusarium wilt of chickpea is highly variable and frequent recurrence of virulent forms have affected chickpea production and exhausted valuable genetic resources. The severity and yield losses of Fusarium wilt differ from place to place owing to existence of physiological races among isolates. Diversity study of fungal population associated with a disease plays a major role in understanding and devising better disease control strategies. The advantages of using molecular markers to understand the distribution of genetic diversity in Foc populations is well understood. The recent development of Diversity Arrays Technology (DArT) offers new possibilities to study the diversity in pathogen population. In this study, we developed DArT markers for Foc population, analysed the genetic diversity existing within and among Foc isolates, compared the genotypic and phenotypic diversity and infer the race scenario of Foc in India.

Results

We report the successful development of DArT markers for Foc and their utility in genotyping of Foc collections representing five chickpea growing agro-ecological zones of India. The DArT arrays revealed a total 1,813 polymorphic markers with an average genotyping call rate of 91.16% and a scoring reproducibility of 100%. Cluster analysis, principal coordinate analysis and population structure indicated that the different isolates of Foc were partially classified based on geographical source. Diversity in Foc population was compared with the phenotypic variability and it was found that DArT markers were able to group the isolates consistent with its virulence group. A number of race-specific unique and rare alleles were also detected.

Conclusion

The present study generated significant information in terms of pathogenic and genetic diversity of Foc which could be used further for development and deployment of region-specific resistant cultivars of chickpea. The DArT markers were proved to be a powerful diagnostic tool to study the genotypic diversity in Foc. The high number of DArT markers allowed a greater resolution of genetic differences among isolates and enabled us to examine the extent of diversity in the Foc population present in India, as well as provided support to know the changing race scenario in Foc population.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-454) contains supplementary material, which is available to authorized users.     

Profiling microRNAs in Eucalyptus grandis reveals no mutual relationship between alterations in miR156 and miR172 expression and adventitious root induction during development

Background

The change from juvenile to mature phase in woody plants is often accompanied by a gradual loss of rooting ability, as well as by reduced microRNA (miR) 156 and increased miR172 expression.

Results

We characterized the population of miRNAs of Eucalyptus grandis and compared the gradual reduction in miR156 and increase in miR172 expression during development to the loss of rooting ability. Forty known and eight novel miRNAs were discovered and their predicted targets are listed. The expression pattern of nine miRNAs was determined during adventitious root formation in juvenile and mature cuttings. While the expression levels of miR156 and miR172 were inverse in juvenile and mature tissues, no mutual relationship was found between high miR156 expression and rooting ability, or high miR172 expression and loss of rooting ability. This is shown both in E. grandis and in E. brachyphylla, in which explants that underwent rejuvenation in tissue culture conditions were also examined.

Conclusions

It is suggested that in these Eucalyptus species, there is no correlation between the switch of miR156 with miR172 expression in the stems and the loss of rooting ability.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-524) contains supplementary material, which is available to authorized users.     

The genetics of domestication of rice bean, Vigna umbellata

Background and Aims

The Asian genus Vigna, to which four cultivated species (rice bean, azuki bean, mung bean and black gram) belong, is suitable for comparative genomics. The aims were to construct a genetic linkage map of rice bean, to identify the genomic regions associated with domestication in rice bean, and to compare these regions with those in azuki bean.

Methods

A genetic linkage map was constructed by using simple sequence repeat and amplified fragment length polymorphism markers in the BC₁F₁ population derived from a cross between cultivated and wild rice bean. Using this map, 31 domestication-related traits were dissected into quantitative trait loci (QTLs). The genetic linkage map and QTLs of rice bean were compared with those of azuki bean.

Key Results

A total of 326 markers converged into 11 linkage groups (LGs), corresponding to the haploid number of rice bean chromosomes. The domestication-related traits in rice bean associated with a few major QTLs distributed as clusters on LGs 2, 4 and 7. A high level of co-linearity in marker order between the rice bean and azuki bean linkage maps was observed. Major QTLs in rice bean were found on LG4, whereas major QTLs in azuki bean were found on LG9.

Conclusions

This is the first report of a genetic linkage map and QTLs for domestication-related traits in rice bean. The inheritance of domestication-related traits was so simple that a few major QTLs explained the phenotypic variation between cultivated and wild rice bean. The high level of genomic synteny between rice bean and azuki bean facilitates QTL comparison between species. These results provide a genetic foundation for improvement of rice bean; interchange of major QTLs between rice bean and azuki bean might be useful for broadening the genetic variation of both species.