Small Number Of Loci Research Articles

Genome divergence and diversification within a geographic mosaic of coevolution.

Despite substantial interest in coevolution's role in diversification, examples of coevolution contributing to speciation have been elusive. Here, we build upon past studies that have shown both coevolution between South Hills crossbills and lodgepole pine (Pinus contorta), and high levels of reproductive isolation between South Hills crossbills and other ecotypes in the North American red crossbill (Loxia curvirostra) complex. We used genotyping by sequencing to generate population genomic data and applied phylogenetic and population genetic analyses to characterize the genetic structure within and among nine of the ecotypes. Although genome-wide divergence was slight between ecotypes (FST =0.011-0.035), we found evidence of relative genetic differentiation (as measured by FST ) between and genetic cohesiveness within many of them. As expected for nomadic and opportunistic breeders, we detected no evidence of isolation by distance. The one sedentary ecotype, the South Hills crossbill, was genetically most distinct because of elevated divergence at a small number of loci rather than pronounced overall genome-wide divergence. These findings suggest that mechanisms related to recent local coevolution between South Hills crossbills and lodgepole pine (e.g. strong resource-based density dependence limiting gene flow) have been associated with genome divergence in the face of gene flow. Our results further characterize a striking example of coevolution driving speciation within perhaps as little as 6000years.

DNA (De)Methylation: The Passive Route to Naïvety?

Mouse pluripotent embryonic stem (ES) cells can exist in distinct yet interchangeable epigenetic states dictated by their culture environment. Previous reports have shown that naïve pluripotent cells grown in the presence of 2i are characterised by global DNA hypomethylation and changes in the abundance and distribution of histone modifications. New research provides insights regarding how this might be achieved.

Geographic cline analysis as a tool for studying genome-wide variation: a case study of pollinator-mediated divergence in a monkeyflower.

A major goal of speciation research is to reveal the genomic signatures that accompany the speciation process. Genome scans are routinely used to explore genome-wide variation and identify highly differentiated loci that may contribute to ecological divergence, but they do not incorporate spatial, phenotypic or environmental data that might enhance outlier detection. Geographic cline analysis provides a potential framework for integrating diverse forms of data in a spatially explicit framework, but has not been used to study genome-wide patterns of divergence. Aided by a first-draft genome assembly, we combined an FCT scan and geographic cline analysis to characterize patterns of genome-wide divergence between divergent pollination ecotypes of Mimulus aurantiacus. FCT analysis of 58872 SNPs generated via RAD-seq revealed little ecotypic differentiation (mean FCT =0.041), although a small number of loci were moderately-to-highly diverged. Consistent with our previous results from the gene MaMyb2, which contributes to differences in flower colour, 130 loci have cline shapes that recapitulate the spatial pattern of trait divergence, suggesting that they may reside in or near the genomic regions that contribute to pollinator isolation. In the narrow hybrid zone between the ecotypes, extensive admixture among individuals and low linkage disequilibrium between markers indicate that most outlier loci are scattered throughout the genome, rather than being restricted to one or a few divergent regions. In addition to revealing the genomic consequences of ecological divergence in this system, we discuss how geographic cline analysis is a powerful but under-utilized framework for studying genome-wide patterns of divergence.

Genetic and epigenetic variation among inbred mouse littermates: identification of inter-individual differentially methylated regions.

BackgroundPhenotypic variability among inbred littermates reared in controlled environments remains poorly understood. Metastable epialleles refer to loci that intrinsically behave in this way and a few examples have been described. They display differential methylation in association with differential expression. For example, inbred mice carrying the agouti viable yellow (Avy) allele show a range of coat colours associated with different DNA methylation states at the locus. The availability of next-generation sequencing, in particular whole genome sequencing of bisulphite converted DNA, allows us, for the first time, to search for metastable epialleles at base pair resolution.ResultsUsing whole genome bisulphite sequencing of DNA from the livers of five mice from the Avy colony, we searched for sites at which DNA methylation differed among the mice. A small number of loci, 356, were detected and we call these inter-individual Differentially Methylated Regions, iiDMRs, 55 of which overlap with endogenous retroviral elements (ERVs). Whole genome resequencing of two mice from the colony identified very few differences and these did not occur at or near the iiDMRs. Further work suggested that the majority of ERV iiDMRs are metastable epialleles; the level of methylation was maintained in tissue from other germ layers and the level of mRNA from the neighbouring gene inversely correlated with methylation state. Most iiDMRs that were not overlapping ERV insertions occurred at tissue-specific DMRs and it cannot be ruled out that these are driven by changes in the ratio of cell types in the tissues analysed.ConclusionsUsing the most thorough genome-wide profiling technologies for differentially methylated regions, we find very few intrinsically epigenetically variable regions that we term iiDMRs. The most robust of these are at retroviral elements and appear to be metastable epialleles. The non-ERV iiDMRs cannot be described as metastable epialleles at this stage but provide a novel class of variably methylated elements for further study.Electronic supplementary materialThe online version of this article (doi:10.1186/s13072-015-0047-z) contains supplementary material, which is available to authorized users.

"Missing" G x E Variation Controls Flowering Time in Arabidopsis thaliana.

Understanding how genetic variation interacts with the environment is essential for understanding adaptation. In particular, the life cycle of plants is tightly coordinated with local environmental signals through complex interactions with the genetic variation (G x E). The mechanistic basis for G x E is almost completely unknown. We collected flowering time data for 173 natural inbred lines of Arabidopsis thaliana from Sweden under two growth temperatures (10°C and 16°C), and observed massive G x E variation. To identify the genetic polymorphisms underlying this variation, we conducted genome-wide scans using both SNPs and local variance components. The SNP-based scan identified several variants that had common effects in both environments, but found no trace of G x E effects, whereas the scan using local variance components found both. Furthermore, the G x E effects appears to be concentrated in a small fraction of the genome (0.5%). Our conclusion is that G x E effects in this study are mostly due to large numbers of allele or haplotypes at a small number of loci, many of which correspond to previously identified flowering time genes.

On the occurrence of false positives in tests of migration under an isolation-with-migration model.

The population genetic study of divergence is often carried out using a Bayesian genealogy sampler, like those implemented in ima2 and related programs, and these analyses frequently include a likelihood ratio test of the null hypothesis of no migration between populations. Cruickshank and Hahn (2014, Molecular Ecology, 23, 3133-3157) recently reported a high rate of false-positive test results with ima2 for data simulated with small numbers of loci under models with no migration and recent splitting times. We confirm these findings and discover that they are caused by a failure of the assumptions underlying likelihood ratio tests that arises when using marginal likelihoods for a subset of model parameters. We also show that for small data sets, with little divergence between samples from two populations, an excellent fit can often be found by a model with a low migration rate and recent splitting time and a model with a high migration rate and a deep splitting time.

On the relative roles of background selection and genetic hitchhiking in shaping human cytomegalovirus genetic diversity.

A central focus of population genetics has been examining the contribution of selective and neutral processes in shaping patterns of intraspecies diversity. In terms of selection specifically, surveys of higher organisms have shown considerable variation in the relative contributions of background selection and genetic hitchhiking in shaping the distribution of polymorphisms, although these analyses have rarely been extended to bacteria and viruses. Here, we study the evolution of a ubiquitous, viral pathogen, human cytomegalovirus (HCMV), by analysing the relationship among intraspecies diversity, interspecies divergence and rates of recombination. We show that there is a strong correlation between diversity and divergence, consistent with expectations of neutral evolution. However, after correcting for divergence, there remains a significant correlation between intraspecies diversity and recombination rates, with additional analyses suggesting that this correlation is largely due to the effects of background selection. In addition, a small number of loci, centred on long noncoding RNAs, also show evidence of selective sweeps. These data suggest that HCMV evolution is dominated by neutral mechanisms as well as background selection, expanding our understanding of linked selection to a novel class of organisms.

Natural hybridization and asymmetric introgression at the distribution margin of two Buddleja species with a large overlap.

BackgroundNatural hybridization in plants is universal and plays an important role in evolution. Based on morphology it has been presumed that hybridization occurred in the genus Buddleja, though genetic studies confirming this assumption have not been conducted to date. The two species B. crispa and B. officinalis overlap in their distributions over a wide range in South-West China, and we aimed to provide genetic evidence for ongoing hybridization in this study.ResultsWe investigated the occurrence of hybrids between the two species at the southern-most edge of the distribution of B. crispa using five nuclear loci and pollination experiments. The genetic data suggest substantial differentiation between the two species as species-specific alleles are separated by at least 7–28 mutations. The natural hybrids found were nearly all F1s (21 of 23), but backcrosses were detected, and some individuals, morphologically indistinguishable from the parental species, showed introgression. Pollen viability test shows that the percentage of viable pollen grains was 50 ± 4 % for B. crispa, and 81 ± 2 % for B. officinalis. This difference is highly significant (t = 7.382, p < 0.0001). Hand cross-pollination experiments showed that B. crispa is not successful as pollen-parent, but B. officinalis is able to pollinate B. crispa to produce viable hybrid seed. Inter-specific seed-set is low (8 seeds per fruit, as opposed to about 65 for intra-specific pollinations), suggesting post-zygotic reproductive barriers. In addition, one of the reference populations also suggests a history of introgression at other localities.ConclusionsThe occurrence of morphologically intermediate individuals between B. crispa and B. officinalis at Xishan Mountain is unequivocally linked to hybridization and almost all examined individuals of the putative hybrids were likely F1s. Despite pollination experiments indicating higher chances for introgression into B. officinalis (hybrids only produced viable seed when crossed with B. officinalis), observed introgression was asymmetrical into B. crispa. This could be due to seeds produced by hybrids not contributing to seedlings, or other factors favoring the establishment of backcrosses towards B. crispa. However, further research will be needed to confirm these observations, as the small number of plants used for the pollination experiments could have introduced an artifact, for example if used individuals were more or less compatible than the species average, and also the small number of loci used could convey a picture of introgression that is not representative for the whole genome.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-015-0539-9) contains supplementary material, which is available to authorized users.

Cell-composition effects in the analysis of DNA methylation array data: a mathematical perspective.

BackgroundThe impact of cell-composition effects in analysis of DNA methylation data is now widely appreciated. With the availability of a reference data set consisting of DNA methylation measurements on isolated cell types, it is possible to impute cell proportions and adjust for them, but there is increasing interest in methods that adjust for cell composition effects when reference sets are incomplete or unavailable.ResultsIn this article we present a theoretical basis for one such method, showing that the total effect of a phenotype on DNA methylation can be decomposed into orthogonal components, one representing the effect of phenotype on proportions of major cell types, the other representing either subtle effects in composition or global effects at focused loci, and that it is possible to separate these two types of effects in a finite data set. We demonstrate this principle empirically on nine DNA methylation data sets, showing that the first few principal components generally contain a majority of the information on cell-type present in the data, but that later principal components nevertheless contain information about a small number of loci that may represent more focused associations. We also present a new method for determining the number of linear terms to interpret as cell-mixture effects and demonstrate robustness to the choice of this parameter.ConclusionsTaken together, our work demonstrates that reference-free algorithms for cell-mixture adjustment can produce biologically valid results, separating cell-mediated epigenetic effects (i.e. apparent effects arising from differences in cell composition) from those that are not cell mediated, and that in general the interpretation of associations evident from DNA methylation should be carefully considered.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-015-0527-y) contains supplementary material, which is available to authorized users.

Genetic perturbation of the maize methylome.

DNA methylation can play important roles in the regulation of transposable elements and genes. A collection of mutant alleles for 11 maize (Zea mays) genes predicted to play roles in controlling DNA methylation were isolated through forward- or reverse-genetic approaches. Low-coverage whole-genome bisulfite sequencing and high-coverage sequence-capture bisulfite sequencing were applied to mutant lines to determine context- and locus-specific effects of these mutations on DNA methylation profiles. Plants containing mutant alleles for components of the RNA-directed DNA methylation pathway exhibit loss of CHH methylation at many loci as well as CG and CHG methylation at a small number of loci. Plants containing loss-of-function alleles for chromomethylase (CMT) genes exhibit strong genome-wide reductions in CHG methylation and some locus-specific loss of CHH methylation. In an attempt to identify stocks with stronger reductions in DNA methylation levels than provided by single gene mutations, we performed crosses to create double mutants for the maize CMT3 orthologs, Zmet2 and Zmet5, and for the maize DDM1 orthologs, Chr101 and Chr106. While loss-of-function alleles are viable as single gene mutants, the double mutants were not recovered, suggesting that severe perturbations of the maize methylome may have stronger deleterious phenotypic effects than in Arabidopsis thaliana.

Establishment of regions of genomic activity during the Drosophila maternal to zygotic transition.

We describe the genome-wide distributions and temporal dynamics of nucleosomes and post-translational histone modifications throughout the maternal-to-zygotic transition in embryos of Drosophila melanogaster. At mitotic cycle 8, when few zygotic genes are being transcribed, embryonic chromatin is in a relatively simple state: there are few nucleosome free regions, undetectable levels of the histone methylation marks characteristic of mature chromatin, and low levels of histone acetylation at a relatively small number of loci. Histone acetylation increases by cycle 12, but it is not until cycle 14 that nucleosome free regions and domains of histone methylation become widespread. Early histone acetylation is strongly associated with regions that we have previously shown to be bound in early embryos by the maternally deposited transcription factor Zelda, suggesting that Zelda triggers a cascade of events, including the accumulation of specific histone modifications, that plays a role in the subsequent activation of these sequences.

BBCA: Improving the scalability of *BEAST using random binning.

BackgroundSpecies tree estimation can be challenging in the presence of gene tree conflict due to incomplete lineage sorting (ILS), which can occur when the time between speciation events is short relative to the population size. Of the many methods that have been developed to estimate species trees in the presence of ILS, *BEAST, a Bayesian method that co-estimates the species tree and gene trees given sequence alignments on multiple loci, has generally been shown to have the best accuracy. However, *BEAST is extremely computationally intensive so that it cannot be used with large numbers of loci; hence, *BEAST is not suitable for genome-scale analyses.ResultsWe present BBCA (boosted binned coalescent-based analysis), a method that can be used with *BEAST (and other such co-estimation methods) to improve scalability. BBCA partitions the loci randomly into subsets, uses *BEAST on each subset to co-estimate the gene trees and species tree for the subset, and then combines the newly estimated gene trees together using MP-EST, a popular coalescent-based summary method. We compare time-restricted versions of BBCA and *BEAST on simulated datasets, and show that BBCA is at least as accurate as *BEAST, and achieves better convergence rates for large numbers of loci.ConclusionsPhylogenomic analysis using *BEAST is currently limited to datasets with a small number of loci, and analyses with even just 100 loci can be computationally challenging. BBCA uses a very simple divide-and-conquer approach that makes it possible to use *BEAST on datasets containing hundreds of loci. This study shows that BBCA provides excellent accuracy and is highly scalable.

Multiparental mapping of plant height and flowering time QTL in partially isogenic sorghum families.

Sorghum varieties suitable for grain production at temperate latitudes show dwarfism and photoperiod insensitivity, both of which are controlled by a small number of loci with large effects. We studied the genetic control of plant height and flowering time in five sorghum families (A–E), each derived from a cross between a tropical line and a partially isogenic line carrying introgressions derived from a common, temperate-adapted donor. A total of 724 F2:3 lines were phenotyped in temperate and tropical environments for plant height and flowering time and scored at 9139 SNPs using genotyping-by-sequencing. Biparental mapping was compared with multiparental mapping in different subsets of families (AB, ABC, ABCD, and ABCDE) using both a GWAS approach, which fit each QTL as a single effect across all families, and using a joint linkage approach, which fit QTL effects as nested within families. GWAS using all families (ABCDE) performed best at the cloned Dw3 locus, whereas joint linkage using all families performed best at the cloned Ma1 locus. Both multiparental approaches yielded apparently synthetic associations due to genetic heterogeneity and were highly dependent on the subset of families used. Comparison of all mapping approaches suggests that a GA2-oxidase underlies Dw1, and that a mir172a gene underlies a Dw1-linked flowering time QTL.

Transmission of epi-alleles with MET1-dependent dense methylation in Arabidopsis thaliana.

DNA methylation in plants targets cytosines in three sequence contexts, CG, CHG and CHH (H representing A, C or T). Each of these patterns has traditionally been associated with distinct DNA methylation pathways with CHH methylation being controlled by the RNA dependent DNA methylation (RdDM) pathway employing small RNAs as a guide for the de novo DOMAINS REARRANGED METHYLTRANSFERASE (DRM2), and maintenance DNA METHYLTRANSFERASE1 (MET1) being responsible for faithful propagation of CG methylation. Here we report an unusual ‘dense methylation’ pattern under the control of MET1, with methylation in all three sequence contexts. We identified epi-alleles of dense methylation at a non coding RNA locus (At4g15242) in Arabidopsis ecotypes, with distinct dense methylation and expression characteristics, which are stably maintained and transmitted in genetic crosses and which can be heritably altered by depletion of MET1. This suggests that, in addition to its classical CG maintenance function, at certain loci MET1 plays a role in creating transcriptional diversity based on the generation of independent epi-alleles. Database inspection identified several other loci with MET1-dependent dense methylation patterns. Arabidopsis ecotypes contain distinct epi-alleles of these loci with expression patterns that inversely correlate with methylation density, predominantly within the transcribed region. In Arabidopsis, dense methylation appears to be an exception as it is only found at a small number of loci. Its presence does, however, highlight the potential for MET1 as a contributor to epigenetic diversity, and it will be interesting to investigate the representation of dense methylation in other plant species.

Genetic Differentiation in Response to Selection for Water-Soluble Carbohydrate Content in Perennial Ryegrass (Lolium perenne L.)

Plant carbohydrates are of increasing interest as renewable feedstocks to replace petrochemicals in the generation of fuels and production of high-value chemicals. Greater understanding of the genetic control of diversity in fructan synthesis and accumulation would facilitate more directed channelling of feedstock to process in a ryegrass biorefinery. Divergent populations produced by phenotypic selection for water-soluble carbohydrate content have been used to investigate relationships between traits, and to identify patterns of genetic differentiation which indicate genomic regions under high and low selection pressure. Selection for high water-soluble carbohydrate content was associated with increased synthesis of large fructan polymers and increased accumulation of above-ground plant biomass, particularly during spring. Three rounds of selection and two rounds of recombination resulted in widespread genetic differentiation across the whole genome, causing reduced allelic richness and increasing homozygosity at some loci. A smaller number of loci were shown to be subject to high selection pressure. Breeding material subjected to many years of selection for water-soluble carbohydrate also showed allelic differences which may reflect the consequences of high selection pressure at some of these same loci. However, some of the loci unaffected in the divergent selection experiment showed similar effects. This might arise from differences in linkage disequilibrium in these two sets of plant materials, but more likely from the different genetic background of the germplasm. This illustrates the complex nature of the water-soluble carbohydrate trait in perennial ryegrass.

Sensory drive speciation and patterns of variation at selectively neutral genes

Speciation by sensory drive can occur if divergent adaptation of sensory systems causes rapid evolution of mating traits and the resulting development of assortative mating. Previous theoretical studies have shown that sensory drive can cause rapid divergent adaptive evolution from one to two phenotypes. In this study, we examined two topics: the possibility of adaptive radiation by sensory drive from one to more than two phenotypes and the relationships of patterns of variation at selectively neutral genes to levels of viability selection, habitat and mating preferences and migration. We conducted individual-based simulations assuming a sensory trait and a mating trait controlled by a small number of loci. We found that adaptive radiation is possible when the number of loci controlling the sensory trait is small; the levels of viability selection, habitat and mating preferences are intermediate; and the emigration rate is high. We also found that emigration rates as well as the levels of habitat and mating preferences are related to F ST values at neutral loci, but F ST proved to be insensitive to a small change in the number of loci controlling the mating trait. This suggests that an estimation of the past population history is possible without an accurate genetic model.

Future directions of forensic DNA databases.

Forensic DNA databases are indispensable tools of the law enforcement system. The purpose of establishing forensic DNA databases was to develop investigative leads for solving crime and usually was the purview of “criminal justice agencies for law enforcement identification purposes” (1). The forensic DNA databases of most countries generally contain two types of profiles: 1) reference profiles from convicted offenders and/or arrestee profiles; these profiles are from known sources; and 2) forensic profiles which derive from crime scenes and are characteristically from unknown sources. In a typical database search, an unknown forensic profile is searched against the convicted offender and arrestee profiles (or can be searched against other unknown forensic profiles) to determine if an association, often called a match or hit, can be found. The hit can be used to develop investigation leads. As of May 2013, China and the United States (US) maintain the two largest forensic DNA databases, containing more than 20 and 12 million profiles and have produced over 410 000 (2) and 185 000 hits (3), respectively. In addition to direct matching between known and unknown sample profiles, profiles from missing persons and their relatives, as well as unidentified human remains, are included in a number of databases (2,3). Missing person identification also is an invaluable module for investigating certain crimes. For example, as of June 2013, China has successfully identified and rescued 2455 trafficked children through the use of its DNA database (2). The US National Missing and Unidentified Persons System (NamUs), which uses other meta-data has successfully solved 3499 missing persons cases as of August 2012 (4). As expected with the great success of the use of forensic DNA databases, new challenges are emerging. The databases are experiencing rapid growth, and thus there is a potential of increased adventitious hits; the power for current and new applications (eg, missing person identification and familial searching) require additional infrastructure support; and there is an increased desire for international data sharing (5-7), which possibly could be hampered if only a relatively small number of loci is shared among laboratories worldwide.

Genome-Wide Analysis of Selection on the Malaria Parasite Plasmodium falciparum in West African Populations of Differing Infection Endemicity

Locally varying selection on pathogens may be due to differences in drug pressure, host immunity, transmission opportunities between hosts, or the intensity of between-genotype competition within hosts. Highly recombining populations of the human malaria parasite Plasmodium falciparum throughout West Africa are closely related, as gene flow is relatively unrestricted in this endemic region, but markedly varying ecology and transmission intensity should cause distinct local selective pressures. Genome-wide analysis of sequence variation was undertaken on a sample of 100 P. falciparum clinical isolates from a highly endemic region of the Republic of Guinea where transmission occurs for most of each year and compared with data from 52 clinical isolates from a previously sampled population from The Gambia, where there is relatively limited seasonal malaria transmission. Paired-end short-read sequences were mapped against the 3D7 P. falciparum reference genome sequence, and data on 136,144 single nucleotide polymorphisms (SNPs) were obtained. Within-population analyses identifying loci showing evidence of recent positive directional selection and balancing selection confirm that antimalarial drugs and host immunity have been major selective agents. Many of the signatures of recent directional selection reflected by standardized integrated haplotype scores were population specific, including differences at drug resistance loci due to historically different antimalarial use between the countries. In contrast, both populations showed a similar set of loci likely to be under balancing selection as indicated by very high Tajima’s D values, including a significant overrepresentation of genes expressed at the merozoite stage that invades erythrocytes and several previously validated targets of acquired immunity. Between-population FST analysis identified exceptional differentiation of allele frequencies at a small number of loci, most markedly for five SNPs covering a 15-kb region within and flanking the gdv1 gene that regulates the early stages of gametocyte development, which is likely related to the extreme differences in mosquito vector abundance and seasonality that determine the transmission opportunities for the sexual stage of the parasite.

Effective combination of xenoprotective transgenes

Abundant expression of a series of xenoprotective transgenes is vital for clinically useful xeno‐donor pigs. Breeding of such multi‐transgenic animals requires that xeno‐transgenes be placed at a single, or a small number of loci in the pig genome to avoid or minimise separation by genetic segregation. We are investigating methods that support transgene expression and also enable transgenes to be “stacked” at a single locus. These include the use of BAC based multi‐transgene vectors, and very recently serial transgene placement at a permissive locus by homologous recombination supported by TALENs [1].BAC vector constructs containing genomic sequences for the complement activation inhibitors CD55, CD46 and CD59 were assessed to determine the minimum size without reduction of expression level. These constructs enable us to express all biologically active CD55 splice variants (membrane bound and soluble forms). Effective xenoprotection also requires ubiquitous expression. We are thus examining endogenous promoter sequences of various lengths as well as the widely used CAGGs (CMV enhancer/chicken β actin) promoter to direct complement inhibitory gene expression levels in all porcine tissues.To investigate means of coexpressing groups of transgenes that provide xeno‐protection at different levels, we developed a series of BAC vector constructs containing the CD55 genomic sequence plus two additional xeno‐transgene cDNAs, including A20 [2], HO1 [3], thrombomodulin [4] and CTLA4‐Ig (LEA29Y) [5].Analysis of porcine adipose MSC stably transfected cell clones revealed expression of all inserted xeno‐transgenes most importantly high expression of complement inhibitory genes (Fig. 1). This demonstrates that xeno‐transgenes can be effectively combined in a single construct, thus minimising the number of transgene loci in the final donor pigs. The huge carrying capacity of BAC vectors enables the addition of further xeno‐transgenes in a single vector construct.We have further examined the system in vivo. Primary cells transfected with up to three BAC constructs were used for nuclear transfer, fetuses explanted and examined for expression with very promising results.

Development of a Hierarchical Variable-Number Tandem Repeat Typing Scheme for Mycobacterium tuberculosis in China

Molecular typing based on variable-number tandem repeats (VNTR) analysis is a promising tool for identifying transmission of Mycobacterium tuberculosis. However, the currently proposed 15- and 24-locus VNTR sets (VNTR-15/24) only have limited resolution and contain too many loci for large-scale typing in high burden countries. To develop an optimal typing scheme in China, we evaluated the resolution and robustness of 25 VNTR loci, using population-based collections of 1362 clinical isolates from six provinces across the country. The resolution of most loci showed considerable variations among regions. By calculating the average resolution of all possible combinations of 20 robust loci, we identified an optimal locus set with a minimum of 9 loci (VNTR-9) that could achieve comparable resolution of the standard VNTR-15. The VNTR-9 had consistently high resolutions in all six regions, and it was highly concordant with VNTR-15 for defining both clustered and unique genotypes. Furthermore, VNTR-9 was phylogenetically informative for classifying lineages/sublineages of M. tuberculosis. Three hypervariable loci (HV-3), VNTR 3232, VNTR 3820 and VNTR 4120, were proved important for further differentiating unrelated clustered strains based on VNTR-9. We propose the optimized VNTR-9 as first-line method and the HV-3 as second-line method for molecular typing of M. tuberculosis in China and surrounding countries. The development of hierarchical VNTR typing methods that can achieve high resolution with a small number of loci could be suitable for molecular epidemiology study in other high burden countries.