Simple identification of mitochondrial lineages in contact zones based on lineage-selective primers

A variety of research projects focus on genetic variation among and within maternal lineages as encompassed by mitochondrial DNA (mtDNA). While mtDNA often differs substantially between species, large differences may also be found within species. The evaluation of such divergent lineages, for example in intraspecific contact zones (hybrid zones), commonly involves sequencing numerous individuals. Large-scale sequencing is both expensive and labour-intensive. Based on sequences from 15 individuals, we devised a simple and quick polymerase chain reaction assay for identification of divergent mtDNA lineages in a secondary contact zone of the side-blotched lizard (Uta stansburiana). The application uses lineage-selective primers to amplify a lineage-diagnostic product, and is based on each group of mtDNA haplotypes being a monophyletic assemblage of haplotypes sharing the same maternal ancestry, deeply divergent from the other group. The assay was tested on a larger sample (n = 147) of specimens from the contact zone, confirming its usefulness in quick and reliable identification of mtDNA lineages. This approach can be modified for other species, provided diagnostic lineage variation is available, and may also be performed in simple laboratory settings while conducting fieldwork.

Multiple interkingdom horizontal gene transfers in Pyrenophora and closely related species and their contributions to phytopathogenic lifestyles

Many studies have reported horizontal gene transfer (HGT) events from eukaryotes, especially fungi. However, only a few investigations summarized multiple interkingdom HGTs involving important phytopathogenic species of Pyrenophora and few have investigated the genetic contributions of HGTs to fungi. We investigated HGT events in P. teres and P. tritici-repentis and discovered that both species harbored 14 HGT genes derived from bacteria and plants, including 12 HGT genes that occurred in both species. One gene coding a leucine-rich repeat protein was present in both species of Pyrenophora and it may have been transferred from a host plant. The transfer of genes from a host plant to pathogenic fungi has been reported rarely and we discovered the first evidence for this transfer in phytopathogenic Pyrenophora. Two HGTs in Pyrenophora underwent subsequent duplications. Some HGT genes had homologs in a few other fungi, indicating relatively ancient transfer events. Functional analyses indicated that half of the HGT genes encoded extracellular proteins and these may have facilitated the infection of plants by Pyrenophora via interference with plant defense-response and the degradation of plant cell walls. Some other HGT genes appeared to participate in carbohydrate metabolism. Together, these functions implied that HGTs may have led to highly efficient mechanisms of infection as well as the utilization of host carbohydrates. Evolutionary analyses indicated that HGT genes experienced amelioration, purifying selection, and accelerated evolution. These appeared to constitute adaptations to the background genome of the recipient. The discovery of multiple interkingdom HGTs in Pyrenophora, their significance to infection, and their adaptive evolution, provided valuable insights into the evolutionary significance of interkingdom HGTs from multiple donors.

Origin of Chinese goldfish and sequential loss of genetic diversity accompanies new breeds

BACKGROUND

Goldfish, Carassius auratus, have experienced strong anthropogenic selection during their evolutionary history, generating a tremendous extent of morphological variation relative to that in native Carassius. To locate the geographic origin of goldfish, we analyzed nucleotide sequences from part of the control region (CR) and the entire cytochrome b (Cytb) mitochondrial DNA genes for 234 goldfish and a large series of native specimens. Four important morphological characteristics used in goldfish taxonomy-body shape, dorsal fin, eye shape, and tailfin-were selected for hypothesis-testing to identify those that better correspond to evolutionary history.

PRINCIPAL FINDING

Haplotypes of goldfish rooted in two sublineages (C5 and C6), which contained the haplotypes of native C. a. auratus from southern China. Values of F(ST) and N(m) revealed a close relationship between goldfish and native C. a. auratus from the lower Yangtze River. An extraordinary, stepwise loss of genetic diversity was detected from native fish to goldfish and from Grass-goldfish relative to other breeds. Significantly negative results for the tests of Tajima's D and Fu and Li's D* and F* were identified in goldfish, including the Grass breed. The results identified eye-shape as being the least informative character for grouping goldfish with respect to their evolutionary history. Fisher's exact test identified matrilineal constraints on domestication.

CONCLUSIONS

Chinese goldfish have a matrilineal origin from native southern Chinese C. a. auratus, especially the lineages from the lower Yangtze River. Anthropogenic selection of the native Carassius eliminated aesthetically unappealing goldfish and this action appeared to be responsible for the stepwise decrease in genetic diversity of domesticated goldfish, a process similar to that reported for the domestication of pigs, rice, and maize. The three-breed taxonomy--Grass-goldfish, Egg-goldfish, and Wen-goldfish--better reflected the history of domestication.

Landscape genetics of leaf-toed geckos in the tropical dry forest of northern Mexico

Habitat fragmentation due to both natural and anthropogenic forces continues to threaten the evolution and maintenance of biological diversity. This is of particular concern in tropical regions that are experiencing elevated rates of habitat loss. Although less well-studied than tropical rain forests, tropical dry forests (TDF) contain an enormous diversity of species and continue to be threatened by anthropogenic activities including grazing and agriculture. However, little is known about the processes that shape genetic connectivity in species inhabiting TDF ecosystems. We adopt a landscape genetic approach to understanding functional connectivity for leaf-toed geckos (Phyllodactylus tuberculosus) at multiple sites near the northernmost limit of this ecosystem at Alamos, Sonora, Mexico. Traditional analyses of population genetics are combined with multivariate GIS-based landscape analyses to test hypotheses on the potential drivers of spatial genetic variation. Moderate levels of within-population diversity and substantial levels of population differentiation are revealed by FST and Dest. Analyses using structure suggest the occurrence of from 2 to 9 genetic clusters depending on the model used. Landscape genetic analysis suggests that forest cover, stream connectivity, undisturbed habitat, slope, and minimum temperature of the coldest period explain more genetic variation than do simple Euclidean distances. Additional landscape genetic studies throughout TDF habitat are required to understand species-specific responses to landscape and climate change and to identify common drivers. We urge researchers interested in using multivariate distance methods to test for, and report, significant correlations among predictor matrices that can impact results, particularly when adopting least-cost path approaches. Further investigation into the use of information theoretic approaches for model selection is also warranted.

Ecological genetics of Chinese rhesus macaque in response to mountain building: all things are not equal

BACKGROUND

Pliocene uplifting of the Qinghai-Tibetan Plateau (QTP) and Quaternary glaciation may have impacted the Asian biota more than any other events. Little is documented with respect to how the geological and climatological events influenced speciation as well as spatial and genetic structuring, especially in vertebrate endotherms. Macaca mulatta is the most widely distributed non-human primate. It may be the most suitable model to test hypotheses regarding the genetic consequences of orogenesis on an endotherm.

METHODOLOGY AND PRINCIPAL FINDINGS

Using a large dataset of maternally inherited mitochondrial DNA gene sequences and nuclear microsatellite DNA data, we discovered two maternal super-haplogroups exist, one in western China and the other in eastern China. M. mulatta formed around 2.31 Ma (1.51-3.15, 95%), and divergence of the two major matrilines was estimated at 1.15 Ma (0.78-1.55, 95%). The western super-haplogroup exhibits significant geographic structure. In contrast, the eastern super-haplogroup has far greater haplotypic variability with little structure based on analyses of six variable microsatellite loci using Structure and Geneland. Analysis using Migrate detected greater gene flow from WEST to EAST than vice versa. We did not detect signals of bottlenecking in most populations.

CONCLUSIONS

Analyses of the nuclear and mitochondrial datasets obtained large differences in genetic patterns for M. mulatta. The difference likely reflects inheritance mechanisms of the maternally inherited mtDNA genome versus nuclear biparentally inherited STRs and male-mediated gene flow. Dramatic environmental changes may be responsible for shaping the matrilineal history of macaques. The timing of events, the formation of M. mulatta, and the divergence of the super-haplogroups, corresponds to both the uplifting of the QTP and Quaternary climatic oscillations. Orogenesis likely drove divergence of western populations in China, and Pleistocene glaciations are likely responsible for genetic structuring in the eastern super-haplogroup via geographic isolation and secondary contact.

Multiple ancient horizontal gene transfers and duplications in lepidopteran species

Eukaryotic horizontal gene transfer (HGT) events are increasingly being discovered yet few reports have summarized multiple occurrences in a wide range of species. We systematically investigated HGT events in the order Lepidoptera by employing a series of filters. Bombyx mori, Danaus plexippus and Heliconius melpomene had 13, 12 and 12 HGTs, respectively, from bacteria and fungi. These HGTs contributed a total of 64 predicted genes: 22 to B. mori, 22 to D. plexippus and 20 to H. melpomene. Several new genes were generated by post-transfer duplications. Post-transfer duplication of a suite of functional HGTs has rarely been reported in higher organisms. The distributional patterns of paralogues for certain genes differed in the three species, indicating potential independent duplication or loss events. All of these HGTs had homologues expressed in some other lepidopterans, indicating ancient transfer events. Most HGTs were involved in the metabolism of sugar and amino acids. These HGTs appeared to have experienced amelioration, purifying selection and accelerated evolution to adapt to the background genome of the recipient. The discovery of ancient, massive HGTs and duplications in lepidopterans and their adaptive evolution provides further insights into the evolutionary significance of the events from donors to multicellular host recipients.

Evolution and expression plasticity of opsin genes in a fig pollinator, Ceratosolen solmsi

Figs and fig pollinators have co-evolved species-specific systems of mutualism. So far, it was unknown how visual opsin genes of pollinators have evolved in the light conditions inside their host figs. We cloned intact full-length mRNA sequences of four opsin genes from a species of fig pollinator, Ceratosolen solmsi, and tested for selective pressure and expressional plasticity of these genes. Molecular evolutionary analysis indicated that the four opsin genes evolved under different selective constraints. Subsets of codons in the two long wavelength sensitive opsin (LW1, LW2) genes were positively selected in ancestral fig pollinators. The ultraviolet sensitive opsin (UV) gene was under strong purifying selection, whereas a relaxation of selective constrains occurred on several amino acids in the blue opsin. RT-qPCR analysis suggested that female and male fig pollinators had different expression patterns possibly due to their distinct lifestyles and different responses to light within the syconia. Co-evolutionary history with figs might have influenced the evolution and expression plasticity of opsin genes in fig pollinators.

Taxonomic position of Eothenomys wardi (Arvicolinae: Cricetidae) based on morphological and molecular analyses with a detailed description of the species

Ward's Red-backed Vole (Eothenomys wardi) is a rodent from the family Cricetidae. This endemic species occurs only in extreme northwestern Yunnan province, China in the Mekong and Salween river divide. It occupies steep cliffs at 2,800 to 4,250 m above sea level on the remote Qinghai-Tibetan plateau. The validity of E. wardi is controversial and no specimens exist apart from the nominal series. In 2010, we collected 38 topotypes of E. wardi from Meri Snow Mountain. The results of our phylogenetic analyses based on nucleotide sequences of the mitochondrial genes cytochrome b (cytb) and cytochrome c oxidase subunit one (COI) suggest that E. wardi is the sister group of E. custos, against its previously presumed sister species or conspecific species E. chinensis. In addition, seven out of 34 morphological characters differentiate E. wardi from other members of the genus Eothenomys. Therefore, we consider E. wardi to be a valid species and we provide its detailed morphological description.

Genealogy and Demographic History of a Widespread Amphibian throughout Indochina

Relatively little is known about spatial patterns of cryptic diversity in tropical species and the processes that generate them. Few studies examine the geographic distribution of genetic lineages in Southeast Asia, an area hypothesized to harbor substantial cryptic diversity. We investigated the evolutionary history of Asian tree frogs of the Polypedates leucomystax complex (n = 172) based on 1800 bp of the mtDNA genes ND1 and cytochrome b and tested hypotheses pertaining to climate, geology, and dispersal patterns. Analyses revealed substantial genetic diversity and lineage divergence throughout the region with evidence for widespread sympatric lineages and a general north versus south clustering. Relaxed molecular clock analysis and tests for demographic expansion identified an initial cladogenesis during the Miocene with subsequent Plio-Pleistocene diversification, with the former corresponding to periods of increased aridity and the onset of monsoonal weather systems. Rates of diversification were relatively constant until the Early Pleistocene when rates increased exponentially. We found equivocal evidence for isolation-by-distance and a potential role of some landscape features as partial barriers to dispersal. Finally, our analyses showed that divergence between insular and mainland populations occurred before Homo sapiens colonized Southeast Asia, suggesting that historical human-mediated dispersal did not drive insular diversification. Our results suggested that demographic expansion in the Late Pleistocene resulted in widespread sympatric lineages in the P. leucomystax complex throughout southern China and Indochina and further clarified the evolutionary history of lineages within P. leucomystax.

River islands, refugia and genetic structuring in the endemic brown frog Rana kukunoris (Anura, Ranidae) of the Qinghai-Tibetan Plateau

Frequently, Pleistocene climatic cycling has been found to be the diver of genetic structuring in populations, even in areas that did not have continental ice sheets, such as on the Qinghai-Tibetan Plateau (QTP). Typically, species distributed on the plateau have been hypothesized to re-treat to south-eastern refugia, especially during the Last Glacial Maximum (LGM). We evaluated sequence variation in the mitochondrial DNA gene Cytb and the nuclear DNA gene RAG-1 in Rana kukunoris, a species endemic to the QTP. Two major lineages, N and S, were identified, and lineage N was further subdivided into N1 and N2. The geographical distribution and genealogical divergences supported the hypothesis of multiple refugia. However, major lineages and sublineages diverged prior to the LGM. Demographical expansion was detected only in lineage S and sublineage N2. Sublineage N1 might have survived several glacial cycles in situ and did not expand after the LGM because of the absence of suitable habitat; it survived in river islands. Genetic analysis and environment modelling suggested that the north-eastern edge of QTP contained a major refugium for R. kukunoris. From here, lineage S dispersed southwards after the LGM. Two microrefugia in northern Qilian Mountains greatly contributed to current level of intraspecific genetic diversity. These results were found to have important implications for the habitat conservation in Northwest China.

Clonal variation in the parthenogenetic rock lizard Lacerta armeniaca

Genetic diversity at 35 allozyme loci was surveyed in seven populations of Lacerta armeniaca. Fixed heterozygotes were present at 16 loci, with homozygotes at 17 loci. Variation occurred at two loci, one in each of two populations, indicating one widespread clone, one restricted clone, and one apparently restricted clone. The low level of variation in this species suggests a recent restricted origin, involving few parental individuals.

Complete genome sequence of a monosense densovirus infecting the cotton bollworm, Helicoverpa armigera

Densoviruses (DNVs) infecting arthropods are members of the family Parvoviridae. Here we report the complete genome sequence of a novel DNV with a monosense genome that infects cotton bollworms (Helicoverpa armigera), named HaDNV-1. Alignment and phylogenetic analysis revealed that HaDNV-1 showed high identity with the genus Iteravirus.

Tissue sampling methods and standards for vertebrate genomics

The recent rise in speed and efficiency of new sequencing technologies have facilitated high-throughput sequencing, assembly and analyses of genomes, advancing ongoing efforts to analyze genetic sequences across major vertebrate groups. Standardized procedures in acquiring high quality DNA and RNA and establishing cell lines from target species will facilitate these initiatives. We provide a legal and methodological guide according to four standards of acquiring and storing tissue for the Genome 10K Project and similar initiatives as follows: four-star (banked tissue/cell cultures, RNA from multiple types of tissue for transcriptomes, and sufficient flash-frozen tissue for 1 mg of DNA, all from a single individual); three-star (RNA as above and frozen tissue for 1 mg of DNA); two-star (frozen tissue for at least 700 μg of DNA); and one-star (ethanol-preserved tissue for 700 μg of DNA or less of mixed quality). At a minimum, all tissues collected for the Genome 10K and other genomic projects should consider each species' natural history and follow institutional and legal requirements. Associated documentation should detail as much information as possible about provenance to ensure representative sampling and subsequent sequencing. Hopefully, the procedures outlined here will not only encourage success in the Genome 10K Project but also inspire the adaptation of standards by other genomic projects, including those involving other biota.