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.

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.

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.

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.

The operations of quantum logic gates with pure and mixed initial states.

The implementations of quantum logic gates realized by the rovibrational states of a C(12)O(16) molecule in the X((1)Σ(+)) electronic ground state are investigated. Optimal laser fields are obtained by using the modified multitarget optimal theory (MTOCT) which combines the maxima of the cost functional and the fidelity for state and quantum process. The projection operator technique together with modified MTOCT is used to get optimal laser fields. If initial states of the quantum gate are pure states, states at target time approach well to ideal target states. However, if the initial states are mixed states, the target states do not approach well to ideal ones. The process fidelity is introduced to investigate the reliability of the quantum gate operation driven by the optimal laser field. We found that the quantum gates operate reliably whether the initial states are pure or mixed.

Elastic properties of self-supported circular thin copper films calculated from equilibrium thermal vibration.

Thermomechanical vibration of ultrathin, self-supported copper films due to thermal fluctuations is studied via the molecular dynamics simulation at room temperature. The elastodynamic theory with pre-stress is adopted to extract the physical properties of the films by comparing with the molecular dynamics data. The edge-clamped circular films consist of several atomic layers of fcc copper with the [100] direction normal to the film surface. From the time-history trajectories of atoms and their Fourier frequency spectrums, it was found that the fundamental resonant frequency non-monotonically varies with the film thickness due to the existence of residual stress in the film. Multiple resonant modes are adopted for modulus calculation and residual stress determination. The value of Young's modulus increases with increasing thickness of the film and the residual stress decreases with increasing thickness. Thicker films exhibit less residual stress, indicating the equilibrium distance between copper atoms changes with the film thickness.

Population structure of Nouelia insignis (Asteraceae), an endangered species in southwestern China, based on chloroplast DNA sequences: recent demographic shrinking.

Nouelia insignis, an endangered species, is distributed in the Jinsha and Nanpan drainage areas in southwestern China. In this study, we examined the genetic diversity and population structure based on the sequences of the cpDNA rpL 16 intron. Low levels of genetic variation were detected within all populations of the endemic species. A gene genealogy of 11 haplotypes recovered two major lineages I and II, with haplotypes H1 and H6 nested as interior nodes, respectively. Haplotype H1 was widespread in all populations, while haplotype H6 was restricted to populations southern of the Jinsha River. Low levels of genetic differentiation were detected, as most F (st) values between populations were zero. This result, however, contradicts previous studies based on allozymes and fingerprinting. Genetic analyses suggested that coancestry due to low evolutionary rates resulted in the lack of geographical subdivision. Molecular dating estimated that the two lineages split about 3.224 MYA (95% CI 1.070-6.089 MYA). Maintenance of ancestral polymorphisms was possibly attributable to a long-standing large effective population size until recently. Postglacial demographic expansion was supported by a unimodal mismatch distribution and star-like phylogenies.

Multilocus analysis of genetic divergence between outcrossing Arabidopsis species: evidence of genome-wide admixture.

• Outcrossing Arabidopsis species that diverged from their inbreeding relative Arabidopsis thaliana 5 million yr ago and display a biogeographical pattern of interspecific sympatry vs intraspecific allopatry provides an ideal model for studying impacts of gene introgression and polyploidization on species diversification. • Flow cytometry analyses detected ploidy polymorphisms of 2× and 4× in Arabidopsis lyrata ssp. kamchatica of Taiwan. Genomic divergence between species/subspecies was estimated based on 98 randomly chosen nuclear genes. Multilocus analyses revealed a mosaic genome in diploid A. l. kamchatica composed of Arabidopsis halleri-like and A. lyrata-like alleles. • Coalescent analyses suggest that the segregation of ancestral polymorphisms alone cannot explain the high inconsistency between gene trees across loci, and that gene introgression via diploid A. l. kamchatica likely distorts the molecular phylogenies of Arabidopsis species. However, not all genes migrated across species freely. Gene ontology analyses suggested that some nonmigrating genes were constrained by natural selection. • High levels of estimated ancestral polymorphisms between A. halleri and A. lyrata suggest that gene flow between these species has not completely ceased since their initial isolation. Polymorphism data of extant populations also imply recent gene flow between the species. Our study reveals that interspecific gene flow affects the genome evolution in Arabidopsis.

Role of hydroxyl groups in the NH(x) (x = 1-3) adsorption on the TiO2 anatase (101) surface determined by a first-principles study.

A spin-polarized density functional theory calculation was carried out to study the adsorption of NH(x) species (x = 1-3) on a TiO2 anatase (101) surface with and without hydroxyl groups by using first-principles calculations. It was found that the present hydroxyl group has the effect of significantly enhancing the adsorption of monodentate adsorbates H2N-Ti(a) compared to that on a bare surface. The nature of the interaction between the adsorbate (NH(x)) and the hydroxylated or bare surface was analyzed by the Mulliken charge and density of states (DOS) calculations. This facilitation of NH2 is caused by the donation of coadsorbed H filling the nonbonding orbital of NH2, resulting in an electron gain in NH2 from the bonding. In addition, the upper valence band, which originally consisted of the mixing of O 2p and Ti 3d orbitals, has been broadened by the two adjacent H 1s and NH2 sigma(y)(b) orbitals joined to the bottom of the original TiO2 valence band. The results are important to understand the OH effect in heterogeneous catalysis.

The nanoindentation of a copper substrate by single-walled carbon nanocone tips: a molecular dynamics study.

This study dealt with deep nanoindentation of a copper substrate with single-walled carbon nanocones (SWCNCs) as the proximal probe tip, using molecular dynamics (MD) simulations. As an important feature, during the indentation the end part of the SWCNC tip will suffer a narrowing effect due to the radial component of resistant compression from the substrate and then forms into a somewhat flat arrowhead-like shape. The effective cross-sectional area of the SWCNC tip inside the substrate that the resistant force is acting on therefore is reduced to lower the normal resistant force on the tip. The narrowing effect is more significant for longer SWCNC tips. Two categories of SWCNCs are therefore classified according to whether the SWCNC tip buckles at its part inside or outside the substrate. SWCNCs of the first category defined in this paper are found able to indent into the substrate up to a desired depth. Further analyses demonstrate that a longer SWCNC tip of the first category will encounter smaller repulsive force during the indentation and thus require less net work to accomplish the indentation process. Raising temperatures will weaken the narrowing effect, so an SWCNC tip of the first category also encounters greater repulsive force and larger net work in the indentation process performed at a higher temperature. Notably, a permanent hollow hole with high aspect ratio will be produced on the copper substrate, while copper atoms in close proximity to the hole are only slightly disordered, especially when the indentation is manipulated at a lower temperature by using a longer SWCNC tip.

Structure-activity relationships of alphaS1-casomorphin using AM1 calculations and molecular dynamics simulations.

This paper investigates the structure-activity relationships of alphaS1-casomorphin (alphaS1-CM) using AM1 calculations and molecular dynamics (MD) simulations. Previous studies have shown that this peptide has remarkable opioid actions, and not only has a high affinity toward all three subtypes (kappa1-kappa3) of the kappa-opioid sites, but also inhibits the proliferation of the T47D human breast cancer cell line. The systematic conformer search performed by the AM1 calculations is based on the torsional angles of the Val2-Pro3 (omega2) and Phe4-Pro5 (omega4) amide bonds. The AM1 results reveal that the alphaS1-CM conformers strongly favor the cis/cis pair of the omega2/omega4 amide bonds in the minimized energy state. Furthermore, the picture of these stable conformers is found to be a strong interaction of the coulomb's force between two terminuses. MD simulations are performed to investigate the features of both the structural stability and pharmacological activity of alphaS1-CM in aqueous solution. The simulation results reveal that the omega2/omega4 amide bonds favor the cis/cis status in the stable state. Furthermore, the pharmacophoric distance between two aromatic rings is found to be 5.0 approximately 5.4A. The chi1 rotamers of the Tyr and Phe residues show a preference for gauche (-) and trans, respectively. The side chain rotamers of alphaS1-CM are competed to those of other opioid ligands with a known potency and selectivity for delta- and mu-opioid receptors. Finally, we address a likely kappa pharmacophore model compared to the delta pharmacophore model.

Electronic properties of capped, finite-length armchair carbon nanotubes in an electric field.

This study investigates the electronic properties of finite-length armchair carbon nanotubes in an electric field (E) using a single-pi-band tight-binding model. Three different tip configurations are considered, namely, open ends with hydrogen terminations (H-terminations), one end capped with half of C60 fullerene and the other end open with H-terminations, and both ends capped with half of C60 fullerene. In general, the electronic states and energy gaps of low-energy electronic structures are highly sensitive to changes in the direction and magnitude of the applied electric field and to the tip configuration. The present results show that the electric field induces a strong modulation of the state energies and energy gaps of the current nanotubes, changes their energy spacings, and prompts the occurrence of semiconductor-metal transitions (SMTs). It is found that the SMTs occur more frequently as the direction of the electric field approaches the symmetry axis or when its magnitude becomes sufficiently large. The present results also indicate that the Fermi levels and energy gaps of the three nanotubes considered in this study are strongly influenced by the cap configuration. Finally, the convergent decay behavior of the energy gap which is observed as the length of the nanotube is increased is also strongly dependent on the tip configuration.

Exponential synchronization of a class of neural networks with time-varying delays.

This paper aims to present a synchronization scheme for a class of delayed neural networks, which covers the Hopfield neural networks and cellular neural networks with time-varying delays. A feedback control gain matrix is derived to achieve the exponential synchronization of the drive-response structure of neural networks by using the Lyapunov stability theory, and its exponential synchronization condition can be verified if a certain Hamiltonian matrix with no eigenvalues on the imaginary axis. This condition can avoid solving an algebraic Riccati equation. Both the cellular neural networks and Hopfield neural networks with time-varying delays are given as examples for illustration.

Globally asymptotic stability of a class of neutral-type neural networks with delays.

Several stability conditions for a class of systems with retarded-type delays are presented in the literature. However, no results have yet been presented for neural networks with neutral-type delays. Accordingly, this correspondence investigates the globally asymptotic stability of a class of neutral-type neural networks with delays. This class of systems includes Hopfield neural networks, cellular neural networks, and Cohen-Grossberg neural networks. Based on the Lyapunov stability method, two delay-independent sufficient stability conditions are derived. These stability conditions are easily checked and can be derived from the connection matrix and the network parameters without the requirement for any assumptions regarding the symmetry of the interconnections. Two illustrative examples are presented to demonstrate the validity of the proposed stability criteria.

Draft genome of Myxosarcina sp. strain GI1, a baeocytous cyanobacterium associated with the marine sponge Terpios hoshinota.

To date, genome sequences (complete or in draft form) from only six baeocytous cyanobacteria in four genera have been reported: Xenococcus, Chroococcidiopsis, Pleurocapsa, and Stanieria. To expand our knowledge on the diversity of baeocytous cyanobacteria, this study sequenced the genome of GI1, which is a Myxosarcina-like baeocytous cyanobacterium. GI1 is of interest not only because of its phylogenetic niche, but also because it is a cyanobiont isolated from the marine cyanobacteriosponge Terpios hoshinota, which has been shown to cause the death of corals. The ~7 Mb draft GI1 genome contains 6,891 protein-coding genes and 62 RNA genes. A comparison of genomes among the sequenced baeocytous cyanobacterial strains revealed the existence or absence of numerous discrete genes involved in nitrogen metabolism. It will be interesting to determine whether these genes are important for cyanobacterial adaptations and interactions between cyanobionts and their marine sponge hosts.

New bioactive chromanes from Litchi chinensis.

Seven new δ-tocotrienols, designated litchtocotrienols A-G (1-7), together with one glorious macrocyclic analogue, macrolitchtocotrienol A (8), and one new meroditerpene chromane, cyclolitchtocotrienol A (9), were isolated from the leaves of Litchi chinensis. Their structures were mainly determined by extensive spectroscopic analysis, and their biological activities were evaluated by cytotoxicity against human gastric adenocarcinoma cell lines (AGS, ATCC CRL-1739) and hepatoma carcinoma cell line (HepG2 2.2.1.5). The structure-activity relationship of the isolated compounds was also discussed.

Long-lived quantum coherence and non-Markovianity of photosynthetic complexes.

Long-lived quantum coherence in photosynthetic pigment-protein complexes has recently been reported at physiological temperature. It has been pointed out that the discrete vibrational modes may be responsible for the long-lived coherence. Here, we propose an analytical non-Markovian model to explain the origin of the long-lived coherence in pigment-protein complexes. We show that the memory effect of the discrete vibrational modes produces a long oscillating tail in the coherence. We further use the recently proposed measure to quantify the non-Markovianity of the system and find out the prolonged coherence is highly correlated to it.

OMACC: an Optical-Map-Assisted Contig Connector for improving de novo genome assembly.

BACKGROUND: Genome sequencing and assembly are essential for revealing the secrets of life hidden in genomes. Because of repeats in most genomes, current programs collate sequencing data into a set of assembled sequences, called contigs, instead of a complete genome. Toward completing a genome, optical mapping is powerful in rendering the relative order of contigs on the genome, which is called scaffolding. However, connecting the neighboring contigs with nucleotide sequences requires further efforts. Nagarajian et al. have recently proposed a software module, FINISH, to close the gaps between contigs with other contig sequences after scaffolding contigs using an optical map. The results, however, are not yet satisfying. RESULTS: To increase the accuracy of contig connections, we develop OMACC, which carefully takes into account length information in optical maps. Specifically, it rescales optical map and applies length constraint for selecting the correct contig sequences for gap closure. In addition, it uses an advanced graph search algorithm to facilitate estimating the number of repeat copies within gaps between contigs. On both simulated and real datasets, OMACC achieves a <10% false gap-closing rate, three times lower than the ~27% false rate by FINISH, while maintaining a similar sensitivity. CONCLUSION: As optical mapping is becoming popular and repeats are the bottleneck of assembly, OMACC should benefit various downstream biological studies via accurately connecting contigs into a more complete genome. AVAILABILITY: http://140.116.235.124/~tliu/omacc.

Effects of GC bias in next-generation-sequencing data on de novo genome assembly.

Next-generation-sequencing (NGS) has revolutionized the field of genome assembly because of its much higher data throughput and much lower cost compared with traditional Sanger sequencing. However, NGS poses new computational challenges to de novo genome assembly. Among the challenges, GC bias in NGS data is known to aggravate genome assembly. However, it is not clear to what extent GC bias affects genome assembly in general. In this work, we conduct a systematic analysis on the effects of GC bias on genome assembly. Our analyses reveal that GC bias only lowers assembly completeness when the degree of GC bias is above a threshold. At a strong GC bias, the assembly fragmentation due to GC bias can be explained by the low coverage of reads in the GC-poor or GC-rich regions of a genome. This effect is observed for all the assemblers under study. Increasing the total amount of NGS data thus rescues the assembly fragmentation because of GC bias. However, the amount of data needed for a full rescue depends on the distribution of GC contents. Both low and high coverage depths due to GC bias lower the accuracy of assembly. These pieces of information provide guidance toward a better de novo genome assembly in the presence of GC bias.