The architecture of the chloroplast psbA-trnH non-coding region in angiosperms

The psbA-trnH intergenic region is among the most variable regions in the angiosperm chloroplast genome. It is a popular tool for plant population genetics and species level phylogenetics and has been proposed as suitable for DNA barcoding studies. This region contains two parts differing in their evolutionary conservation: 1) the psbA 3′UTR (untranslated region) and 2) the psbA-trnH intergenic non-transcribed spacer. We compared the sequence and RNA secondary structure of the psbA 3′ UTR across angiosperms and found consensus motifs corresponding to the stem portions of the RNA stem-loop structures and a consensus TTAGTGTATA box. The psbA-trnH spacer exhibited patterns that can be explained by the independent evolution of large inversions in the psbA 3′UTR and mutational hot spots in the remaining portion of the psbA-trnH spacer. We conclude that a comparison of chloroplast UTRs across angiosperms offer clues to the identity of putative regulatory elements and information about selective constraints imposed on the chloroplast non-coding regions.     

Phylogeny of the Brachytheciaceae (Bryophyta) based on morphology and sequence level data

Brachytheciaceae is often considered a taxonomically difficult group of mosses. For example, morphological variation has led to difficulty in generic delimitation. We used DNA sequence data (chloroplast psbT‐H and trnL‐F and nuclear ITS2) together with morphology (63 characters) to examine the relationships within this family. The combined unaligned length of the DNA sequences used in the phylogenetic analyses varied between 1277 and 1343 bp. For phylogeny reconstruction we performed direct optimization, as implemented in POY. Analyses were performed with three different gap costs and the morphological data partition was weighted both: (1) equal to gap cost, and (2) with a weight of one. The utility of sensitivity analysis has recently been cast into doubt; hence in this study it was performed only to explore the effects of weighting on homology statements and topologies and to enable more detailed comparisons between earlier studies utilizing the direct optimization method. The wide sequence length variation of non‐coding ITS2 sequences resulted in character optimizations (i.e., “alignments”) of very different lengths when various gap costs were applied. Despite this variation, the topologies of equally parsimonious trees remained fairly stable. The inclusion of several outgroups, instead of only one, was observed to increase the congruence between data sets and to slightly increase the resolution. An inversion event in the 9 bp loop region in the chloroplast psbT‐N spacer in mosses has been postulated to include only uninformative variation, thus possibly negatively impacting the phylogeny reconstruction. Despite this inversion, its variation within Brachytheciaceae was clearly congruent with information from other sources, but inclusion of these 9 bp in the analysis had only a minor effect on the phylogenetic results. In the most parsimonious topology, which was obtained with equal weighting of all data, Meteoriaceae and Brachytheciaceae were resolved as monophyletic sister groups, which had recently been suggested based on a few shared morphological characters. Our study revealed some new generic relationships within the Brachytheciaceae, which are discussed in light of the morphological characters traditionally used for generic delimitation.     

Evaluation of six candidate DNA barcoding loci in Ficus (Moraceae) of China

Ficus, with about 755 species, diverse habits and complicated co‐evolutionary history with fig wasps, is a notoriously difficult group in taxonomy. DNA barcoding is expected to bring light to the identification of Ficus but needs evaluation of candidate loci. Based on five plastid loci (rbcL, matK, trnH‐psbA, psbK‐psbI, atpF‐atpH) and a nuclear locus [internal transcribed spacer (ITS)], we calculated genetic distances and DNA barcoding gaps individually and in combination and constructed phylogenetic trees to test their ability to distinguish the species of the genus. A total of 228 samples representing 63 putative species in Ficus (Moraceae) of China were included in this study. The results demonstrated that ITS has the most variable sites, greater intra‐ and inter‐specific divergences, the highest species discrimination rate (72%) and higher primer universality among the single loci. It is followed by psbK‐psbI and trnH‐psbA with moderate variation and considerably lower species discrimination rates (about 19%), whereas matK, rbcL and atpF‐atpH could not effectively separate the species. Among the possible combinations of loci, ITS + trnH‐psbA performed best but only marginally improved species resolution over ITS alone (75% vs. 72%). Therefore, we recommend using ITS as a single DNA barcoding locus in Ficus.     

Solution structures and backbone dynamics of the ribosomal protein S6 and its permutant P54‐55

The ribosomal protein S6 from Thermus thermophilus has served as a model system for the study of protein folding, especially for understanding the effects of circular permutations of secondary structure elements. This study presents the structure of a permutant protein, the 96‐residue P^54‐55, and the structure of its 101‐residue parent protein S6^wt in solution. The data also characterizes the effects of circular permutation on the backbone dynamics of S6. Consistent with crystallographic data on S6^wt, the overall solution structures of both P^54‐55 and S6^wt show a β‐sheet of four antiparallel β‐strands with two α‐helices packed on one side of the sheet. In clear contrast to the crystal data, however, the solution structure of S6^wt reveals a disordered loop in the region between β‐strands 2 and 3 (Leu43‐Phe60) instead of a well‐ordered stretch and associated hydrophobic mini‐core observed in the crystal structure. Moreover, the data for P^54‐55 show that the joined wild‐type N‐ and C‐terminals form a dynamically robust stretch with a hairpin structure that complies with the in silico design. Taken together, the results explain why the loop region of the S6^wt structure is relatively insensitive to mutational perturbations, and why P^54‐55 is more stable than S6^wt: the permutant incision at Lys54‐Asp55 is energetically neutral by being located in an already disordered loop whereas the new hairpin between the wild‐type N‐ and C‐termini is stabilizing.     

Contrasting population genetic structures using allozymes and the inversion polymorphism in Drosophila buzzatii

Second chromosome inversion and genotypic frequencies at seven allozyme loci, differentially associated with inversions, were determined in seven natural populations of Drosophila buzzatii. The patterns of variation of allozymes and the inversion polymorphisms were significantly different, indicating the role of adaptive differentiation for the latter. Moreover, the patterns of population structure varied among allozyme loci, suggesting the operation of diversifying selection for certain loci. Differentiation was negligible for Leucyl‐amino peptidase (Lap) and Peptidase‐2 (Pep‐2), low to moderate for Aldehyde oxidase (Aldox), Peptidase‐1 (Pep‐1) and Esterase‐1 (Est‐1) and high for Esterase‐2 (Est‐2) and Xanthine dehydrogenase (Xdh). Significant linkage disequilibria were detected between inversions and Aldox, Est‐1, Est‐2 and Xdh. Multiple regression analyses of inversion and allele frequencies on environmental variables revealed the existence of clines for inversions, Est‐1, Est‐2, Xdh and Aldox along altitudinal, latitudinal and/or climatic gradients. Tests using conditional allele frequencies showed that Est‐1 and Aldox clines could be accounted for by hitchhiking with inversions, whereas natural selection should be invoked to explain the clines observed for Est‐2 and Xdh.     

Molecular,genetic and evolutionary analysis of a paracentric inversion in Arabidopsis thaliana

Chromosomal inversions can provide windows onto the cytogenetic, molecular, evolutionary and demographic histories of a species. Here we investigate a paracentric 1.17‐Mb inversion on chromosome 4 of Arabidopsis thaliana with nucleotide precision of its borders. The inversion is created by Vandal transposon activity, splitting an F‐box and relocating a pericentric heterochromatin segment in juxtaposition with euchromatin without affecting the epigenetic landscape. Examination of the RegMap panel and the 1001 Arabidopsis genomes revealed more than 170 inversion accessions in Europe and North America. The SNP patterns revealed historical recombinations from which we infer diverse haplotype patterns, ancient introgression events and phylogenetic relationships. We find a robust association between the inversion and fecundity under drought. We also find linkage disequilibrium between the inverted region and the early flowering Col‐FRIGIDA allele. Finally, SNP analysis elucidates the origin of the inversion to South‐Eastern Europe approximately 5000 years ago and the FRI‐Col allele to North‐West Europe, and reveals the spreading of a single haplotype to North America during the 17th to 19th century. The ‘American haplotype’ was identified from several European localities, potentially due to return migration.     

The use of a non-coding region of chloroplast DNA in phylogenetic studies of the subtribeSonchinae (Asteraceae:Lactuceae)

The systematic utility of sequences from a non-coding region of chloroplast DNA (cpDNA) betweenpsbA andtrnH^(GUG) was examined by assessing phylogenetic relationships in subtribeSonchinae (Asteraceae:Lactuceae). Primers constructed against highly conserved regions of tRNA genes were used for PCR amplification and sequencing. ThepsbA-trnH intergenic spacer contains several insertions and deletions (indels) inSonchinae with the length varying from 385 to 450 bp. Sequence divergence ranges from 0.00% to 7.54% withinSonchinae, with an average of 2.4%. Average sequence divergence inSonchus subg.Sonchus is 2.0%, while the mean for subg.Dendrosonchus and its close relatives in Macaronesia (the woodySonchus alliance) is 1.0%. Our results suggest that this region does not evolve rapidly enough to resolve relationships among closely related genera or insular endemics in theAsteraceae. The phylogenetic utility ofpsbA-trnH sequences of the non-coding cpDNA was compared to sequences from the ITS region of nuclear ribosomal DNA. The results suggest that ITS sequences evolve nearly four times faster thanpsbA-trnH intergenic spacer sequences. Furthermore, the ITS sequences provide more variable and phylogenetically informative sites and generate more highly resolved trees with more strongly supported clades, and thus are more suitable for phylogenetic comparisons at lower taxonomic levels than thepsbA-trnH intergenic chloroplast sequences.     

Speciation and chromosomal evolution of the Baikalian endemic chironomids of the genus Sergentia Kief. (Diptera, Chironomidae): Karyotype divergence and chromosomal polymorphism in the populations of deep-water species Sergentia nebulosa Linevitsh et al. and Sergentia assimilis Proviz V. et Proviz L.

Specific karyotype structure and chromosomal polymorphism was investigated in the populations of Sergentia nebulosa Linevitsh et al., 1984 and Sergentia assimilis Proviz V. et Proviz L., 1999, the deep-water endemic chironomid species (Diptera, Chironomidae) from the Baikal Lake. The distinguishing feature of the karyotypes of these species, compared to the other Baikalian Sergentia, is well-developed nucleolus in region 6 of arm C. Both species display the presence of interspecific population polymorphism, determined by the structure of this arm. In some populations, chromosome regions from 4 to 6 contain a homozygous inversion, which is absent in the other populations. The distinguishing karyotype feature of S. assimilis, which shares fluctuating homozygous inversions with the other species, is the presence of two species-specific homozygous inversions. These are the secondary overlapping inversion in arm A, regions 2 to 7, and the inversion in regions 4 to 10 of arm G. Both species of interest contain nucleolus organizer in region 10 of arm G. In populations of S. nebulosa, six heterozygous inversions localized in arms A, B, C, F, and G were discovered. The highest number of heterozygotes for inversions (71%) was observed in the population from Southern Baikal. In arm B of S. assimilis, one heterozygous inversion and heterozygosity for nucleolus organizer in the chromosome region 16 was detected. Chromosomal evolution of Baikalian Sergentia, and the role of inversion polymorphism in the population adaptation is discussed. Original Russian Text ? V.I. Proviz, 2008, published in Genetika, 2008, Vol. 44, No. 12, pp. 1627–1637.     

Efficient induction of heritable inversions in plant genomes using the CRISPR/Cas system

During the evolution of plant genomes, sequence inversions occurred repeatedly making the respective regions inaccessible for meiotic recombination and thus for breeding. Therefore, it is important to develop technologies that allow the induction of inversions within chromosomes in a directed and efficient manner. Using the Cas9 nuclease from Staphylococcus aureus (SaCas9), we were able to obtain scarless heritable inversions with high efficiency in the model plant Arabidopsis thaliana. Via deep sequencing, we defined the patterns of junction formation in wild‐type and in the non‐homologous end‐joining (NHEJ) mutant ku70‐1. Surprisingly, in plants deficient of KU70, inversion induction is enhanced, indicating that KU70 is required for tethering the local broken ends together during repair. However, in contrast to wild‐type, most junctions are formed by microhomology‐mediated NHEJ and thus are imperfect with mainly deletions, making this approach unsuitable for practical applications. Using egg‐cell‐specific expression of Cas9, we were able to induce heritable inversions at different genomic loci and at intervals between 3 and 18 kb, in the percentage range, in the T1 generation. By screening individual lines, inversion frequencies of up to the 10% range were found in T2. Most of these inversions had scarless junctions and were without any sequence change within the inverted region, making the technology attractive for use in crop plants. Applying our approach, it should be possible to reverse natural inversions and induce artificial ones to break or fix linkages between traits at will.     

STRUCTURE AND POPULATION GENETICS OF THE BREAKPOINTS OF A POLYMORPHIC INVERSION IN DROSOPHILA SUBOBSCURA

Drosophila subobscura is a paleartic species of the obscura group with a rich chromosomal polymorphism. To further our understanding on the origin of inversions and on how they regain variation, we have identified and sequenced the two breakpoints of a polymorphic inversion of D. subobscura—inversion 3 of the O chromosome—in a population sample. The breakpoints could be identified as two rather short fragments (～300 bp and 60 bp long) with no similarity to any known transposable element family or repetitive sequence. The presence of the ～300‐bp fragment at the two breakpoints of inverted chromosomes implies its duplication, an indication of the inversion origin via staggered double‐strand breaks. Present results and previous findings support that the mode of origin of inversions is neither related to the inversion age nor species‐group specific. The breakpoint regions do not consistently exhibit the lower level of variation within and stronger genetic differentiation between arrangements than more internal regions that would be expected, even in moderately small inversions, if gene conversion were greatly restricted at inversion breakpoints. Comparison of the proximal breakpoint region in species of the obscura group shows that this breakpoint lies in a small high‐turnover fragment within a long collinear region (～300 kb).     

Limits of the distal inversion in the t complex of the house mouse: Evidence from linkage disequilibria

The suppression of crossing-over and the consequent linkage disequilibrium of genetic markers within the t complex of the house mouse is caused by two large and two short inversions. The inversions encompass a region that is some 15 centiMorgans (cM) long in the homologous wild-type chromosome. The limits of the proximal inversions are reasonably welldefined, those of the distal inversions much less so. We have recently obtained seven new DNA markers (D17Tu) which in wild-type chromosomes map into the region presumably involved in the distal inversions of the t chromosomes. To find out whether the corresponding loci do indeed reside within the inversions, we have determined their variability among 26 complete and 12 partial t haplotypes. In addition, we also tested the same collection of t haplotypes for their variability at five D17Leh, Hba-ps4, Pim-1, and Crya-1 loci. The results suggest that the distal end of the most distal inversion lies between the loci D17Leh467 and D17Tu26. The proximal end of the large distal inversion was mapped to the region between the D17Tu43 and Hba-ps4 loci, but this assignment is rather ambiguous. The loci Pim-1, Crya-1, and the H-2 complex, which have been mapped between the Hba-sp4 and Grr within the large distal inversion, behave as if they recombine from time to time with their wildtype homologs.     

Molecular phylogeny of Solms‐laubachia (Brassicaceae) s.l., based on multiple nuclear and plastid DNA sequences,and its biogeographic implications

Abstract The Hengduan Mountains region of south‐west China is a noted biodiversity hotspot, but the geographic origins and historical assembly of its rich endemic flora, including the sky‐island species of Solms‐laubachia Muschl. (Brassicaceae), have been little studied. Previous molecular studies on the phylogeny of Solms‐laubachia showed it to be paraphyletic, leading to considerable expansion not only of its taxonomic limits, but also its geographic range, with the inclusion of taxa from outside the Hengduan region. However, these studies provided little resolution of interspecific relationships, preventing inferences about historical biogeography within the clade. In the present study, new sequence data from two nuclear genes (LEAFY and G3pdh) and two chloroplast intergenic spacers (petN–psbM and psbM–trnD) were combined with existing markers to increase phylogenetic signals. Phaeonychium villosum (Maxim.) Al‐Shehbaz was found to be nested within Solms‐laubachia s.l. In general, phylogenetic relationships appear to be a good predictor of geography, with the Hengduan Mountain endemics embedded in a paraphyletic grade of species from the western Himalayas and central Asia, but they also imply morphological homoplasy. Incongruence was detected between the nuclear and chloroplast gene trees, perhaps resulting from incomplete lineage sorting of ancestral polymorphisms. The crown age of Solms‐laubachia s.l. was estimated to be approximately 1.42–3.68 mya, using Bayesian relaxed molecular clock analysis. Historical biogeographic analysis using a parametric dispersal–extinction–cladogenesis model inferred central Asia and the western Himalayas as most probable ancestral range of Solms‐laubachia s.l., and estimated higher rates of eastward expansion than westward during the diversification of descendant lineages. In summary, our results suggest that Solms‐laubachia s.l. originated during the Pliocene in central Asia, and subsequently migrated eastward into the Hengduan Mountains, colonizing sky‐island, alpine scree‐slope habitats that may have provided novel ecological opportunity and accelerated speciation, ultimately establishing this region as the present center of diversity of the genus.     

Widespread occurrence of small inversions in the chloroplast genomes of land plants

Large inversions are well characterized in the chloroplast genomes of land plants. In contrast, reports of small inversions are rare and involve limited plant groups. In this study, we report the widespread occurrence of small inversions ranging from 5 to 50 bp in fully and partially sequenced chloroplast genomes of both monocots and dicots. We found that small inversions were much more common than large inversions. The small inversions were scattered over the chloroplast genome including the IR, SSC, and LSC regions. Several small inversions were uncovered in chloroplast genomes even though they shared the same overall gene order. The majority of these small inversions were located within 100 bp downstream of the 3' ends of genes. All had inverted repeat sequences, ranging from 11 to 24 bp, at their ends. Such small inversions form stem-loop hairpin structures that usually have the function of stabilizing the corresponding mRNA molecules. Intra-molecular recombination between the inverted sequences in the stem-forming regions are responsible for generating flip-flop orientations of the loops. The presence of two different orientations of the stem-loop in the trnL-F noncoding region of a single species of Jasminum elegans suggests that a short inversion can be generated within a short period of time. Small inversions of non-coding sequences may influence sequence alignment and character interpretation in phylogeny reconstructions, as shown in nine species of Jasminum. Many small inversions may have been generated by parallel or back mutation events during chloroplast genome evolution. Our data indicate that caution is needed when using chloroplast non-coding sequences for phylogenetic analysis.     

Evolutionary history of Gymnocarpos (Caryophyllaceae) in the arid regions from North Africa to Central Asia

Gymnocarpos has only about ten species distributed in the arid regions of Asia and Africa, but it exhibits a geographical disjunction between eastern Central Asia and western North Africa and Minor Asia. We sampled eight species of the genus and sequenced two chloroplast regions (rps16 and psbB–psbH), and the nuclear rDNA (ITS) to study the phylogeny and biogeography. The results of the phylogenetic analyses corroborated that Gymnocarpos is monophyletic, in the phylogenetic tree two well supported clades are recognized: clade 1 includes Gymnocarpos sclerocephalus and G. decandrus, mainly the North African group, whereas clade 2 comprises the remaining species, mainly in the Southern Arabian Peninsula. Molecular dating analysis revealed that the divergence age of Gymnocarpos was c. 31.33 Mya near the Eocene and Oligocene transition boundary, the initial diversification within Gymnocarpos dated to c. 6.69 Mya in the late Miocene, and the intraspecific diversification mostly occurred during the Quaternary climate oscillations. Ancestral area reconstruction suggested that the Southern Arabian Peninsula was the ancestral area for Gymnocarpos. Our conclusions revealed that the aridification since mid‐late Miocene significantly affected the diversification of the genus in these areas.     

Genomic evidence for role of inversion 3RP of Drosophila melanogaster in facilitating climate change adaptation

Chromosomal inversion polymorphisms are common in animals and plants, and recent models suggest that alternative arrangements spread by capturing different combinations of alleles acting additively or epistatically to favour local adaptation. It is also thought that inversions typically maintain favoured combinations for a long time by suppressing recombination between alternative chromosomal arrangements. Here, we consider patterns of linkage disequilibrium and genetic divergence in an old inversion polymorphism in Drosophila melanogaster (In(3R)Payne) known to be associated with climate change adaptation and a recent invasion event into Australia. We extracted, karyotyped and sequenced whole chromosomes from two Australian populations, so that changes in the arrangement of the alleles between geographically separated tropical and temperate areas could be compared. Chromosome‐wide linkage disequilibrium (LD) analysis revealed strong LD within the region spanned by In(3R)Payne. This genomic region also showed strong differentiation between the tropical and the temperate populations, but no differentiation between different karyotypes from the same population, after controlling for chromosomal arrangement. Patterns of differentiation across the chromosome arm and in gene ontologies were enhanced by the presence of the inversion. These data support the notion that inversions are strongly selected by bringing together combinations of genes, but it is still not clear if such combinations act additively or epistatically. Our data suggest that climatic adaptation through inversions can be dynamic, reflecting changes in the relative abundance of different forms of an inversion and ongoing evolution of allelic content within an inversion.     

Molecular Evolution and Phylogenetic Utility of the Chloroplastrpl16Intron inChusqueaand the Bambusoideae (Poaceae)

Phylogenetic relationships withinChusquea,a diverse genus of neotropical woody bamboos, and among selected members of the Bambusoideae were explored usingrpl16intron sequence data from the chloroplast genome. Mechanisms of mutation, including slipped-strand mispairing, secondary structure, minute inversions, and base substitutions, were examined within therpl16intron, and their effects on sequence alignment and phylogenetic analysis were investigated. Thirty-five bamboo sequences were generated and two separate matrices were analyzed using maximum parsimony. In the first, 23 sequences fromChusquea,1 ofNeurolepis,and 3 outgroups were included.Neurolepiswas supported as sister toChusquea, Chusqueawas strongly supported as a monophyletic lineage, and three species ofChusqueasubg.Rettbergiawere resolved as the most basal clade within the genus. In the second analysis, 15 sequences, 14 from across the subfamily and 1 outgroup, were included. A Bambusoideae clade was recovered with the Olyreae/Parianeae (herbaceous bamboos) and the Bambuseae (woody bamboos) each supported as monophyletic. Two clades corresponding to temperate and tropical woody bamboos were derived within the Bambuseae and the tropical taxa were further split into New World and Old World clades. Therpl16intron in bamboos was found to be susceptible to frequent length mutations of multiple origins, nonindependent character evolution, and regions of high mutability, all of which created difficulties in alignment and phylogenetic analysis; nonetheless therpl16intron is phylogenetically informative at the inter- and intrageneric levels in bamboos.