Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus

Determining evolutionary origin and phylogenetic relationships of mallard-like ducks of Oceania, greater Indonesia, and the Philippines with ddRAD-seq data

AIM

We aim to determine the evolutionary origins and population genetics of mallard-like ducks of Oceania, greater Indonesia, and the Philippines.

LOCATION

Oceania, greater Indonesia, and the Philippines.

TAXON

Mallard (Anas platyrhynchos), Pacific black duck (A. superciliosa spp.), and Philippine duck (A. luzonica) METHODS: Thousands of nuclear ddRAD-seq loci and the mitochondrial DNA control region were assayed across individuals representative of each species' range. We assessed population structure and phylogenetic relationships, as well as estimated demographic histories to reconstruct the biogeographical history of each species.

RESULTS

Philippine and Pacific black ducks represent unique genetic lineages that diverged from the mallard 1-2 million years ago. We find no support for the Philippine duck representing a hybrid species as once posited; however, their low levels of genetic diversity requires further attention. We find a lack of substructure among Philippine ducks. However, we found pronounced differentiation between subspecies of Pacific black ducks, especially between A. s. superciliosa from New Zealand and A. s. rogersi from Australia, Papua New Guinea, and Timor-Leste, Indonesia. Anas superciliosa pelewensis gave mixed results; individuals from the Solomon Islands were differentiated from the other subspecies, but those from the island of Aunu'u, American Samoa, were genetically more similar to A. s. rogersi than A. s. pelewensis samples from the Solomon Islands. Finally, we find limited evidence of interspecific gene flow at evolutionary scales, and mallard introgression among contemporary samples.

MAIN CONCLUSIONS

Mallard-like ducks radiated across Oceania, greater Indonesia, and the Philippines within the last 2 million years. Only the Pacific black duck showed unique sub-structuring that largely followed known sub-species ranges, except for A. s. pelewensis. We posit that the high interrelatedness among Solomon Island samples suggests that their genetic distinctiveness may simply be the result of high levels of genetic drift. In contrast, we conclude that mainland Australian Pacific black ducks were the most likely source for the recent colonization of American Samoa. As a result, our findings suggest that either the A. s. pelewensis subspecies designations and/or its geographical range may require re-evaluation. Continued re-evaluation of evolutionary and taxonomic relationships is necessary when attempting to reconstruct and understand biogeographical histories, with important implications towards any attempts to implement conservation strategies.

Where east meets west: Phylogeography of the high Arctic North American brant goose

Genetic variation in Arctic species is often influenced by vicariance during the Pleistocene, as ice sheets fragmented the landscape and displaced populations to low- and high-latitude refugia. The formation of secondary contact or suture zones during periods of ice sheet retraction has important consequences on genetic diversity by facilitating genetic connectivity between formerly isolated populations. Brant geese (Branta bernicla) are a maritime migratory waterfowl (Anseriformes) species that almost exclusively uses coastal habitats. Within North America, brant geese are characterized by two phenotypically distinct subspecies that utilize disjunct breeding and wintering areas in the northern Pacific and Atlantic. In the Western High Arctic of Canada, brant geese consist of individuals with an intermediate phenotype that are rarely observed nesting outside this region. We examined the genetic structure of brant geese populations from each subspecies and areas consisting of intermediate phenotypes using mitochondrial DNA (mtDNA) control region sequence data and microsatellite loci. We found a strong east-west partition in both marker types consistent with refugial populations. Within subspecies, structure was also observed at mtDNA while microsatellite data suggested the presence of only two distinct genetic clusters. The Western High Arctic (WHA) appears to be a secondary contact zone for both Atlantic and Pacific lineages as mtDNA and nuclear genotypes were assigned to both subspecies, and admixed individuals were observed in this region. The mtDNA sequence data outside WHA suggests no or very restricted intermixing between Atlantic and Pacific wintering populations which is consistent with published banding and telemetry data. Our study indicates that, although brant geese in the WHA are not a genetically distinct lineage, this region may act as a reservoir of genetic diversity and may be an area of high conservation value given the potential of low reproductive output in this species.

Numerous insertions of mitochondrial DNA in the genome of the northern mole vole, Ellobius talpinus

BACKGROUND

Ellobius talpinus is a subterranean rodent representing an attractive model in population ecology studies due to its highly special lifestyle and sociality. In such studies, mitochondrial DNA (mtDNA) is widely used. However, if nuclear copies of mtDNA, aka NUMTs, are present, they may co-amplify with the target mtDNA fragment, generating misleading results. The aim of this study was to determine whether NUMTs are present in E. talpinus.

METHODS AND RESULTS

PCR amplification of the putative mtDNA CytB-D-loop fragment using 'universal' primers from 56 E. talpinus samples produced multiple double peaks in 90% of the sequencing chromatograms. To reveal NUMTs, molecular cloning and sequencing of PCR products of three specimens was conducted, followed by phylogenetic analysis. The pseudogene nature of three out of the seven detected haplotypes was confirmed by their basal positions in relation to other Ellobius haplotypes in the phylogenetic tree. Additionally, 'haplotype B' was basal in relation to other E. talpinus haplotypes and found present in very distant sampling sites. BLASTN search revealed 195 NUMTs in the E. talpinus nuclear genome, including fragments of all four PCR amplified pseudogenes. Although the majority of the NUMTs studied were short, the entire mtDNA had copies in the nuclear genome. The most numerous NUMTs were found for rrnL, COXI, and D-loop.

CONCLUSIONS

Numerous NUMTs are present in E. talpinus and can be difficult to discriminate against mtDNA sequences. Thus, in future population or phylogenetic studies in E. talpinus, the possibility of cryptic NUMTs amplification should always be taken into account.

The meaning of wild: Genetic and adaptive consequences from large-scale releases of domestic mallards

The translocation of individuals around the world is leading to rising incidences of anthropogenic hybridization, particularly between domestic and wild congeners. We apply a landscape genomics approach for thousands of mallard (Anas platyrhynchos) samples across continental and island populations to determine the result of over a century of supplementation practices. We establish that a single domestic game-farm mallard breed is the source for contemporary release programs in Eurasia and North America, as well as for established feral populations in New Zealand and Hawaii. In particular, we identify central Europe and eastern North America as epicenters of ongoing anthropogenic hybridization, and conclude that the release of game-farm mallards continues to affect the genetic integrity of wild mallards. Conversely, self-sustaining feral populations in New Zealand and Hawaii not only show strong differentiation from their original stock, but also signatures of local adaptation occurring in less than a half-century since game-farm mallard releases have ceased. We conclude that 'wild' is not singular, and that even feral populations are capable of responding to natural processes. Although considered paradoxical to biological conservation, understanding the capacity for wildness among feral and feral admixed populations in human landscapes is critical as such interactions increase in the Anthropocene.

Frequency and types of alternative breeding strategies employed by nesting American black ducks in North Carolina

Although most birds are considered to be at least partially monogamous, molecular evidence continues to uncover that many species can have multiple sexual mates. Many species of Waterfowl (Order Anseriformes) consistently deploy alternative breeding strategies, and although cavity nesting species have been well studied, few attempts to understand rates of alternative breeding strategies exist in the Anatini tribe. Here, we assay mitochondrial DNA and thousands of nuclear markers across 20 broods of American black ducks (Anas rubripes; "black duck") that included 19 females and 172 offspring to study population structure as well as types and rates of secondary breeding strategies in coastal North Carolina. First, we report high levels of relatedness among nesting black ducks and offspring and while 17 (of 19) females were of pure black duck descent, three were found to be black duck x mallard (A. platyrhynchos) hybrids. Next, we evaluated for mismatched mitochondrial DNA and paternity identities across each female's clutch to determine types and frequency of alternative or secondary breeding strategies. Although we report that nest parasitism occurred in two nests, 37% (7 of 19) of the sampled nests were multi-paternal as a result of extra-pair copulation. In addition to being part of a mix of strategies used to increase fecundity by successfully breeding females, we posit nest densities providing easier alternative mate access for males also explains high rates of extra-pair copulation among our sampled black ducks. Ultimately, however, while some proportion of females of many species engage in forms of secondary breeding strategies, we conclude that the decision to do so appears to be seasonally flexible for each individual.

Interspecific forced copulations generate most hybrids in broadly sympatric ducks

Although rare, hybrids are more common in broadly sympatric waterfowl than in any other avian family; yet, the behavioral ecology explaining their generation has remained controversial. Leading hypotheses are forced interspecific copulations, mis-imprinting caused by mixed broods, and scarcity of conspecific mates. Using a large sample of hybrid ducks solicited from North American hunters we evaluated these hypotheses by genetically determining the mother and father species of F₁ hybrids. Based on abundances in areas where their breeding ranges overlap, the frequency of hybrids varied greatly from expectations, with hybrids between species within recently derived clades being much more frequent than those between more divergent clades. Forced copulations, as measured by large phallus-length asymmetries between parentals, strongly predicted the father species of most F₁ hybrids. Thus, most Anas acuta x A. platyrhynchos (Northern Pintail x Mallard) F₁s were sired by A. acuta, and most A. platyrhynchos x Mareca strepera (Mallard x Gadwall) F₁s were sired by A. platyrhynchos. Siring asymmetries were consistent with phallus length asymmetries in five additional parental combinations, but none had samples large enough to be individually statistically significant. The exception to this trend was our sample of nine A. platyrhynchos x Mareca americana (Mallard x Gadwall) F₁s, for which a large phallus asymmetry failed to predict the father species. Hybrids were rare in brood parasitic species, suggesting mis-imprinting to be an unlikely cause of most hybrids; however, our samples of hybrids from regular brood parasites were inadequate to strongly address this hypothesis. We could test the scarcity of mates hypothesis for only a single hybrid combination and it contradicted our prediction: most F₁M. Penelope x M. americana (Eurasian x American Wigeon) were sired by M. penelope, strongly contradicting our prediction that female M. penelope wintering in enormous flocks of M. americana (American Wigeon) on the west coast of North America would have difficulty finding conspecific mates. In general, our results support interspecific forced copulations as the predominant behavioral mechanism generating hybrids in North temperate waterfowl.

Genomic and morphological data shed light on the complexities of shared ancestry between closely related duck species

Causes for genomic and morphological similarities among recently radiated species are often multifaceted and are further convoluted among species that readily interbreed. Here, we couple genomic and morphological trait comparisons to test the extent that ancestry and gene flow explain the retention of mallard-like traits within a sister species, the Mexican duck. First, we confirm that these taxa remain genetically structured, and that Mexican ducks exhibit an isolation-by-distance pattern. Despite the assumption of wide-spread hybridization, we found only a few late-stage hybrids, all from the southwestern USA. Next, assessing 23 morphological traits, we developed a genetically-vetted morphological key that is > 97% accurate in distinguishing across sex-age cohorts of Mexican ducks, mallards, and hybrids. During key development, we determined that 25% of genetically pure, immature male Mexican ducks of the northern population naturally displayed mallard-like traits in their formative plumage. In fact, applying this key to 55 museum specimens, we identified that only four of the 14 specimens originally classified as phenotypic hybrids were truly hybrids. We discuss how genomic and morphological comparisons shed light into the mechanism(s) underlying the evolution of complex phenotypic traits in recent radiations, and how misunderstanding the true morphological diversity within Mexican ducks resulted in taxonomic revisions that hindered conservation efforts.

Nest parasitism, promiscuity, and relatedness among wood ducks

Nest parasitism is a common reproductive strategy used by many species of cavity nesting birds. Among these, the wood duck (Aix sponsa) is known to have evolved very specific strategies of when and whom to parasitize that is often based on population and/or environmental queues. Here, we investigated the genetic relationship of two female wood ducks competing over an artificial nesting box in Delaware, including the continued incubation of one female despite the death and body remains of the other female throughout the incubation process. We test whether such an extreme case of nest parasitism can be explained by relatedness, egg lineage composition, or a combination of other factors. To do so, we extracted genomic DNA from blood and tissue of the females, as well as chorioallantoic membranes of all viable and inviable eggs. Subsequently, we assessed relatedness among females and eggs based on hundreds of nuclear loci and the mitochondrial control region. We concluded that (1) the two incubating females were entirely unrelated, (2) the single clutch is in fact represented by a minimum of four unrelated females, and (3) a single female can lay eggs sired by different males. The latter finding is the first direct evidence for successful extra-pair copulation in wood ducks. With decreasing costs and increasing effectiveness, genomic methods have the potential to provide important insights into more complex ecological and evolutionary tactics of such populations.

Detection of heteroplasmy and nuclear mitochondrial pseudogenes in the Japanese spiny lobster Panulirus japonicus

Partial mtDNA cytochrome oxidase subunit I (COI) fragments and near entire stretch of 12S rDNA (12S) and control region (Dloop) of the Japanese spiny lobster (Panulirus japonicus) (n = 3) were amplified by PCR and used for direct nucleotide sequencing and for clone library-based nucleotide sequence analysis. Nucleotide sequences of a total of 75 clones in COI, 77 in 12S and 92 in Dloop were determined. Haplotypes of the clones matched with those obtained by direct sequencing were determined to be genuine mtDNA sequence of the individual. Phylogenetic analysis revealed several distinct groups of haplotypes in all three regions. Genuine mtDNA sequences were observed to form a group with their closely related variables, and most of these variables may be due to amplification error but a few to be heteroplasmy. Haplotypes determined as nuclear mitochondrial pseudogenes (NUMTs) formed distinct groups. Nucleotide sequence divergence (K2P distance) between genuine haplotypes and NUMTs were substantial (7.169-23.880% for COI, 1.336-23.434% for 12S, and 7.897-71.862% for Dloop). These values were comparable to or smaller than those between species of the genus Panulirus, indicating that integration of mtDNA into the nuclear genome is a continuous and dynamic process throughout pre- and post-speciation events. Double peaks in electropherograms obtained by direct nucleotide sequencing were attributed to common nucleotides shared by multiple NUMTs. Information on the heteroplasmy and NUMTs would be very important for addressing their impact on direct nucleotide sequencing and for quality control of nucleotide sequences obtained.

Sharp Increase of Problematic Mitogenomes of Birds: Causes, Consequences, and Remedies

Authentic DNA sequences are crucial for reliable evolutionary inference. Concerns about the identification of DNA sequences have been voiced several times in the past but few quantitative studies exist. Mitogenomes play important roles in phylogenetics, phylogeography, population genetics, and DNA identification. However, the large number of mitogenomes being published routinely, often in brief data papers, has raised questions about their authenticity. In this study, we quantify problematic mitogenomes of birds and their reusage in other papers. Of 1,876 complete or partial mitogenomes of birds published until January 1, 2020, the authenticity of 1,559 could be assessed with sequences of conspecifics. Of these, 78 (5.0%) were found to be problematic, including 45 curated reference sequences. Problems were due to misidentification (33), chimeras of two or three species (23), sequencing errors/numts (18), incorrect sequence assembly (1), mislabeling at GenBank but not in the final paper (2), or vice versa (1). The number of problematic mitogenomes has increased sharply since 2012. Worryingly, these problematic sequences have been reused 436 times in other papers, including 385 times in phylogenies. No less than 53% of all mitogenomic phylogenies/networks published until January 1, 2020 included at least one problematic mitogenome. Problematic mitogenomes have resulted in incorrect phylogenetic hypotheses and proposals for unwarranted taxonomic revision, and may have compromised comparative analyses and measurements of divergence times. Our results indicate that a major upgrade of quality control measures is warranted. We propose a comprehensive set of measures that may serve as a new standard for publishing mitogenome sequences.

Genetic structure and diversity in Neochen jubata (Aves: Anatidae) from the Araguaia River, GO, Brazil

Abstract: The Orinoco Goose (Neochen jubata) is a few-known and endemic Anatidae to South America, inhabiting sandy beaches along medium and large rivers, with a well-developed riparian forest and in swamp savannas and large freshwater baths. Recent data indicate the presence of longitudinal migratory behavior, and despite them, there are no records on the genetic profile of this species. The Araguaia River region, in the municipality of Luiz Alves, Goiás, receives an undetermined number of ducks seasonally, and there is little information about the individuals who visit this place, constituting the ideal scenario for a study able to offer a genetic overview perspective of this species and to understand the relationship between these individuals better. For this, we genetically characterized 61 individuals sampled in three distinct years of collection using microsatellite molecular markers and mitochondrial DNA. Genetic diversity analyses revealed low levels of heterozygosity for all sampled groups. However, they are within the equilibrium proposed by Hardy-Weinberg (HWE), as inbreeding or drift are not acting in these groups. The parentage analysis supports it, showing a high number of unrelated individuals over the years. AMOVA showed a significant difference among groups. These results may reflect the structure of this migratory species in that region, with the paired differentiation test of individuals from 2013 and 2014 being more similar to each other than those from other years, indicating a possible genetic structure diagnosed by the years of capture. However, there is a high allelic sharing among the three sampled groups, suggesting that these individuals are a population that connects over time and that they have a philopatric relationship with the location. The results found in this study constitute an initial milestone for the genetic knowledge of the mallard duck that should be raised in many other genetic studies.

Demographic history of the fragmented yellowthroated bulbul (Pycnonotus xantholaemus) population in the Deccan Peninsula, India

The yellow-throated bulbul (YTB) is an endemic passerine restricted to scrub forests along hill slopes with exposed rocky outcrops in the Deccan Peninsula, India. It is found in small, discontinuous populations and is vulnerable to extinction due to ongoing habitat loss and subsequent population decline. To assess the genetic connectivity and past demography, we sequenced 1050 nucleotide base pairs of the mitochondrial control region of 60 individuals that represent distinct populations in the geographic range of the species. We recovered 39 haplotypes defined by 81 variable sites. Haplotype diversity was high with low nucleotide diversity, suggesting rapid population growth from a founder population with a small effective population size. The negative values of Tajima’sDand Fu’sFsand small positive value of Ramos-Onsins and Rozas’R2suggest deviation from neutrality and population expansion. The haplotype network and demographic expansion parameters further suggest historical population expansion. Mismatch analysis statistics and Bayesian skyline plots estimate population expansion during the late Pleistocene. Although the species presently occurs in small, disconnected we found no structuring of the population. Dispersal events are the most likely explanation for the absence of genetic structuring in the YTB population. These results represent important data for the design of a conservation plan for this endemic and globally threatened species.

Biogeographical patterns and processes in the genus group Scotussae (Acrididae: Melanoplinae): an integrative approach

A biogeographical study of the genus group Scotussae, a clade of grasshoppers endemic to the subtropical temperate region of the La Plata Basin, South America, was performed within a phylogenetic context to test whether wing reduction reflects evolutionary and ecological processes within the clade. We used an integrative biogeographical approach to determine the role of geohistorical events, geography, ecology and phylogenetic niche conservatism on the distribution and diversification processes of the group. We performed a total evidence phylogenetic analysis and tested the phylogenetic signal of ecological niche traits (niche optimum and niche breadth). We also assessed the degree to which phylogenetic distance is correlated with geographical and ecological niche traits and we used BioGeoBEARS to estimate ancestral ranges. The results provided evidence for phylogenetic niche conservatism as well as a significant association between phylogeny and both geographical and, more strongly, ecological traits. Two main clades were clearly associated with wing development, and evidence points to the evolutionary and ecological processes within these two groups being different. The Brachypterous clade shows evidence that allopatric speciation was the main source of diversification, while for the Macropterous clade sympatric speciation seems more likely.

Convergent changes in muscle metabolism depend on duration of high-altitude ancestry across Andean waterfowl

High-altitude environments require that animals meet the metabolic O2 demands for locomotion and thermogenesis in O2-thin air, but the degree to which convergent metabolic changes have arisen across independent high-altitude lineages or the speed at which such changes arise is unclear. We examined seven high-altitude waterfowl that have inhabited the Andes (3812-4806 m elevation) over varying evolutionary time scales, to elucidate changes in biochemical pathways of energy metabolism in flight muscle relative to low-altitude sister taxa. Convergent changes across high-altitude taxa included increased hydroxyacyl-coA dehydrogenase and succinate dehydrogenase activities, decreased lactate dehydrogenase, pyruvate kinase, creatine kinase, and cytochrome c oxidase activities, and increased myoglobin content. ATP synthase activity increased in only the longest established high-altitude taxa, whereas hexokinase activity increased in only newly established taxa. Therefore, changes in pathways of lipid oxidation, glycolysis, and mitochondrial oxidative phosphorylation are common strategies to cope with high-altitude hypoxia, but some changes require longer evolutionary time to arise.

Assessing changes in genomic divergence following a century of human-mediated secondary contact among wild and captive-bred ducks

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human-mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD-seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game-farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game-farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait-capture array targeting thousands of loci in century-old (1842-1915) and contemporary (2009-2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game-farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game-farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens.

Holarctic phylogeographic structure of Eurasian wigeon (Mareca penelope)

The Eurasian wigeon (Mareca penelope) is one of the most numerous migrant species of waterfowl in the Palearctic. Annually, significant part of the world’s wigeon population makes seasonal flights over distances from tens to thousands or more kilometers. According to different estimates based on banding data, five geographic populations of the species were described in the Palearctic. However, distinct borders between the populations have not been identified. At the same time, no phylogeographic studies have been carried out for the complete native range of wigeon so far. In addition to the fundamental importance of such a study, knowledge of the genetic structure of populations is necessary for the development of measures to increase the number of and preserve this valuable game species. The aim of our work was a phylogeographic analysis of the wigeon across its vast native range in the Palearctic including ducks wintering in North America. We examined genetic diversity and differentiation of wigeon populations identified with banding data, phylogenetic relationships of mtDNA haplotypes and demographic history of populations and species as a whole by sequencing a 661 base-pair 5’-fragment of the mitochondrial control region from 195 individual ducks collected throughout the Palearctic and Nearctic. Genetic diversity was high in all studied populations. A reconstruction of haplotypes phylogeny revealed the absence of geographic structure in the data. Nonetheless, analysis of molecular variance (AMOVA) identified two groups of populations: EuropeanSiberian and East Asian. The former included wigeons from Europe, Siberia and the Atlantic coast of North America, and the latter comprised ducks from Russian Far East, Kamchatka Peninsula, Chukotka Autonomous District, the Aleutian Islands, Alaska, and the Pacific coast of North America.

Identifying hybrids & the genomics of hybridization: Mallards & American black ducks of Eastern North America

Resolving evolutionary relationships and establishing population structure depends on molecular diagnosability that is often limited for closely related taxa. Here, we use 3,200 ddRAD-seq loci across 290 mallards, American black ducks, and putative hybrids to establish population structure and estimate hybridization rates. We test between traditional assignment probability and accumulated recombination events based analyses to assign hybrids to generational classes. For hybrid identification, we report the distribution of recombination events complements ADMIXTURE simulation by extending resolution past F4 hybrid status; however, caution against hybrid assignment based on accumulated recombination events due to an inability to resolve F1 hybrids. Nevertheless, both analyses suggest that there are relatively few backcrossed stages before a lineage's hybrid ancestry is lost and the offspring are effectively parental again. We conclude that despite high rates of observed interspecific hybridization between mallards and black ducks in the middle part of the 20th century, our results do not support the predicted hybrid swarm. Conversely, we report that mallard samples genetically assigned to western and non-western clusters. We indicate that these non-western mallards likely originated from game-farm stock, suggesting landscape level gene flow between domestic and wild conspecifics.

The ectoparasites of hybrid ducks in New Zealand (Mallard x Grey Duck)

We studied the population genetics of one population sample of hybrid Mallard x Grey Ducks and their lice in New Zealand. We aimed to document the relationship between ectoparasite load and host phenotype, and test for an association between the mtDNA diversity of the lice and their hosts, which is predicted based on maternal care. We found three feather lice species previously described for these hosts: Anaticola crassicornis (wing louse), Anatoecus dentatus (head louse), and Trinoton querquedulae (body louse). No new or rare lice species were uncovered. Most ducks in our sample were more Mallard-like than Grey Duck-like hybrids for the five colour and plumage traits examined. We confirm that based solely on phenotypic characters it is difficult to distinguish between Mallards, hybrids and Grey Ducks. We detected no association between the number of lice and host phenotype for two of the three louse species (while controlling for bird size). However, the Grey Duck-like hybrids had fewer head lice (A. dentatus) than their Mallard-like counterparts. Only three of the 40 hosts had mtDNA haplotypes that characterise Grey Ducks. We present the first genetic data of Anaticola crassicornis, Anatoecus dentatus and Trinoton querquedulae from New Zealand waterfowl. We found that the lice mtDNA had greater sequence diversity than the homologous gene for the ducks. A mitochondrial phylogeny for A. crassicornis collected from hosts worldwide has been previously published, and we added our novel data to infer evolutionary relationships among worldwide populations of this louse. None of the three lice species showed a close association of parasite and host mtDNA lineage despite lack of paternal care in these duck species.

•

We studied population genetics of hybrid Mallard x Grey Ducks and their lice in New Zealand.

•

We found three feather lice species previously described for these hosts.

•

Most ducks were more Mallard-like than Grey Duck-like hybrids for the traits examined.

•

Only three of the forty hosts had mtDNA haplotypes that characterise Grey Ducks.

•

We present the first genetic data for the three lice species from New Zealand.

Flyway structure in the circumpolar greater white-fronted goose

Dispersal and migratory behavior are influential factors in determining how genetic diversity is distributed across the landscape. In migratory species, genetic structure can be promoted via several mechanisms including fidelity to distinct migratory routes. Particularly within North America, waterfowl management units have been delineated according to distinct longitudinal migratory flyways supported by banding data and other direct evidence. The greater white-fronted goose (Anser albifrons) is a migratory waterfowl species with a largely circumpolar distribution consisting of up to six subspecies roughly corresponding to phenotypic variation. We examined the rangewide population genetic structure of greater white-fronted geese using mtDNA control region sequence data and microsatellite loci from 23 locales across North America and Eurasia. We found significant differentiation in mtDNA between sampling locales with flyway delineation explaining a significant portion of the observed genetic variation (~12%). This is concordant with band recovery data which shows little interflyway or intercontinental movements. However, microsatellite loci revealed little genetic structure suggesting a panmictic population across most of the Arctic. As with many high-latitude species, Beringia appears to have played a role in the diversification of this species. A common Beringian origin of North America and Asian populations and a recent divergence could at least partly explain the general lack of structure at nuclear markers. Further, our results do not provide strong support for the various taxonomic proposals for this species except for supporting the distinctness of two isolated breeding populations within Cook Inlet, Alaska (A. a. elgasi) and Greenland (A. a. flavirostris), consistent with their subspecies status.

Influence of invasion history on rapid morphological divergence across island populations of an exotic bird

There is increasing evidence that exotic populations may rapidly differentiate from those in their native range and that differences also arise among populations within the exotic range. Using morphological and DNA-based analyses, we document the extent of trait divergence among native North American and exotic Hawaiian populations of northern cardinal (Cardinalis cardinalis). Furthermore, using a combination of historical records and DNA-based analyses, we evaluate the role of founder effects in producing observed trait differences. We measured and compared key morphological traits across northern cardinal populations in the native and exotic ranges to assess whether trait divergence across the Hawaiian Islands, where this species was introduced between 1929 and 1931, reflected observed variation across native phylogeographic clades in its native North America. We used and added to prior phylogenetic analyses based on a mitochondrial locus to identify the most likely native source clade(s) for the Hawaiian cardinal populations. We then used Approximate Bayesian Computation (ABC) to evaluate the role of founder effects in producing the observed differences in body size and bill morphology across native and exotic populations. We found cardinal populations on the Hawaiian Islands had morphological traits that diverged substantially across islands and overlapped the trait space of all measured native North American clades. The phylogeographic analysis identified the eastern North American clade (C. cardinalis cardinalis) as the most likely and sole native source for all the Hawaiian cardinal populations. The ABC analyses supported written accounts of the cardinal's introduction that indicate the original 300 cardinals shipped to Hawaii were simultaneously and evenly released across Hawaii, Kauai, and Oahu. Populations on each island likely experienced bottlenecks followed by expansion, with cardinals from the island of Hawaii eventually colonizing Maui unaided. Overall, our results suggest that founder effects had limited impact on morphological trait divergence of exotic cardinal populations in the Hawaiian archipelago, which instead reflect postintroduction events.

The little brown bat nuclear genome contains an entire mitochondrial genome: Real or artifact?

Nuclear mitochondrial DNA sequences (NUMTs) have been documented in almost all eukaryotic genomes studied. Recently, with the number of sequenced genomes increasing, extremely large NUMTs, even a nearly entire mitochondrial genome, have been reported in some plants and animals. However, few such studies provided strong experimental evidences for these important discoveries. In this study using a computer-based search method an entire mitochondrial genome (NUMT-1) was found in the nuclear genome of a bat species (Myotis lucifugus). This super-large NUMT shared a same scaffold with a 754bp nuclear genomic sequence and a second NUMT (NUMT-2, 3292bp). If NUMT-1 was real, it will be the largest NUMT found in animals and this finding will provide valuable insights into the mode of generation of NUMTs in the nuclear genome. Unfortunately, although the initial sequencing technology of the published M. lucifugus genome makes the possibility of artifact less likely, our results from both the PCR amplification followed by Sanger sequencing and mapping method based on the whole-genome resequencing datasets suggested that the scaffold containing the entire mitochondrial genome was artifact possibly due to a misassembly of mitochondrial and the nuclear DNA sequences. Our current study highlights the necessity to validate the authenticity of extremely large NUMTs identified in previous searches on whole-genome sequence assemblies.

Genetic structure among greater white-fronted goose populations of the Pacific Flyway

An understanding of the genetic structure of populations in the wild is essential for long‐term conservation and stewardship in the face of environmental change. Knowledge of the present‐day distribution of genetic lineages (phylogeography) of a species is especially important for organisms that are exploited or utilize habitats that may be jeopardized by human intervention, including climate change. Here, we describe mitochondrial (mtDNA) and nuclear genetic (microsatellite) diversity among three populations of a migratory bird, the greater white‐fronted goose (Anser albifrons), which breeds discontinuously in western and southwestern Alaska and winters in the Pacific Flyway of North America. Significant genetic structure was evident at both marker types. All three populations were differentiated for mtDNA, whereas microsatellite analysis only differentiated geese from the Cook Inlet Basin. In sexual reproducing species, nonrandom mate selection, when occurring in concert with fine‐scale resource partitioning, can lead to phenotypic and genetic divergence as we observed in our study. If mate selection does not occur at the time of reproduction, which is not uncommon in long‐lived organisms, then mechanisms influencing the true availability of potential mates may be obscured, and the degree of genetic and phenotypic diversity may appear incongruous with presumed patterns of gene flow. Previous investigations revealed population‐specific behavioral, temporal, and spatial mechanisms that likely influence the amount of gene flow measured among greater white‐fronted goose populations. The degree of observed genetic structuring aligns well with our current understanding of population differences pertaining to seasonal movements, social structure, pairing behavior, and resource partitioning.

Incipient speciation with gene flow on a continental island: Species delimitation of the Hainan Hwamei (Leucodioptron canorum owstoni, Passeriformes, Aves)

Because of their isolation, continental islands (e.g., Madagascar) are often thought of as ideal systems to study allopatric speciation. However, many such islands have been connected intermittently to their neighboring continent during recent periods of glaciation, which may cause frequent contact between the diverging populations on the island and continent. As a result, the speciation processes on continental islands may not meet the prerequisites for strictly allopatric speciation. We used multiple lines of evidence to re-evaluate the taxonomic status of the Hainan Hwamei (Leucodioptron canorum owstoni), which is endemic to Hainan, the largest continental island in the South China Sea. Our analysis of mitochondrial DNA and twelve nuclear loci suggests that the Hainan Hwamei can be regarded as an independent species (L. owstoni); the morphological traits of the Hainan Hwamei also showed significant divergence from those of their mainland sister taxon, the Chinese Hwamei (L. canorum). We also inferred the divergence history of the Hainan and Chinese Hwamei to see whether their divergence was consistent with a strictly allopatric model. Our results suggest that the two Hwameis split only 0.2 million years ago with limited asymmetrical post-divergence gene flow. This implies that the Hainan Hwamei is an incipient species and that speciation occurred through ecologically divergent selection and/or assortative mating rather than a strictly allopatric process.

Genetic Diversity and Phylogenetic Analysis of South-East Asian Duck Populations Based on the mtDNA D-loop Sequences

The maternally inherited mitochondrial DNA (mtDNA) D–loop region is widely used for exploring genetic relationships and for investigating the origin of various animal species. Currently, domestic ducks play an important role in animal protein supply. In this study, partial mtDNA D–loop sequences were obtained from 145 samples belonging to six South-East Asian duck populations and commercial duck population. All these populations were closely related to the mallard duck (Anas platyrhynchos), as indicated by their mean overall genetic distance. Sixteen nucleotide substitutions were identified in sequence analyses allowing the distinction of 28 haplotypes. Around 42.76% of the duck sequences were classified as Hap_02, which completely matched with Anas platyrhynchos duck species. The neighbor-joining phylogenetic tree also revealed that South-East Asian duck populations were closely related to Anas platyrhynchos. Network profiles were also traced using the 28 haplotypes. Overall, results showed that those duck populations D-loop haplotypes were shared between several duck breeds from Korea and Bangladesh sub continental regions. Therefore, these results confirmed that South-East Asian domestic duck populations have been domesticated from Anas platyrhynchos duck as the maternal origins.

Heteroplasmy and ancient translocation of mitochondrial DNA to the nucleus in the Chinese Horseshoe Bat (Rhinolophus sinicus) complex

The utility and reliability of mitochondrial DNA sequences in phylogenetic and phylogeographic studies may be compromised by widespread and undetected nuclear mitochondrial copies (numts) as well as heteroplasmy within individuals. Both numts and heteroplasmy are likely to be common across diverse taxa yet few studies have characterised their frequencies and variation at the intra-specific level. Here we report the presence of both numts and heteroplasmy in the mitochondrial control region of the Chinese horseshoe bat Rhinolophus sinicus. In total we generated 123 sequences from 18 bats, which contained two different numt clades (i.e. Numt-1 and Numt-2) and one mtDNA clade. The sequence divergence between Numt-1 and Numt-2 was 16.8% and each numt type was found in all four R. sinicus taxa, suggesting either two ancient translocations of mitochondrial DNA into the nucleus from the same source taxon, or a single translocation from different source taxa that occurred before the split of R. sinicus into different lineages. Within the mtDNA clade, phylogenetic relationships among the four taxa of R. sinicus were similar to those seen in previous results. Based on PCR comparisons, heteroplasmy was inferred between almost all individuals of R. sinicus with respect to sequence variation. Consistent with introgression of mtDNA between Central sinicus and septentrionalis, individuals from these two taxa exhibited similar signatures of repeated sequences in the control region. Our study highlights the importance of testing for the presence of numts and heteroplasmy when applying mtDNA markers to phylogenetic studies.

Association of 5-methylcytosine and 5-hydroxymethylcytosine with mitochondrial DNA content and clinical and biochemical parameters in hepatocellular carcinoma

Increasing epidemiological evidence has indicated that inherited variations of mitochondrial DNA (mtDNA) copy number affect the genetic susceptibility of many malignancies in a tumour-specific manner and that DNA methylation also plays an important role in controlling gene expression during the differentiation and development of hepatocellular carcinoma (HCC). Our previous study demonstrated that HCC tissues showed a lower 5-hydroxymethylcytosine (5-hmC) content when compared to tumour-adjacent tissues, but the relationship among 5-hmC, 5-methylcytosine (5-mC) and mtDNA content in HCC patients is still unknown. This study aimed to clarify the correlation among mtDNA content, 5-mC and 5-hmC by quantitative real-time PCR and liquid chromatography tandem mass spectrometry analysis. We demonstrated that 5-hmC correlated with tumour size [odds ratio (OR) 0.847, 95% confidence interval (CI) 0.746–0.962, P = 0.011], and HCC patients with a tumour size ≥5.0 cm showed a lower 5-hmC content and higher levels of fasting plasma aspartate aminotransferase, the ratio of alanine amiotransferase to aspartate aminotransferase, γ-glutamyltransferase, alpha-fetoprotein than those with a tumour size <5 cm (all P<0.05). We further revealed that the mtDNA content of HCC tumour tissues was 225.97(105.42, 430.54) [median (25th Percentile, 75th Percentile)] and was negatively correlated with 5-mC content (P = 0.035), but not 5-hmC content, in genomic DNA from HCC tumour tissues.

Computation vs. cloning: evaluation of two methods for haplotype determination

Nuclear sequence data, often from multiple loci, are increasingly being employed in analyses of population structure and history, yet there has been relatively little evaluation of methods for accurately and efficiently separating the alleles or haplotypes in heterozygous individuals. We compared the performance of a computational method of haplotype reconstruction and standard cloning methods using a highly variable intron (ornithine decarboxylase, intron 6) in three closely related species of dabbling ducks (genus Anas). Cloned sequences from 32 individuals were compared to results obtained from phase 2.1.1 . phase correctly identified haplotypes in 28 of 30 heterozygous individuals when the underlying model assumed no recombination. Haplotypes of the remaining two individuals were also inferred correctly except for unique polymorphisms, the phase of which was appropriately indicated as uncertain (phase probability = 0.5). For a larger set of 232 individuals, results were essentially identical regardless of the recombination model used and haplotypes for all 30 of the tested heterozygotes were correctly inferred, with the exception of uncertain phase for unique polymorphisms in one individual. In contrast, initial sequences of one clone per sample yielded accurate haplotype determination in only 26 of 30 individuals; polymerase chain reaction (PCR)/cloning errors resulting from misincorporation of individual nucleotides could be recognized and avoided by comparison to direct sequences, but errors due to PCR recombination resulted in incorrect haplotype reconstruction in four individuals. The accuracy of haplotypes reconstructed by phase, even when dealing with a relatively small number of samples and numerous variable sites, suggests broad utility of computational approaches for reducing the cost and improving the efficiency of data collection from nuclear sequence loci.

Range-wide genetic population structure of common pochard (Aythya ferina): a potentially important vector of highly pathogenic avian influenza viruses

An understanding of the distribution and spatial structure of the natural vectors of zoonothic pathogens is of interest for effective disease control and prevention. Here, we investigate the range-wide population genetic structure of common pochard (Aythya ferina), a long-distance migratory duck and potential vector of highly pathogenic avian influenza. We collected several hundred samples from breeding and wintering grounds across Eurasia including some H5N1-positive individuals and generated partial sequences of the mitochondrial control region and multilocus microsatellite genotypes. Genetic differentiation among breeding populations was significant for both marker types but higher for maternally inherited mtDNA than for biparentally inherited nuclear markers. There was only weak genetic divergence between ducks sampled in Europe and East Asia, and genetic differentiation between populations was not generally associated with geographical distance. No evidence of genetic substructure was detected for ducks sampled on the European wintering grounds. Our results suggest limited breeding-site fidelity, especially in females, but extensive population admixture on the wintering grounds. The specific role of pochards as natural vectors of zoonotic pathogens and in particular H5N1 remains to be clarified but our results point to wintering grounds as potential hotspots for disease transmission.

Genetic and phenotypic divergence between low- and high-altitude populations of two recently diverged cinnamon teal subspecies

Spatial variation in the environment can lead to divergent selection between populations occupying different parts of a species' range, and ultimately lead to population divergence. The colonization of new areas can thus facilitate divergence in beneficial traits, yet with little differentiation at neutral genetic markers. We investigated genetic and phenotypic patterns of divergence between low- and high-altitude populations of cinnamon teal inhabiting normoxic and hypoxic regions in the Andes and adjacent lowlands of South America. Cinnamon teal showed strong divergence in body size (PC1; P(ST) = 0.56) and exhibited significant frequency differences in a single nonsynonymous α-hemoglobin amino acid polymorphism (Asn/Ser-α9; F(ST) = 0.60) between environmental extremes, despite considerable admixture of mtDNA and intron loci (F(ST) = 0.004-0.168). Inferences of strong population segregation were further supported by the observation of few mismatched individuals in either environmental extreme. Coalescent analyses indicated that the highlands were most likely colonized from lowland regions but following divergence, gene flow has been asymmetric from the highlands into the lowlands. Multiple selection pressures associated with high-altitude habitats, including cold and hypoxia, have likely shaped morphological and genetic divergence within South American cinnamon teal populations.

Breeding site fidelity and winter admixture in a long-distance migrant, the tufted duck (Aythya fuligula)

Long-distance migrants are, by definition, highly mobile but it is poorly understood if this leads to high rates of gene flow and an essentially panmictic global population structure. Genetic divergence in migratory species could be promoted, for example, by fidelity to distinct migratory pathways. In this study, we investigate the population genetic structure of tufted duck (Aythya fuligula), a long-distance migrant with a largely continuous breeding distribution across Eurasia. Distinct, longitudinally oriented flyways have been postulated based on geographically disjunct wintering areas and are supported by evidence from ringing data. We generated sequences of the mitochondrial control region and multi-locus microsatellite genotypes for several hundreds of samples from the European and Asian breeding and wintering grounds including some individuals infected with highly pathogenic avian influenza virus H5N1. Significant differentiation between breeding sites was observed for both marker types, but F(ST) values were approximately 10 times higher for maternally inherited mitochondrial DNA than for biparentally transmitted nuclear markers. The genetic differentiation between the postulated European and Asian flyways was similar to that observed within continents and, in general, genetic divergence was not associated with geographic distance. Neither marker type showed evidence of genetic substructure among aggregations on the European wintering grounds. Our results suggest some breeding site fidelity, especially in females, but extensive population admixture on the wintering grounds. Several scenarios may explain the observed lack of genetic divergence between Europe and Asia including non-equilibrium conditions following a recent range expansion or contemporary gene flow across the postulated migratory divides.

Ancient DNA reveals genetic stability despite demographic decline: 3,000 years of population history in the endemic Hawaiian petrel

In the Hawaiian Islands, human colonization, which began approximately 1,200 to 800 years ago, marks the beginning of a period in which nearly 75% of the endemic avifauna became extinct and the population size and range of many additional species declined. It remains unclear why some species persisted whereas others did not. The endemic Hawaiian petrel (Pterodroma sandwichensis) has escaped extinction, but colonies on two islands have been extirpated and populations on remaining islands have contracted. We obtained mitochondrial DNA sequences from 100 subfossil bones, 28 museum specimens, and 289 modern samples to investigate patterns of gene flow and temporal changes in the genetic diversity of this endangered species over the last 3,000 years, as Polynesians and then Europeans colonized the Hawaiian Islands. Genetic differentiation was found to be high between both modern and ancient petrel populations. However, gene flow was substantial between the extirpated colonies on Oahu and Molokai and modern birds from the island of Lanai. No significant reductions in genetic diversity occurred over this period, despite fears in the mid-1900s that this species may have been extinct. Simulations show that even a decline to a stable effective population size of 100 individuals would result in the loss of only 5% of the expected heterozygosity. Simulations also show that high levels of genetic diversity may be retained due to the long generation time of this species. Such decoupling between population size and genetic diversity in long-lived species can have important conservation implications. It appears that a pattern of dispersal from declining colonies, in addition to long generation time, may have allowed the Hawaiian petrel to escape a severe genetic bottleneck, and the associated extinction vortex, and persist despite a large population decline after human colonization.