Genome-wide expression analysis of vegetative organs during developmental and herbicide-induced whole plant senescence in Arabidopsis thaliana.

BACKGROUND: Whole plant senescence represents the final stage in the life cycle of annual plants, characterized by the decomposition of aging organs and transfer of nutrients to seeds, thereby ensuring the survival of next generation. However, the transcriptomic profile of vegetative organs during this death process remains to be fully elucidated, especially regarding the distinctions between natural programmed death and artificial sudden death induced by herbicide. RESULTS: Differential genes expression analysis using RNA-seq in leaves and roots of Arabidopsis thaliana revealed that natural senescence commenced in leaves at 45-52 days after planting, followed by roots initiated at 52-60 days. Additionally, both organs exhibited similarities with artificially induced senescence by glyphosate. Transcription factors Rap2.6L and WKRY75 appeared to serve as central mediators of regulatory changes during natural senescence, as indicated by co-expression networks. Furthermore, the upregulation of RRTF1, exclusively observed during natural death, suggested its role as a regulator of jasmonic acid and reactive oxygen species (ROS) responses, potentially triggering nitrogen recycling in leaves, such as the glutamate dehydrogenase (GDH) shunt. Root senescence was characterized by the activation of AMT2;1 and GLN1;3, facilitating ammonium availability for root-to-shoot translocation, likely under the regulation of PDF2.1. CONCLUSIONS: Our study offers valuable insights into the transcriptomic interplay between phytohormones and ROS during whole plant senescence. We observed distinct regulatory networks governing nitrogen utilization in leaf and root senescence processes. Furthermore, the efficient allocation of energy from vegetative organs to seeds emerges as a critical determinant of population sustainability of annual Arabidopsis.

Microbial Communities Along 2,3,7,8-tetrachlorodibenzodioxin Concentration Gradient in Soils Polluted with Agent Orange Based on Metagenomic Analyses.

The 2,3,7,8-tetrachlorodibenzodioxin (TCDD), a contaminant in Agent Orange released during the US-Vietnam War, led to a severe environmental crisis. Approximately, 50 years have passed since the end of this war, and vegetation has gradually recovered from the pollution. Soil bacterial communities were investigated by 16S metagenomics in habitats with different vegetation physiognomies in Central Vietnam, namely, forests (S0), barren land (S1), grassland (S2), and developing woods (S3). Vegetation complexity was negatively associated with TCDD concentrations, revealing the reasoning behind the utilization of vegetation physiognomy as an indicator for ecological succession along the gradient of pollutants. Stark changes in bacterial composition were detected between S0 and S1, with an increase in Firmicutes and a decrease in Acidobacteria and Bacteroidetes. Notably, dioxin digesters Arthrobacter, Rhodococcus, Comamonadaceae, and Bacialles were detected in highly contaminated soil (S1). Along the TCDD gradients, following the dioxin decay from S1 to S2, the abundance of Firmicutes and Actinobacteria decreased, while that of Acidobacteria increased; slight changes occurred at the phylum level from S2 to S3. Although metagenomics analyses disclosed a trend toward bacterial communities before contamination with vegetation recovery, non-metric multidimensional scaling analysis unveiled a new trajectory deviating from the native state. Recovery of the bacterial community may have been hindered, as indicated by lower bacterial diversity in S3 compared to S0 due to a significant loss of bacterial taxa and recruitment of fewer colonizers. The results indicate that dioxins significantly altered the soil microbiomes into a state of disorder with a deviating trajectory in restoration.

Behavioral and brain- transcriptomic synchronization between the two opponents of a fighting pair of the fish Betta splendens.

Conspecific male animals fight for resources such as food and mating opportunities but typically stop fighting after assessing their relative fighting abilities to avoid serious injuries. Physiologically, how the fighting behavior is controlled remains unknown. Using the fighting fish Betta splendens, we studied behavioral and brain-transcriptomic changes during the fight between the two opponents. At the behavioral level, surface-breathing, and biting/striking occurred only during intervals between mouth-locking. Eventually, the behaviors of the two opponents became synchronized, with each pair showing a unique behavioral pattern. At the physiological level, we examined the expression patterns of 23,306 brain transcripts using RNA-sequencing data from brains of fighting pairs after a 20-min (D20) and a 60-min (D60) fight. The two opponents in each D60 fighting pair showed a strong gene expression correlation, whereas those in D20 fighting pairs showed a weak correlation. Moreover, each fighting pair in the D60 group showed pair-specific gene expression patterns in a grade of membership analysis (GoM) and were grouped as a pair in the heatmap clustering. The observed pair-specific individualization in brain-transcriptomic synchronization (PIBS) suggested that this synchronization provides a physiological basis for the behavioral synchronization. An analysis using the synchronized genes in fighting pairs of the D60 group found genes enriched for ion transport, synaptic function, and learning and memory. Brain-transcriptomic synchronization could be a general phenomenon and may provide a new cornerstone with which to investigate coordinating and sustaining social interactions between two interacting partners of vertebrates.

Endophytic Microbiome of Biofuel Plant Miscanthus sinensis (Poaceae) Interacts with Environmental Gradients.

Miscanthus in Taiwan occupies a cline along altitude and adapts to diverse environments, e.g., habitats of high salinity and volcanoes. Rhizospheric and endophytic bacteria may help Miscanthus acclimate to those stresses. The relative contributions of rhizosphere vs. endosphere compartments to the adaptation remain unknown. Here, we used targeted metagenomics to compare the microbial communities in the rhizosphere and endosphere among ecotypes of M. sinensis that dwell habitats under different stresses. Proteobacteria and Actinobacteria predominated in the endosphere. Diverse phyla constituted the rhizosphere microbiome, including a core microbiome found consistently across habitats. In endosphere, the predominance of the bacteria colonizing from the surrounding soil suggests that soil recruitment must have subsequently determined the endophytic microbiome in Miscanthus roots. In endosphere, the bacterial diversity decreased with the altitude, likely corresponding to rising limitation to microorganisms according to the species-energy theory. Specific endophytes were associated with different environmental stresses, e.g., Pseudomonas spp. for alpine and Agrobacterium spp. for coastal habitats. This suggests Miscanthus actively recruits an endosphere microbiome from the rhizosphere it influences.

Revegetation on abandoned salt ponds relieves the seasonal fluctuation of soil microbiomes.

BACKGROUND: Salt pond restoration aims to recover the environmental damages that accumulated over the long history of salt production. Of the restoration strategies, phytoremediation that utilizes salt-tolerant plants and soil microorganisms to reduce the salt concentrations is believed to be environmentally-friendly. However, little is known about the change of bacterial community during salt pond restoration in the context of phytoremediation. In the present study, we used 16S metagenomics to compare seasonal changes of bacterial communities between the revegetated and barren salterns at Sicao, Taiwan. RESULTS: In both saltern types, Proteobacteria, Planctomycetes, Chloroflexi, and Bacteroidetes were predominant at the phylum level. In the revegetated salterns, the soil microbiomes displayed high species diversities and underwent a stepwise transition across seasons. In the barren salterns, the soil microbiomes fluctuated greatly, indicating that mangroves tended to stabilize the soil microorganism communities over the succession. Bacteria in the order Halanaerobiaceae and archaea in the family Halobacteriaceae that were adapted to high salinity exclusively occurred in the barren salterns. Among the 441 persistent operational taxonomic units detected in the revegetated salterns, 387 (87.5%) were present as transient species in the barren salterns. Only 32 persistent bacteria were exclusively detected in the revegetated salterns. Possibly, salt-tolerant plants provided shelters for those new colonizers. CONCLUSIONS: The collective data indicate that revegetation tended to stabilize the microbiome across seasons and enriched the microbial diversity in the salterns, especially species of Planctomycetes and Acidobacteria.

Common Stress Transcriptome Analysis Reveals Functional and Genomic Architecture Differences Between Early and Delayed Response Genes.

To identify the similarities among responses to diverse environmental stresses, we analyzed the transcriptome response of rice roots to three rhizotoxic perturbations (chromium, ferulic acid and mercury) and identified common early-transient, early-constant and delayed gene inductions. Common early response genes were mostly associated with signal transduction and hormones, and delayed response genes with lipid metabolism. Network component analysis revealed complicated interactions among common genes, the most highly connected signaling hubs being PP2C68, MPK5, LRR-RLK and NPR1. Gene architecture studies revealed different conserved promoter motifs and a different ratio of CpG island distribution between early and delayed genes. In addition, early-transient genes had more exons and a shorter first exon. IMEter was used to calculate the transcription regulation effects of introns, with greater effects for the first introns of early-transient than delayed genes. The higher Ka/Ks (non-synonymous/synonymous mutation) ratio of early-constant genes than early-transient, delayed and the genome median demonstrates the rapid evolution of early-constant genes. Our results suggest that finely tuned transcriptional control in response to environmental stress in rice depends on genomic architecture and signal intensity and duration.

Comparative transcriptome analysis of Gastrodia elata (Orchidaceae) in response to fungus symbiosis to identify gastrodin biosynthesis-related genes.

BACKGROUND: Gastrodia elata Blume (Orchidaceae) is an important Chinese medicine with several functional components. In the life cycle of G. elata, the orchid develops a symbiotic relationship with two compatible mycorrhizal fungi Mycena spp. and Armillaria mellea during seed germination to form vegetative propagation corm and vegetative growth to develop tubers, respectively. Gastrodin (p-hydroxymethylphenol-beta-D-glucoside) is the most important functional component in G. elata, and gastrodin significantly increases from vegetative propagation corms to tubers. To address the gene regulation mechanism in gastrodin biosynthesis in G. elata, a comparative analysis of de novo transcriptome sequencing among the vegetative propagation corms and tubers of G. elata and A. mellea was conducted using deep sequencing. RESULTS: Transcriptome comparison between the vegetative propagation corms and juvenile tubers of G. elata revealed 703 differentially expressed unigenes, of which 298 and 405 unigenes were, respectively up-regulated (fold-change ≥ 2, q-value < 0.05, the trimmed mean of M-values (TMM)-normalized fragments per kilobase of transcript per Million mapped reads (FPKM) > 10) and down-regulated (fold-change ≤ 0.5, q-value <0.05, TMM-normalized FPKM > 10) in juvenile tubers. After Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 112 up-regulated unigenes with KEGG Ortholog identifiers (KOids) or enzyme commission (EC) numbers were assigned to 159 isogroups involved in seventy-eight different pathways, and 132 down-regulated unigenes with KOids or EC numbers were assigned to 168 isogroups, involved in eighty different pathways. The analysis of the isogroup genes from all pathways revealed that the two unigenes TRINITY_DN54282_c0_g1 (putative monooxygenases) and TRINITY_DN50323_c0_g1 (putative glycosyltransferases) might participate in hydroxylation and glucosylation in the gastrodin biosynthetic pathway. CONCLUSIONS: The gene expression of the two unique unigenes encoding monooxygenase and glycosyltransferase significantly increases from vegetative propagation corms to tubers, and the molecular basis of gastrodin biosynthesis in the tubers of G. elata is proposed.

Cross-Species, Amplifiable EST-SSR Markers for Amentotaxus Species Obtained by Next-Generation Sequencing.

Amentotaxus, a genus of Taxaceae, is an ancient lineage with six relic and endangered species. Four Amentotaxus species, namely A. argotaenia, A. formosana, A. yunnanensis, and A. poilanei, are considered a species complex because of their morphological similarities. Small populations of these species are allopatrically distributed in Asian forests. However, only a few codominant markers have been developed and applied to study population genetic structure of these endangered species. In this study, we developed and characterized polymorphic expressed sequence tag-simple sequence repeats (EST-SSRs) from the transcriptome of A. formosana. We identified 4955 putative EST-SSRs from 68,281 unigenes as potential molecular markers. Twenty-six EST-SSRs were selected for estimating polymorphism and transferability among Amentotaxus species, of which 23 EST-SSRs were polymorphic within Amentotaxus species. Among these, the number of alleles ranged from 1-4, the polymorphism information content ranged from 0.000-0.692, and the observed and expected heterozygosity were 0.000-1.000 and 0.080-0.740, respectively. Population genetic structure analyses confirmed that A. argotaenia and A. formosana were separate species and A. yunnanensis and A. poilanei were the same species. These novel EST-SSRs can facilitate further population genetic structure research of Amentotaxus species.

Adaptive divergence with gene flow in incipient speciation of Miscanthus floridulus/sinensis complex (Poaceae).

Young incipient species provide ideal materials for untangling the process of ecological speciation in the presence of gene flow. The Miscanthus floridulus/sinensis complex exhibits diverse phenotypic and ecological differences despite recent divergence (approximately 1.59 million years ago). To elucidate the process of genetic differentiation during early stages of ecological speciation, we analyzed genomic divergence in the Miscanthus complex using 72 randomly selected genes from a newly assembled transcriptome. In this study, rampant gene flow was detected between species, estimated as M = 3.36 × 10(-9) to 1.20 × 10(-6) , resulting in contradicting phylogenies across loci. Nevertheless, beast analyses revealed the species identity and the effects of extrinsic cohesive forces that counteracted the non-stop introgression. As expected, early in speciation with gene flow, only 3-13 loci were highly diverged; two to five outliers (approximately 2.78-6.94% of the genome) were characterized by strong linkage disequilibrium, and asymmetrically distributed among ecotypes, indicating footprints of diversifying selection. In conclusion, ecological speciation of incipient species of Miscanthus probably followed the parapatric model, whereas allopatric speciation cannot be completely ruled out, especially between the geographically isolated northern and southern M. sinensis, for which no significant gene flow across oceanic barriers was detected. Divergence between local ecotypes in early-stage speciation began at a few genomic regions under the influence of natural selection and divergence hitchhiking that overcame gene flow.

Highly diversified fungi are associated with the achlorophyllous orchid Gastrodia flavilabella.

BACKGROUND: Mycoheterotrophic orchids are achlorophyllous plants that obtain carbon and nutrients from their mycorrhizal fungi. They often show strong preferential association with certain fungi and may obtain nutrients from surrounding photosynthetic plants through ectomycorrhizal fungi. Gastrodia is a large genus of mycoheterotrophic orchids in Asia, but Gastrodia species' association with fungi has not been well studied. We asked two questions: (1) whether certain fungi were preferentially associated with G. flavilabella, which is an orchid in Taiwan and (2) whether fungal associations of G. flavilabella were affected by the composition of fungi in the environment. RESULTS: Using next-generation sequencing, we studied the fungal communities in the tubers of Gastrodia flavilabella and the surrounding soil. We found (1) highly diversified fungi in the G. flavilabella tubers, (2) that Mycena species were the predominant fungi in the tubers but minor in the surrounding soil, and (3) the fungal communities in the G. flavilabella tubers were clearly distinct from those in the surrounding soil. We also found that the fungal composition in soil can change quickly with distance. CONCLUSIONS: G. flavilabella was associated with many more fungi than previously thought. Among the fungi in the tuber of G. flavilabella, Mycena species were predominant, different from the previous finding that adult G. elata depends on Armillaria species for nutritional supply. Moreover, the preferential fungus association of G. flavilabella was not significantly influenced by the composition of fungi in the environment.

Ecological genomics in Xanthomonas: the nature of genetic adaptation with homologous recombination and host shifts.

BACKGROUND: Comparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants. RESULTS: Using whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting. CONCLUSION: Altogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification.

Biogeography of the Phalaenopsis amabilis species complex inferred from nuclear and plastid DNAs.

BACKGROUND: Phalaenopsis is one of the important commercial orchids in the world. Members of the P. amabilis species complex represent invaluable germplasm for the breeding program. However, the phylogeny of the P. amabilis species complex is still uncertain. The Phalaenopsis amabilis species complex (Orchidaceae) consists of subspecies amabilis, moluccana, and rosenstromii of P. amabilis, as well as P. aphrodite ssp. aphrodite, P. ap. ssp. formosana, and P. sanderiana. The aims of this study were to reconstruct the phylogeny and biogeographcial patterns of the species complex using Neighbor Joining (NJ), Maxinum Parsimony (MP), Bayesian Evolutionary Analysis Sampling Trees (BEAST) and Reconstruct Ancestral State in Phylogenies (RASP) analyses based on sequences of internal transcribed spacers 1 and 2 from the nuclear ribosomal DNA and the trnH-psbA spacer from the plastid DNA. RESULTS: A pattern of vicariance, dispersal, and vicariance + dispersal among disjunctly distributed taxa was uncovered based on RASP analysis. Although two subspecies of P. aphrodite could not be differentiated from each other in dispersal state, they were distinct from P. amabilis and P. sanderiana. Within P. amabilis, three subspecies were separated phylogenetically, in agreement with the vicariance or vicariance + dispersal scenario, with geographic subdivision along Huxley's, Wallace's and Lydekker's Lines. Molecular dating revealed such subdivisions among taxa of P. amabilis complex dating back to the late Pleistocene. Population-dynamic analyses using a Bayesian skyline plot suggested that the species complex experienced an in situ range expansion and population concentration during the late Last Glacial Maximum (LGM). CONCLUSIONS: Taxa of the P. amabilis complex with disjunct distributions were differentiated due to vicariance or vicariance + dispersal, with events likely occurring in the late Pleistocene. Demographic growth associated with the climatic oscillations in the Würm glacial period followed the species splits. Nevertheless, a subsequent population slowdown occurred in the late LGM due to extinction of regional populations. The reduction of suitable habitats resulted in geographic fragmenttation of the remaining taxa.

Taming the wild: resolving the gene pools of non-model Arabidopsis lineages.

BACKGROUND: Wild relatives in the genus Arabidopsis are recognized as useful model systems to study traits and evolutionary processes in outcrossing species, which are often difficult or even impossible to investigate in the selfing and annual Arabidopsis thaliana. However, Arabidopsis as a genus is littered with sub-species and ecotypes which make realizing the potential of these non-model Arabidopsis lineages problematic. There are relatively few evolutionary studies which comprehensively characterize the gene pools across all of the Arabidopsis supra-groups and hypothesized evolutionary lineages and none include sampling at a world-wide scale. Here we explore the gene pools of these various taxa using various molecular markers and cytological analyses. RESULTS: Based on ITS, microsatellite, chloroplast and nuclear DNA content data we demonstrate the presence of three major evolutionary groups broadly characterized as A. lyrata group, A. halleri group and A. arenosa group. All are composed of further species and sub-species forming larger aggregates. Depending on the resolution of the marker, a few closely related taxa such as A. pedemontana, A. cebennensis and A. croatica are also clearly distinct evolutionary lineages. ITS sequences and a population-based screen based on microsatellites were highly concordant. The major gene pools identified by ITS sequences were also significantly differentiated by their homoploid nuclear DNA content estimated by flow cytometry. The chloroplast genome provided less resolution than the nuclear data, and it remains unclear whether the extensive haplotype sharing apparent between taxa results from gene flow or incomplete lineage sorting in this relatively young group of species with Pleistocene origins. CONCLUSIONS: Our study provides a comprehensive overview of the genetic variation within and among the various taxa of the genus Arabidopsis. The resolved gene pools and evolutionary lineages will set the framework for future comparative studies on genetic diversity. Extensive population-based phylogeographic studies will also be required, however, in particular for A. arenosa and their affiliated taxa and cytotypes.

Footprints of natural and artificial selection for photoperiod pathway genes in Oryza.

Asian rice, Oryza sativa, consists of two major subspecies, indica and japonica, which are physiologically differentiated and adapted to different latitudes. Genes for photoperiod sensitivity are likely targets of selection along latitude. We examined the footprints of natural and artificial selections for four major genes of the photoperiod pathway, namely PHYTOCHROME B (PhyB), HEADING DATE 1 (Hd1), HEADING DATE 3a (Hd3a), and EARLY HEADING DATE 1 (Ehd1), by investigation of the patterns of nucleotide polymorphisms in cultivated and wild rice. Geographical subdivision between tropical and subtropical O. rufipogon was found for all of the photoperiod genes in plants divided by the Tropic of Cancer (TOC). All of these genes, except for PhyB, were characterized by the existence of clades that split a long time ago and that corresponded to latitudinal subdivisions, and revealed a likely diversifying selection. Ssp. indica showed close affinity to tropical O. rufipogon for all genes, while ssp. japonica, which has a much wider range of distribution, displayed complex patterns of differentiation from O. rufipogon, which reflected various agricultural needs in relation to crop yield. In japonica, all genes, except Hd3a, were genetically differentiated at the TOC, while geographical subdivision occurred at 31°N in Hd3a, probably the result of varying photoperiods. Many other features of the photoperiod genes revealed domestication signatures, which included high linkage disequilibrium (LD) within genes, the occurrence of frequent and recurrent non-functional Hd1 mutants in cultivated rice, crossovers between subtropical and tropical alleles of Hd1, and significant LD between Hd1 and Hd3a in japonica and indica.

Characterization of 42 microsatellite markers from poison ivy, Toxicodendron radicans (Anacardiaceae).

Poison ivy, Toxicodendron radicans, and poison oaks, T. diversilobum and T. pubescens, are perennial woody species of the Anacardiaceae and are poisonous, containing strong allergens named urushiols that cause allergic contact dermatitis. Poison ivy is a species distributed from North America to East Asia, while T. diversilobum and T. pubescens are distributed in western and eastern North America, respectively. Phylogreography and population structure of these species remain unclear. Here, we developed microsatellite markers, via constructing a magnetic enriched microsatellite library, from poison ivy. We designed 51 primer pairs, 42 of which successfully yielded products that were subsequently tested for polymorphism in poison oak, and three subspecies of poison ivy. Among the 42 loci, 38 are polymorphic, while 4 are monomorphic. The number of alleles and the expected heterozygosity ranged from 1 to 12 and from 0.10 to 0.87, respectively, in poison ivy, while varied from 2 to 8 and, from 0.26 to 0.83, respectively in poison oak. Genetic analysis revealed distinct differentiation between poison ivy and poison oak, whereas slight genetic differentiation was detected among three subspecies of poison ivy. These highly polymorphic microsatellite fingerprints enable biologists to explore the population genetics, phylogeography, and speciation in Toxicodendron.

New insights into polyploid evolution and dynamic nature of Ludwigia section Isnardia (Onagraceae).

BACKGROUND: While polyploids are common in plants, the evolutionary history and natural dynamics of most polyploid groups are still unclear. Owing to plentiful earlier systematic studies, Ludwigia sect. Isnardia (comprising 22 wetland taxa) is an ideal allopolyploid complex to investigate polyploid evolution and natural dynamics within and among taxa. With a considerable sampling, we concentrated on revisiting earlier phylogenies of Isnardia, reevaluating the earlier estimated age of the most recent common ancestor (TMRCA), exploring the correlation between infraspecific genetic diversity and ploidy levels, and inspecting interspecific gene flows among taxa. RESULTS: Phylogenetic trees and network concurred with earlier phylogenies and hypothesized genomes by incorporating 192 atpB-rbcL and ITS sequences representing 91% of Isnardia taxa. Moreover, we detected three multi-origin taxa. Our findings on L. repens and L. sphaerocarpa were consistent with earlier studies; L. arcuata was reported as a multi-origin taxon here, and an additional evolutionary scenario of L. sphaerocarpa was uncovered, both for the first time. Furthermore, estimated Isnardia TMRCA ages based on our data (5.9 or 8.9 million years ago) are in accordance with earlier estimates, although younger than fossil dates (Middle Miocene). Surprisingly, infraspecific genetic variations of Isnardia taxa did not increase with ploidy levels as anticipated from many other polyploid groups. In addition, the exuberant, low, and asymmetrical gene flows among Isnardia taxa indicated that the reproductive barriers may be weakened owing to allopolyploidization, which has rarely been reported. CONCLUSIONS: The present research gives new perceptions of the reticulate evolution and dynamic nature of Isnardia and points to gaps in current knowledge about allopolyploid evolution.

Deciphering microbial community dynamics along the fermentation course of soy sauce under different temperatures using metagenomic analysis.

Fermented soy sauce consists of microorganisms that exert beneficial effects. However, the microbial community dynamics during the fermentation course is poorly characterized. Soy sauce production is classified into the stages of mash fermentation with koji (S0), brine addition (S1), microbial transformation (S2), flavor creation (S3), and fermentation completion (S4). In this study, microbial succession was investigated across stages at different temperatures using metagenomics analyses. During mash fermentation, Aspergillus dominated the fungal microbiota in all stages, while the bacterial composition was dominated by Bacillus at room temperature and by a diverse composition of enriched lactic acid bacteria (LAB) at a controlled temperature. Compared with a stable fungal composition, bacterial dynamics were mostly attributable to fluctuations of LAB, which break down carbohydrates into lactic acid. After adding brine, increased levels of Enterococcus and decreased levels of Lactococcus from S1 to S4 may reflect differences in salinity tolerance. Staphylococcus, as a fermentation starter at S0, stayed predominant throughout fermentation and hydrolyzed soybean proteins. Meanwhile, Rhizopus and Penicillium may improve the flavor. The acidification of soy sauce was likely attributable to production of organic acids by Bacillus and LAB under room temperature and controlled temperature conditions, respectively. Metagenomic analysis revealed that microbial succession was associated with the fermentation efficiency and flavor enhancement. Controlled temperature nurture more LAB than uncontrolled temperatures and may ensure the production of lactic acid for the development of soy sauce flavor.

Isolation and characterization of polymorphic microsatellite loci from Metapenaeopsis barbata using PCR-Based Isolation of Microsatellite Arrays (PIMA).

The red-spot prawn, Metapenaeopsis barbata, is a commercially important, widely distributed demersal species in the Indo-West Pacific Ocean. Overfishing has made its populations decline in the past decade. To study conservation genetics, eight polymorphic microsatellite loci were isolated. Genetic characteristics of the SSR (simple sequence repeat) fingerprints were estimated in 61 individuals from adjacent seas of Taiwan and China. The number of alleles, ranging from 2 to 4, as well as observed and expected heterozygosities in populations, ranging from 0.048 to 0.538, and 0.048 and 0.654, respectively, were detected. No deviation from Hardy-Weinberg expectations was detected at either locus. No significant linkage disequilibrium was detected in locus pairs. The polymorphic microsatellite loci will be useful for investigations of the genetic variation, population structure, and conservation genetics of this species.

Genetic population structure of the alpine species Rhododendron pseudochrysanthum sensu lato (Ericaceae) inferred from chloroplast and nuclear DNA.

BACKGROUND: A complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence. We examined the phylogeography and the evolutionary history of the Rhododendron pseudochrysanthum s. l. RESULTS: Systematic inconsistency was detected between gene genealogies of the cpDNA and nrDNA. Rooted at R. hyperythrum and R. formosana, both trees lacked reciprocal monophyly for all members of the complex. For R. pseudochrysanthum s.l., the spatial distribution of the cpDNA had a noteworthy pattern showing high genetic differentiation (FST=0.56-0.72) between populations in the Yushan Mountain Range and populations of the other mountain ranges. CONCLUSION: Both incomplete lineage sorting and interspecific hybridization/introgression may have contributed to the lack of monophyly among R. hyperythrum, R. formosana and R. pseudochrysanthum s.l. Independent colonizations, plus low capabilities of seed dispersal in current environments, may have resulted in the genetic differentiation between populations of different mountain ranges. At the population level, the populations of Central, and Sheishan Mountains may have undergone postglacial demographic expansion, while populations of the Yushan Mountain Range are likely to have remained stable ever since the colonization. In contrast, the single population of the Alishan Mountain Range with a fixed cpDNA haplotype may have experienced bottleneck/founder's events.

Conservation genetics and phylogeography of endangered and endemic shrub Tetraena mongolica (Zygophyllaceae) in Inner Mongolia, China.

BACKGROUND: Tetraena mongolica (Zygophyllaceae), an endangered endemic species in western Inner Mongolia, China. For endemic species with a limited geographical range and declining populations, historical patterns of demography and hierarchical genetic structure are important for determining population structure, and also provide information for developing effective and sustainable management plans. In this study, we assess genetic variation, population structure, and phylogeography of T. mongolica from eight populations. Furthermore, we evaluate the conservation and management units to provide the information for conservation. RESULTS: Sequence variation and spatial apportionment of the atpB-rbcL noncoding spacer region of the chloroplast DNA were used to reconstruct the phylogeography of T. mongolica. A total of 880 bp was sequenced from eight extant populations throughout the whole range of its distribution. At the cpDNA locus, high levels of genetic differentiation among populations and low levels of genetic variation within populations were detected, indicating that most seed dispersal was restricted within populations. CONCLUSIONS: Demographic fluctuations, which led to random losses of genetic polymorphisms from populations, due to frequent flooding of the Yellow River and human disturbance were indicated by the analysis of BEAST skyline plot. Nested clade analysis revealed that restricted gene flow with isolation by distance plus occasional long distance dispersal is the main evolutionary factor affecting the phylogeography and population structure of T. mongolica. For setting a conservation management plan, each population of T. mongolica should be recognized as a conservation unit.