Chromosomal evolution, environmental heterogeneity, and migration drive spatial patterns of species richness in Calochortus (Liliaceae).

We used nuclear genomic data and statistical models to evaluate the ecological and evolutionary processes shaping spatial variation in species richness in Calochortus (Liliaceae, 74 spp.). Calochortus occupies diverse habitats in the western United States and Mexico and has a center of diversity in the California Floristic Province, marked by multiple orogenies, winter rainfall, and highly divergent climates and substrates (including serpentine). We used sequences of 294 low-copy nuclear loci to produce a time-calibrated phylogeny, estimate historical biogeography, and test hypotheses regarding drivers of present-day spatial patterns in species number. Speciation and species coexistence require reproductive isolation and ecological divergence, so we examined the roles of chromosome number, environmental heterogeneity, and migration in shaping local species richness. Six major clades-inhabiting different geographic/climatic areas, and often marked by different base chromosome numbers (n = 6 to 10)-began diverging from each other ~10.3 Mya. As predicted, local species number increased significantly with local heterogeneity in chromosome number, elevation, soil characteristics, and serpentine presence. Species richness is greatest in the Transverse/Peninsular Ranges where clades with different chromosome numbers overlap, topographic complexity provides diverse conditions over short distances, and several physiographic provinces meet allowing immigration by several clades. Recently diverged sister-species pairs generally have peri-patric distributions, and maximum geographic overlap between species increases over the first million years since divergence, suggesting that chromosomal evolution, genetic divergence leading to gametic isolation or hybrid inviability/sterility, and/or ecological divergence over small spatial scales may permit species co-occurrence.

Global analysis of Poales diversification - parallel evolution in space and time into open and closed habitats.

Poales are one of the most species-rich, ecologically and economically important orders of plants and often characterise open habitats, enabled by unique suites of traits. We test six hypotheses regarding the evolution and assembly of Poales in open and closed habitats throughout the world, and examine whether diversification patterns demonstrate parallel evolution. We sampled 42% of Poales species and obtained taxonomic and biogeographic data from the World Checklist of Vascular Plants database, which was combined with open/closed habitat data scored by taxonomic experts. A dated supertree of Poales was constructed. We integrated spatial phylogenetics with regionalisation analyses, historical biogeography and ancestral state estimations. Diversification in Poales and assembly of open and closed habitats result from dynamic evolutionary processes that vary across lineages, time and space, most prominently in tropical and southern latitudes. Our results reveal parallel and recurrent patterns of habitat and trait transitions in the species-rich families Poaceae and Cyperaceae. Smaller families display unique and often divergent evolutionary trajectories. The Poales have achieved global dominance via parallel evolution in open habitats, with notable, spatially and phylogenetically restricted divergences into strictly closed habitats.

Respirable Metals, Bacteria, and Fungi during a Saharan-Sahelian Dust Event in Houston, Texas.

Although airborne bacteria and fungi can impact human, animal, plant, and ecosystem health, very few studies have investigated the possible impact of their long-range transport in the context of more commonly measured aerosol species, especially those present in an urban environment. We report first-of-kind simultaneous measurements of the elemental and microbial composition of North American respirable airborne particulate matter concurrent with a Saharan-Sahelian dust episode. Comprehensive taxonomic and phylogenetic profiles of microbial communities obtained by 16S/18S/ITS rDNA sequencing identified hundreds of bacteria and fungi, including several cataloged in the World Health Organization's lists of global priority human pathogens along with numerous other animal and plant pathogens and (poly)extremophiles. While elemental analysis sensitively tracked long-range transported Saharan dust and its mixing with locally emitted aerosols, microbial diversity, phylogeny, composition, and abundance did not well correlate with the apportioned African dust mass. Bacterial/fungal diversity, phylogenetic signal, and community turnover were strongly correlated to apportioned sources (especially vehicular emissions and construction activities) and elemental composition (especially calcium). Bacterial communities were substantially more dissimilar from each other across sampling days than were fungal communities. Generalized dissimilarity modeling revealed that daily compositional turnover in both communities was linked to calcium concentrations and aerosols from local vehicles and Saharan dust. Because African dust is known to impact large areas in northern South America, the Caribbean Basin, and the southern United States, the microbiological impacts of such long-range transport should be assessed in these regions.

Adaptive associations among life history, reproductive traits, environment, and origin in the Wisconsin angiosperm flora.

PREMISE: We tested 25 classic and novel hypotheses regarding trait-origin, trait-trait, and trait-environment relationships to account for flora-wide variation in life history, habit, and especially reproductive traits using a plastid DNA phylogeny of most native (96.6%, or 1494/1547 species) and introduced (87.5%, or 690/789 species) angiosperms in Wisconsin, USA. METHODS: We assembled data on life history, habit, flowering, dispersal, mating system, and occurrence across open/closed/mixed habitats across species in the state phylogeny. We used phylogenetically structured analyses to assess the strength and statistical significance of associations predicted by our models. RESULTS: Introduced species are more likely to be annual herbs, occupy open habitats, have large, visually conspicuous, hermaphroditic flowers, and bear passively dispersed seeds. Among native species, hermaphroditism is associated with larger, more conspicuous flowers; monoecy is associated with small, inconspicuous flowers and passive seed dispersal; and dioecy is associated with small, inconspicuous flowers and fleshy fruits. Larger flowers with more conspicuous colors are more common in open habitats, and in understory species flowering under open (spring) canopies; fleshy fruits are more common in closed habitats. Wind pollination may help favor dioecy in open habitats. CONCLUSIONS: These findings support predictions regarding how breeding systems depend on flower size, flower color, and fruit type, and how those traits depend on habitat. This study is the first to combine flora-wide phylogenies with complete trait databases and phylogenetically structured analyses to provide powerful tests of evolutionary hypotheses about reproductive traits and their variation with geographic source, each other, and environmental conditions.

Phylogeography, population genetics and distribution modelling reveal vulnerability of Scirpus longii (Cyperaceae) and the Atlantic Coastal Plain Flora to climate change.

A proactive approach to conservation must be predictive, anticipating how habitats will change and which species are likely to decline or prosper. We use composite species distribution modelling to identify suitable habitats for 18 members of the North American Atlantic Coastal Plain Flora (ACPF) since the Last Glacial Maximum and project these into the future. We then use Scirpus longii (Cyperaceae), a globally imperiled ACPF sedge with many of the characteristics of extinction vulnerability, as a case study. We integrate phylogeographical and population genetic analyses and species distribution modelling to develop a broad view of its current condition and prognosis for conservation. We use genotyping-by-sequencing to characterize the genomes of 142 S. longii individuals from 20 populations distributed throughout its range (New Jersey to Nova Scotia). We measure the distribution of genetic diversity in the species and reconstruct its phylogeographical history using the snapp and rase models. Extant populations of S. longii originated from a single refugium south of the Laurentide ice sheet around 25 ka. The genetic diversity of S. longii is exceedingly low, populations exhibit little genetic structure and the species is slightly inbred. Projected climate scenarios indicate that nearly half of extant populations of S. longii will be exposed to unsuitable climate by 2070. Similar changes in suitable habitat will occur for many other northern ACPF species-centres of diversity will shift northward and Nova Scotia may become the last refuges for those species not extinguished.

Inferring hypothesis-based transitions in clade-specific models of chromosome number evolution in sedges (Cyperaceae).

Large-scale changes in chromosome number have been associated with diversification rate shifts in many lineages of plants. For instance, several ancient rounds of polyploidization events have been inferred to promote genomic differentiation and/or isolation and, consequently, angiosperm diversification. Dysploidy, although less studied, has been suggested to also play an important role in angiosperm diversification. In this article, we aim to elucidate the role of chromosomal rearrangements on lineage diversification by analyzing a new comprehensive sedge (Cyperaceae) phylogenetic tree. Our null hypothesis is that the mode and tempo of chromosome evolution are to be homogeneous across the complete phylogeny. In order to discern patterns of diversification shifts and chromosome number changes within the family tree, we tested clade-specific chromosome evolution models for several subtrees according to previously reported increments of diversification rates. Results show that a complex, heterogeneous model composed of different clade-specific chromosome evolution transitions are significantly supported against the null hypothesis of a model with no chromosome number model transition events along the phylogeny. This could suggest a link between diversification and changes in chromosome number evolution although other possibilities are not discarded. Our methodological approach may allow identifying different patterns of chromosome evolution, as found for Cyperaceae, for other lineages at different evolutionary levels.

Role of cytochrome P450 2B sequence variation and gene copy number in facilitating dietary specialization in mammalian herbivores.

Theory postulates that dietary specialization in mammalian herbivores is enabled by a specialized set of liver enzymes that process the high concentrations of similar plant secondary metabolites (PSMs) in the diets of specialists. To investigate whether qualitative and quantitative differences in detoxification mechanisms distinguish dietary specialists from generalists, we compared the sequence diversity and gene copy number of detoxification enzymes in two woodrat species: a generalist, the white-throated woodrat (Neotoma albigula) and a juniper specialist, Stephens' woodrat (N. stephensi). We focused on enzymes in the cytochrome P450 subfamily 2B (CYP2B), because previous research suggests this subfamily plays a key role in the processing of PSMs. For both woodrat species, we obtained and sequenced CYP2B cDNA, generated CYP2B phylogenies, estimated CYP2B gene copy number and created a homology model of the active site. We found that the specialist possessed on average ~5 more CYP2B gene copies than the generalist, but the specialist's CYP2B sequences were less diverse. Phylogenetic analysis of putative CYP2B homologs resolved woodrat species as reciprocally monophyletic and suggested evolutionary convergence of distinct homologs on similar key amino acid residues in both species. Homology modelling of the CYP2B enzyme suggests that interspecific differences in substrate preference and function likely result from amino acid differences in the enzyme active site. The characteristics of CYP2B in the specialist, that is greater gene copy number coupled with less sequence variation, are consistent with specialization to a narrow range of dietary toxins.

Comparative transcriptomics and genomic patterns of discordance in Capsiceae (Solanaceae).

The integration of genomics and phylogenetics allows new insight into the structure of gene tree discordance, the relationships among gene position, gene history, and rate of evolution, as well as the correspondence of gene function, positive selection, and gene ontology enrichment across lineages. We explore these issues using the tribe Capsiceae (Solanaceae), which is comprised of the genera Lycianthes and Capsicum (peppers). In combining the annotated genomes of Capsicum with newly sequenced transcriptomes of four species of Lycianthes and Capsicum, we develop phylogenies for 6747 genes, and construct a backbone species tree using both concordance and explicit phylogenetic network approaches. We quantify phylogenetic discordance among individual gene trees, measure their rates of synonymous and nonsynonymous substitution, and test whether they were positively selected along any branch of the phylogeny. We then map these genes onto the annotated Capsicum genome and test whether rates of evolution, gene history, and gene ontology vary significantly with gene position. We observed substantial discordance among gene trees. A bifurcating species tree placing Capsicum within a paraphyletic Lycianthes was supported over all phylogenetic networks. Rates of synonymous and nonsynonymous substitution varied 41-fold and 130-fold among genes, respectively, and were significantly lower in pericentromeric regions. We found that results of concordance tree analyses vary depending on the subset of genes used, and that genes within the pericentromeric regions only capture a portion of the observed discordance. We identified 787 genes that have been positively selected throughout the diversification history of Capsiceae, and discuss the importance of these genes as targets for investigation of economically important traits in the domesticated peppers.

Spatial phylogenetics reveals evolutionary constraints on the assembly of a large regional flora.

PREMISE OF THE STUDY: We used spatial phylogenetics to analyze the assembly of the Wisconsin flora, linking processes of dispersal and niche evolution to spatial patterns of floristic and phylogenetic diversity and testing whether phylogenetic niche conservatism can account for these patterns. METHODS: We used digitized records and a new molecular phylogeny for 93% of vascular plants in Wisconsin to estimate spatial variation in species richness and phylogenetic α and ß diversity in a native flora shaped mainly by postglacial dispersal and response to environmental gradients. We developed distribution models for all species and used these to infer fine-scale variation in potential diversity, phylogenetic distance, and interspecific range overlaps. We identified 11 bioregions based on floristic composition, mapped areas of neo- and paleo-endemism to establish new conservation priorities and predict how community-assembly patterns should shift with climatic change. KEY RESULTS: Spatial phylogenetic turnover most strongly reflects differences in temperature and spatial distance. For all vascular plants, assemblages shift from phylogenetically clustered to overdispersed northward, contrary to most other studies. This pattern is lost for angiosperms alone, illustrating the importance of phylogenetic scale. CONCLUSIONS: Species ranges and assemblage composition appear driven primarily by phylogenetic niche conservatism. Closely related species are ecologically similar and occupy similar territories. The average level and geographic structure of plant phylogenetic diversity within Wisconsin are expected to greatly decline over the next half century, while potential species richness will increase throughout the state. Our methods can be applied to allochthonous communities throughout the world.

Monocot plastid phylogenomics, timeline, net rates of species diversification, the power of multi-gene analyses, and a functional model for the origin of monocots.

PREMISE OF THE STUDY: We present the first plastome phylogeny encompassing all 77 monocot families, estimate branch support, and infer monocot-wide divergence times and rates of species diversification. METHODS: We conducted maximum likelihood analyses of phylogeny and BAMM studies of diversification rates based on 77 plastid genes across 545 monocots and 22 outgroups. We quantified how branch support and ascertainment vary with gene number, branch length, and branch depth. KEY RESULTS: Phylogenomic analyses shift the placement of 16 families in relation to earlier studies based on four plastid genes, add seven families, date the divergence between monocots and eudicots+Ceratophyllum at 136 Mya, successfully place all mycoheterotrophic taxa examined, and support recognizing Taccaceae and Thismiaceae as separate families and Arecales and Dasypogonales as separate orders. Only 45% of interfamilial divergences occurred after the Cretaceous. Net species diversification underwent four large-scale accelerations in PACMAD-BOP Poaceae, Asparagales sister to Doryanthaceae, Orchidoideae-Epidendroideae, and Araceae sister to Lemnoideae, each associated with specific ecological/morphological shifts. Branch ascertainment and support across monocots increase with gene number and branch length, and decrease with relative branch depth. Analysis of entire plastomes in Zingiberales quantifies the importance of non-coding regions in identifying and supporting short, deep branches. CONCLUSIONS: We provide the first resolved, well-supported monocot phylogeny and timeline spanning all families, and quantify the significant contribution of plastome-scale data to resolving short, deep branches. We outline a new functional model for the evolution of monocots and their diagnostic morphological traits from submersed aquatic ancestors, supported by convergent evolution of many of these traits in aquatic Hydatellaceae (Nymphaeales).

Divergence times, historical biogeography, and shifts in speciation rates of Myrtales.

We examine the eudicot order Myrtales, a clade with strong Gondwanan representation for most of its families. Although previous phylogenetic studies greatly improved our understanding of intergeneric and interspecific relationships within the order, our understanding of inter-familial relationships still remains unresolved; hence, we also lack a robust time-calibrated chronogram to address hypotheses (e.g., biogeography and diversification rates) that have implicit time assumptions. Six loci (rbcL, ndhF, matK, matR, 18S, and 26S) were amplified and sequenced for 102 taxa across Myrtales for phylogenetic reconstruction and ten fossil priors were utilized to produce a chronogram in BEAST. Combretaceae is identified as the sister clade to all remaining families with moderate support, and within the latter clade, two strongly supported groups are seen: (1) Onagraceae+Lythraceae, and (2) Melastomataceae+the Crypteroniaceae, Alzateaceae, Penaeaceae clade along with Myrtaceae+Vochysiaceae. Divergence time estimates suggest Myrtales diverged from Geraniales ∼124Mya during the Aptian of the Early Cretaceous. The crown date for Myrtales is estimated at ∼116Mya (Albian-Aptian). BioGeoBEARS showed significant improvement in the likelihood score when the "jump dispersal" parameter was added. South America and/or Africa are implicated as important ancestral areas in all deeper nodes. BAMM analyses indicate that the best configuration included three significant shifts in diversification rates within Myrtales: near the crown of Melastomataceae (∼67-64Mya), along the stem of subfamily Myrtoideae (Myrtaceae; ∼75Mya), and along the stem of tribe Combreteae (Combretaceae; ∼50-45Mya). Issues with conducting diversification analyses more generally are examined in the context of scale, taxon sampling, and larger sets of phylogenetic trees.

Evolution of geographical place and niche space: Patterns of diversification in the North American sedge (Cyperaceae) flora.

The role of geography and ecology in speciation are often discussed in the context of phylogenetic niche conservatism (PNC), the propensity of lineages to retain ancestral niche related traits. However, a recent paradigm shift focuses instead on measuring divergence of these traits in conjunction with patterns of speciation. Under this framework, we analyzed the diversification of North America's third most diverse family, Cyperaceae ("sedges"), using a modified Parsimony Analysis of Endemicity approach to identify floristic regions and ordination statistics to quantify species distribution in a continuous manner. Utilizing over 200,000 georeferenced specimens, we characterized the geographical distribution and climatic and edaphic niche space occupied by each species. We constructed a supermatrix phylogeny of the North American sedge flora, aided in part by the sequencing of all sedges of Wisconsin, and employed a multifaceted approach to assess the role of geographical and ecological divergence on lineage diversification. In addition to measuring phylogenetic signal for these traits, we also measured pairwise phylogenetic distance of species within floristic regions, calculated rates of speciation, and tested for correlations of speciation rate to tempo of geographical and ecological evolution. Our analyses consistently show that evolutionarily related species tend to be geographically unrelated. Rates of geographical and ecological diversification are closely linked to tempo of speciation, and exploration of geographical place coincides with divergence in ecological niche space. We highlight the benefits of treating geography in a continuous manner, and stress the importance of employing a diverse suite of analytical approaches in testing hypotheses regarding the evolution of range and niche.

Phylogenomics and historical biogeography of the monocot order Liliales: out of Australia and through Antarctica.

We present the first phylogenomic analysis of relationships among all ten families of Liliales, based on 75 plastid genes from 35 species in 29 genera, and 97 additional plastomes stratified across angiosperm lineages. We used a supermatrix approach to extend our analysis to 58 of 64 genera of Liliales, and calibrated the resulting phylogeny against 17 fossil dates to produce a new timeline for monocot evolution. Liliales diverged from other monocots 124 Mya and began splitting into separate families 113 Mya. Our data support an Australian origin for Liliales, with close relationships between three pairs of lineages (Corsiaceae/Campynemataceae, Philesiaceae/Ripogonaceae, tribes Alstroemerieae/Luzuriageae) in South America and Australia or New Zealand reflecting teleconnections of these areas via Antarctica. Long-distance dispersal (LDD) across the Pacific and Tasman Sea led to re-invasion of New Zealand by two lineages (Luzuriaga, Ripogonum); LDD allowed Campynemanthe to colonize New Caledonia after its submergence until 37 Mya. LDD permitted Colchicaceae to invade East Asia and Africa from Australia, and re-invade Africa from Australia. Periodic desert greening permitted Gloriosa and Iphigenia to colonize Southeast Asia overland from Africa, and Androcymbium-Colchicum to invade the Mediterranean from South Africa. Melanthiaceae and Liliaceae crossed the Bering land-bridge several times from the Miocene to the Pleistocene.

Orchid phylogenomics and multiple drivers of their extraordinary diversification.

Orchids are the most diverse family of angiosperms, with over 25 000 species,more than mammals, birds and reptiles combined. Tests of hypotheses to account for such diversity have been stymied by the lack of a fully resolved broad-scale phylogeny. Here,we provide such a phylogeny, based on 75 chloroplast genes for 39 species representing all orchid subfamilies and 16 of 17 tribes, time-calibrated against 17 angiosperm fossils. Asupermatrix analysis places an additional 144 species based on three plastid genes. Orchids appear to have arisen roughly 112 million years ago (Mya); the subfamilies Orchidoideae and Epidendroideae diverged from each other at the end of the Cretaceous; and the eight tribes and three previously unplaced subtribes of the upper epidendroids diverged rapidly from each other between 37.9 and 30.8 Mya. Orchids appear to have undergone one significant acceleration of net species diversification in the orchidoids, and two accelerations and one deceleration in the upper epidendroids. Consistent with theory, such accelerations were correlated with the evolution of pollinia, the epiphytic habit, CAM photosynthesis, tropical distribution (especially in extensive cordilleras),and pollination via Lepidoptera or euglossine bees. Deceit pollination appears to have elevated the number of orchid species by one-half but not via acceleration of the rate of net diversification. The highest rate of net species diversification within the orchids (0.382 sp sp(-1) My(-1)) is 6.8 times that at the Asparagales crown.

Spatial phylogenetics of two topographic extremes of the Hengduan Mountains in southwestern China and its implications for biodiversity conservation.

Previous attempts to elucidate the drivers of speciation mechanisms and spatial distribution patterns of biodiversity in mountain regions have treated different floras within a single geological region as one flora, ignoring the potential contributions of high habitat/ecosystem heterogeneity. Furthermore, current conservation strategies largely focus on forest ecosystems and/or specific flagship species, ignoring marginal ecosystems, leaving species in these ecosystems at risk. Here, we compared the spatial patterns of biodiversity and the potential drivers of these patterns in the river valley and subnival ecosystems of the Hengduan Mountains region (HDM) in southwestern China. Specifically, we compared spatial patterns of diversity, endemism, and threatened species in these ecosystems based on both traditional measurements and recent phylogenetic approaches. We then examined how those patterns were related to environmental factors and human activity in these same regions. We found that the middle-southern HDM supports the highest diversity and endemism for the river valley and subnival ecosystems; however, the distribution patterns of neo- and paleo-endemism in these two ecosystems differ. Regression models indicate that habitat diversity and paleo-climatic fluctuation are important drivers of diversity and endemism for these two ecosystems. Temperature and precipitation, however, showed different influences on the spatial patterns in different ecosystems. Categorical analysis of neo- and paleo-endemism (CANAPE) indicated that most endemism centers are not covered by current nature reserves. Moreover, the intensity of human activity is highest in the southern and southeastern HDM, which coincides with the distribution patterns of diversity, mixed-endemism and high-priority (and threatened) species. These findings suggest that different floras within a single geographic/floristic region respond differently to environmental factors and show different spatial phylogenetic patterns. We, therefore, recommend that future research into the drivers of biodiversity consider the contributions of various ecosystem types within a single geological region. This study also provides a theoretical basis for protecting habitat diversity. Our work confirms that current conservation efforts are insufficient to protect ecosystem diversity in the river valley and subnival ecosystems of the Hengduan Mountains. Therefore, we recommend the establishment of nature reserves in the regions identified in this study; furthermore, we strongly recommend improving current and establishing new management policies for biodiversity conservation in this region.

Genomic Polymorphism Associated with the Emergence of Virulent Isolates of Mycobacterium bovis in the Nile Delta.

Mycobacterium bovis is responsible for bovine tuberculosis in both animals and humans. Despite being one of the most important global zoonotic disease, data related to the ecology and pathogenicity of bovine tuberculosis is scarce, especially in developing countries. In this report, we examined the dynamics of M. bovis transmission among dairy cattle in the Nile Delta of Egypt. Animals belonging to 27 herds from 7 governorates were tested by the Single Intradermal Comparative Skin Tuberculin (SICST), as a preliminary screen for the presence of bovine tuberculosis. Positive SICST reactors were identified in 3% of the animals spread among 40% of the examined herds. Post-mortem examination of slaughtered reactors confirmed the presence of both pulmonary and/or digestive forms of tuberculosis in > 50% of the examined animals. Targeted and whole-genome analysis of M. bovis isolates indicated the emergences of a predominant spoligotype (SB0268) between 2013-2015, suggesting a recent clonal spread of this isolate within the Nile Delta. Surprisingly, 2 isolates belonged to M. bovis BCG group, which are not allowed for animal vaccination in Egypt, while the rest of isolates belonged to the virulent M. bovis clonal complex European 2 present in Latin America and several European countries. Analysis of strain virulence in the murine model of tuberculosis indicated the emergence of a more virulent strain (MBE4) with a specific genotype. More analysis is needed to understand the molecular basis for successful spread of virulent isolates of bovine tuberculosis among animals and to establish genotype/phenotype association.

Publisher Correction: Genomic analysis of the emergence of drug-resistant strains of Mycobacterium tuberculosis in the Middle East.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genomic analysis of the emergence of drug-resistant strains of Mycobacterium tuberculosis in the Middle East.

Tuberculosis (TB) represents a significant challenge to public health authorities, especially with the emergence of drug-resistant (DR) and multidrug-resistant (MDR) isolates of Mycobacterium tuberculosis. We sought to examine the genomic variations among recently isolated strains of M. tuberculosis in two closely related countries with different population demography in the Middle East. Clinical isolates of M. tuberculosis from both Egypt and Saudi Arabia were subjected to phenotypic and genotypic analysis on gene and genome-wide levels. Isolates with MDR phenotypes were highly prevalent in Egypt (up to 35%) despite its relatively stable population structure (sympatric pattern). MDR-TB isolates were not identified in the isolates from Saudi Arabia despite its active guest worker program (allopatric pattern). However, tuberculosis isolates from Saudi Arabia, where lineage 4 was more prevalent (>65%), showed more diversity than isolates from Egypt, where lineage 3 was the most prevalent (>75%). Phylogenetic and molecular dating analyses indicated that lineages from Egypt were recently diverged (~78 years), whereas those from Saudi Arabia were diverged by over 200 years. Interestingly, DR isolates did not appear to cluster together or spread more widely than drug-sensitive isolates, suggesting poor treatment as the main cause for emergence of drug resistance rather than more virulence or more capacity to persist.

Divergence and isolation of cryptic sympatric taxa within the annual legume Amphicarpaea bracteata.

The amphicarpic annual legume Amphicarpaea bracteata is unusual in producing aerial and subterranean cleistogamous flowers that always self-fertilize and, less commonly, aerial chasmogamous flowers that outcross. Although both morphologic and genetic variants are known in this highly selfing species, debate continues over whether this variation is continuous, reflecting the segregation of standing genetic variation, or discontinuous, reflecting distinct taxa that rarely intercross. We characterized SNP variation in 128 individuals in southern Wisconsin to assess within- and among-population variation at 3928 SNPs. We also assessed genotype and leaf morphology in an additional 76 individuals to connect phenotypic variation with genetic variation. Genetic variation maps onto three strongly divergent and highly inbred genetic groups showing little relation to site location. Each group has a distinct phenotype, but the divergence of these groups differs from the varietal divisions previously identified based on morphological characters. Like previous authors, we argue that the taxonomy of this species should be revised. Despite extensive sympatry, estimates of among-group migration rates are low, and hybrid individuals were at low frequency (<2%) in our dataset. Restricted gene flow likely results from high selfing rates and partial reproductive incompatibility as evidenced by the U-shaped distribution of pairwise F ST values reflecting "islands" of genomic divergence. These islands may be associated with hybrid incompatibility loci that arose in allopatry. The coexistence of lineages within sites may reflect density-dependent attack by species-specific strains of pathogenic fungi and/or root-nodulating bacteria specializing on distinct genotypes.