Sequencing of 15 622 gene‐bearing BACs clarifies the gene‐dense regions of the barley genome

Barley (Hordeum vulgare L.) possesses a large and highly repetitive genome of 5.1 Gb that has hindered the development of a complete sequence. In 2012, the International Barley Sequencing Consortium released a resource integrating whole‐genome shotgun sequences with a physical and genetic framework. However, because only 6278 bacterial artificial chromosome (BACs) in the physical map were sequenced, fine structure was limited. To gain access to the gene‐containing portion of the barley genome at high resolution, we identified and sequenced 15 622 BACs representing the minimal tiling path of 72 052 physical‐mapped gene‐bearing BACs. This generated ~1.7 Gb of genomic sequence containing an estimated 2/3 of all Morex barley genes. Exploration of these sequenced BACs revealed that although distal ends of chromosomes contain most of the gene‐enriched BACs and are characterized by high recombination rates, there are also gene‐dense regions with suppressed recombination. We made use of published map‐anchored sequence data from Aegilops tauschii to develop a synteny viewer between barley and the ancestor of the wheat D‐genome. Except for some notable inversions, there is a high level of collinearity between the two species. The software HarvEST:Barley provides facile access to BAC sequences and their annotations, along with the barley–Ae. tauschii synteny viewer. These BAC sequences constitute a resource to improve the efficiency of marker development, map‐based cloning, and comparative genomics in barley and related crops. Additional knowledge about regions of the barley genome that are gene‐dense but low recombination is particularly relevant.     

Sex‐biased transcriptome divergence along a latitudinal gradient

Sex‐dependent gene expression is likely an important genomic mechanism that allows sex‐specific adaptation to environmental changes. Among Drosophila species, sex‐biased genes display remarkably consistent evolutionary patterns; male‐biased genes evolve faster than unbiased genes in both coding sequence and expression level, suggesting sex differences in selection through time. However, comparatively little is known of the evolutionary process shaping sex‐biased expression within species. Latitudinal clines offer an opportunity to examine how changes in key ecological parameters also influence sex‐specific selection and the evolution of sex‐biased gene expression. We assayed male and female gene expression in Drosophila serrata along a latitudinal gradient in eastern Australia spanning most of its endemic distribution. Analysis of 11 631 genes across eight populations revealed strong sex differences in the frequency, mode and strength of divergence. Divergence was far stronger in males than females and while latitudinal clines were evident in both sexes, male divergence was often population specific, suggesting responses to localized selection pressures that do not covary predictably with latitude. While divergence was enriched for male‐biased genes, there was no overrepresentation of X‐linked genes in males. By contrast, X‐linked divergence was elevated in females, especially for female‐biased genes. Many genes that diverged in D. serrata have homologs also showing latitudinal divergence in Drosophila simulans and Drosophila melanogaster on other continents, likely indicating parallel adaptation in these distantly related species. Our results suggest that sex differences in selection play an important role in shaping the evolution of gene expression over macro‐ and micro‐ecological spatial scales.     

Expression of tomato reference genes using established primer sets: Stability across experimental set‐ups

Reliable reference genes are critical for relative quantification using quantitative real‐time PCR (qPCR). Ten tomato genes (Solanum lycopersicum) and their respective primer sets, which have been used over the last 6 years as references in expression studies, were evaluated for their performance using leaf tissue samples grown under semi‐controlled conditions and infected with grey mould (Botrytis cinerea) or late blight (Phytophthora infestans). The target genes coding for U6 snRNA‐associated Sm‐like protein LSm7, calcineurin B‐like protein and V‐type proton ATPase were the most stable expressed of all the genes tested in three experimental repetitions. Evaluation of candidate reference genes with geNorm and NormFinder softwares yielded the lowest mean values for their respective primer sets LSM7, SlCBL1 and SlATPase, suggesting stable expression. However, SlATPase primer set revealed a comparably high intra‐group variation and was thus not considered further. In follow‐up experiments with P. infestans, the geNorm and NormFinder values of primer sets LSM7 and SlCBL1 were even lower, indicating the stability of their expression also under these conditions. Primer efficiency differed by ‐18 to +5 percentage points from values presented in the literature. Our findings show that a reference primer set which delivers the best results in one system may be outperformed by another under different experimental conditions, thus recommending a reassessment of both expression stability and qPCR efficiency whenever the biological or technical experimental set‐up is changed. On the basis of our results, we recommend the use of LSM7 and SlCBL1 as reference primer sets for gene expression studies on plant tissue derived from open or semi‐controlled conditions.     

Molecular Characterization of Aster Yellows Subgroup 16SrI‐B Phytoplasma in Verbena × hybrida

Symptoms resembling phytoplasma disease were observed on Verbena × hybrida in Alanya, Turkey, during October 2013. Infected plants were growing as perennials in a flower border and showed symptoms of discoloured flowers, poor flower clusters, inflorescences with a small number of developed flowers and thickened fruit stalks. Electron microscopy examination of the ultra‐thin sections revealed polymorphic bodies in the phloem tissue of leaf midribs. The phytoplasma aetiology of this disease was confirmed by polymerase chain reaction of the 16S rRNA gene, the 16–23S rRNA intergenic spacer region and the start of the 23S rRNA gene using universal phytoplasma‐specific primer pair P1A/P7A, two ribosomal protein (rp) genes (rpl22 and rps3) (the group‐specific primer pair rp(I)F1A/rp(I)R1A) and the Tuf gene (group‐specific fTufAy/rTufAy primers) generating amplicons of 1.8 kbp, 1.2 kbp and 940 bp, respectively. Comparison of the amplified sequences with those available in GenBank allowed classification of the phytoplasma into aster yellows subgroups 16SrI‐B, rpI‐B and tufI‐B. This is the first report about molecular detection and identification of natural infection of the genus Verbena by phytoplasma and occurrence of the aster yellows group phytoplasma on an ornamental plant in Turkey.     

GABI‐DUPLO: a collection of double mutants to overcome genetic redundancy in Arabidopsis thaliana

Owing to duplication events in its progenitor, more than 90% of the genes in the Arabidopsis thaliana genome are members of multigene families. A set of 2108 gene families, each consisting of precisely two unlinked paralogous genes, was identified in the nuclear genome of A. thaliana on the basis of sequence similarity. A systematic method for the creation of double knock‐out lines for such gene pairs, designated as DUPLO lines, was established and 200 lines are now publicly available. Their initial phenotypic characterisation led to the identification of seven lines with defects that emerge only in the adult stage. A further six lines display seedling lethality and 23 lines were lethal before germination. Another 14 lines are known to show phenotypes under non‐standard conditions or at the molecular level. Knock‐out of gene pairs with very similar coding sequences or expression profiles is more likely to produce a mutant phenotype than inactivation of gene pairs with dissimilar profiles or sequences. High coding sequence similarity and highly similar expression profiles are only weakly correlated, implying that promoter and coding regions of these gene pairs display different degrees of diversification.     

Single‐copy gene‐based chromosome painting in cucumber and its application for chromosome rearrangement analysis in Cucumis

Chromosome painting based on fluorescence in situ hybridization (FISH) has played an important role in chromosome identification and research into chromosome rearrangements, diagnosis of chromosome abnormalities and evolution in human and animal species. However, it has not been applied widely in plants due to the large amounts of dispersed repetitive sequences in chromosomes. In the present work, a chromosome painting method for single‐copy gene pools in Cucumis sativus was successfully developed. Gene probes with sizes above 2 kb were detected consistently. A cucumber karyotype was constructed based on FISH using a cocktail containing chromosome‐specific gene probes. This single‐copy gene‐based chromosome painting (ScgCP) technique was performed by PCR amplification, purification, pooling, labeling and hybridization onto chromosome spreads. Gene pools containing sequential genes with an interval less than 300 kb yielded painting patterns on pachytene chromosomes. Seven gene pools corresponding to individual chromosomes unambiguously painted each chromosome pair of C. sativus. Three mis‐aligned regions on chromosome 4 were identified by the painting patterns. A probe pool comprising 133 genes covering the 8 Mb distal end of chromosome 4 was used to evaluate the potential utility of the ScgCP technique for chromosome rearrangement research through cross‐species FISH in the Cucumis genus. Distinct painting patterns of this region were observed in C. sativus, C. melo and C. metuliferus species. A comparative chromosome map of this region was constructed between cucumber and melon. With increasing sequence resources, this ScgCP technique may be applied on any other sequenced species for chromosome painting research.     

Development of Streptococcus gordonii‐specific quantitative real‐time polymerase chain reaction primers based on the nucleotide sequence of rpoB

In this study, Streptococcus gordonii‐specific quantitative real‐time polymerase chain reaction (qPCR) primers, RTSgo‐F2/RTSgo‐R2, were developed based on the nucleotide sequences of RNA polymerase β‐subunit gene (rpoB). The specificity of the RTSgo‐F2/RTSgo‐R2 primers was assessed by conventional PCR on 99 strains comprising 63 oral bacterial species, including the type strain and eight clinical isolates of S. gordonii. PCR products were amplified from the genomic DNAs of only S. gordonii strains. The qPCR primers were able to detect as little as 40 fg of S. gordonii genomic DNA at a cycle threshold value of 33. These findings suggest that these qPCR primers detect S. gordonii with high specificity and sensitivity.     

Mitochondrial phylogeography and subspecies of the wide‐ranging sub‐Saharan leopard tortoise Stigmochelys pardalis (Testudines: Testudinidae) – a case study for the pitfalls of pseudogenes and GenBank sequences

The leopard tortoise (Stigmochelys pardalis) is the most widely distributed sub‐Saharan tortoise species, with a range extending from the Horn of Africa all over eastern Africa to the Republic of South Africa, Namibia and southernmost Angola. Using 1938 bp of mitochondrial DNA (cyt b gene, partial ND4 gene plus adjacent tRNA genes) from a nearly range‐wide sampling, we examined its phylogeographic structure and compared our findings with previously published GenBank sequences. We identified seven major clades that are largely parapatrically distributed. A few records of distinct haplotypes at the same locality or in close proximity could be the result of translocation of tortoises by man. The greatest diversity occurs in the south of the species’ range, with five out of the seven clades. Testing for isolation‐by‐distance suggests that the observed phylogeographic structure is the result of restricted geographical gene flow and not of historical vicariance. This is in sharp contrast to wide‐ranging thermophilic reptiles from the western Palaearctic, whose phylogeographic structure was significantly shaped by Pleistocene range interruptions, but also by earlier dispersal and vicariant events. Most cyt b sequences of S. pardalis from GenBank turned out to be nuclear pseudogenes, or to be of chimerical origin from such pseudogenes and authentic mitochondrial sequences, which argues for caution regarding uncritical usage of GenBank sequences. The recent revalidation of the two subspecies of S. pardalis was based on such a chimerical sequence that was erroneously identified with the subspecies S. p. babcocki. Furthermore, according to our data, the distribution of mitochondrial clades does match neither the traditional subspecies ranges nor the pronounced geographical size variation of leopard tortoises. We conclude that there is no rationale for recognizing subspecies within S. pardalis.     

Recombinant production of omega‐3 fatty acids in Escherichia coli using a gene cluster isolated from Shewanella baltica MAC1

Aim: To isolate eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) genes from Shewanella baltica MAC1 and to examine recombinant production of EPA and DHA in E. coli to investigate cost‐effective, sustainable and convenient alternative sources for fish oils. Methods and Results: A fosmid library was prepared from the genomic DNA of S. baltica MAC1 and was screened for EPA and DHA genes by colony hybridization using a partial fragment of the S. baltica MAC1 pfaA and pfaD genes as probes. Analysis of total fatty acids isolated from transgenic E. coli positive for pfaA and pfaD genes by gas chromatography and gas chromatography‐mass spectrometry indicated recombinant production of both EPA and DHA. Analysis of the complete nucleotide sequence for the isolated gene cluster showed 16 putative open reading frames (ORFs). Among those, four ORFs showed homology with pfaA, pfaB, pfaC and pfaD genes of the EPA and/or DHA biosynthesis gene clusters; however, the protein domains of these genes were different from other EPA/DHA biosynthesis genes. Conclusions: The EPA and DHA gene cluster was cloned successfully. The transgenic E. coli strain carrying the omega‐3 gene cluster was able to produce both EPA and DHA. The isolated gene cluster contained all the genes required for the recombinant production of both EPA and DHA in E. coli. Significance and Impact of the Study: These findings have implications for any future use of the EPA and DHA gene cluster in other micro‐organisms, notably those being used for fermentation. Recombinant production of both EPA and DHA by E. coli or any other micro‐organism has great potential to add economic value to a variety of industrial and agricultural products.     

Detection and Molecular Characterization of a 16SrII‐A* Phytoplasma in Grapefruit (Citrus paradisi) with Huanglongbing‐like Symptoms in China

In the year 2010, in a survey in Guangxi Province, China, to detect and characterize phytoplasmas in a huanglongbing (HLB)‐infected grapefruit (Citrus paradisi) orchard, 87 leaf samples with symptoms of blotchy mottle were collected from symptomatic grapefruit trees, and 320 leaf samples from symptomless trees adjacent to the symptomatic trees. Nested polymerase chain reaction (PCR) using universal phytoplasma primer set P1/P7 followed by primer set fU5/rU3 identified 7 (8.0%) positive samples from symptomatic samples but none from symptomless samples. Of the 87 symptomatic samples, 77 (88.5%) were positive for ‘Candidatus Liberibacter asiaticus’ and 5 for both phytoplasma and ‘Ca. L. asiaticus’. Sequence analysis indicated that seven 881‐bp amplicons, amplified by nested phytoplasma primer sets P1/P7 and fU5/rU3, shared 100.0% sequence identity with each other. Genome walking was then performed based on the 881 bp known sequences, and 5111 bp of upstream and downstream sequences were obtained. The total 5992 bp sequences contained a complete rRNA operon, composed of a 16S rRNA gene, a tRNA^Ile gene, a 23S rRNA gene and a 5S rRNA gene followed by eight tRNA genes. Phylogenetic analysis and virtual restriction fragment length polymorphism analysis confirmed the phytoplasma was a variant (16SrII‐A*) of phytoplasma subgroup 16SrII‐A. As phytoplasmas were only detected in blotchy‐mottle leaves, the 16SrII‐A* phytoplasma identified was related to HLB‐like symptoms.     

Identification of ‘Candidatus Phytoplasma asteris' (16SrI‐B) Causing Witches' Broom Disease of Mallotus japonicus in Korea

Mallotus japonicus with witches' broom disease were observed in Jeollabuk‐do, Korea. A phytoplasma from the infected leaves was identified, based on the 16S rDNA, 16S‐23S intergenic spacer region, and fragment of rp operon and tuf gene sequences. The 16S rDNA sequences exhibited maximum (99.7%) similarity with Iranian lettuce phytoplasma, the rp operon sequences exhibited 100% similarity with Goldenrain stunt phytoplasma, and the tuf gene sequences exhibited 99.8% similarity with Japanese spurge yellows phytoplasma. Results of the sequence analysis and phylogenetic studies confirmed that the phytoplasma associated with M. japonicus in Korea was an isolate of Aster Yellows group (subgroup16SrI‐B).     

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast‐Z effect in Heliconius

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex‐ and tissue‐specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male‐biased than female‐biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic‐specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast‐Z effect where hemizygosity influences the rate of adaptive substitutions.