USGS Science⌕ Search

SEARCH · USGS Science

Results for “Genes”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58Linked to original sources
USGS2026

Cnidarian–algal partnerships structure bacterial communities during strobilation in Cassiopea xamachana

Cnidarian–algal (Symbiodiniaceae) symbioses rely on complex interactions among the cnidarian host, algal symbionts, and associated bacterial communities. In the upside-down jellyfish Cassiopea xamachana , the polyp-to-medusa transition (strobilation) requires the establishment of symbiosis with Symbiodiniaceae algal partners, yet bacterial community dynamics during this developmental process remain unknown. Here, we experimentally induced symbiosis in aposymbiotic polyps using four algal treatments: xenic Symbiodinium microadriaticum (native symbiont), xenic Breviolum minutum , antibiotic-treated B. minutum , and a photosynthetically impaired B. minutum mutant. We combined 16S rRNA gene sequencing with measurements of photosynthetic efficiency, asexual budding, and algal surface N-glycan profiles to characterize holobiont assembly during symbiosis onset and strobilation. Algal treatment structured bacterial communities in both algal cultures and polyp tissues. Our analyses identified a set of amplicon sequence variants that consistently distinguished strobilating polyps from non-strobilating aposymbiotic and mutant polyps, in addition to potential bacterial biomarkers associated with successful metamorphosis. Strobilation was associated with the enrichment of bacterial communities putatively involved in sulfur and nitrogen cycling, whereas non-strobilating aposymbiotic and mutant polyps were characterized by opportunistic bacteria and increased community variability. Together, these results reveal coordinated changes in algal physiology, surface glycan profiles, and bacterial community structure associated with successful strobilation in C. xamachana and support a model in which tripartite host–alga–bacteria interactions influence cnidarian life stage transitions.

ISME Communications↗
USGS2024

Transglobal spread of an ecologically relevant sea urchin parasite

Mass mortality of the dominant coral reef herbivore Diadema antillarum in the Caribbean in the early 1980s contributed to a persistent phase shift from coral- to algal-dominated reefs. In 2022, a scuticociliate most closely related to Philaster apodigitiformis caused further mass mortality of D. antillarum across the Caribbean, leading to >95% mortality at affected sites. Mortality was also reported in the related species Diadema setosum in the Mediterranean in 2022, though the causative agent of the Mediterranean outbreak has not yet been determined. In April 2023, mass mortality of Diadema setosum occurred along the Sultanate of Oman's coastline. Urchins displayed signs compatible with scuticociliatosis including abnormal behavior, drooping and loss of spines, followed by tissue necrosis and death. Here we report the detection of an 18S rRNA gene sequence in abnormal urchins from Muscat, Oman, that is identical to the Philaster strain responsible for D. antillarum mass mortality in the Caribbean. We also show that scuticociliatosis signs can be elicited in Diadema setosum by experimental challenge with the cultivated Philaster strain associated with Caribbean scuticociliatosis. These results demonstrate the Philaster sp. associated with D. antillarum mass mortality has rapidly spread to geographically distant coral reefs, compelling global-scale awareness and monitoring for this devastating condition through field surveys, microscopy, and molecular microbiological approaches, and prompting investigation of long-range transmission mechanisms.

ISME Journal↗
USGS2025

Manipulation of the Symbiodiniaceae microbiome confers multigenerational impacts on symbioses and reproductive ecology of its Exaiptasia diaphana host

Symbiodiniaceae-associated microbiota strongly influence cnidarian symbioses. We systematically reduced the bacterial and fungal communities associated with Symbiodiniaceae to study potential effects on the cnidarian holobiont Exaiptasia diaphana (Aiptasia). Clonal anemones were inoculated with xenic Breviolum minutum (SSB01) and microbiome-manipulated cultures after antibacterial or antifungal treatment. The asexual reproduction of pedal laceration allowed for three generations of clonal aposymbiotic Aiptasia to be utilised in this study, from the initial adult generation (G0), to the first (G1), and second (G2) generation. We inoculated small and large G1 Aiptasia with SSB01 algae and monitored onset of symbiosis, rate of algal proliferation, and holobiont characteristics. Sequencing the 16S and 18S rRNA gene regions identified significant differences in the bacterial and fungal communities of the G0 and G1 generations, alongside differences between the size classes of small and large G1 anemones. The microbiome of larger G1 individuals was distinct to the smaller G1 anemones, suggesting a microbiome maturation process. Control Breviolum minutum cultures exhibited a significantly greater proliferation rate in large G1 anemones when compared to antibacterial or antifungal treated cultures, whereas the opposite trend was documented in the small G1 anemones. Although no differences were observed between algal photochemical parameters, or the growth and polyp activity of G1 juveniles, we observed a significant influence in the production of G2 clones between treatments. Overall, we provide strong ecological implications of manipulating Symbiodiniaceae microbiome, not for the algae themselves, but for the maturation of the host Aiptasia, as well as for the cnidarian holobiont over multiple generations.

The ISME Journal↗
USGS2025

Sulfide stress tolerance as a controller of methane production in temperate wetlands

Wetlands are a major source of methane emissions and contribute to the observed increase in atmospheric methane over the last 20 years. Methane production in wetlands is the final step of carbon decomposition performed by anaerobic archaea. Although hydrogen/carbon dioxide and acetate are the substrates most often attributed to methanogenesis, other substrates—such as methylated compounds—may additionally play important roles in driving methane production in wetland systems. Here we conducted mesocosm experiments combined with genome-resolved metatranscriptomics to investigate the impact of diverse methanogenic substrate amendment on methanogenesis in two high methane-emitting wetlands with distinct geochemistry, termed P7 and P8. Methanol amendment resulted in high methane production at both sites, whereas acetate and formate amendment only stimulated methanogenesis in P7 mesocosms, where aqueous sulfide concentrations were lower. In P7 sediments, formate amendment fueled acetogenic microbes that produced acetate, which was subsequently utilized by acetoclastic methanogens. In contrast to expression profiles in P7 mesocosms, active methylotrophic methanogen genomes from P8 showed increased expression of genes related to membrane remodeling and DNA damage repair, indicative of stress tolerance mechanisms to counter sulfide toxicity. Methylotrophic methanogenesis generates higher free energy yields than acetoclastic methanogenesis, which likely enables allocation of more energy toward stress responses. These findings contribute to the growing body of literature highlighting methylotrophic methanogenesis as an important methane production pathway in wetlands. By using less competitive substrates like methanol that provide greater energy yields, methylotrophic methanogens may invest in physiological strategies that provide competitive advantages across a range of environmental stresses.

North Dakota↗
USGS2026

Diverse cyanopeptides follow distinct temporal succession patterns in freshwater harmful algal blooms

Toxic cyanobacterial harmful algal blooms (cyanoHABs) threaten freshwater resources globally and are intensifying with increasing eutrophication. Bloom toxicity is strongly influenced by intraspecific variation in the biosynthetic repertoires of toxic cyanobacteria, yet few studies examine the diversity of cyanobacterial cyanopeptides beyond hepatotoxic microcystins . To understand the dynamics and drivers of cyanopeptide diversity in cyanoHABs, we analyzed temporal patterns of cyanobacteria, metabolites, and their biosynthetic gene clusters (BGCs) in western Lake Erie using a 7-year time series (2016–2022) of metagenomic and metabolomic data. Our findings demonstrate that shifts from Microcystis to Dolichospermum occur later in the bloom season, coinciding with lower temperatures. Modules of co-varying BGCs (biosynthesis modules) from these genera were identified with hierarchical clustering, with uncharacterized BGCs among the most abundant. Biosynthesis modules rich in nonribosomal peptide synthetases (NRPS) peaked in early August, coinciding with elevated levels of inorganic nitrogen, warmer temperatures, and high Microcystis abundance. In contrast, modules rich in polyketide synthases (PKS) and ribosomally synthesized and post-translationally modified peptides (RiPPs) peaked following the Microcystis maximum in mid-August. Metabolomic analyses confirmed that metabolites followed shared seasonal patterns with their associated biosynthesis modules, forming three phases characterized by (i) microcystins, (ii) anabaenopeptins and aeruginosins, and (iii) aerucyclamides. These phases co-varied with bottom-up and top-down pressures, with later phases coinciding with increased microbially processed organic nitrogen and reduced detection of grazers. This study demonstrates consistent seasonal patterns of cyanobacterial metabolite succession and co-occurrence beyond microcystins, suggesting tradeoffs between biosynthetic resource demands and ecological controls.

Ohio↗
USGS2025

Identical sequence types of Yersinia ruckeri associated with lethal disease in wild-caught invasive Blue Catfish and cultured hybrid catfish (Channel Catfish ♀ × Blue Catfish ♂) from disparate aquatic ecosystems

Objective The Blue Catfish Ictalurus furcatus is commonly raised in warmwater aquaculture in the United States to produce Channel Catfish I. punctatus × Blue Catfish hybrids. It is also a prominent aquatic invasive species of concern in the mid-Atlantic region of the United States. Here, Yersina ruckeri was isolated from moribund Blue Catfish and hybrid catfish from disparate regions of the USA. The goal of the research here was to compare these Y. ruckeri strains to each other and other known strains for which adequate sequence data was available. In addition, we sought to determine if the strain from Blue Catfish was pathogenic to Rainbow Trout Oncorhynchus mykiss . Methods Moribund hybrid catfish from culture ponds in Mississippi were processed for diagnostic evaluation in March 2016. In April 2022, a moribund Blue Catfish specimen was collected from a tributary of the Nanticoke River in Maryland. Bacterial isolates were identified and characterized using biochemical tests, antimicrobial sensitivity screening, serotyping, and complete or partial genome sequencing. Disease pathology was described via histology. The isolate from Blue Catfish was used in challenge experiments to determine if it was pathogenic to Rainbow Trout. Multilocus sequencing typing was conducted using the PubMLST database. Results Biochemical testing was consistent with Y. ruckeri . A draft genome of the Y. ruckeri isolate was assembled based on Oxford Nanopore Technology sequencing and identified a single genomic replicon (3,791,418 bp) consistent in size to other Y. ruckeri genomes and a pLT plasmid (60, 933 bp). The challenge study demonstrated no significant virulence of this isolate for Rainbow Trout ( Y. ruckeri ). This isolate was most similar to other strains isolated from ictalurids. Notably, the gyrase B gene from this isolate was identical to that of archived strains isolated from moribund Mississippi hybrid catfish aquaculture during 2016 and these isolates share identical PubMLST sequence type profiles. Similarly, they shared a pLT plasmid that differed by only 6 bp. This plasmid has never been reported from trout isolates and appears to be unique to ictalurids. Conclusions Analyses here provide preliminary genetic evidence that geographically distant (Maryland and Mississippi, USA) isolates of Y. ruckeri from ictalurids are genetically similar to each other and Y. ruckeri (strain SC09) that infects ictalurids in China. This strain is not a biothreat to Rainbow Trout at typical culture temperatures.

Journal of Aquatic Animal Health↗
USGS2001

Subspecific affinity of black bears in the White River National Wildlife Refuge

The black bear population of the White River National Wildlife Refuge (NWR) is adjacent to populations of black bear in Louisiana ( Urusus americanus luteolus ) which are listed as threatened under the U.S. Endangered Species Act. Wildlife management plans can pose restrictions on bear harvests and timber extraction; therefore the management plan for the White River NWR is sensitive to subspecific classification of its bear population. The objective of this study was to analyze genetic variation in the White River NWR and seven adjacent populations of black bears to assess the subspecific affinity of the White River NWR population. Here we report the variation at seven microsatellite DNA loci among eight black bear populations. The patterns of genetic variation gave strong support for distinguishing a southern group of black bears comprised of the White River, Arkansas; Tensas River, Louisiana; Upper Atchafalaya, Louisiana; Lower Atchafalaya, Louisiana; and Alabama/Mississippi populations. Phylogenetic analysis of individual variation suggested that historical black bear introductions into Arkansas and Louisiana affected gene pools of certain southern receiving populations, but did not significantly change interpopulation relatedness. Phylogenetic inferences at both the population and individual levels support the hypothesis that the White River NWR population of black bears belongs to the U. a. luteolus subspecies.

Arkansas↗
USGS2022

Give and take: Effects of genetic admixture on mutation load in endangered Florida panthers

Genetic admixture is a biological event inherent to genetic rescue programs aimed at the long-term conservation of endangered wildlife. Although the success of such programs can be measured by the increase in genetic diversity and fitness of subsequent admixed individuals, predictions supporting admixture costs to fitness due to the introduction of novel deleterious alleles are necessary. Here, we analyzed nonsynonymous variation from conserved genes to quantify and compare levels of mutation load (i.e. proportion of deleterious alleles and genotypes carrying these alleles) among endangered Florida panthers and non-endangered Texas pumas. Specifically, we used canonical (i.e. non-admixed) Florida panthers, Texas pumas, and F 1 (canonical Florida × Texas) panthers dating from a genetic rescue program and Everglades National Park panthers with Central American ancestry resulting from an earlier admixture event. We found neither genetic drift nor selection significantly reduced overall proportions of deleterious alleles in the severely bottlenecked canonical Florida panthers. Nevertheless, the deleterious alleles identified were distributed into a disproportionately high number of homozygous genotypes due to close inbreeding in this group. Conversely, admixed Florida panthers (either with Texas or Central American ancestry) presented reduced levels of homozygous genotypes carrying deleterious alleles but increased levels of heterozygous genotypes carrying these variants relative to canonical Florida panthers. Although admixture is likely to alleviate the load of standing deleterious variation present in homozygous genotypes, our results suggest that introduced novel deleterious alleles (temporarily present in heterozygous state) in genetically rescued populations could potentially be expressed in subsequent generations if their effective sizes remain small.

Journal of Heredity↗
USGS2023

Assembly of the largest squamate reference genome to date: The western fence lizard, Sceloporus occidentalis

Spiny lizards (genus Sceloporus ) have long served as important systems for studies of behavior, thermal physiology, dietary ecology, vector biology, speciation, and biogeography. The western fence lizard, Sceloporus occidentalis , is found across most of the major biogeographical regions in the western United States and northern Baja California, Mexico, inhabiting a wide range of habitats, from grassland to chaparral to open woodlands. As small ectotherms, Sceloporus lizards are particularly vulnerable to climate change, and S. occidentalis has also become an important system for studying the impacts of land use change and urbanization on small vertebrates. Here, we report a new reference genome assembly for S. occidentalis , as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomics strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 608 scaffolds spanning 2,856 Mb, has a contig N50 of 18.9 Mb, a scaffold N50 of 98.4 Mb, and BUSCO completeness score of 98.1% based on the tetrapod gene set. This reference genome will be valuable for understanding ecological and evolutionary dynamics in S. occidentalis , the species status of the California endemic island fence lizard ( S. becki ), and the spectacular radiation of Sceloporus lizards.

California↗
USGS2023

Genomic insights into isolation of the threatened Florida crested caracara (Caracara plancus)

We conducted a population genomic study of the crested caracara ( Caracara plancus ) using samples ( n = 290) collected from individuals in Florida, Texas, and Arizona, United States. Crested caracaras are non-migratory raptors ranging from the southern tip of South America to the southern United States, including a federally protected relict population in Florida long thought to have been isolated since the last ice age. Our objectives were to evaluate genetic diversity and population structure of Florida’s apparently isolated population and to evaluate taxonomic relationships of crested caracaras at the northern edge of their range. Using DNA purified from blood samples, we conducted double-digest restriction site associated DNA sequencing and sequenced the mitochondrial ND2 gene. Analyses of population structure using over 9,000 SNPs suggest that two major clusters are best supported, one cluster including only Florida individuals and the other cluster including Arizona and Texas individuals. Both SNPs and mitochondrial haplotypes reveal the Florida population to be highly differentiated genetically from Arizona and Texas populations, whereas, Arizona and Texas populations are moderately differentiated from each other. The Florida population’s mitochondrial haplotypes form a separate monophyletic group, while Arizona and Texas populations share mitochondrial haplotypes. Results of this study provide substantial genetic evidence that Florida’s crested caracaras have experienced long-term isolation from caracaras in Arizona and Texas and thus, represent a distinct evolutionary lineage possibly warranting distinction as an Evolutionarily Significant Unit (ESU) or subspecies. This study will inform conservation strategies focused on long-term survival of Florida’s distinct, panmictic population.

Journal of Heredity↗
USGS2026

A chromosome-level genome assembly of a vernal pool specialist amphibian, the Western Spadefoot, Spea hammondii

We assembled and annotated a chromosome-level genome for the Western Spadefoot, Spea hammondii (Anura, Scaphiopodidae) representing one of only three amphibians included in the California Conservation Genomics Project (CCGP). Spea hammondii is a vernal pool breeding anuran native to California and northwestern Baja California which has undergone both range contractions and local extirpations across its distribution, primarily due to habitat loss and degradation and drought. The species is recognized by the state of California as a Species of Special Concern and is proposed for listing under the United States Endangered Species Act. Using the established CCGP pipeline, this S. hammondii genome was produced using Pacific Biosciences HiFi long-reads and Omni-C proximity ligation, resulting in a de novo genome assembly 1.14 Gb in length, distributed across 479 scaffolds (scaffold N50 = 120.8 Mb; largest scaffold = 183.6 Mb) with a BUSCO completeness score of 90.9% using a conserved tetrapod ortholog set. Our assembly shows high base accuracy (quality value [QV] = 63.7) and low frameshift error in coding regions (QV 50.42). Annotation of this genome yielded 20,434 genes with a BUSCO completeness score of 94.7%. This genome assembly, in combination with range-wide resequencing data from CCGP, will facilitate statewide population genomic assessments to delineate conservation units, quantify inbreeding and genomic load, and test for adaptive variation associated with vernal pool hydrology and drought tolerance, all of which are important considerations in the proposed federal listing.

Journal of Heredity↗
USGS2012

Genetic population substructure in bison at Yellowstone National Park

The Yellowstone National Park bison herd is 1 of only 2 populations known to have continually persisted on their current landscape since pre-Columbian times. Over the last century, the census size of this herd has fluctuated from around 100 individuals to over 3000 animals. Previous studies involving radiotelemetry, tooth wear, and parturition timing provide evidence of at least 2 distinct groups of bison within Yellowstone National Park. To better understand the biology of Yellowstone bison, we investigated the potential for limited gene flow across this population using multilocus Bayesian clustering analysis. Two genetically distinct and clearly defined subpopulations were identified based on both genotypic diversity and allelic distributions. Genetic cluster assignments were highly correlated with sampling locations for a subgroup of live capture individuals. Furthermore, a comparison of the cluster assignments to the 2 principle winter cull sites revealed critical differences in migration patterns across years. The 2 Yellowstone subpopulations display levels of differentiation that are only slightly less than that between populations which have been geographically and reproductively isolated for over 40 years. The identification of cryptic population subdivision and genetic differentiation of this magnitude highlights the importance of this biological phenomenon in the management of wildlife species.

Journal of Heredity↗
USGS2019

Elevated heterozygosity in adults relative to juveniles provides evidence of viability selection on eagles and falcons

Viability selection yields adult populations that are more genetically variable than those of juveniles, producing a positive correlation between heterozygosity and survival. Viability selection could be the result of decreased heterozygosity across many loci in inbred individuals and a subsequent decrease in survivorship resulting from the expression of the deleterious alleles. Alternatively, locus-specific differences in genetic variability between adults and juveniles may be driven by forms of balancing selection, including heterozygote advantage, frequency-dependent selection, or selection across temporal and spatial scales. We use a pooled-sequencing approach to compare genome-wide and locus-specific genetic variability between 74 golden eagle ( Aquila chrysaetos ), 62 imperial eagle ( Aquila heliaca ), and 69 prairie falcon ( Falco mexicanus ) juveniles and adults. Although genome-wide genetic variability is comparable between juvenile and adult golden eagles and prairie falcons, imperial eagle adults are significantly more heterozygous than juveniles. This evidence of viability selection may stem from a relatively smaller imperial eagle effective population size and potentially greater genetic load. We additionally identify ~2000 single-nucleotide polymorphisms across the 3 species with extreme differences in heterozygosity between juveniles and adults. Many of these markers are associated with genes implicated in immune function or olfaction. These loci represent potential targets for studies of how heterozygote advantage, frequency-dependent selection, and selection over spatial and temporal scales influence survivorship in avian species. Overall, our genome-wide data extend previous studies that used allozyme or microsatellite markers and indicate that viability selection may be a more common evolutionary phenomenon than often appreciated.

Journal of Heredity↗
USGS2021

Dynamic landscapes in northwestern North America structured populations of wolverines (Gulo gulo)

Cyclic climatic and glacial fluctuations of the Late Quaternary produced a dynamic biogeographic history for high latitudes. To refine our understanding of this history in northwestern North America, we explored geographic structure in a wide-ranging carnivore, the wolverine ( Gulo gulo ). We examined genetic variation in populations across mainland Alaska, coastal Southeast Alaska, and mainland western Canada using nuclear microsatellite genotypes and sequence data from the mitochondrial DNA (mtDNA) control region and Cytochrome b ( Cytb ) gene. Data from maternally inherited mtDNA reflect stable populations in Northwest Alaska, suggesting the region harbored wolverine populations since at least the Last Glacial Maximum (LGM; 21 Kya), consistent with their persistence in the fossil record of Beringia. Populations in Southeast Alaska are characterized by minimal divergence, with no genetic signature of long-term refugial persistence (consistent with the lack of pre-Holocene fossil records there). The Kenai Peninsula population exhibits mixed signatures depending on marker type: mtDNA data indicate stability (i.e., historical persistence) and include a private haplotype, whereas biparentally inherited microsatellites exhibit relatively low variation and a lack of private alleles consistent with a more recent Holocene colonization of the peninsula. Our genetic work is largely consistent with the early 20 th century taxonomic hypothesis that wolverines on the Kenai Peninsula belong to a distinct subspecies. Our finding of significant genetic differentiation of wolverines inhabiting the Kenai Peninsula, coupled with the peninsula’s burgeoning human population and the wolverine’s known sensitivity to anthropogenic impacts, provides valuable foundational data that can be used to inform conservation and management prescriptions for wolverines inhabiting these landscapes.

Alaska, British Columbia, Northwest Territories, N↗
USGS1995

A complex alloantigen system in Florida sandhill cranes, Grus canadensis pratensis: Evidence for the major histocompatibility (B) system

The B blood group system constitutes the major histocompatibility complex ( Mhc ) in birds. The Mhc is a cluster of genes largely devoted to the processing and presentation of antigen. The Mhc is highly polymorphic in many species and, thus, useful in the evaluation of genetic diversity for fitness traits within populations of a variety of animals. Correlations found between particular Mhc haplotypes and resistance to certain diseases emphasize the importance of understanding the functional significance of diversity of the Mhc , particularly in species threatened with extinction. As part of studies focused on genetic diversity in wild birds, serological techniques were used to define a highly polymorphic alloantigen system in seven families of Florida sandhill cranes ( Grus canadensis pratensis ). The results of analyses with antisera produced within the crane families and with chicken Mhc antigen-specific reagents revealed a single major alloantigen system that is likely the Mhc of the Florida sandhill crane. Preliminary experiments indicate that these crane alloantisera will provide a means of defining the Mhc in other species of cranes.

Journal of Heredity↗
USGS1996

Phylogenetic relationships within the Alcidae (Charadriiformes: Aves) inferred from total molecular evidence

The Alcidae is a unique assemblage of Northern Hemisphere seabirds that forage by "flying" underwater. Despite obvious affinities among the species, their evolutionary relationships are unclear. We analyzed nucleotide sequences of 1,045 base pairs of the mitochondrial cytochrome b gene and allelic profiles for 37 allozyme loci in all 22 extant species. Trees were constructed on independent and combined data sets using maximum parsimony and distance methods that correct for superimposed changes. Alternative methods of analysis produced only minor differences in relationships that were supported strongly by bootstrapping or standard error tests. Combining sequence and allozyme data into a single analysis provided the greatest number of relationships receiving strong support. Addition of published morphological and ecological data did not improve support for any additional relationship. All analyses grouped species into six distinct lineages: (1) the dovekie ( Alle alle ) and auks, (2) guillemots, (3) brachyramphine murrelets, (4) synthliboramphine murrelets, (5) true auklets, and (6) the rhinoceros auklet ( Cerorhinca monocerata ) and puffins. The two murres (genus Uria ) were sister taxa, and the black guillemot ( Cepphus grylle ) was basal to the other guillemots. The Asian subspecies of the marbled murrelet ( Brachyramphus marmoratus perdix ) was the most divergent brachyramphine murrelet, and two distinct lineages occurred within the synthliboramphine murrelets. Cassin's auklet ( Ptychoramphus aleuticus ) and the rhinoceros auklet were basal to the other auklets and puffins, respectively, and the Atlantic ( Fratercula arctica ) and horned ( Fratercula corniculata ) puffins were sister taxa. Several relationships among tribes, among the dovekie and auks, and among the auklets could not be resolved but resembled "star" phylogenies indicative of adaptive radiations at different depths within the trees.

Molecular Biology and Evolution↗
USGS2021

Integrating sequence capture and restriction-site associated DNA sequencing to resolve recent radiations of pelagic seabirds

The diversification of modern birds has been shaped by a number of radiations. Rapid diversification events make reconstructing the evolutionary relationships among taxa challenging due to the convoluted effects of incomplete lineage sorting (ILS) and introgression. Phylogenomic data sets have the potential to detect patterns of phylogenetic incongruence, and to address their causes. However, the footprints of ILS and introgression on sequence data can vary between different phylogenomic markers at different phylogenetic scales depending on factors such as their evolutionary rates or their selection pressures. We show that combining phylogenomic markers that evolve at different rates, such as paired-end double-digest restriction site-associated DNA (PE-ddRAD) and ultraconserved elements (UCEs), allows a comprehensive exploration of the causes of phylogenetic discordance associated with short internodes at different timescales. We used thousands of UCE and PE-ddRAD markers to produce the first well-resolved phylogeny of shearwaters, a group of medium-sized pelagic seabirds that are among the most phylogenetically controversial and endangered bird groups. We found that phylogenomic conflict was mainly derived from high levels of ILS due to rapid speciation events. We also documented a case of introgression, despite the high philopatry of shearwaters to their breeding sites, which typically limits gene flow. We integrated state-of-the-art concatenated and coalescent-based approaches to expand on previous comparisons of UCE and RAD-Seq data sets for phylogenetics, divergence time estimation, and inference of introgression, and we propose a strategy to optimize RAD-Seq data for phylogenetic analyses. Our results highlight the usefulness of combining phylogenomic markers evolving at different rates to understand the causes of phylogenetic discordance at different timescales.

Systematic Biology↗
USGS2016

Barrier displacement on a neutral landscape: Towards a theory of continental biogeography

Macroevolutionary theory posits three processes leading to lineage diversification and the formation of regional biotas: dispersal (species geographic range expansion), speciation (species lineage splitting), and extinction (species lineage termination). The Theory of Island Biogeography (TIB) predicts species richness values using just two of these processes; dispersal and extinction. Yet most species on Earth live on continents or continental shelves, and the dynamics of evolutionary diversification at regional and continental scales are qualitatively different from those that govern the formation of species richness on biogeographic islands. Certain geomorphological processes operating perennially on continental platforms displace barriers to gene flow and organismal dispersal, and affect all three terms of macroevolutionary diversification. For example, uplift of a dissected landscape and river capture both merge and separate portions of adjacent areas, allowing dispersal and larger geographic ranges, vicariant speciation and smaller geographic ranges, and extinction when range sizes are subdivided below a minimum persistence threshold. The TIB also does not predict many biogeographic and phylogenetic patterns widely observed in continentally distributed taxa, including: 1, power function-like species-area relationships; 2, log-normal distribution of species geographic range sizes, in which most species have restricted ranges (are endemic) and few species have broad ranges (are cosmopolitan); 3, mid-domain effects with more species towards the geographic center, and more early-branching, species-poor clades towards the geographic periphery; 4, exponential rates of net diversification with log-linear accumulation of lineages through geological time; and 5, power function-like relationships between species-richness and clade diversity, in which most clades are species-poor and few clades are species-rich. Current theory does not provide a robust mechanistic framework to connect these seemingly disparate patterns. Here we present SEAMLESS (Spatially-Explicit Area Model of Landscape Evolution by SimulationS) that generates clade diversification by moving geographic barriers on a continuous, neutral landscape. SEAMLESS is a neutral Landscape Evolution Model (LEM) that treats species and barriers as functionally equivalent with respect to model parameters. SEAMLESS differs from other model-based biogeographic methods (e.g. Lagrange, GeoSSE, BayArea, BioGeoBEARS) by modeling properties of dispersal barriers rather than areas, and by modeling the evolution of species lineages on a continuous landscape, rather than the evolution of geographic ranges along branches of a phylogeny. SEAMLESS shows how dispersal is required to maintain species richness and avoid clade-wide extinction, demonstrates that ancestral range size does not predict species richness, and provides a unified explanation for the suite of commonly observed biogeographic and phylogenetic patterns listed above. SEAMLESS explains how a simple barrier-displacement mechanism affects lineage diversification under neutral conditions, and is advanced here towards the formulation of a general theory of continental biogeography.

Systematic Biology↗