Natural variation, functional divergence, and local adaptation of nucleotide binding site sequences in Rhododendron (Ericaceae)

To further understand natural variation and local adaptation in the evolution of plant defense, we analyzed polymorphism data of nucleotide-binding site (NBS) sequences of Rhododendron at both the species and population levels. Multiple duplication events were found in NBS sequence evolution in Rhododendron genomes, which resulted in six clades: A–F. Our results of several NBS clade pair comparisons showed significant evolutionary rate changes based on differences in substitution rates between NBS-encoding protein clades (type I functional divergence). Pairwise comparisons of NBS clades further revealed that many amino acids displayed radical biochemical property changes causing a shift in amino acid preferences between NBS-encoding protein clades (type II functional divergence). Such divergent evolution of NBSs is likely a consequence of positive selection related to differentiation of recognition signals in response to different pathogens. Primers specific to clades B and C, which differed in the number of radical amino acid changes causing type II functional divergence and levels of nucleotide diversities, were further used to amplify population clades B and C NBS sequences of Rhododendron formosanum populations. Higher levels of net nucleotide divergences (measured by D _a) between R. formosanum populations were found based on NBS sequences of population clade B compared to population clade C, suggesting local adaptation of population clade B NBS sequences. Local adaptation can be further inferred for R. formosanum population clade B NBS sequences because of significant Φ _ST based on variation in nonsynonymous substitutions. Furthermore, local adaptation was also suggested by no significant correlation of population pairwise F _ST between population clades B and C in R. formosanum.     

Glucosinolate biosynthetic genes in Brassica rapa

Glucosinolates (GS) are a group of amino acid-derived secondary metabolites found throughout the Cruciferae family. Glucosinolates and their degradation products play important roles in pathogen and insect interactions, as well as in human health. In order to elucidate the glucosinolate biosynthetic pathway in Brassica rapa, we conducted comparative genomic analyses of Arabidopsis thaliana and B. rapa on a genome-wide level. We identified 102 putative genes in B. rapa as the orthologs of 52 GS genes in A. thaliana. All but one gene was successfully mapped on 10 chromosomes. Most GS genes exist in more than one copy in B. rapa. A high co-linearity in the glucosinolate biosynthetic pathway between A. thaliana and B. rapa was also established. The homologous GS genes in B. rapa and A. thaliana share 59-91% nucleotide sequence identity and 93% of the GS genes exhibit synteny between B. rapa and A. thaliana. Moreover, the structure and arrangement of the B. rapa GS (BrGS) genes correspond with the known evolutionary divergence of B. rapa, and may help explain the profiles and accumulation of GS in B. rapa.     

Genome-wide analysis and stress-responsive expression of CCCH zinc finger family genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Background

Ubiquitous CCCH nucleic acid-binding motif is found in a wide-variety of organisms. CCCH genes are involved in plant developmental processes and biotic and abiotic stress responses. Brassica rapa is a vital economic crop and classical model plant of polyploidy evolution, but the functions of CCCH genes in B. rapa are unclear.

Results

In this study, 103 CCCH genes in B. rapa were identified. A comparative analysis of the chromosomal position, gene structure, domain organization and duplication event between B. rapa and Arabidopsis thaliana were performed. Results showed that CCCH genes could be divided into 18 subfamilies, and segmental duplication might mainly contribute to this family expansion. C-X_7/8-C-X₅-C₃-H was the most commonly found motif, but some novel CCCH motifs were also found, along with some loses of typical CCCH motifs widespread in other plant species. The multifarious gene structures and domain organizations implicated functional diversity of CCCH genes in B. rapa. Evidence also suggested functional redundancy in at least one subfamily due to high conservation between members. Finally, the expression profiles of subfamily-IX genes indicated that they are likely involved in various stress responses.

Conclusion

This study provides the first genome-wide characterization of the CCCH genes in B. rapa. The results suggest that B. rapa CCCH genes are likely functionally divergent, but mostly involved in plant development and stress response. These results are expected to facilitate future functional characterization of this potential RNA-binding protein family in Brassica crops.

     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

Spatial evolutionary and ecological vicariance analysis (SEEVA), a novel approach to biogeography and speciation research,with an example from Brazilian Gentianaceae

Aim Spatial evolutionary and ecological vicariance analysis (SEEVA) is a simple analytical method that evaluates environmental or ecological divergence associated with evolutionary splits. It integrates evolutionary hypotheses, phylogenetic data, and spatial, temporal, environmental and geographical information to elucidate patterns. Using a phylogeny of Prepusa Mart. and Senaea Taub. (Angiospermae: Gentianaceae), SEEVA is used to describe the radiation and ecological patterns of this basal gentian group across south‐eastern Brazil. Location Latin America, global. Methods Environmental data for 151 geolocated botanical collections, associated with specimens from seven species, were compiled with Arc GIS, and were matched with geolocated base layers of eight climatological variables, as well as one each of geological, soil type, elevational and vegetation variables. Sister groups were defined on the basis of the six nested nodes that defined the phylogenetic tree of these two genera. A (0, 1)‐scaled divergence index (D) was defined and tested for each of 12 environmental and for each of the six phylogenetic nodes, by means of contingency analyses. We contrast divergence indices of nested clades, allopatric and sympatric sister clades. Results The level of ecological divergence between sister clades/species, defined in terms of D measures, was substantial for five of six nodes, with 21 of 72 environmental comparisons having D > 0.75. Soil types and geological age of bedrock were strongly divergent only for basal nodes in the phylogeny, by contrast with temperature and precipitation, which exhibited strong divergence at all nodes. There has been strong divergence and progressive occupation of wetter and colder habitats throughout the history of Prepusa. Nodes separating allopatric sister clades exhibited larger niche divergence than did those separating sympatric sister clades. Main conclusions SEEVA provides a multi‐source, direct analysis method for correlating field collections, phylogenetic hypotheses, species distributions and georeferenced environmental data. Using SEEVA, it was possible to quantify and test the divergence between sister lineages, illustrating both niche conservatism and ecological specialization. SEEVA permits elucidation of historical and ecological vicariance for evolutionary lineages, and is amenable to wide application, taxonomically, geographically and ecologically.     

Comparative transcriptional profiling and evolutionary analysis of the GRAM domain family in eukaryotes

The GRAM domain was found in glucosyltransferases, myotubularins and other membrane-associated proteins. So far, functions for majority of these proteins are yet to be uncovered. In order to address the evolutionary and functional significance of this family members, we have performed a comprehensive investigation on their genome-wide identification, phylogenetic relationship and expression divergence in five different organisms representing monocot/dicot plants, vertebrate/invertebrate animals and yeast, namely, Oryza sativa, Arabidopsis thaliana, Mus musculus, Drosophila melanogaster and Saccharomyces cerevisiae, respectively. We have identified 65 members of GRAM domain family from these organisms. Our data revealed that this family was an ancient group and various organisms had evolved into different family sizes. Large-scale genome duplication and divergence in both expression patterns and functions were significantly contributed to the expansion and retention of this family. Mouse and Drosophila members showed higher divergences in their proteins as indicated by higher Ka/Ks ratios and possessed multiple domains in various combinations. However, in plants, their protein functions were possibly retained with a relatively low divergence as signified by lower Ka/Ks ratios and only one additional domain was combined during evolution. On the other hand, this family in all five organisms exhibited high divergence in their expression patterns both at tissue level and under various biotic and abiotic stresses. These highly divergent expression patterns unraveled the complexity of functions of GRAM domain family. Each member may play specialized roles in a specific tissue or stress condition and may function as regulators of environmental and hormonal signaling.     

Extensive Functional Diversification of the Populus Glutathione S-Transferase Supergene Family

Identifying how genes and their functions evolve after duplication is central to understanding gene family radiation. In this study, we systematically examined the functional diversification of the glutathione S-transferase (GST) gene family in Populus trichocarpa by integrating phylogeny, expression, substrate specificity, and enzyme kinetic data. GSTs are ubiquitous proteins in plants that play important roles in stress tolerance and detoxification metabolism. Genome annotation identified 81 GST genes in Populus that were divided into eight classes with distinct divergence in their evolutionary rate, gene structure, expression responses to abiotic stressors, and enzymatic properties of encoded proteins. In addition, when all the functional parameters were examined, clear divergence was observed within tandem clusters and between paralogous gene pairs, suggesting that subfunctionalization has taken place among duplicate genes. The two domains of GST proteins appear to have evolved under differential selective pressures. The C-terminal domain seems to have been subject to more relaxed functional constraints or divergent directional selection, which may have allowed rapid changes in substrate specificity, affinity, and activity, while maintaining the primary function of the enzyme. Our findings shed light on mechanisms that facilitate the retention of duplicate genes, which can result in a large gene family with a broad substrate spectrum and a wide range of reactivity toward different substrates.     

Phylogenetic and Structural Relationships of the PR5 Gene Family Reveal an Ancient Multigene Family Conserved in Plants and Select Animal Taxa

Pathogenesis-related group 5 (PR5) plant proteins include thaumatin, osmotin, and related proteins, many of which have antimicrobial activity. The recent discovery of PR5-like (PR5-L) sequences in nematodes and insects raises questions about their evolutionary relationships. Using complete plant genome data and discovery of multiple insect PR5-L sequences, phylogenetic comparisons among plants and animals were performed. All PR5/PR5-L protein sequences were mined from genome data of a member of each of two main angiosperm groups—the eudicots (Arabidoposis thaliana) and the monocots (Oryza sativa)—and from the Caenorhabditis nematode (C. elegans and C. briggsase). Insect PR5-L sequences were mined from EST databases and GenBank submissions from four insect orders: Coleoptera (Diaprepes abbreviatus and Biphyllus lunatus), Orthoptera (Schistocerca gregaria), Hymenoptera (Lysiphlebus testaceipes), and Hemiptera (Toxoptera citricida). Parsimony and Bayesian phylogenetic analyses showed that the PR5 family is paraphyletic in plants, likely arising from 10 genes in a common ancestor to monocots and eudicots. After evolutionary divergence of monocots and eudicots, PR5 genes increased asymmetrically among the 10 clades. Insects and nematodes contain multiple sequences (seven PR5-Ls in nematodes and at least three in some insects) all related to the same plant clade, with nematode and insect sequences separating as two clades. Protein structural homology modeling showed strong similarity among animal and plant PR5/PR5-Ls, with divergence only in surface-exposed loops. Sequence and structural conservation among PR5/PR5-Ls suggests an important and conserved role throughout the evolutionary divergence of the diverse organisms from which they reside. [Reviewing Editor: Dr. Rafael Zardoya]     

Genome-wide identification,classification, and expression analysis of the phytocyanin gene family in <Emphasis Type="Italic">Phalaenopsis equestris</Emphasis>

Phytocyanins (PCs) are ancient blue copper-binding proteins in plants that bind to single type I copper atoms and function as electron transporters. PCs play an important role in plant development and stress resistance. Many PCs are considered to be chimeric arabinogalactan proteins (AGPs). Previously, 38, 62, and 84 PC genes were identified in Arabidopsis thaliana, Oryza sativa, and Brassica rapa, respectively. In this study, we identified 30 putative PC genes in the orchid Phalaenopsis equestris through comprehensive bioinformatics analysis. Based on phylogeny and motif constitution, the P. equestris phytocyanins (PePCs) were divided into five subclasses: 10 early nodulin-like proteins, 10 uclacyanin-like proteins, five stellacyanin-like proteins, four plantacyanin-like proteins, and one unknown protein. Structural and glycosylation predictions suggested that 16 PePCs were glycosylphosphatidylinositol-anchored proteins localized to the plasma membrane, 22 PePCs contain N-glycosylation sites, and 14 are chimeric AGPs. Phylogenetic analysis indicated that each subfamily was derived from a common ancestor before the divergence of monocot and dicot lineages and that the expansion of the PC subfamilies occurred after the divergence of orchids and Arabidopsis. The number of exons in PC genes was conserved. Expression analysis in four tissues revealed that nine PC genes were highly expressed in flowers, stems, and roots, suggesting that these genes play important roles in growth and development in P. equestris. The results of this study lay the foundation for further analysis of the functions of this gene family in plants.     

Cytochrome oxidase I sequences reveal possible cryptic diversity in the cosmopolitan symbiotic copepod Nesippus orientalis Heller, 1868 (Pandaridae: Siphonostomatoida) on elasmobranch hosts from the KwaZulu-Natal coast of South Africa

Over the past decade, numerous molecular phylogenetic studies uncovered cryptic diversity within the Copepoda, yet very few investigations focused on symbiotic copepods. Here we report mitochondrial DNA cytochrome oxidase I diversity in the cosmopolitan elasmobranch symbiont Nesippus orientalis off the KwaZulu-Natal coast of South Africa. Analysis of partial COI sequences of copepods sampled from a diversity of shark hosts, revealed the presence of two divergent clades. Diversity within the clades does not appear to be structured based on host species, host individual, geographic locality or time of sampling. However, divergence between the two clades seems to be related to host species. Phylogenetic analyses of representatives from the two clades, along with Nesippus spp., Caligus spp. and Lepeophtheirus spp. outgroups, further supports the distinction between the two clades. Future molecular phylogenetic investigations of widespread copepod symbionts most likely will reveal far greater levels of biodiversity than currently recognized.     

Genome-wide acceleration of protein evolution in flies (Diptera)

Background

The rate of molecular evolution varies widely between proteins, both within and among lineages. To what extent is this variation influenced by genome-wide, lineage-specific effects? To answer this question, we assess the rate variation between insect lineages for a large number of orthologous genes.

Results

When compared to the beetle Tribolium castaneum, we find that the stem lineage of flies and mosquitoes (Diptera) has experienced on average a 3-fold increase in the rate of evolution. Pairwise gene comparisons between Drosophila and Tribolium show a high correlation between evolutionary rates of orthologous proteins.

Conclusion

Gene specific divergence rates remain roughly constant over long evolutionary times, modulated by genome-wide, lineage-specific effects. Among the insects analysed so far, it appears that the Tribolium genes show the lowest rates of divergence. This has the practical consequence that homology searches for human genes yield significantly better matches in Tribolium than in Drosophila. We therefore suggest that Tribolium is better suited for comparisons between phyla than the widely employed dipterans.     

Loss/retention and evolution of NBS-encoding genes upon whole genome triplication of Brassica rapa

A genome triplication took place in the ancestor of Brassiceae species after the split of the Arabidopsis lineage. The postfragmentation and shuffling of the genome turned the ancestral hexaploid back to diploids and caused the radiation of Brassiceae species. The course of speciation was accompanied by the loss of duplicate genes and also influenced the evolution of retained genes. Of all the genes, those encoding NBS domains are typical R genes that confer resistance to invading pathogens. In this study, using the genome of Arabidopsis thaliana as a reference, we examined the loss/retention of orthologous NBS-encoding loci in the tripled Brassica rapa genome and discovered differential loss/retention frequencies. Further analysis indicated that loci of different retention ratios showed different evolutionary patterns. The loci of classesII and III (maintaining two and three syntenic loci, respectively, multi-loci) show sharper expansions by tandem duplications, have faster evolutionary rates and have more potential to be associated with novel gene functions. On the other hand, the loci that are retained at the minimal rate (keeping only one locus, class I, single locus) showed opposite patterns. Phylogenetic analysis indicated that recombination and translocation events were common among multi-loci in B. rapa, and differential evolutionary patterns between multi- and single-loci are likely the consequence of recombination. Investigations towards other gene families demonstrated different evolutionary characteristics between different gene families. The evolution of genes is more likely determined by the property of each gene family, and the whole genome triplication provided only a specific condition.     

Allopatric divergence and secondary contact without genetic admixture for Arichanna perimelaina (Lepidoptera: Geometridae), an alpine moth endemic to the Hengduan Mountains

Mountain systems, especially at high altitudes, are an excellent model for determining the mechanisms underlying high species diversity and endemism. Herein, we elucidate the evolutionary history of the alpine moth Arichanna perimelaina (Wehrli), which is endemic to the Hengduan Mountains (HM) region in southwest China, based on three mitochondrial genes and two nuclear genes. Our results revealed six deeply divergent clades that corresponded to populations in different mountain systems in the HM region. Bayesian divergence time estimations suggested a mid‐ to late Pleistocene genetic divergence. The results also showed that the Mt Yulong (YL) region was a refugium and valley corridors established by glaciation during the Pleistocene allowed populations on the separate mountains to migrate. The reproductive isolation among the different clades on contact zone in the YL region may be associated with the asynchronous mating rhythms and/or the divergent mate recognition caused by the ecological source of divergent selection. Allopatric divergence associated with complex topographies and climatic oscillations, regional dispersal via valley corridors and the suitable refugium of the YL region shaped the genetic divergence and distribution pattern of A. perimelaina in the HM region. These findings highlight the essential role of complex terrain and climatic fluctuations in shaping the unique phylogeographic history of a narrow alpine moth, and provide insights into the mechanisms underlying high species richness and endemism in the HM region.     

Environmental versus Anthropogenic Effects on Population Adaptive Divergence in the Freshwater Snail Lymnaea stagnalis

Repeated pesticide contaminations of lentic freshwater systems located within agricultural landscapes may affect population evolution in non-target organisms, especially in species with a fully aquatic life cycle and low dispersal ability. The issue of evolutionary impact of pollutants is therefore conceptually important for ecotoxicologists. The impact of historical exposure to pesticides on genetic divergence was investigated in the freshwater gastropod Lymnaea stagnalis, using a set of 14 populations from contrasted environments in terms of pesticide and other anthropogenic pressures. The hypothesis of population adaptive divergence was tested on 11 life-history traits, using Q _ST -F _ST comparisons. Despite strong neutral differentiation (mean F _ST = 0.291), five adult traits or parameters were found to be under divergent selection. Conversely, two early expressed traits showed a pattern consistent with uniform selection or trait canalization, and four adult traits appeared to evolve neutrally. Divergent selection patterns were mostly consistent with a habitat effect, opposing pond to ditch and channel populations. Comparatively, pesticide and other human pressures had little correspondence with evolutionary patterns, despite hatching rate impairment associated with global anthropogenic pressure. Globally, analyses revealed high genetic variation both at neutral markers and fitness-related traits in a species used as model in ecotoxicology, providing empirical support for the need to account for genetic and evolutionary components of population response in ecological risk assessment.     

Evolutionary history of mitogen-activated protein kinase (MAPK) genes in Lotus,Medicago, and Phaseolus

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that mediate various signaling pathways associated with biotic and abiotic stress responses in eukaryotes. The MAPK genes form a 3-tier signal transduction cascade between cellular stimuli and physiological responses. Recent identification of soybean MAPKs and availability of genome sequences from other legume species allowed us to identify their MAPK genes. The main objectives of this study were to identify MAPKs in 3 legume species, Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, and to assess their phylogenetic relationships. We used approaches in comparative genomics for MAPK gene identification and named the newly identified genes following Arabidopsis MAPK nomenclature model. We identified 19, 18, and 15 MAPKs and 7, 4, and 9 MAPKKs in the genome of Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, respectively. Within clade placement of MAPKs and MAPKKs in the 3 legume species were consistent with those in soybean and Arabidopsis. Among 5 clades of MAPKs, 4 founder clades were consistent to MAPKs of other plant species and orthologs of MAPK genes in the fifth clade-"Clade E" were consistent with those in soybean. Our results also indicated that some gene duplication events might have occurred prior to eudicot-monocot divergence. Highly diversified MAPKs in soybean relative to those in 3 other legume species are attributable to the polyploidization events in soybean. The identification of the MAPK genes in the legume species is important for the legume crop improvement; and evolutionary relationships and functional divergence of these gene members provide insights into plant genome evolution.     

Preferential retention,expression profile and potential functional diversity analysis of HD-Zip gene family in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Homeodomain-Leu zipper (HD-Zip) gene family performs important biological functions related to organ development, photomorphogenesis and abiotic stress response in higher plants. However, systematic analysis of HD-Zip genes in Brassica rapa has not been performed. In the present study, a bioinformatics approach was used to identify and characterize the BraHD-Zip gene family in B. rapa. A total of 88 members were identified. All putative BraHD-Zip proteins contained a clear HD and LZ combined domain. Eighty-seven BraHD-Zips were non-randomly located on ten chromosomes. This gene family was mainly expanded following the whole genome triplication event and was preferentially over-retained relative to its neighboring genes in B. rapa. On phylogenetic analysis, the BraHD-Zips could be categorized into four distinct major groups (I–IV). Each group exhibited variant gene structures and motif distributions. Some syntenic orthologous gene pairs presented diverse expression profiles, which indicate that these gene pairs may be involved in the development of new functions during evolution. In summary, our analysis provided genome-wide insights into the expansion, preferential retention, expression profiles and functional diversity of BraHD-Zip genes following whole genome triplication in B. rapa.     

Genome-Wide Survey and Expression Analysis of the Putative Non-Specific Lipid Transfer Proteins in Brassica rapa L

Background

Plant non-specific lipid transfer proteins (nsLtps) are small, basic proteins encoded by multigene families and have reported functions in many physiological processes such as mediating phospholipid transfer, defense reactions against phytopathogens, the adaptation of plants to various environmental conditions, and sexual reproduction. To date, no genome-wide overview of the Brassica rapa nsLtp (BrnsLtp) gene family has been performed. Therefore, as the first step and as a helpful strategy to elucidate the functions of BrnsLtps, a genome-wide study for this gene family is necessary.

Methodology/Principal Finding

In this study, a total of 63 putative BrnsLtp genes were identified through a comprehensive in silico analysis of the whole genome of B. rapa. Based on the sequence similarities, these BrnsLtps was grouped into nine types (I, II, III, IV, V, VI, VIII, IX, and XI). There is no type VII nsLtps in B. rapa, and a new type, XI nsLtps, was identified in B. rapa. Furthermore, nine type II AtLtps have no homologous genes in B. rapa. Gene duplication analysis demonstrated that the conserved collinear block of each BrnsLtp is highly identical to those in Arabidopsis and that both segmental duplications and tandem duplications seem to play equal roles in the diversification of this gene family. Expression analysis indicated that 29 out of the 63 BrnsLtps showed specific expression patterns. After careful comparison and analysis, we hypothesize that some of the type I BrnsLtps may function like Arabidopsis pathogenesis-related-14 (PR-14) proteins to protect the plant from phytopathogen attack. Eleven BrnsLtps with inflorescence-specific expression may play important roles in sexual reproduction. Additionally, BrnsLtpI.3 may have functions similar to Arabidopsis LTP1.

Conclusions/Significance

The genome-wide identification, bioinformatic analysis and expression analysis of BrnsLtp genes should facilitate research of this gene family and polyploidy evolution and provide new insight towards elucidating their biological functions in plants.     

Clade-specific chromosomal rearrangements and loss of subtelomeric adhesins in Candida auris

Candida auris is an emerging fungal pathogen of rising concern due to global spread, the ability to cause healthcare-associated outbreaks, and antifungal resistance. Genomic analyses revealed that early contemporaneously detected cases of C. auris were geographically stratified into four major clades. While Clades I, III, and IV are responsible for ongoing outbreaks of invasive and multidrug-resistant infections, Clade II, also termed the East Asian clade, consists primarily of cases of ear infection, is often susceptible to all antifungal drugs, and has not been associated with outbreaks. Here, we generate chromosome-level assemblies of twelve isolates representing the phylogenetic breadth of these four clades and the only isolate described to date from Clade V. This Clade V genome is highly syntenic with those of Clades I, III, and IV, although the sequence is highly divergent from the other clades. Clade II genomes appear highly rearranged, with translocations occurring near GC-poor regions, and large subtelomeric deletions in most chromosomes, resulting in a substantially different karyotype. Rearrangements and deletion lengths vary across Clade II isolates, including two from a single patient, supporting ongoing genome instability. Deleted subtelomeric regions are enriched in Hyr/Iff-like cell-surface proteins, novel candidate cell wall proteins, and an ALS-like adhesin. Cell wall proteins from these families and other drug-related genes show clade-specific signatures of selection in Clades I, III, and IV. Subtelomeric dynamics and the conservation of cell surface proteins in the clades responsible for global outbreaks causing invasive infections suggest an explanation for the different phenotypes observed between clades.