Uncloaking a cryptic, threatened rail with molecular markers: origins, connectivity and demography of a recently-discovered population

Reference genome of the Virginia rail, Rallus limicola

The Virginia rail, Rallus limicola, is a member of the family Rallidae, which also includes many other species of secretive and poorly studied wetland birds. It is recognized as a single species throughout its broad distribution in North America where it is exploited as a game bird, often with generous harvest limits, despite a lack of systematic population surveys and evidence of declines in many areas due to wetland loss and degradation. To help advance understanding of the phylogeography, biology, and ecology of this elusive species, we report the first reference genome assembly for the Virginia rail, produced as part of the California Conservation Genomics Project (CCGP). We produced a de novo genome assembly using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology with an estimated sequencing error rate of 0.191%. The assembly consists of 1,102 scaffolds spanning 1.39 Gb, with a contig N50 of 11.0 Mb, scaffold N50 of 25.3 Mb, largest contig of 45 Mb, and largest scaffold of 128.4 Mb. It has a high BUSCO completeness score of 96.9% and represents the first genome assembly available for the genus Rallus. This genome assembly will help resolve questions about the complex evolutionary history of rails and evaluate the potential of rails for adaptive evolution in the face of growing threats from climate change and habitat loss and fragmentation. It will also provide a valuable resource for rail conservation efforts by quantifying Virginia rail vagility, population connectivity, and effective population sizes.

Reference genome of the black rail, Laterallus jamaicensis

The black rail, Laterallus jamaicensis, is one of the most secretive and poorly understood birds in the Americas. Two of its five subspecies breed in North America: the Eastern black rail (L. j. jamaicensis), found primarily in the southern and mid-Atlantic states, and the California black rail (L. j. coturniculus), inhabiting California and Arizona, are recognized across the highly disjunct distribution. Population declines, due primarily to wetland loss and degradation, have resulted in conservation status listings for both subspecies. To help advance understanding of the phylogeography, biology, and ecology of this elusive species, we report the first reference genome assembly for the black rail, produced as part of the California Conservation Genomics Project (CCGP). We produced a de novo genome assembly using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology with an estimated sequencing error rate of 0.182%. The assembly consists of 964 scaffolds spanning 1.39 Gb, with a contig N50 of 7.4 Mb, scaffold N50 of 21.4 Mb, largest contig of 44.8 Mb, and largest scaffold of 101.2 Mb. The assembly has a high BUSCO completeness score of 96.8% and represents the first genome assembly available for the genus Laterallus. This genome assembly can help resolve questions about the complex evolutionary history of rails, assess black rail vagility and population connectivity, estimate effective population sizes, and evaluate the potential of rails for adaptive evolution in the face of growing threats from climate change, habitat loss and fragmentation, and disease.

Stability of patch-turnover relationships under equilibrium and nonequilibrium metapopulation dynamics driven by biogeography

Two controversial tenets of metapopulation biology are whether patch quality and the surrounding matrix are more important to turnover (colonisation and extinction) than biogeography (patch area and isolation) and whether factors governing turnover during equilibrium also dominate nonequilibrium dynamics. We tested both tenets using 18 years of surveys for two secretive wetland birds, black and Virginia rails, during (1) a period of equilibrium with stable occupancy and (2) after drought and arrival of West Nile Virus (WNV), which resulted in WNV infections in rails, increased extinction and decreased colonisation probabilities modified by WNV, nonequilibrium dynamics for both species and occupancy decline for black rails. Area (primarily) and isolation (secondarily) drove turnover during both stable and unstable metapopulation dynamics, greatly exceeding the effects of patch quality and matrix conditions. Moreover, slopes between turnover and patch characteristics changed little between equilibrium and nonequilibrium, confirming the overriding influences of biogeographic factors on turnover.

The taxonomic basis of subspecies listed as threatened and endangered under the endangered species act

More than 170 subspecies are listed as threatened or endangered under the US Endangered Species Act. Most of these subspecies were described decades ago on the basis of geographical variation in morphology using relatively primitive taxonomic methods. The US Fish and Wildlife Service defaults to subspecies descriptions by taxonomists working with specific groups of organisms, but there is no single definition of subspecies across plants and animals. Valid tests today usually entail molecular analyses of variation within and among populations, although there is no reason that behavioral, ecological or molecular characters could not be used, and include tests for significant differences between samples of the putative endangered subspecies and its nearest geographic relatives. We evaluated data gathered since subspecies listed under the ESA were described finding about one-third are valid (distinct evolutionary taxa), one-third are not, and one-third have not been tested. Therefore, it should not be assumed that because a subspecies occurs in a checklist, it is taxonomically valid. If the US Fish and Wildlife Service intends to continue listing subspecies, we suggest that they convene taxonomic experts representing various groups of organisms to provide a minimal set of criteria for a subspecies to be listed under the ESA.

Evolutionary History of the Galápagos Rail Revealed by Ancient Mitogenomes and Modern Samples

The biotas of the Galápagos Islands are one of the best studied island systems and have provided a broad model for insular species’ origins and evolution. Nevertheless, some locally endemic taxa, such as the Galápagos Rail Laterallus spilonota, remain poorly characterized. Owing to its elusive behavior, cryptic plumage, and restricted distribution, the Galápagos Rail is one of the least studied endemic vertebrates of the Galapagos Islands. To date, there is no genetic data for this species, leaving its origins, relationships to other taxa, and levels of genetic diversity uncharacterized. This lack of information is critical given the adverse fate of island rail species around the world in the recent past. Here, we examine the genetics of Galápagos Rails using a combination of mitogenome de novo assembly with multilocus nuclear and mitochondrial sequencing from both modern and historical samples. We show that the Galápagos Rail is part of the “American black rail clade”, sister to the Black Rail L. jamaicensis, with a colonization of Galápagos dated to 1.2 million years ago. A separate analysis of one nuclear and two mitochondrial markers in the larger population samples demonstrates a shallow population structure across the islands, possibly due to elevated island connectivity. Additionally, birds from the island Pinta possessed the lowest levels of genetic diversity, possibly reflecting past population bottlenecks associated with overgrazing of their habitat by invasive goats. The modern and historical data presented here highlight the low genetic diversity in this endemic rail species and provide useful information to guide conservation efforts.

Inferring the timing of long-distance dispersal between Rail metapopulations using genetic and isotopic assignments

The stochastic and infrequent nature of long-distance dispersal often makes it difficult to detect. We quantified the frequency, distance, and timing of long-distance dispersal in a nonmigratory, secretive wetland bird, the California Black Rail (Laterallus jamaicensis coturniculus), between an inland and a coastal metapopulation separated by greater than 100 km. Using 15 microsatellites in conjunction with stable carbon, nitrogen, and sulfur isotopes, we classified Rails as residents of their capture population, recent migrants that dispersed to their capture population less than one year before capture, established migrants that dispersed to their capture population more than one year before capture, and seasonal migrants that dispersed away from their capture population to forage, but returned the next season. Most Rails (195 of 204, or 95.6%) were classified as residents, but we detected two established migrants that had moved >100 km more than a year before capture. Seven Rails appeared to be seasonal migrants, but comparisons of feather isotope values with isotope values from wetland soils indicated that the isotope values in the feathers of these Rails likely resulted from natural environmental variation (e.g., source element effects) rather than long-distance dispersal of individuals. Thus, these seven Rails were most likely misassigned by isotopic population assignments due to small-scale variation in the isoscape. Using genetic data in conjunction with isotopic data allowed us to not only infer the timing of long-distance dispersal events, but to successfully track long-distance movements of nonmigratory Rails between metapopulations even when environmental variation of isotopes occurred across small spatial scales.

Marshes as "Mountain Tops": Genetic Analyses of the Critically Endangered São Paulo Marsh Antwren (Aves: Thamnophilidae)

Small populations of endangered species can be impacted by genetic processes such as drift and inbreeding that reduce population viability. As such, conservation genetic analyses that assess population levels of genetic variation and levels of gene flow can provide important information for managing threatened species. The São Paulo Marsh Antwren (Formicivora paludicola) is a recently-described and critically endangered bird from São Paulo State (Brazil) whose total estimated population is around 250-300 individuals, distributed in only 15 isolated marshes around São Paulo metropolitan region. We used microsatellite DNA markers to estimate the population genetic characteristics of the three largest remaining populations of this species all within 60 km of each other. We detected a high and significant genetic structure between all populations (overall FST = 0.103) which is comparable to the highest levels of differentiation ever documented for birds, (e.g., endangered birds found in isolated populations on the tops of African mountains), but also evidence for first-generation immigrants, likely from small local unsampled populations. Effective population sizes were small (between 28.8-99.9 individuals) yet there are high levels of genetic variability within populations and no evidence for inbreeding. Conservation implications of this work are that the high levels of genetic structure suggests that translocations between populations need to be carefully considered in light of possible local adaptation and that remaining populations of these birds should be managed as conservation units that contain both main populations studied here but also small outlying populations which may be a source of immigrants.

Impact of cattle grazing on the occupancy of a cryptic, threatened rail

Impacts of livestock grazing in arid and semiarid environments are often concentrated in and around wetlands where animals congregate for water, cooler temperatures, and green forage. We assessed the impacts of winter-spring (November-May) cattle grazing on marsh vegetation cover and occupancy of a highly secretive marsh bird that relies on dense vegetation cover, the California Black Rail (Laterallus jamaicensis coturniculus), in the northern Sierra Nevada foothills of California, U.S.A. Using detection-nondetection data collected during repeated call playback surveys at grazed vs. ungrazed marshes and a "random changes in occupancy" parameterization of a multi-season occupancy model, we examined relationships between occupancy and habitat covariates, while accounting for imperfect detection. Marsh vegetation cover was significantly lower at grazed marshes than at ungrazed marshes during the grazing season in 2007 but not in 2008. Winter-spring grazing had little effect on Black Rail occupancy at irrigated marshes. However, at nonirrigated marshes fed by natural springs and streams, grazed sites had lower occupancy than ungrazed sites. Black Rail occupancy was positively associated with marsh area, irrigation as a water source, and summer vegetation cover, and negatively associated with marsh isolation. Residual dry matter (RDM), a commonly used metric of grazing intensity, was significantly associated with summer marsh vegetation cover at grazed sites but not spring cover. Direct monitoring of marsh vegetation cover, particularly at natural spring- or stream-fed marshes, is recommended to prevent negative impacts to rails from overgrazing.