Species-specific Gene Families Research Articles

Cicer super-pangenome provides insights into species evolution and agronomic trait loci for crop improvement in chickpea.

Chickpea (Cicer arietinum L.)-an important legume crop cultivated in arid and semiarid regions-has limited genetic diversity. Efforts are being undertaken to broaden its diversity by utilizing its wild relatives, which remain largely unexplored. Here, we present the Cicer super-pangenome based on the de novo genome assemblies of eight annual Cicer wild species. We identified 24,827 gene families, including 14,748 core, 2,958 softcore, 6,212 dispensable and 909 species-specific gene families. The dispensable genome was enriched for genes related to key agronomic traits. Structural variations between cultivated and wild genomes were used to construct a graph-based genome, revealing variations in genes affecting traits such as flowering time, vernalization and disease resistance. These variations will facilitate the transfer of valuable traits from wild Cicer species into elite chickpea varieties through marker-assisted selection or gene-editing. This study offers valuable insights into the genetic diversity and potential avenues for crop improvement in chickpea.

Genomic and Transcriptomic Insights into the Genetic Basis of Foam Secretion in Rice Spittlebug Callitettix versicolor.

Many animal species produce protective foams, the majority of which exhibit evolutionary adaptability. Although the function and composition of foams have been widely studied, the genetic basis of foam secretion remains unknown. Unlike most species that produce foam under specific situations, spittlebugs continuously secrete foams throughout all nymphal stages. Here, we capitalize on the rice spittlebug (Callitettix versicolor) to explore the genetic basis of foam secretion through genomic and transcriptomic approaches. Our comparative genomic analysis for C. versicolor and eight other insect species reveals 606 species-specific gene families and 66 expanded gene families, associated with carbohydrate and lipid metabolism. These functions are in accordance with the composition of foams secreted by spittlebugs. Transcriptomic analyses of malpighian tubules across developmental stages detected 3192 differentially expressed genes. Enrichment analysis of these genes highlights functions also revealed by our comparative genomic analysis and aligns with previous histochemical and morphological observations of foam secretion. This consistency suggests the important roles of these candidate genes in foam production. Our study not only provides novel insights into the genetic basis of foam secretion in rice spittlebugs but also contributes valuable knowledge for future evolutionary studies of spittlebugs and the development of pest control strategies for C. versicolor.

Genome and pan-genome assembly of asparagus bean (Vigna unguiculata ssp. sesquipedialis) reveal the genetic basis of cold adaptation.

Asparagus bean (Vigna unguiculata ssp. sesquipedialis) is an important cowpea subspecies. We assembled the genomes of Ningjiang 3 (NJ, 550.31 Mb) and Dubai bean (DB, 564.12 Mb) for comparative genomics analysis. The whole-genome duplication events of DB and NJ occurred at 64.55 and 64.81 Mya, respectively, while the divergence between soybean and Vigna occurred in the Paleogene period. NJ genes underwent positive selection and amplification in response to temperature and abiotic stress. In species-specific gene families, NJ is mainly enriched in response to abiotic stress, while DB is primarily enriched in respiration and photosynthesis. We established the pan-genomes of four accessions (NJ, DB, IT97K-499-35 and Xiabao II) and identified 20,336 (70.5%) core genes present in all the accessions, 6,507 (55.56%) variable genes in two individuals, and 2,004 (6.95%) unique genes. The final pan genome is 616.35 Mb, and the core genome is 399.78 Mb. The variable genes are manifested mainly in stress response functions, ABC transporters, seed storage, and dormancy control. In the pan-genome sequence variation analysis, genes affected by presence/absence variants were enriched in biological processes associated with defense responses, immune system processes, signal transduction, and agronomic traits. The results of the present study provide genetic data that could facilitate efficient asparagus bean genetic improvement, especially in producing cold-adapted asparagus bean.

Chromosome-Level Genome Assembly and Multi-Omics Dataset Provide Insights into Isoflavone and Puerarin Biosynthesis in Pueraria lobata (Wild.) Ohwi.

Pueraria lobata (wild.) Ohwi is a leguminous plant and one of the traditional Chinese herbal medicines. Its puerarin extract is widely used in the pharmaceutical industry. This study reported a chromosome-level genome assembly for P. lobata and its characteristics. The genome size was ~939.2 Mb, with a contig N50 of 29.51 Mbp. Approximately 97.82% of the assembled sequences were represented by 11 pseudochromosomes. We identified that the repetitive sequences accounted for 63.50% of the P. lobata genome. A total of 33,171 coding genes were predicted, of which 97.34% could predict the function. Compared with other species, P. lobata had 757 species-specific gene families, including 1874 genes. The genome evolution analysis revealed that P. lobata was most closely related to Glycine max and underwent two whole-genome duplication (WGD) events. One was in a gamma event shared by the core dicotyledons at around 65 million years ago, and another was in the common ancestor shared by legume species at around 25 million years ago. The collinearity analysis showed that 61.45% of the genes (54,579 gene pairs) in G. max and P. lobata had collinearity. In this study, six unique PlUGT43 homologous genes were retrieved from the genome of P. lobata, and no 2-hydroxyisoflavanone 8-C-glucoside was found in the metabolites. This also revealed that the puerarin synthesis was mainly from the glycation of daidzein. The combined transcriptome and metabolome analysis suggested that two bHLHs, six MYBs and four WRKYs were involved in the expression regulation of puerarin synthesis structural genes. The genetic information obtained in this study provided novel insights into the biological evolution of P. lobata and leguminous species, and it laid the foundation for further exploring the regulatory mechanism of puerarin synthesis.

Comparative genomic analysis of Babesia duncani responsible for human babesiosis

BackgroundHuman babesiosis, caused by parasites of the genus Babesia, is an emerging and re-emerging tick-borne disease that is mainly transmitted by tick bites and infected blood transfusion. Babesia duncani has caused majority of human babesiosis in Canada; however, limited data are available to correlate its genomic information and biological features.ResultsWe generated a B. duncani reference genome using Oxford Nanopore Technology (ONT) and Illumina sequencing technology and uncovered its biological features and phylogenetic relationship with other Apicomplexa parasites. Phylogenetic analyses revealed that B. duncani form a clade distinct from B. microti, Babesia spp. infective to bovine and ovine species, and Theileria spp. infective to bovines. We identified the largest species-specific gene family that could be applied as diagnostic markers for this pathogen. In addition, two gene families show signals of significant expansion and several genes that present signatures of positive selection in B. duncani, suggesting their possible roles in the capability of this parasite to infect humans or tick vectors.ConclusionsUsing ONT sequencing and Illumina sequencing technologies, we provide the first B. duncani reference genome and confirm that B. duncani forms a phylogenetically distinct clade from other Piroplasm parasites. Comparative genomic analyses show that two gene families are significantly expanded in B. duncani and may play important roles in host cell invasion and virulence of B. duncani. Our study provides basic information for further exploring B. duncani features, such as host-parasite and tick-parasite interactions.

Insights from C. elegans into Microsporidia Biology and Host-Pathogen Relationships.

Microsporidia are poorly understood, ubiquitous eukaryotic parasites that are completely dependent on their hosts for replication. With the discovery of microsporidia species naturally infecting the genetically tractable transparent nematode C. elegans, this host has been used to explore multiple areas of microsporidia biology. Here we review results about microsporidia infections in C. elegans, which began with the discovery of the intestinal-infecting species Nematocida parisii. Recent findings include new species identification in the Nematocida genus, with more intestinal-infecting species, and also a species with broader tissue tropism, the epidermal and muscle-infecting species Nematocida displodere. This species has a longer polar tube infection apparatus, which may enable its wider tissue range. After invasion, multiple Nematocida species appear to fuse host cells, which likely promotes their dissemination within host organs. Localized proteomics identified Nematocida proteins that have direct contact with the C. elegans intestinal cytosol and nucleus, and many of these host-exposed proteins belong to expanded, species-specific gene families. On the host side, forward genetic screens have identified regulators of the Intracellular Pathogen Response (IPR), which is a transcriptional response induced by both microsporidia and the Orsay virus, which is also a natural, obligate intracellular pathogen of the C. elegans intestine. The IPR constitutes a novel immune/stress response that promotes resistance against microsporidia, virus, and heat shock. Overall, the Nematocida/C. elegans system has provided insights about strategies for microsporidia pathogenesis, as well as innate defense pathways against these parasites.

Zanthoxylum-specific whole genome duplication and recent activity of transposable elements in the highly repetitive paleotetraploid Z. bungeanum genome

Zanthoxylum bungeanum is an important spice and medicinal plant that is unique for its accumulation of abundant secondary metabolites, which create a characteristic aroma and tingling sensation in the mouth. Owing to the high proportion of repetitive sequences, high heterozygosity, and increased chromosome number of Z. bungeanum, the assembly of its chromosomal pseudomolecules is extremely challenging. Here, we present a genome sequence for Z. bungeanum, with a dramatically expanded size of 4.23 Gb, assembled into 68 chromosomes. This genome is approximately tenfold larger than that of its close relative Citrus sinensis. After the divergence of Zanthoxylum and Citrus, the lineage-specific whole-genome duplication event η-WGD approximately 26.8 million years ago (MYA) and the recent transposable element (TE) burst ~6.41 MYA account for the substantial genome expansion in Z. bungeanum. The independent Zanthoxylum-specific WGD event was followed by numerous fusion/fission events that shaped the genomic architecture. Integrative genomic and transcriptomic analyses suggested that prominent species-specific gene family expansions and changes in gene expression have shaped the biosynthesis of sanshools, terpenoids, and anthocyanins, which contribute to the special flavor and appearance of Z. bungeanum. In summary, the reference genome provides a valuable model for studying the impact of WGDs with recent TE activity on gene gain and loss and genome reconstruction and provides resources to accelerate Zanthoxylum improvement.

Insights into the Host Specificity of a New Oomycete Root Pathogen, Pythium brassicum P1: Whole Genome Sequencing and Comparative Analysis Reveals Contracted Regulation of Metabolism, Protein Families, and Distinct Pathogenicity Repertoire.

Pythium brassicum P1 Stanghellini, Mohammadi, Förster, and Adaskaveg is an oomycete root pathogen that has recently been characterized. It only attacks plant species belonging to Brassicaceae family, causing root necrosis, stunting, and yield loss. Since P. brassicum P1 is limited in its host range, this prompted us to sequence its whole genome and compare it to those of broad host range Pythium spp. such as P. aphanidermatum and P. ultimum var. ultimum. A genomic DNA library was constructed with a total of 374 million reads. The sequencing data were assembled using SOAPdenovo2, yielding a total genome size of 50.3 Mb contained in 5434 scaffolds, N50 of 30.2 Kb, 61.2% G+C content, and 13,232 putative protein-coding genes. Pythium brassicum P1 had 175 species-specific gene families, which is slightly below the normal average. Like P. ultimum, P. brassicum P1 genome did not encode any classical RxLR effectors or cutinases, suggesting a significant difference in virulence mechanisms compared to other oomycetes. Pythium brassicum P1 had a much smaller proportions of the YxSL sequence motif in both secreted and non-secreted proteins, relative to other Pythium species. Similarly, P. brassicum P1 had the fewest Crinkler (CRN) effectors of all the Pythium species. There were 633 proteins predicted to be secreted in the P. brassicum P1 genome, which is, again, slightly below average among Pythium genomes. Pythium brassicum P1 had only one cadherin gene with calcium ion-binding LDRE and DxND motifs, compared to Pythium ultimum having four copies. Pythium brassicum P1 had a reduced number of proteins falling under carbohydrate binding module and hydrolytic enzymes. Pythium brassicum P1 had a reduced complement of cellulase and pectinase genes in contrast to P. ultimum and was deficient in xylan degrading enzymes. The contraction in ABC transporter families in P. brassicum P1 is suggested to be the result of a lack of diversity in nutrient uptake and therefore host range.

Genome assembly of Scorias spongiosa and comparative genomics provide insights into ecological adaptation of honeydew-dependent sooty mould fungi

Sooty moulds are fungi of economic importance and with unique lifestyle mainly growing on insect honeydew. However, the lack of genomic data hinders investigation of genetic mechanisms underlying their ecological adaptation. With long-read sequencing technology, we generated the genome of Scorias spongiosa, an extraordinary sooty mould fungus associated with honeydew of colony aphids and producing large fruiting bodies. A 24.21 Mb high-quality genome assembly with a N50 length of 3.37 Mb was obtained. The genome contained 7758 protein coding genes, 97.13% of which were homologous to known genes, and approximately 0.29 Mb repeat sequences. Comparative genomics showed S. spongiosa lost relatively more gene families and contained fewer species-specific genes and gene families, with many CAZyme families and sugar transporters reduced or absent. This study not only promotes understanding of the ecological adaptation of sooty moulds, but also provides valuable genomic data resource for future comparative genomic and genetic studies.

Additional description and genome analyses of Caenorhabditis auriculariae representing the basal lineage of genus Caenorhabditis

Caenorhabditis auriculariae, which was morphologically described in 1999, was re-isolated from a Platydema mushroom-associated beetle. Based on the re-isolated materials, some morphological characteristics were re-examined and ascribed to the species. In addition, to clarify phylogenetic relationships with other Caenorhabditis species and biological features of the nematode, the whole genome was sequenced and assembled into 109.5 Mb with 16,279 predicted protein-coding genes. Molecular phylogenetic analyses based on ribosomal RNA and 269 single-copy genes revealed the species is closely related to C. sonorae and C. monodelphis placing them at the most basal clade of the genus. C. auriculariae has morphological characteristics clearly differed from those two species and harbours a number of species-specific gene families, indicating its usefulness as a new outgroup species for Caenorhabditis evolutionary studies. A comparison of carbohydrate-active enzyme (CAZy) repertoires in genomes, which we found useful to speculate about the lifestyle of Caenorhabditis nematodes, suggested that C. auriculariae likely has a life-cycle with tight-association with insects.

Eukaryote-Conserved Methylarginine Is Absent in Diplomonads and Functionally Compensated in Giardia.

Methylation is a common posttranslational modification of arginine and lysine in eukaryotic proteins. Methylproteomes are best characterized for higher eukaryotes, where they are functionally expanded and evolved complex regulation. However, this is not the case for protist species evolved from the earliest eukaryotic lineages. Here, we integrated bioinformatic, proteomic, and drug-screening data sets to comprehensively explore the methylproteome of Giardia duodenalis—a deeply branching parasitic protist. We demonstrate that Giardia and related diplomonads lack arginine-methyltransferases and have remodeled conserved RGG/RG motifs targeted by these enzymes. We also provide experimental evidence for methylarginine absence in proteomes of Giardia but readily detect methyllysine. We bioinformatically infer 11 lysine-methyltransferases in Giardia, including highly diverged Su(var)3-9, Enhancer-of-zeste and Trithorax proteins with reduced domain architectures, and novel annotations demonstrating conserved methyllysine regulation of eukaryotic elongation factor 1 alpha. Using mass spectrometry, we identify more than 200 methyllysine sites in Giardia, including in species-specific gene families involved in cytoskeletal regulation, enriched in coiled-coil features. Finally, we use known methylation inhibitors to show that methylation plays key roles in replication and cyst formation in this parasite. This study highlights reduced methylation enzymes, sites, and functions early in eukaryote evolution, including absent methylarginine networks in the Diplomonadida. These results challenge the view that arginine methylation is eukaryote conserved and demonstrate that functional compensation of methylarginine was possible preceding expansion and diversification of these key networks in higher eukaryotes.

Draft Genome Assembly of Floccularia luteovirens, an Edible and Symbiotic Mushroom on Qinghai-Tibet Plateau.

Floccularia luteovirens, also known as “Yellow mushroom”, is an edible ectomycorrhizal fungus widely distributed in the Qinghai-Tibet Plateau alpine meadow. So far, little genomic information is known about F. luteovirens, which is not conductive to the protection and utilization of it. In this manuscript, we present a first draft genome assembly and annotation of F. luteovirens. The fruiting body of F. luteovirens was sequenced with PacBio Sequel and Illumina Hiseq 2500 system. The assembled genome size was 28.8 Mb, and comprising 183 contigs with a N50 contig size of 571 kb. A total of 8,333 protein-coding genes were predicted and 7,999 genes were further assigned to different public protein databases. Besides, 400 CAZymes were identified in F. luteovirens. Phylogenetic analysis suggested that F. luteovirens should belong to the Agaricaceae family. Time tree result showed that the speciation of F. luteovirens happened approximately 170 Million years ago. Furthermore, 357 species-specific gene families were annotated against KEGG and GO database. This genome assembly and annotation should be an essential genomic foundation for understanding the phylogenetic, metabolic and symbiotic traits of F. luteovirens.

Three New Genome Assemblies Support a Rapid Radiation in Musa acuminata (Wild Banana).

Edible bananas result from interspecific hybridization between Musa acuminata and Musa balbisiana, as well as among subspecies in M. acuminata. Four particular M. acuminata subspecies have been proposed as the main contributors of edible bananas, all of which radiated in a short period of time in southeastern Asia. Clarifying the evolution of these lineages at a whole-genome scale is therefore an important step toward understanding the domestication and diversification of this crop. This study reports the de novo genome assembly and gene annotation of a representative genotype from three different subspecies of M. acuminata. These data are combined with the previously published genome of the fourth subspecies to investigate phylogenetic relationships. Analyses of shared and unique gene families reveal that the four subspecies are quite homogenous, with a core genome representing at least 50% of all genes and very few M. acuminata species-specific gene families. Multiple alignments indicate high sequence identity between homologous single copy-genes, supporting the close relationships of these lineages. Interestingly, phylogenomic analyses demonstrate high levels of gene tree discordance, due to both incomplete lineage sorting and introgression. This pattern suggests rapid radiation within Musa acuminata subspecies that occurred after the divergence with M. balbisiana. Introgression between M. a. ssp. malaccensis and M. a. ssp. burmannica was detected across the genome, though multiple approaches to resolve the subspecies tree converged on the same topology. To support evolutionary and functional analyses, we introduce the PanMusa database, which enables researchers to exploration of individual gene families and trees.

Molecular footprints of inshore aquatic adaptation in Indo-Pacific humpback dolphin (Sousa chinensis)

The Indo-Pacific humpback dolphin, Sousa chinensis, being a member of cetaceans, had fully adapted to inshore waters. As a threatened marine mammal, little molecular information available for understanding the genetic basis of ecological adaptation. We firstly sequenced and obtained the draft genome map of S. chinensis. Phylogenetic analysis in this study, based on the single copy orthologous genes of the draft genome, is consistent with traditional phylogenetic classification. The comparative genomic analysis indicated that S. chinensis had 494 species-specific gene families, which involved immune, DNA repair and sensory systems associated with the potential adaption mechanism. We also identified the expansion and positive selection genes in S. chinensis lineage to investigate the potential adaptation mechanism. Our study provided the potential insight into the molecular bases of ecological adaptation in Indo-Pacific humpback dolphin and will be also valuable for future understanding the ecological adaptation and evolution of cetaceans at the genomic level.

Whole genome and transcriptome analysis reveal adaptive strategies and pathogenesis of Calonectria pseudoreteaudii to Eucalyptus

BackgroundLeaf blight caused by Calonectria spp. is one of the most destructive diseases to affect Eucalyptus nurseries and plantations. These pathogens mainly attack Eucalyptus, a tree with a diversity of secondary metabolites employed as defense-related phytoalexins. To unravel the fungal adaptive mechanisms to various phytoalexins, we examined the genome of C. pseudoreteaudii, which is one of the most aggressive pathogens in southeast Asia.ResultsA 63.7 Mb genome with 14,355 coding genes of C. pseudoreteaudii were assembled. Genomic comparisons identified 1785 species-specific gene families in C. pseudoreteaudii. Most of them were not annotated and those annotated genes were enriched in peptidase activity, pathogenesis, oxidoreductase activity, etc. RNA-seq showed that 4425 genes were differentially expressed on the eucalyptus(the resistant cultivar E. grandis×E.camaldulensis M1) tissue induced medium. The annotation of GO term and KEGG pathway indicated that some of the differential expression genes were involved in detoxification and transportation, such as genes encoding ABC transporters, degrading enzymes of aromatic compounds and so on.ConclusionsPotential genomic determinants of phytoalexin detoxification were identified in C. pseudoreteaudii by comparison with 13 other fungi. This pathogen seems to employ membrane transporters and degradation enzymes to detoxify Eucalyptus phytoalexins. Remarkably, the Calonectria genome possesses a surprising number of secondary metabolism backbone enzyme genes involving toxin biosynthesis. It is also especially suited for cutin and lignin degradation. This indicates that toxin and cell wall degrading enzymes may act important roles in the establishment of Calonectria leaf blight. This study provides further understanding on the mechanism of pathogenesis in Calonectria.

Genome Sequence of the Freshwater Yangtze Finless Porpoise.

The Yangtze finless porpoise (Neophocaena asiaeorientalis ssp. asiaeorientalis) is a subspecies of the narrow-ridged finless porpoise (N. asiaeorientalis). In total, 714.28 gigabases (Gb) of raw reads were generated by whole-genome sequencing of the Yangtze finless porpoise, using an Illumina HiSeq 2000 platform. After filtering the low-quality and duplicated reads, we assembled a draft genome of 2.22 Gb, with contig N50 and scaffold N50 values of 46.69 kilobases (kb) and 1.71 megabases (Mb), respectively. We identified 887.63 Mb of repetitive sequences and predicted 18,479 protein-coding genes in the assembled genome. The phylogenetic tree showed a relationship between the Yangtze finless porpoise and the Yangtze River dolphin, which diverged approximately 20.84 million years ago. In comparisons with the genomes of 10 other mammals, we detected 44 species-specific gene families, 164 expanded gene families, and 313 positively selected genes in the Yangtze finless porpoise genome. The assembled genome sequence and underlying sequence data are available at the National Center for Biotechnology Information under BioProject accession number PRJNA433603.

Insights into the molecular evolution of peptidase inhibitors in arthropods.

Peptidase inhibitors are key proteins involved in the control of peptidases. In arthropods, peptidase inhibitors modulate the activity of peptidases involved in endogenous physiological processes and peptidases of the organisms with which they interact. Exploring available arthropod genomic sequences is a powerful way to obtain the repertoire of peptidase inhibitors in every arthropod species and to understand the evolutionary mechanisms involved in the diversification of this kind of proteins. A genomic comparative analysis of peptidase inhibitors in species belonging to different arthropod taxonomic groups was performed. The results point out: i) species or clade-specific presence is shown for several families of peptidase inhibitors; ii) multidomain peptidase inhibitors are commonly found in many peptidase inhibitor families; iii) several families have a wide range of members in different arthropod species; iv) several peptidase inhibitor families show species-specific (or clade-specific) gene family expansions; v) functional divergence may be assumed for particular clades; vi) passive expansions may be used by natural selection to fix adaptations. In conclusion, conservation and divergence of duplicated genes and the potential recruitment as peptidase inhibitors of proteins from other families are the main mechanisms used by arthropods to fix diversity. This diversity would be associated to the control of target peptidases and, as consequence, to adapt to specific environments.

Analysis of gene gain and loss in the evolution of predatory bacteria

Predatory bacteria are ubiquitously distributed in nature in including in aquatic environments, sewage, intestinal tracts of animals and humans, rhizophere and, soils. However, our understanding of their evolutionary history is limited. Results of recent studies have shown that acquiring novel genes is a major force driving bacterial evolution. Therefore, to gain a better understanding of the impact of gene gain and loss in the evolution of bacterial predators, this study employed comparative genomic approaches to identify core-set gene families and species-specific gene families, and model gene gain and loss events among 11 genomes that represented diverse lineages. In total, 1977 gene families were classified. Of these 509 (pattern 11111111111) were present all of the 11 species. Among the non-core set gene families, 52 were present only in saltwater bacteria predators and had no ortholog in the other genomes. Similarly 109 and 44 were present only in the genomes of Micavibrio spp. and Bdellovibrio spp., respectively. In this study, the gain loss mapping engine GLOOME was selected to analyze and estimate the expectations and probabilities of both gain and loss events in the predatory bacteria. In total, 354 gene families were involved in significant gene gain events, and 407 gene families were classified into gene loss events with high supported value. Moreover, 18 families from the core set gene family were identified as putative genes under positive selection. The results of this study suggest that acquisition of particular genes that encode functional proteins in metabolism and cellular processes and signaling, especially ABC systems, may help bacterial predators adapt to surrounding environmental changes and present different predation strategies for survival in their habitats.

Microsynteny and phylogenetic analysis of tandemly organised miRNA families across five members of Brassicaceae reveals complex retention and loss history

Plant genomes are characterized by the presence of large miRNA gene families which are few in number. The expansion of miRNA families is thought to be driven by gene and genome duplication. Some members of these miRNA gene families are tandemly arranged and their analysis is of interest because such organisation may indicate origin through tandem duplication and also to investigate whether some such tandem clusters have similar expression patterns, and whether these are regulated through a common set of cis-regulatory elements (eg. promoters and enhancers). As a first step, we undertake a comprehensive study using micro-synteny analyses of tandemly organised miRNA families across the Brassicaceae spanning an evolutionary time scale of ca. 45 million years, among Arabidopsis, Capsella, Brassica and Thellungiella species, to address the following questions: Are most miRNA gene families present as tandem clusters? To what extent are these tandem patterns retained? To what extent can family sizes be ascribed to genome duplication? Our analysis of thirteen tandemly organised miRNA families revealed that synteny is largely conserved among Arabidopsis thaliana, A. lyrata and Capsella rubella, which form a clade spanning approximately between 6.2–9.8 my (Acarkan et al., 2000) [1]. On the other hand, comparison of sequences from these species with Brassica rapa, B. oleracea and Thellungiella halophila, which form a separate clade spanning 31 my (Franzke et al., 2011) [2] reveals many differences. The latter clade reveals several paralogous duplications that probably resulted from whole genome duplication, as well as disrupted synteny. Phylogenetic analyses of precursor sequences generally support the history inferred from synteny analysis. Synteny and phylogenetic analysis of six members of the tandemly organised miR169 family suggest that the Brassicaceae ancestral state consisted of a “dimer as a unit” which may have undergone direct local duplication to retain the transcriptional orientation followed by lineage specific changes. MiR169, to the best of our knowledge, is one of the largest tandemly organised miRNA gene family across plant kingdom and further analysis should reveal the generality of this pattern of evolution. The conserved organisation of miR395A-B-C and miR395 D-E-F as two clusters on same chromosome/scaffold across A. thaliana, B. rapa and salsuginea demonstrates retention of the large chromosomal segment across the two lineages. MiRNA family miR845 was detected only in Arabidopsis species and Thellungiella indicating a complex loss and retention history. MiR447A-B family was only found in A. thaliana indicating that it is a species-specific gene family of recent origin.

High quality reference genome of drumstick tree (Moringa oleifera Lam.), a potential perennial crop.

The drumstick tree (Moringa oleifera Lam.) is a perennial crop that has gained popularity in certain developing countries for its high-nutrition content and adaptability to arid and semi-arid environments. Here we report a high-quality draft genome sequence of M. oleifera. This assembly represents 91.78% of the estimated genome size and contains 19,465 protein-coding genes. Comparative genomic analysis between M. oleifera and related woody plant genomes helps clarify the general evolution of this species, while the identification of several species-specific gene families and positively selected genes in M. oleifera may help identify genes related to M. oleifera's high protein content, fast-growth, heat and stress tolerance. This reference genome greatly extends the basic research on M. oleifera, and may further promote applying genomics to enhanced breeding and improvement of M. oleifera.